#include <darcy_flow_mh.hh>
Definition at line 406 of file darcy_flow_mh.hh.
void P1_CouplingAssembler::add_sides |
( |
const Element * |
ele, |
|
|
unsigned int |
shift, |
|
|
vector< int > & |
dofs, |
|
|
vector< double > & |
dirichlet |
|
) |
| |
void P1_CouplingAssembler::assembly |
( |
LinSys & |
ls | ) |
|
P1 coonection of different dimensions
- demonstrated convergence, but still major open questions: ? in all test cases the error on the fracture is less on the left and greater on the right with incresing trend tried:
- changed order of 1d elements (no change)
- higher precision of ngh output and linear solver (no change) ? seems that there should be some factor 6.0 in the communication term ? in the case of infinite k2 -> simplest 1d-constant communication, the biggest difference on borders, numerical solution greater then analytical
TODO:
- full implementation of Dirichlet BC ( also for 2d sides)
Definition at line 721 of file darcy_flow_mh.cc.
vector<double> P1_CouplingAssembler::dirichlet |
|
private |
vector<int> P1_CouplingAssembler::dofs |
|
private |
arma::vec P1_CouplingAssembler::rhs |
|
private |
The documentation for this class was generated from the following files: