19 #include <dune/common/parallel/mpihelper.hh>
20 #include <dune/common/parametertreeparser.hh>
21 #include <dune/common/shared_ptr.hh>
22 #include <dune/common/classname.hh>
23 #include <dune/common/exceptions.hh>
24 #include <dune/common/fvector.hh>
26 #include <dune/geometry/type.hh>
27 #include <dune/geometry/quadraturerules.hh>
29 #include <dune/grid/sgrid.hh>
30 #include <dune/grid/onedgrid.hh>
31 #include <dune/grid/io/file/vtk/subsamplingvtkwriter.hh>
32 #include <dune/grid/yaspgrid.hh>
34 #include <dune/grid/uggrid.hh>
37 #include<dune/grid/albertagrid.hh>
38 #include <dune/grid/albertagrid/dgfparser.hh>
41 #include<dune/grid/uggrid.hh>
44 #include<dune/grid/alugrid.hh>
45 #include<dune/grid/io/file/dgfparser/dgfalu.hh>
46 #include<dune/grid/io/file/dgfparser/dgfparser.hh>
48 #include <dune/grid/utility/structuredgridfactory.hh>
49 #include <dune/grid/io/file/gmshreader.hh>
51 #include <dune/istl/bvector.hh>
52 #include <dune/istl/operators.hh>
53 #include <dune/istl/solvers.hh>
54 #include <dune/istl/solvercategory.hh>
55 #include <dune/istl/preconditioners.hh>
56 #include <dune/istl/io.hh>
58 #include <dune/istl/paamg/amg.hh>
96 typedef typename T::ctype
ctype;
97 static const int dim = T::dimension;
103 FieldVector<ctype,dimworld> lowerLeft(0.0);
104 FieldVector<ctype,dimworld> upperRight(1.0);
105 array<unsigned int,dim> elements; elements.fill(cells);
107 StructuredGridFactory<T> factory;
109 if (meshtype==Dune::GeometryType::cube)
110 gridp = factory.createCubeGrid(lowerLeft,upperRight,elements);
111 else if (meshtype==Dune::GeometryType::simplex)
112 gridp = factory.createSimplexGrid(lowerLeft,upperRight,elements);
115 DUNE_THROW(GridError, className<StructuredGrid>()
116 <<
"::StructuredGrid(): grid type must be simplex or cube ");
122 array<double,dimworld> lower_left, array<double,dimworld> upper_right,
123 array<unsigned int,dim> cells)
125 FieldVector<ctype,dimworld> lowerLeft;
126 FieldVector<ctype,dimworld> upperRight;
127 array<unsigned int,dim> elements;
132 lowerLeft[i] = lower_left[i];
133 upperRight[i] = upper_right[i];
135 for (
size_t i=0; i<
dim; i++)
137 elements[i] = cells[i];
140 StructuredGridFactory<T> factory;
142 if (meshtype==Dune::GeometryType::cube)
143 gridp = factory.createCubeGrid(lowerLeft,upperRight,elements);
144 else if (meshtype==Dune::GeometryType::simplex)
145 gridp = factory.createSimplexGrid(lowerLeft,upperRight,elements);
148 DUNE_THROW(GridError, className<StructuredGrid>()
149 <<
"::StructuredGrid(): grid type must be simplex or cube ");
178 return gridp.operator->();
188 return gridp.operator->();
193 Dune::shared_ptr<T> gridp;
208 StructuredGrid (Dune::GeometryType::BasicType meshtype,
unsigned int cells,
int overlap=1)
211 if (meshtype!=Dune::GeometryType::cube)
212 std::cout <<
"StructuredGrid(): element type " << meshtype <<
" is ignored" << std::endl;
215 Dune::FieldVector<double,dimworld> L(1.0);
216 Dune::array<int,dimworld> N(Dune::fill_array<int,dimworld>(cells));
217 std::bitset<dimworld> B(
false);
221 gridp = shared_ptr<Grid>(
new Grid(Dune::MPIHelper::getCommunicator(),L,N,B,overlap));
223 gridp = shared_ptr<Grid>(
new Grid(L,N,B,overlap));
229 array<double,dimworld> lower_left, array<double,dimworld> upper_right,
230 array<unsigned int,dim> cells,
int overlap=1)
234 if(std::abs(lower_left[d]) > std::abs(upper_right[d])*1
e-10)
235 DUNE_THROW(GridError, className<StructuredGrid>()
236 <<
"::createCubeGrid(): The lower coordinates "
237 "must be at the origin for YaspGrid.");
240 if (meshtype!=Dune::GeometryType::cube)
241 std::cout <<
"StructuredGrid(): element type " << meshtype <<
" is ignored" << std::endl;
244 Dune::FieldVector<double,dimworld> L;
245 Dune::array<int,dimworld> N;
246 std::bitset<dimworld> B(
false);
249 L[i] = upper_right[i];
255 gridp = shared_ptr<Grid>(
new Grid(Dune::MPIHelper::getCommunicator(),L,N,B,overlap));
257 gridp = shared_ptr<Grid>(
new Grid(L,N,B,overlap));
263 array<double,dimworld> lower_left, array<double,dimworld> upper_right,
264 array<unsigned int,dim> cells, array<bool,dim> periodic,
int overlap=1)
268 if(std::abs(lower_left[d]) > std::abs(upper_right[d])*1
e-10)
269 DUNE_THROW(GridError, className<StructuredGrid>()
270 <<
"::createCubeGrid(): The lower coordinates "
271 "must be at the origin for YaspGrid.");
274 if (meshtype!=Dune::GeometryType::cube)
275 std::cout <<
"StructuredGrid(): element type " << meshtype <<
" is ignored" << std::endl;
278 Dune::FieldVector<double,dimworld> L;
279 Dune::array<int,dimworld> N;
280 std::bitset<dimworld> B(
false);
283 L[i] = upper_right[i];
290 gridp = shared_ptr<Grid>(
new Grid(Dune::MPIHelper::getCommunicator(),L,N,B,overlap));
292 gridp = shared_ptr<Grid>(
new Grid(L,N,B,overlap));
321 return gridp.operator->();
331 return gridp.operator->();
335 Dune::shared_ptr<Grid> gridp;
346 static const int dim = T::dimension;
350 UnstructuredGrid (std::string filename,
bool verbose =
true,
bool insert_boundary_segments=
true)
352 Dune::GridFactory<T> factory;
353 Dune::GmshReader<T>::read(factory,filename,verbose,insert_boundary_segments);
354 gridp = shared_ptr<T>(factory.createGrid());
382 return gridp.operator->();
392 return gridp.operator->();
396 Dune::shared_ptr<T> gridp;
405 template<
typename GV,
typename C,
typename R,
unsigned int degree,
unsigned int dim, Dune::GeometryType::BasicType gt>
409 template<
typename GV,
typename C,
typename R,
unsigned int degree,
unsigned int dim>
417 femp = shared_ptr<FEM>(
new FEM(gridview));
424 shared_ptr<FEM> femp;
427 template<
typename GV,
typename C,
typename R,
unsigned int degree,
unsigned int dim>
435 femp = shared_ptr<FEM>(
new FEM(gridview));
442 shared_ptr<FEM> femp;
454 template<
typename Gr
id,
unsigned int degree, Dune::GeometryType::BasicType gt, MeshType mt, SolverCategory::Category st,
typename BCType,
typename GV =
typename Gr
id::LeafGr
idView>
458 template<
typename Gr
id,
typename BCType,
typename GV>
464 CGCONBase (Grid& grid,
const BCType& bctype,
const GV& gv)
466 conp = shared_ptr<CON>(
new CON(grid,
true,bctype));
471 conp = shared_ptr<CON>(
new CON(grid,
true,bctype));
474 template<
typename GFS>
478 template<
typename GFS,
typename DOF>
481 shared_ptr<CON> conp;
484 template<
typename Gr
id,
typename BCType,
typename GV>
490 CGCONBase (Grid& grid,
const BCType& bctype,
const GV& gv)
492 conp = shared_ptr<CON>(
new CON(grid,
true,bctype));
497 conp = shared_ptr<CON>(
new CON(grid,
true,bctype));
500 template<
typename GFS>
504 template<
typename GFS,
typename DOF>
507 shared_ptr<CON> conp;
510 template<
typename Gr
id,
unsigned int degree, Dune::GeometryType::BasicType gt,
typename BCType,
typename GV>
516 CGCONBase (Grid& grid,
const BCType& bctype,
const GV& gv)
518 conp = shared_ptr<CON>(
new CON());
523 conp = shared_ptr<CON>(
new CON());
526 template<
typename GFS>
530 template<
typename GFS,
typename DOF>
533 shared_ptr<CON> conp;
536 template<
typename Gr
id,
unsigned int degree, Dune::GeometryType::BasicType gt,
typename BCType,
typename GV>
542 CGCONBase (Grid& grid,
const BCType& bctype,
const GV& gv)
544 conp = shared_ptr<CON>(
new CON());
549 conp = shared_ptr<CON>(
new CON());
552 template<
typename GFS>
556 template<
typename GFS,
typename DOF>
563 if (gfs.gridView().comm().size()>1)
564 gfs.gridView().communicate(adddh,Dune::InteriorBorder_All_Interface,Dune::ForwardCommunication);
567 shared_ptr<CON> conp;
570 template<
typename Gr
id,
unsigned int degree, Dune::GeometryType::BasicType gt,
typename BCType,
typename GV>
576 CGCONBase (Grid& grid,
const BCType& bctype,
const GV& gv)
578 conp = shared_ptr<CON>(
new CON(gv));
583 conp = shared_ptr<CON>(
new CON(grid.leafGridView()));
586 template<
typename GFS>
590 template<
typename GFS,
typename DOF>
593 shared_ptr<CON> conp;
598 template<
typename T,
typename N,
unsigned int degree,
typename BCType,
599 Dune::GeometryType::BasicType gt,
MeshType mt, SolverCategory::Category st = SolverCategory::sequential,
600 typename VBET=ISTLVectorBackend<> >
606 typedef typename T::LeafGridView
GV;
608 static const int dim = T::dimension;
614 typedef typename FEMB::FEM
FEM;
615 typedef typename CONB::CON
CON;
623 typedef typename GFS::template ConstraintsContainer<N>::Type
CC;
628 : gv(grid.leafGridView()), femb(gv), conb(grid,bctype)
630 gfsp = shared_ptr<GFS>(
new GFS(gv,femb.getFEM(),conb.getCON()));
631 gfsp->name(
"cgspace");
634 conb.postGFSHook(*gfsp);
635 ccp = shared_ptr<CC>(
new CC());
640 return femb.getFEM();
645 return femb.getFEM();
686 conb.make_consistent(*gfsp,x);
692 conb.make_consistent(*gfsp,x);
698 conb.make_consistent(*gfsp,xout);
704 conb.make_consistent(*gfsp,xout);
711 shared_ptr<GFS> gfsp;
721 template<SolverCategory::Category st>
732 conp = shared_ptr<CON>(
new CON());
736 template<
typename GFS,
typename DOF>
739 shared_ptr<CON> conp;
749 conp = shared_ptr<CON>(
new CON());
753 template<
typename GFS,
typename DOF>
756 shared_ptr<CON> conp;
766 conp = shared_ptr<CON>(
new CON());
770 template<
typename GFS,
typename DOF>
777 if (gfs.gridView().comm().size()>1)
778 gfs.gridView().communicate(adddh,Dune::InteriorBorder_All_Interface,Dune::ForwardCommunication);
781 shared_ptr<CON> conp;
786 template<
typename T,
typename N,
unsigned int degree,
787 Dune::GeometryType::BasicType gt, SolverCategory::Category st = SolverCategory::sequential,
795 typedef typename T::LeafGridView
GV;
797 static const int dim = T::dimension;
806 typedef typename CONB::CON
CON;
811 typedef typename GFS::template ConstraintsContainer<N>::Type
CC;
817 femp = shared_ptr<FEM>(
new FEM());
818 gfsp = shared_ptr<GFS>(
new GFS(gv,*femp));
821 ccp = shared_ptr<CC>(
new CC());
839 template<
class BCTYPE>
854 conb.make_consistent(*gfsp,x);
860 conb.make_consistent(*gfsp,x);
866 conb.make_consistent(*gfsp,xout);
872 conb.make_consistent(*gfsp,xout);
878 shared_ptr<FEM> femp;
879 shared_ptr<GFS> gfsp;
885 template<
typename T,
typename N,
unsigned int degree,
886 Dune::GeometryType::BasicType gt, SolverCategory::Category st = SolverCategory::sequential,
888 typename VBET=ISTLVectorBackend<> >
895 typedef typename T::LeafGridView
GV;
897 static const int dim = T::dimension;
906 typedef typename CONB::CON
CON;
911 typedef typename GFS::template ConstraintsContainer<N>::Type
CC;
917 femp = shared_ptr<FEM>(
new FEM());
918 gfsp = shared_ptr<GFS>(
new GFS(gv,*femp));
921 ccp = shared_ptr<CC>(
new CC());
939 template<
class BCTYPE>
954 conb.make_consistent(*gfsp,x);
960 conb.make_consistent(*gfsp,x);
966 conb.make_consistent(*gfsp,xout);
972 conb.make_consistent(*gfsp,xout);
978 shared_ptr<FEM> femp;
979 shared_ptr<GFS> gfsp;
985 template<
typename T,
typename N,
unsigned int degree,
986 Dune::GeometryType::BasicType gt, SolverCategory::Category st = SolverCategory::sequential,
994 typedef typename T::LeafGridView
GV;
996 static const int dim = T::dimension;
1001 typedef typename CONB::CON
CON;
1006 typedef typename GFS::template ConstraintsContainer<N>::Type
CC;
1012 femp = shared_ptr<FEM>(
new FEM());
1013 gfsp = shared_ptr<GFS>(
new GFS(gv,*femp));
1016 ccp = shared_ptr<CC>(
new CC());
1034 template<
class BCTYPE>
1049 conb.make_consistent(*gfsp,x);
1055 conb.make_consistent(*gfsp,x);
1061 conb.make_consistent(*gfsp,xout);
1067 conb.make_consistent(*gfsp,xout);
1073 shared_ptr<FEM> femp;
1074 shared_ptr<GFS> gfsp;
1080 template<
typename T,
typename N,
unsigned int degree,
1081 Dune::GeometryType::BasicType gt, SolverCategory::Category st = SolverCategory::sequential,
1083 typename VBET=ISTLVectorBackend<> >
1090 typedef typename T::LeafGridView
GV;
1092 static const int dim = T::dimension;
1097 typedef typename CONB::CON
CON;
1102 typedef typename GFS::template ConstraintsContainer<N>::Type
CC;
1108 femp = shared_ptr<FEM>(
new FEM());
1109 gfsp = shared_ptr<GFS>(
new GFS(gv,*femp));
1112 ccp = shared_ptr<CC>(
new CC());
1130 template<
class BCTYPE>
1145 conb.make_consistent(*gfsp,x);
1151 conb.make_consistent(*gfsp,x);
1157 conb.make_consistent(*gfsp,xout);
1163 conb.make_consistent(*gfsp,xout);
1169 shared_ptr<FEM> femp;
1170 shared_ptr<GFS> gfsp;
1177 template<
typename T,
typename N,
1178 Dune::GeometryType::BasicType gt, SolverCategory::Category st = SolverCategory::sequential,
1179 typename VBET=ISTLVectorBackend<> >
1186 typedef typename T::LeafGridView
GV;
1188 static const int dim = T::dimension;
1193 typedef typename CONB::CON
CON;
1198 typedef typename GFS::template ConstraintsContainer<N>::Type
CC;
1204 femp = shared_ptr<FEM>(
new FEM(Dune::GeometryType(gt,
dim)));
1205 gfsp = shared_ptr<GFS>(
new GFS(gv,*femp));
1208 ccp = shared_ptr<CC>(
new CC());
1226 template<
class BCTYPE>
1241 conb.make_consistent(*gfsp,x);
1247 conb.make_consistent(*gfsp,x);
1253 conb.make_consistent(*gfsp,xout);
1259 conb.make_consistent(*gfsp,xout);
1265 shared_ptr<FEM> femp;
1266 shared_ptr<GFS> gfsp;
1273 template<
typename FS,
typename Functor>
1275 :
public GridFunctionBase<GridFunctionTraits<typename FS::GV, typename FS::NT,
1276 1,FieldVector<typename FS::NT,1> >
1277 ,UserFunction<FS,Functor> >
1294 std::vector<double> x__(x.size());
1295 for (
size_t i=0; i<x.size(); ++i) x__[i]=x_[i];
1301 return fs.getGFS().gridView();
1310 template<
typename FS,
typename LOP, SolverCategory::Category st = SolverCategory::sequential>
1317 typename FS::NT,
typename FS::NT,
typename FS::NT,
1318 typename FS::CC,
typename FS::CC>
GO;
1323 gop = shared_ptr<GO>(
new GO(fs.getGFS(),fs.getCC(),fs.getGFS(),fs.getCC(),lop));
1345 return gop.operator->();
1355 return gop.operator->();
1363 template<
typename FS,
typename LOP>
1370 typename FS::NT,
typename FS::NT,
typename FS::NT,
1371 typename FS::CC,
typename FS::CC,
true>
GO;
1376 gop = shared_ptr<GO>(
new GO(fs.getGFS(),fs.getCC(),fs.getGFS(),fs.getCC(),lop));
1398 return gop.operator->();
1408 return gop.operator->();
1416 template<
typename FS,
typename LOP, SolverCategory::Category st = SolverCategory::sequential>
1423 typename FS::NT,
typename FS::NT,
typename FS::NT,
1424 typename FS::CC,
typename FS::CC>
GO;
1429 gop = shared_ptr<GO>(
new GO(fs.getGFS(),fs.getCC(),fs.getGFS(),fs.getCC(),lop,mbe));
1451 return gop.operator->();
1461 return gop.operator->();
1469 template<
typename FS,
typename LOP>
1476 typename FS::NT,
typename FS::NT,
typename FS::NT,
1477 typename FS::CC,
typename FS::CC,
true>
GO;
1482 gop = shared_ptr<GO>(
new GO(fs.getGFS(),fs.getCC(),fs.getGFS(),fs.getCC(),lop,mbe));
1504 return gop.operator->();
1514 return gop.operator->();
1525 template<
typename FSU,
typename FSV,
typename LOP, SolverCategory::Category st>
1532 typename FSU::NT,
typename FSU::NT,
typename FSU::NT,
1533 typename FSU::CC,
typename FSV::CC>
GO;
1538 gop = shared_ptr<GO>(
new GO(fsu.getGFS(),fsu.getCC(),fsv.getGFS(),fsv.getCC(),lop));
1560 return gop.operator->();
1570 return gop.operator->();
1578 template<
typename FSU,
typename FSV,
typename LOP>
1585 typename FSU::NT,
typename FSU::NT,
typename FSU::NT,
1586 typename FSU::CC,
typename FSV::CC,
true>
GO;
1591 gop = shared_ptr<GO>(
new GO(fsu.getGFS(),fsu.getCC(),fsv.getGFS(),fsv.getCC(),lop));
1613 return gop.operator->();
1623 return gop.operator->();
1632 template<
typename GO1,
typename GO2,
bool implicit = true>
1643 gop = shared_ptr<GO>(
new GO(*go1,*go2));
1665 return gop.operator->();
1675 return gop.operator->();
1684 template<
typename FS,
typename ASS, SolverCategory::Category st = SolverCategory::sequential>
1692 int verbose_=1,
bool reuse_=
false,
bool usesuperlu_=
true)
1694 lsp = shared_ptr<LS>(
new LS(maxiter_,verbose_,reuse_,usesuperlu_));
1709 template<
typename FS,
typename ASS>
1717 int verbose_=1,
bool reuse_=
false,
bool usesuperlu_=
true)
1719 lsp = shared_ptr<LS>(
new LS(fs.getGFS(),maxiter_,verbose_,reuse_,usesuperlu_));
1734 template<
typename FS,
typename ASS>
1742 int verbose_=1,
bool reuse_=
false,
bool usesuperlu_=
true)
1744 lsp = shared_ptr<LS>(
new LS(fs.getGFS(),maxiter_,verbose_,reuse_,usesuperlu_));
1759 template<
typename FS,
typename ASS, SolverCategory::Category st = SolverCategory::sequential>
1767 int steps_=5,
int verbose_=1)
1769 lsp = shared_ptr<LS>(
new LS(maxiter_,verbose_));
1784 template<
typename FS,
typename ASS>
1792 int steps_=5,
int verbose_=1)
1794 lsp = shared_ptr<LS>(
new LS(fs.getGFS(),maxiter_,steps_,verbose_));
1809 template<
typename FS,
typename ASS>
1817 int steps_=5,
int verbose_=1)
1819 lsp = shared_ptr<LS>(
new LS(fs.getGFS(),fs.getCC(),maxiter_,steps_,verbose_));
1836 template<
typename FS,
typename ASS, SolverCategory::Category st = SolverCategory::sequential>
1845 lsp = shared_ptr<LS>(
new LS(maxiter_,verbose_));
1860 template<
typename FS,
typename ASS>
1869 lsp = shared_ptr<LS>(
new LS(fs.getGFS(),maxiter_,verbose_));
1884 template<
typename FS,
typename ASS>
1893 lsp = shared_ptr<LS>(
new LS(fs.getGFS(),fs.getCC(),maxiter_,3,verbose_));
1909 template<
typename FS,
typename ASS, SolverCategory::Category st = SolverCategory::sequential>
1918 lsp = shared_ptr<LS>(
new LS());
1934 template<
typename FS,
typename ASS>
1943 lsp = shared_ptr<LS>(
new LS(fs.getGFS()));
1959 template<
typename FS,
typename ASS>
1968 lsp = shared_ptr<LS>(
new LS(fs.getGFS()));
ISTLSolverBackend_CG_SSOR(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int steps_=5, int verbose_=1)
Definition: pdelab.hh:1766
DGPkSpace(const GV &gridview)
Definition: pdelab.hh:815
void setNonConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:689
const FEM & getFEM() const
Definition: pdelab.hh:421
const T * operator->() const
Definition: pdelab.hh:390
void clearConstraints()
Definition: pdelab.hh:1233
const LS & getLS() const
Definition: pdelab.hh:1698
void maskForeignDOFs(X &x) const
Mask out all DOFs not owned by the current process with 0.
Definition: parallelhelper.hh:106
const T & operator*() const
Definition: pdelab.hh:181
void make_consistent(const GFS &gfs, DOF &x) const
Definition: pdelab.hh:754
const FS::GV & getGridView() const
Definition: pdelab.hh:1299
T::LeafGridView GV
Definition: pdelab.hh:895
T::ctype ctype
Definition: pdelab.hh:896
T::ctype ctype
Definition: pdelab.hh:1091
const LS * operator->() const
Definition: pdelab.hh:1853
void make_consistent(const GFS &gfs, DOF &x) const
Definition: pdelab.hh:505
const GO & operator*() const
Definition: pdelab.hh:1507
const CON & getCON() const
Definition: pdelab.hh:752
LS & getLS()
Definition: pdelab.hh:1747
HangingNodesDirichletConstraints< Grid, HangingNodesConstraintsAssemblers::CubeGridQ1Assembler, BCType > CON
Definition: pdelab.hh:488
GO * operator->()
Definition: pdelab.hh:1502
GO::Jacobian MAT
Definition: pdelab.hh:1478
Definition: pdelab.hh:455
Dune::PDELab::OneStepGridOperator< typename GO1::GO, typename GO2::GO, implicit > GO
Definition: pdelab.hh:1638
VBET VBE
Definition: pdelab.hh:807
void evaluate(const typename Traits::ElementType &e, const typename Traits::DomainType &x, typename Traits::RangeType &y) const
Evaluate the GridFunction at given position.
Definition: pdelab.hh:1289
const CON & getCON() const
Definition: pdelab.hh:735
DGCONBase< st > CONB
Definition: pdelab.hh:1096
CON & getCON()
Definition: pdelab.hh:476
CGCONBase(Grid &grid, const BCType &bctype)
Definition: pdelab.hh:495
DGCONBase< st > CONB
Definition: pdelab.hh:1000
void setConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:851
Definition: noconstraints.hh:16
T * operator->()
Definition: pdelab.hh:380
LS & operator*()
Definition: pdelab.hh:1923
const GO & operator*() const
Definition: pdelab.hh:1348
StructuredGrid(Dune::GeometryType::BasicType meshtype, array< double, dimworld > lower_left, array< double, dimworld > upper_right, array< unsigned int, dim > cells)
Definition: pdelab.hh:121
OPBLocalFiniteElementMap< ctype, NT, degree, dim, gt > FEM
Definition: pdelab.hh:803
static const int dimworld
Definition: pdelab.hh:609
BackendVectorSelector< GFS, N >::Type DOF
Definition: pdelab.hh:909
OneStepGlobalAssembler(GO1 &go1, GO2 &go2)
Definition: pdelab.hh:1641
const LS * operator->() const
Definition: pdelab.hh:1901
Dune::PDELab::GridOperator< typename FS::GFS, typename FS::GFS, LOP, MBE, typename FS::NT, typename FS::NT, typename FS::NT, typename FS::CC, typename FS::CC > GO
Definition: pdelab.hh:1424
CON & getCON()
Definition: pdelab.hh:751
void make_consistent(const GFS &gfs, DOF &x) const
Definition: pdelab.hh:557
GFS & getGFS()
Definition: pdelab.hh:1215
LS & operator*()
Definition: pdelab.hh:1724
GO & getGO()
Definition: pdelab.hh:1327
void setNonConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:857
const FEM & getFEM() const
Definition: pdelab.hh:439
Definition: l2orthonormal.hh:247
const LS & getLS() const
Definition: pdelab.hh:1972
GalerkinGlobalAssembler(const FS &fs, LOP &lop)
Definition: pdelab.hh:1374
void setConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:1238
const GO * operator->() const
Definition: pdelab.hh:1512
CC & getCC()
Definition: pdelab.hh:661
DGQkOPBSpace(const GV &gridview)
Definition: pdelab.hh:915
T Grid
Definition: pdelab.hh:794
const CON & getCON() const
Definition: pdelab.hh:555
void setNonConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:1148
LS & operator*()
Definition: pdelab.hh:1948
CGFEMBase(const GV &gridview)
Definition: pdelab.hh:415
Solver to be used for explicit time-steppers with (block-)diagonal mass matrix.
Definition: novlpistlsolverbackend.hh:629
const LS & operator*() const
Definition: pdelab.hh:1751
GridFunctionTraits< typename FS::GV, typename FS::NT, 1, FieldVector< typename FS::NT, 1 > > Traits
Definition: pdelab.hh:1281
CON & getCON()
Definition: pdelab.hh:528
CONB::CON CON
Definition: pdelab.hh:806
void copyConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:1250
CC & getCC()
Definition: pdelab.hh:1125
DGCONBase< st > CONB
Definition: pdelab.hh:905
const GO & getGO() const
Definition: pdelab.hh:1439
LS * operator->()
Definition: pdelab.hh:1949
const GFS & getGFS() const
Definition: pdelab.hh:1122
const T & operator*() const
Definition: pdelab.hh:385
const CON & getCON() const
Definition: pdelab.hh:477
void postGFSHook(const GFS &gfs)
Definition: pdelab.hh:587
const CON & getCON() const
Definition: pdelab.hh:529
T::LeafGridView GV
Definition: pdelab.hh:795
LS & getLS()
Definition: pdelab.hh:1797
T Grid
Definition: pdelab.hh:894
Parallel P0 constraints for nonoverlapping grids with ghosts.
Definition: p0ghost.hh:16
N NT
Definition: pdelab.hh:998
CGCONBase(Grid &grid, const BCType &bctype)
Definition: pdelab.hh:547
static const int dimworld
Definition: pdelab.hh:347
void assembleConstraints(const BCTYPE &bctype)
Definition: pdelab.hh:1227
LS * operator->()
Definition: pdelab.hh:1924
GFS::template ConstraintsContainer< N >::Type CC
Definition: pdelab.hh:911
void clearConstraints()
Definition: pdelab.hh:846
DGCONBase()
Definition: pdelab.hh:730
void setConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:951
MeshType
Definition: pdelab.hh:448
Definition: pdelab.hh:1526
CONB::CON CON
Definition: pdelab.hh:906
const LS & operator*() const
Definition: pdelab.hh:1975
BackendVectorSelector< GFS, N >::Type DOF
Definition: pdelab.hh:1004
const FEM & getFEM() const
Definition: pdelab.hh:1020
GO * operator->()
Definition: pdelab.hh:1663
LS & getLS()
Definition: pdelab.hh:1896
GO::Jacobian MAT
Definition: pdelab.hh:1534
LS * operator->()
Definition: pdelab.hh:1700
GV::Traits::template Codim< 0 >::Entity ElementType
codim 0 entity
Definition: function.hh:118
StructuredGrid(Dune::GeometryType::BasicType meshtype, unsigned int cells)
Definition: pdelab.hh:101
LS & operator*()
Definition: pdelab.hh:1824
LS & getLS()
Definition: pdelab.hh:1946
const Grid * operator->() const
Definition: pdelab.hh:329
LS * operator->()
Definition: pdelab.hh:1775
void assembleConstraints(const BCTYPE &bctype)
Definition: pdelab.hh:1035
static const int dim
Definition: pdelab.hh:346
PkLocalFiniteElementMap< GV, C, R, degree > FEM
Definition: pdelab.hh:413
Nonoverlapping parallel BiCGStab solver with Jacobi preconditioner.
Definition: novlpistlsolverbackend.hh:554
const Grid & operator*() const
Definition: pdelab.hh:324
LS & operator*()
Definition: pdelab.hh:1973
void setNonConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:957
T Grid
Definition: pdelab.hh:993
Definition: pdelab.hh:1633
const LS & operator*() const
Definition: pdelab.hh:1950
T::ctype ctype
Definition: pdelab.hh:1187
HangingNodesDirichletConstraints< Grid, HangingNodesConstraintsAssemblers::SimplexGridP1Assembler, BCType > CON
Definition: pdelab.hh:462
void copyConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:1154
ISTLBackend_SEQ_CG_SSOR LS
Definition: pdelab.hh:1764
OPBLocalFiniteElementMap< ctype, NT, degree, dim, gt, N, Dune::PB::BasisType::Qk > FEM
Definition: pdelab.hh:903
P0ParallelConstraints CON
Definition: pdelab.hh:763
LS & getLS()
Definition: pdelab.hh:1697
const CC & getCC() const
Definition: pdelab.hh:837
const CC & getCC() const
Definition: pdelab.hh:1224
void postGFSHook(const GFS &gfs)
Definition: pdelab.hh:553
const LS * operator->() const
Definition: pdelab.hh:1827
Overlapping parallel CGS solver with SSOR preconditioner.
Definition: ovlpistlsolverbackend.hh:637
LS * operator->()
Definition: pdelab.hh:1800
DGCONBase< st > CONB
Definition: pdelab.hh:1192
ConformingDirichletConstraints CON
Definition: pdelab.hh:514
GridFunctionSpace< GV, FEM, CON, VBE > GFS
Definition: pdelab.hh:908
void set_constrained_dofs(const CG &cg, typename XG::ElementType x, XG &xg)
construct constraints from given boundary condition function
Definition: common/constraints.hh:772
const GFS & getGFS() const
Definition: pdelab.hh:931
void setNonConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:1244
const GO & operator*() const
Definition: pdelab.hh:1454
Backend using ISTL matrices.
Definition: istl/descriptors.hh:69
CON & getCON()
Definition: pdelab.hh:768
void clearConstraints()
Definition: pdelab.hh:946
LS & operator*()
Definition: pdelab.hh:1874
GFS::template ConstraintsContainer< N >::Type CC
Definition: pdelab.hh:623
void copyNonConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:701
void postGFSHook(const GFS &gfs)
Definition: pdelab.hh:527
Dune::PDELab::GridOperator< typename FS::GFS, typename FS::GFS, LOP, MBE, typename FS::NT, typename FS::NT, typename FS::NT, typename FS::CC, typename FS::CC > GO
Definition: pdelab.hh:1318
ISTLBackend_OVLP_BCGS_SSORk< typename FS::GFS, typename FS::CC > LS
Definition: pdelab.hh:1889
Dune::shared_ptr< T > getSharedPtr()
Definition: pdelab.hh:358
const GO & getGO() const
Definition: pdelab.hh:1601
ISTLSolverBackend_ExplicitDiagonal(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1)
Definition: pdelab.hh:1941
Definition: pdelab.hh:1760
GO::Jacobian MAT
Definition: pdelab.hh:1319
CONB::CON CON
Definition: pdelab.hh:615
const LS & operator*() const
Definition: pdelab.hh:1726
T Grid
Definition: pdelab.hh:95
Backend using (possibly nested) ISTL BCRSMatrices.
Definition: bcrsmatrixbackend.hh:187
Dune::PDELab::istl::BCRSMatrixBackend MBE
Definition: pdelab.hh:1474
GridFunctionSpace< GV, FEM, CON, VBE > GFS
Definition: pdelab.hh:1099
Definition: pdelab.hh:889
GO & getGO()
Definition: pdelab.hh:1647
void copy_nonconstrained_dofs(const CG &cg, const XG &xgin, XG &xgout)
Definition: common/constraints.hh:976
GO & getGO()
Definition: pdelab.hh:1486
Dune::PDELab::ISTLBackend_CG_AMG_SSOR< typename ASS::GO > LS
Definition: pdelab.hh:1739
const LS & getLS() const
Definition: pdelab.hh:1823
VTKGridFunctionAdapter< DGF > VTKF
Definition: pdelab.hh:812
static const int dim
Definition: adaptivity.hh:82
Grid * operator->()
Definition: pdelab.hh:319
Backend for sequential conjugate gradient solver with SSOR preconditioner.
Definition: seqistlsolverbackend.hh:333
static const int dim
Definition: pdelab.hh:797
Nonoverlapping parallel CG solver preconditioned by block SSOR.
Definition: novlpistlsolverbackend.hh:842
void copyNonConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:1064
static const int dimworld
Definition: pdelab.hh:1093
BackendVectorSelector< GFS, N >::Type DOF
Definition: pdelab.hh:1100
void clearConstraints()
Definition: pdelab.hh:1137
VBET VBE
Definition: pdelab.hh:907
Standard grid operator implementation.
Definition: gridoperator.hh:34
LS & getLS()
Definition: pdelab.hh:1722
GFS & getGFS()
Definition: pdelab.hh:1023
const GFS & getGFS() const
Definition: pdelab.hh:1026
StructuredGrid(Dune::GeometryType::BasicType meshtype, array< double, dimworld > lower_left, array< double, dimworld > upper_right, array< unsigned int, dim > cells, int overlap=1)
Definition: pdelab.hh:228
FEMB::FEM FEM
Definition: pdelab.hh:614
BackendVectorSelector< GFS, N >::Type DOF
Definition: pdelab.hh:1196
const GFS & getGFS() const
Definition: pdelab.hh:1218
Dune::FieldVector< GV::Grid::ctype, GV::dimension > DomainType
domain type in dim-size coordinates
Definition: function.hh:49
Definition: pdelab.hh:1910
N NT
Definition: pdelab.hh:1094
Grid & operator*()
Definition: pdelab.hh:314
const CC & getCC() const
Definition: pdelab.hh:937
ISTLSolverBackend_ExplicitDiagonal(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1)
Definition: pdelab.hh:1966
Dune::PDELab::DiscreteGridFunction< GFS, DOF > DGF
Definition: pdelab.hh:810
const LS * operator->() const
Definition: pdelab.hh:1976
VTKGridFunctionAdapter< DGF > VTKF
Definition: pdelab.hh:624
LS & operator*()
Definition: pdelab.hh:1898
GlobalAssembler(const FSU &fsu, const FSV &fsv, LOP &lop)
Definition: pdelab.hh:1589
CC & getCC()
Definition: pdelab.hh:1221
LS * operator->()
Definition: pdelab.hh:1899
const LS & getLS() const
Definition: pdelab.hh:1897
CGCONBase(Grid &grid, const BCType &bctype)
Definition: pdelab.hh:469
QkLocalFiniteElementMap< GV, C, R, degree > FEM
Definition: pdelab.hh:431
const GFS & getGFS() const
Definition: pdelab.hh:655
ISTLSolverBackend_CG_SSOR(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int steps_=5, int verbose_=1)
Definition: pdelab.hh:1791
CON & getCON()
Definition: pdelab.hh:734
static const int dimworld
Definition: pdelab.hh:798
const GO * operator->() const
Definition: pdelab.hh:1621
Dune::PDELab::istl::BCRSMatrixBackend MBE
Definition: pdelab.hh:1421
GO::Jacobian MAT
Definition: pdelab.hh:1372
Definition: pdelab.hh:1417
Nonoverlapping parallel CG solver preconditioned with AMG smoothed by SSOR.
Definition: novlpistlsolverbackend.hh:1045
ISTLMatrixBackend MBE
Definition: pdelab.hh:1315
UserFunction(const FS &fs_, const Functor &f_)
constructor
Definition: pdelab.hh:1284
LS & operator*()
Definition: pdelab.hh:1774
V & raw(V &v)
Returns the raw ISTL object associated with v, or v itself it is already an ISTL object.
Definition: backend/istl/utility.hh:26
BackendVectorSelector< GFS, N >::Type DOF
Definition: pdelab.hh:809
Definition: pdelab.hh:1837
const GO & operator*() const
Definition: pdelab.hh:1616
VTKGridFunctionAdapter< DGF > VTKF
Definition: pdelab.hh:1103
VBET VBE
Definition: pdelab.hh:617
Definition: pdelab.hh:1311
GO & getGO()
Definition: pdelab.hh:1542
const LS & operator*() const
Definition: pdelab.hh:1900
CGCONBase(Grid &grid, const BCType &bctype)
Definition: pdelab.hh:521
Definition: qkdggl.hh:353
void setConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:1142
LS & getLS()
Definition: pdelab.hh:1872
ISTLSolverBackend_ExplicitDiagonal(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1)
Definition: pdelab.hh:1916
GO & operator*()
Definition: pdelab.hh:1338
void setConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:683
VTKGridFunctionAdapter< DGF > VTKF
Definition: pdelab.hh:912
StructuredGrid(Dune::GeometryType::BasicType meshtype, unsigned int cells, int overlap=1)
Definition: pdelab.hh:208
FEM & getFEM()
Definition: pdelab.hh:420
CON & getCON()
Definition: pdelab.hh:588
N NT
Definition: pdelab.hh:799
void make_consistent(const GFS &gfs, DOF &x) const
Definition: pdelab.hh:479
CONB::CON CON
Definition: pdelab.hh:1001
GO * operator->()
Definition: pdelab.hh:1343
CGSpace(Grid &grid, const BCType &bctype)
Definition: pdelab.hh:627
Dune::shared_ptr< T > getSharedPtr()
Definition: pdelab.hh:154
GFS & getGFS()
Definition: pdelab.hh:828
ISTLSolverBackend_CG_AMG_SSOR(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1, bool reuse_=false, bool usesuperlu_=true)
Definition: pdelab.hh:1741
void clearConstraints()
Definition: pdelab.hh:678
T::ctype ctype
Definition: pdelab.hh:345
GO::Jacobian MAT
Definition: pdelab.hh:1639
CGCONBase< Grid, degree, gt, mt, st, BCType > CONB
Definition: pdelab.hh:612
GO * operator->()
Definition: pdelab.hh:1611
ISTLBackend_OVLP_CG_SSORk< typename FS::GFS, typename FS::CC > LS
Definition: pdelab.hh:1814
const FEM & getFEM() const
Definition: pdelab.hh:925
CON & getCON()
Definition: pdelab.hh:554
const LS & getLS() const
Definition: pdelab.hh:1922
LS & getLS()
Definition: pdelab.hh:1822
const CON & getCON() const
Definition: pdelab.hh:503
void copyConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:1058
const FEM & getFEM() const
Definition: pdelab.hh:1116
N NT
Definition: pdelab.hh:899
T::LeafGridView GV
Definition: pdelab.hh:1186
GO & operator*()
Definition: pdelab.hh:1658
T Grid
Definition: pdelab.hh:1089
FEM & getFEM()
Definition: pdelab.hh:638
GO & operator*()
Definition: pdelab.hh:1497
P0Space(const GV &gridview)
Definition: pdelab.hh:1202
const GO * operator->() const
Definition: pdelab.hh:1459
GO & operator*()
Definition: pdelab.hh:1553
Parallel P0 constraints for overlapping grids.
Definition: p0.hh:15
void constraints(const GFS &gfs, CG &cg, const bool verbose=false)
construct constraints
Definition: common/constraints.hh:723
LS & getLS()
Definition: pdelab.hh:1772
const LS & operator*() const
Definition: pdelab.hh:1776
Definition: pdelab.hh:988
T::ctype ctype
Definition: pdelab.hh:995
R RangeType
range type
Definition: function.hh:61
VBET VBE
Definition: pdelab.hh:1098
CONB::CON CON
Definition: pdelab.hh:1097
Definition: pdelab.hh:1685
leaf of a function tree
Definition: function.hh:577
void make_consistent(const GFS &gfs, DOF &x) const
Definition: pdelab.hh:771
Definition: pdelab.hh:1084
Definition: pdelab.hh:340
const CC & getCC() const
Definition: pdelab.hh:1128
GlobalAssembler(const FSU &fsu, const FSV &fsv, LOP &lop)
Definition: pdelab.hh:1536
static const int dim
Definition: pdelab.hh:996
T::ctype ctype
Definition: pdelab.hh:796
LS & operator*()
Definition: pdelab.hh:1699
GO & getGO()
Definition: pdelab.hh:1433
void copyConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:963
const LS & operator*() const
Definition: pdelab.hh:1826
const LS * operator->() const
Definition: pdelab.hh:1777
Definition: parallelhelper.hh:45
void copyConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:695
Dune::PDELab::P0LocalFiniteElementMap< ctype, NT, dim > FEM
Definition: pdelab.hh:1191
const LS & getLS() const
Definition: pdelab.hh:1849
const LS * operator->() const
Definition: pdelab.hh:1877
T * operator->()
Definition: pdelab.hh:176
const LS & getLS() const
Definition: pdelab.hh:1723
CGCONBase(Grid &grid, const BCType &bctype, const GV &gv)
Definition: pdelab.hh:490
Overlapping parallel conjugate gradient solver preconditioned with AMG smoothed by SSOR...
Definition: ovlpistlsolverbackend.hh:1066
const GO & getGO() const
Definition: pdelab.hh:1548
N NT
Definition: pdelab.hh:620
void assembleConstraints(const BCTYPE &bctype)
Definition: pdelab.hh:840
const GO & operator*() const
Definition: pdelab.hh:1563
LS * operator->()
Definition: pdelab.hh:1875
Definition: genericdatahandle.hh:623
void copyNonConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:969
static const int dim
Definition: pdelab.hh:97
GalerkinGlobalAssemblerNewBackend(const FS &fs, LOP &lop, const MBE &mbe)
Definition: pdelab.hh:1427
static const int dimworld
Definition: pdelab.hh:898
Dune::PDELab::ISTLBackend_SEQ_CG_AMG_SSOR< typename ASS::GO > LS
Definition: pdelab.hh:1689
LS * operator->()
Definition: pdelab.hh:1750
GFS::template ConstraintsContainer< N >::Type CC
Definition: pdelab.hh:1102
BackendVectorSelector< GFS, N >::Type DOF
Definition: pdelab.hh:621
StructuredGrid(Dune::GeometryType::BasicType meshtype, array< double, dimworld > lower_left, array< double, dimworld > upper_right, array< unsigned int, dim > cells, array< bool, dim > periodic, int overlap=1)
Definition: pdelab.hh:262
GO & operator*()
Definition: pdelab.hh:1606
T::LeafGridView GV
Definition: pdelab.hh:606
ISTLMatrixBackend MBE
Definition: pdelab.hh:1583
const T & getGrid() const
Definition: pdelab.hh:166
const LS * operator->() const
Definition: pdelab.hh:1752
P0ParallelGhostConstraints CON
Definition: pdelab.hh:746
Sequential conjugate gradient solver preconditioned with AMG smoothed by SSOR.
Definition: seqistlsolverbackend.hh:613
Dune::PDELab::ISTLBackend_OVLP_ExplicitDiagonal< typename FS::GFS > LS
Definition: pdelab.hh:1939
Traits::Jacobian Jacobian
Definition: gridoperator/onestep.hh:53
const GO & operator*() const
Definition: pdelab.hh:1401
static const unsigned int value
Definition: gridfunctionspace/tags.hh:175
VBET VBE
Definition: pdelab.hh:1002
Dune::PDELab::GridOperator< typename FSU::GFS, typename FSV::GFS, LOP, MBE, typename FSU::NT, typename FSU::NT, typename FSU::NT, typename FSU::CC, typename FSV::CC, true > GO
Definition: pdelab.hh:1586
Definition: l2orthonormal.hh:247
LS * operator->()
Definition: pdelab.hh:1725
const GO & getGO() const
Definition: pdelab.hh:1333
static const int dimworld
Definition: pdelab.hh:98
Dune::PDELab::ISTLBackend_NOVLP_ExplicitDiagonal< typename FS::GFS > LS
Definition: pdelab.hh:1964
Dune::PDELab::DiscreteGridFunction< GFS, DOF > DGF
Definition: pdelab.hh:910
GalerkinGlobalAssemblerNewBackend(const FS &fs, LOP &lop, const MBE &mbe)
Definition: pdelab.hh:1480
void copyNonConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:1160
Dune::PDELab::ISTLBackend_NOVLP_CG_AMG_SSOR< typename ASS::GO > LS
Definition: pdelab.hh:1714
BackendVectorSelectorHelper< Backend, GridFunctionSpace, FieldType >::Type Type
Definition: backendselector.hh:14
DGCONBase< st > CONB
Definition: pdelab.hh:805
ISTLSolverBackend_CG_AMG_SSOR(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1, bool reuse_=false, bool usesuperlu_=true)
Definition: pdelab.hh:1716
const LS * operator->() const
Definition: pdelab.hh:1926
Definition: pdelab.hh:450
const T & getGrid() const
Definition: pdelab.hh:370
DGCONBase()
Definition: pdelab.hh:764
ISTLMatrixBackend MBE
Definition: pdelab.hh:1368
T & operator*()
Definition: pdelab.hh:171
const GFS & getGFS() const
Definition: pdelab.hh:831
LS * operator->()
Definition: pdelab.hh:1974
Solver to be used for explicit time-steppers with (block-)diagonal mass matrix.
Definition: ovlpistlsolverbackend.hh:831
const LS & getLS() const
Definition: pdelab.hh:1798
void setNonConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:1052
const CC & getCC() const
Definition: pdelab.hh:1032
Dune::PDELab::GridOperator< typename FS::GFS, typename FS::GFS, LOP, MBE, typename FS::NT, typename FS::NT, typename FS::NT, typename FS::CC, typename FS::CC, true > GO
Definition: pdelab.hh:1477
NonoverlappingConformingDirichletConstraints< GV > CON
Definition: pdelab.hh:574
const FEM & getFEM() const
Definition: pdelab.hh:643
LS * operator->()
Definition: pdelab.hh:1825
UnstructuredGrid(std::string filename, bool verbose=true, bool insert_boundary_segments=true)
Definition: pdelab.hh:350
ISTLSolverBackend_IterativeDefault(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1)
Definition: pdelab.hh:1867
CGFEMBase(const GV &gridview)
Definition: pdelab.hh:433
T Grid
Definition: pdelab.hh:605
GO & getGO()
Definition: pdelab.hh:1595
static const int dimworld
Definition: pdelab.hh:1189
GridFunctionSpace< GV, FEM, CON, VBE > GFS
Definition: pdelab.hh:1195
T::LeafGridView GV
Definition: pdelab.hh:1090
Dune::PDELab::DiscreteGridFunction< GFS, DOF > DGF
Definition: pdelab.hh:1101
traits class holding the function signature, same as in local function
Definition: function.hh:176
ISTLBackend_NOVLP_CG_SSORk< typename ASS::GO > LS
Definition: pdelab.hh:1789
GO & operator*()
Definition: pdelab.hh:1391
Definition: pdelab.hh:722
GalerkinGlobalAssembler(const FS &fs, LOP &lop)
Definition: pdelab.hh:1321
const LS * operator->() const
Definition: pdelab.hh:1951
FEM & getFEM()
Definition: pdelab.hh:1019
GridFunctionSpace< GV, FEM, CON, VBE > GFS
Definition: pdelab.hh:808
const LS * operator->() const
Definition: pdelab.hh:1802
FEM & getFEM()
Definition: pdelab.hh:924
T::ctype ctype
Definition: pdelab.hh:96
DGQkSpace(const GV &gridview)
Definition: pdelab.hh:1010
void make_consistent(const GFS &gfs, DOF &x) const
Definition: pdelab.hh:737
A grid function space.
Definition: gridfunctionspace.hh:109
const GO * operator->() const
Definition: pdelab.hh:1673
const LS & getLS() const
Definition: pdelab.hh:1873
ISTLMatrixBackend MBE
Definition: pdelab.hh:1637
FEM & getFEM()
Definition: pdelab.hh:1211
GO * operator->()
Definition: pdelab.hh:1449
LS & getLS()
Definition: pdelab.hh:1971
YaspGrid< dim > Grid
Definition: pdelab.hh:203
CC & getCC()
Definition: pdelab.hh:1029
Definition: pdelab.hh:406
LS & getLS()
Definition: pdelab.hh:1921
Dune::PDELab::DiscreteGridFunction< GFS, DOF > DGF
Definition: pdelab.hh:622
void assembleConstraints(const BCType &bctype)
Definition: pdelab.hh:672
const LS & getLS() const
Definition: pdelab.hh:1947
CC & getCC()
Definition: pdelab.hh:934
const LS * operator->() const
Definition: pdelab.hh:1727
GFS & getGFS()
Definition: pdelab.hh:928
wrap a GridFunction so it can be used with the VTKWriter from dune-grid.
Definition: common/vtkexport.hh:22
CGCONBase(Grid &grid, const BCType &bctype, const GV &gv)
Definition: pdelab.hh:516
Definition: pdelab.hh:449
extend conforming constraints class by processor boundary
Definition: conforming.hh:152
Dune::PDELab::DiscreteGridFunction< GFS, DOF > DGF
Definition: pdelab.hh:1005
void copy_constrained_dofs(const CG &cg, const XG &xgin, XG &xgout)
Definition: common/constraints.hh:917
const GO & getGO() const
Definition: pdelab.hh:1653
void make_consistent(const GFS &gfs, DOF &x) const
Definition: pdelab.hh:591
GO::Jacobian MAT
Definition: pdelab.hh:1587
T Grid
Definition: pdelab.hh:1185
Hanging Node constraints construction.
Definition: hangingnode.hh:320
static const int dim
Definition: pdelab.hh:1092
CGFEMBase< GV, ctype, N, degree, dim, gt > FEMB
Definition: pdelab.hh:611
VTKGridFunctionAdapter< DGF > VTKF
Definition: pdelab.hh:1007
const CON & getCON() const
Definition: pdelab.hh:769
const CON & getCON() const
Definition: pdelab.hh:589
GFS::template ConstraintsContainer< N >::Type CC
Definition: pdelab.hh:1198
T Grid
Definition: pdelab.hh:344
VTKGridFunctionAdapter< DGF > VTKF
Definition: pdelab.hh:1199
T::ctype ctype
Definition: pdelab.hh:607
T::LeafGridView GV
Definition: pdelab.hh:994
GFS & getGFS()
Definition: pdelab.hh:1119
const GO * operator->() const
Definition: pdelab.hh:1568
void make_consistent(const GFS &gfs, DOF &x) const
Definition: pdelab.hh:531
Definition: pdelab.hh:601
const FEM & getFEM() const
Definition: pdelab.hh:1212
LS * operator->()
Definition: pdelab.hh:1851
static const int dim
Definition: pdelab.hh:608
ISTLSolverBackend_CG_AMG_SSOR(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1, bool reuse_=false, bool usesuperlu_=true)
Definition: pdelab.hh:1691
const LS & operator*() const
Definition: pdelab.hh:1876
void postGFSHook(const GFS &gfs)
Definition: pdelab.hh:501
LS & operator*()
Definition: pdelab.hh:1850
CGCONBase(Grid &grid, const BCType &bctype, const GV &gv)
Definition: pdelab.hh:542
VBET VBE
Definition: pdelab.hh:1194
const Grid & getGrid() const
Definition: pdelab.hh:309
OverlappingConformingDirichletConstraints CON
Definition: pdelab.hh:540
GO * operator->()
Definition: pdelab.hh:1396
LS & operator*()
Definition: pdelab.hh:1749
const GO * operator->() const
Definition: pdelab.hh:1353
Backend for sequential BiCGSTAB solver with SSOR preconditioner.
Definition: seqistlsolverbackend.hh:246
const LS & operator*() const
Definition: pdelab.hh:1701
N NT
Definition: pdelab.hh:1190
ISTLSolverBackend_IterativeDefault(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1)
Definition: pdelab.hh:1843
QkDGGLLocalFiniteElementMap< ctype, NT, degree, dim > FEM
Definition: pdelab.hh:1095
FEM & getFEM()
Definition: pdelab.hh:1115
Definition: gridoperator/onestep.hh:14
LS & getLS()
Definition: pdelab.hh:1848
ISTLSolverBackend_CG_SSOR(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int steps_=5, int verbose_=1)
Definition: pdelab.hh:1816
const GO & getGO() const
Definition: pdelab.hh:1386
GO & operator*()
Definition: pdelab.hh:1444
NoConstraints CON
Definition: pdelab.hh:729
void copyConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:863
T & getGrid()
Definition: pdelab.hh:160
extend conforming constraints class by processor boundary
Definition: conforming.hh:101
FEM & getFEM()
Definition: pdelab.hh:824
const LS & getLS() const
Definition: pdelab.hh:1773
const GO & operator*() const
Definition: pdelab.hh:1668
GridFunctionSpace< GV, FEM, CON, VBE > GFS
Definition: pdelab.hh:618
Definition: pdelab.hh:789
GridFunctionSpace< GV, FEM, CON, VBE > GFS
Definition: pdelab.hh:1003
GO::Jacobian MAT
Definition: pdelab.hh:1425
Dune::shared_ptr< Grid > getSharedPtr()
Definition: pdelab.hh:297
const LS & getLS() const
Definition: pdelab.hh:1748
const E & e
Definition: interpolate.hh:172
void copyNonConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:869
Grid::ctype ctype
Definition: pdelab.hh:204
CONB::CON CON
Definition: pdelab.hh:1193
DGQkGLSpace(const GV &gridview)
Definition: pdelab.hh:1106
const FEM & getFEM() const
Definition: pdelab.hh:825
const LS & operator*() const
Definition: pdelab.hh:1801
Dune::PDELab::GridOperator< typename FSU::GFS, typename FSV::GFS, LOP, MBE, typename FSU::NT, typename FSU::NT, typename FSU::NT, typename FSU::CC, typename FSV::CC > GO
Definition: pdelab.hh:1533
GFS & getGFS()
Definition: pdelab.hh:649
Dune::PDELab::DiscreteGridFunction< GFS, DOF > DGF
Definition: pdelab.hh:1197
void assembleConstraints(const BCTYPE &bctype)
Definition: pdelab.hh:940
T & operator*()
Definition: pdelab.hh:375
void assembleConstraints(const BCTYPE &bctype)
Definition: pdelab.hh:1131
CON & getCON()
Definition: pdelab.hh:502
const LS & operator*() const
Definition: pdelab.hh:1925
void copyNonConstrainedDOFS(const DOF &xin, DOF &xout) const
Definition: pdelab.hh:1256
convert a grid function space and a coefficient vector into a grid function
Definition: gridfunctionspaceutilities.hh:54
GO * operator->()
Definition: pdelab.hh:1558
void postGFSHook(const GFS &gfs)
Definition: pdelab.hh:475
CGCONBase(Grid &grid, const BCType &bctype, const GV &gv)
Definition: pdelab.hh:576
CGCONBase(Grid &grid, const BCType &bctype, const GV &gv)
Definition: pdelab.hh:464
ISTLBackend_SEQ_BCGS_SSOR LS
Definition: pdelab.hh:1841
Grid & getGrid()
Definition: pdelab.hh:303
static const int dimworld
Definition: pdelab.hh:997
const CC & getCC() const
Definition: pdelab.hh:667
Definition: pdelab.hh:1180
static const int dim
Definition: pdelab.hh:897
GFS::template ConstraintsContainer< N >::Type CC
Definition: pdelab.hh:811
const T * operator->() const
Definition: pdelab.hh:186
QkDGLocalFiniteElementMap< ctype, NT, degree, dim > FEM
Definition: pdelab.hh:999
DGCONBase()
Definition: pdelab.hh:747
Dune::PDELab::BackendMatrixSelector< MB, Domain, Range, JF >::Type Jacobian
The type of the jacobian.
Definition: gridoperator.hh:46
ISTLSolverBackend_IterativeDefault(const FS &fs, const ASS &ass, unsigned maxiter_=5000, int verbose_=1)
Definition: pdelab.hh:1891
ISTLMatrixBackend MBE
Definition: pdelab.hh:1530
void clearConstraints()
Definition: pdelab.hh:1041
const LS * operator->() const
Definition: pdelab.hh:1702
Overlapping parallel BiCGStab solver with SSOR preconditioner.
Definition: ovlpistlsolverbackend.hh:571
void setConstrainedDOFS(DOF &x, NT nt) const
Definition: pdelab.hh:1046
GFS::template ConstraintsContainer< N >::Type CC
Definition: pdelab.hh:1006
GO & getGO()
Definition: pdelab.hh:1380
Dune::PDELab::GridOperator< typename FS::GFS, typename FS::GFS, LOP, MBE, typename FS::NT, typename FS::NT, typename FS::NT, typename FS::CC, typename FS::CC, true > GO
Definition: pdelab.hh:1371
void set_nonconstrained_dofs(const CG &cg, typename XG::ElementType x, XG &xg)
Definition: common/constraints.hh:943
CC & getCC()
Definition: pdelab.hh:834
CGCONBase(Grid &grid, const BCType &bctype)
Definition: pdelab.hh:581
const GO & getGO() const
Definition: pdelab.hh:1492
T & getGrid()
Definition: pdelab.hh:364
Definition: pdelab.hh:1274
LS & operator*()
Definition: pdelab.hh:1799
static const int dim
Definition: pdelab.hh:1188
Solver to be used for explicit time-steppers with (block-)diagonal mass matrix.
Definition: seqistlsolverbackend.hh:433
ISTLBackend_NOVLP_BCGS_Jacobi< typename FS::GFS > LS
Definition: pdelab.hh:1865
Dune::PDELab::ISTLBackend_SEQ_ExplicitDiagonal LS
Definition: pdelab.hh:1914
const GO * operator->() const
Definition: pdelab.hh:1406
const LS & operator*() const
Definition: pdelab.hh:1852
FEM & getFEM()
Definition: pdelab.hh:438
Dirichlet Constraints construction.
Definition: conforming.hh:36