You are here: Concepts>Concepts Web>Class documentation

Class documentation of Concepts

Loading...
Searching...
No Matches
concepts.pypp.hh
1#ifndef concepts_pypp_hh
2#define concepts_pypp_hh
3
4#include "basics/typedefs.hh"
5#include "basics/defines.hh"
6#include "basics/vectorsMatrices.hh"
7#include "basics/outputOperator.hh"
8#include "basics/cloneable.hh"
9#include "basics/level.hh"
10#include "basics/parallel.hh" // isParallelRunning()
11#include "basics/Zm.hh"
12
13#include "toolbox/array.hh"
14#include "toolbox/set.hh"
15#include "toolbox/inputOutput.hh"
16#include "toolbox/sequence.hh"
17#include "toolbox/elementGraphics.hh" //ElementGraphicsBase, ElementGraphics
18#include "toolbox/stringFunc.hh"
19#include "toolbox/resultsTable.hh"
20
21#include "integration/quadRule.hh"
22#include "integration/shapefunction.hh"
23
24#include "formula/elementFormula.hh"
25#include "formula/formula.hh"
26#include "formula/elementFormulaContainer.hh"
27#include "formula/parsedFormula.hh"
28#include "formula/boundary.hh"
29#include "formula/piecewiseFormula.hh"
30#include "formula/frmE_sum.hh"
31#include "formula/frmE_product.hh"
32#include "formula/vectorElementFormula.hh"
33#include "formula/bessel.hh"
34#include "formula/zylindricalFormula.hh"
35
36#include "geometry/connector.hh"
37#include "geometry/connectorSet.hh"
38#include "geometry/mesh.hh"
39#include "geometry/meshImport2D.hh"
40#include "geometry/meshImport2Dgeneral.hh"
41#include "geometry/meshImport2Dgmsh.hh"
42#include "geometry/boundaryConditions.hh"
43#include "geometry/square.hh"
44#include "geometry/semantics.hh"
45#include "geometry/cell.hh"
46#include "geometry/cell2D.hh"
47#include "geometry/topology.hh" //Quad
48#include "geometry/subdivision.hh" //::Subdivision
49#include "geometry/normalVector.hh"
50#include "geometry/elementMaps.hh" //concepts::Map2d
51#include "geometry/mutableMesh.hh" //mutableMesh, buildEdgeMesh
52#include "geometry/meshRelations.hh" // CellToCellMapping
53#include "geometry/normalVector.hh"
54#include "geometry/formula.hh"
55#include "geometry/frmE_normalVector.hh"
56#include "geometry/diameter.hh"
57
58#include "space/space.hh"
59#include "space/element.hh"
60#include "space/elementPairs.hh"
61#include "space/function.hh"
62#include "space/tmatrix.hh"
63#include "space/spacePreBuilder.hh"
64#include "space/spaceSet.hh" //concepts::Set< concepts:: IndexRange >
65#include "space/postProcess.hh" //concepts::CellPostprocess<>
66#include "space/dirichlet.hh" //Dirichlet<>
67#include "space/hpMethod.hh"
68#include "space/integral.hh" //l2 norms
69#include "space/spaceGraph.hh" //SpaceGraph
70
71#include "function/vector.hh"
72#include "function/basis.hh"
73#include "function/elementFormula.hh"
74#include "function/linearForm.hh"
75
76#include "operator/matrix.hh"
77#include "operator/hashedSMatrixIterator.hh"
78#include "operator/sparseMatrix.hh"
79#include "operator/bilinearForm.hh"
80#include "operator/cg.hh"
81#include "operator/compositions.hh"
82#include "operator/solverFabric.hh"
83#ifdef HAS_MUMPS
84#include "operator/mumps.hh"
85#endif
86#ifdef HAS_SuperLU
87#include "operator/superLU.hh"
88#include "operator/superLUfabric.hh"
89#endif
90#include "operator/cgFabric.hh"
91#include "operator/diagonalFabric.hh"
92#include "operator/gmresFabric.hh"
93#include "operator/denseMatrix.hh"
94
95#include "hp1D/element.hh"
96#include "hp2D/hpAdaptiveSpaceH1.hh"
97
98#include "matfile/matfileIO.hh"
99
100#include "estimator/estimator.hh"
101#include "estimator/explicitResidual.hh"
102
103inline void export_concepts_templates(){
104
105#ifdef GEN_PYTHON
106
107 //basics
108
109 std::cout << sizeof(concepts::Point<concepts::Real,1>);
110 std::cout << sizeof(concepts::Point<concepts::Cmplx,1>);
111 std::cout << sizeof(concepts::UnitNd<1>);
112 std::cout << sizeof(concepts::Point<concepts::Real,2>);
113 std::cout << sizeof(concepts::Point<concepts::Cmplx,2>);
114 std::cout << sizeof(concepts::UnitNd<2>);
115 std::cout << sizeof(concepts::Point<concepts::Real,3>);
116 std::cout << sizeof(concepts::UnitNd<3>);
117 std::cout << sizeof(concepts::Level<1>);
118 std::cout << sizeof(concepts::Level<2>);
119 std::cout << sizeof(concepts::Level<3>);
120 std::cout << sizeof(concepts::Mapping<concepts::Real,2>);
121 concepts::Point<concepts::Real,2> p2(1.,2.);
122 concepts::Point<concepts::Real,2> p2_copy(p2);
123 std::cout << p2_copy;
124
125 concepts::Mapping<concepts::Real,2> map2(1.,2.,3.,4.);
126 std::cout << map2;
127
128 std::cout << concepts::isParallelRunning();
129
130 //toolbox
131 std::cout << sizeof(concepts::Array<unsigned int>);
132 std::cout << sizeof(concepts::Array<int>);
133 concepts::Array<unsigned int> arr(2);
134 std::cout <<arr<<std::endl;
135
136 std::cout << sizeof(concepts::Array<concepts::IndexRange>);
137 std::cout << sizeof(concepts::Array<concepts::Attribute>);
138 std::cout << sizeof(concepts::Array<concepts::Real>);
139
140 std::cout << sizeof(concepts::BaseSet<unsigned int>);
141 std::cout << sizeof(concepts::Set<unsigned int>);
142
143 std::cout << sizeof(concepts::BaseSet< concepts::IndexRange > );
144 std::cout << sizeof(concepts::Set< concepts::IndexRange > );
145
146 std::cout << sizeof(concepts::BaseSet< concepts::Attribute > );
147 std::cout << sizeof(concepts::Set<concepts::Attribute>);
148
149 std::cout << sizeof(concepts::Set<concepts::Real>);
150
151 std::cout << sizeof(concepts::Scan< concepts::Cell> );
152 std::cout << sizeof(concepts::Scan< concepts::Cell1> );
153 std::cout << sizeof(concepts::Scan< concepts::Cell2> );
154 std::cout << sizeof(concepts::Scan< concepts::Cell3> );
155 std::cout << sizeof(concepts::Scan< concepts::Element<concepts::Real> > );
156 std::cout << sizeof(concepts::Scan< concepts::ElementWithCell<concepts::Real> > );
157 std::cout << sizeof(concepts::Scan< hp1D::BaseElement<concepts::Real> > );
158
159
160 std::cout << sizeof(concepts::Sequence<concepts::Edge2d*>);
161 std::cout << sizeof(concepts::Sequence<concepts::MappingEdge2d*>);
162 std::cout << sizeof(concepts::Sequence<concepts::Cell1*>);
163 std::cout << sizeof(concepts::Sequence<concepts::Set<unsigned int> >);
164 std::cout << sizeof(concepts::Sequence<concepts::ElementWithCell<concepts::Real> *>);
165 std::cout << sizeof(concepts::Sequence<concepts::Attribute>);
166
167
168 std::cout << sizeof(concepts::ElementGraphics<concepts::Real>);
169 std::cout << sizeof(concepts::ElementGraphics<concepts::Cmplx>);
170
171 //integration
172 std::cout << sizeof(concepts::ShapeFunction1D<concepts::Real>);
173
174 //formula
175 std::cout << sizeof(concepts::Formula<concepts::Real>);
176 std::cout << sizeof(concepts::Formula<concepts::Cmplx>);
177 std::cout << sizeof(concepts::Formula<concepts::Real2d>);
178 std::cout << sizeof(concepts::Formula<concepts::Cmplx2d>);
179 std::cout << sizeof(concepts::ElementFormula<concepts::Real,concepts::Real>);
180 std::cout << sizeof(concepts::ElementFormula<concepts::Cmplx,concepts::Real>);
181 std::cout << sizeof(concepts::ElementFormula<concepts::Point<concepts::Real,2>,concepts::Real>);
182 std::cout << sizeof(concepts::ElementFormula<concepts::Point<concepts::Cmplx, 2>,concepts::Real>);
183 std::cout << sizeof(concepts::ElementFormula<concepts::Mapping<concepts::Real,2>,concepts::Real>);
184
185 std::cout << sizeof(concepts::ElementFormulaContainer<concepts::Real,concepts::Real>);
186 std::cout << sizeof(concepts::ElementFormulaContainer<concepts::Cmplx,concepts::Real>);
187 std::cout << sizeof(concepts::ElementFormulaContainer<concepts::Point<concepts::Real,2>,concepts::Real>);
188 std::cout << sizeof(concepts::ElementFormulaContainer<concepts::Point<concepts::Cmplx,2>,concepts::Real>);
189 std::cout << sizeof(concepts::ElementFormulaContainer<concepts::Mapping<concepts::Real,2>,concepts::Real>);
190 //concepts::ElementFormulaContainer<concepts::Mapping<concepts::Real,2>,concepts::Real> f();
191
192 // concepts::ElementFormulaContainer<concepts::Real,concepts::Real> fr1_r_r;
193 // concepts::ElementFormulaContainer<cfoncepts::Cmplx,concepts::Real> fr1_c_r;
194 // concepts::ElementFormulaContainer<concepts::Real,concepts::Real> fr2_r_r;
195 // concepts::ElementFormulaContainer<concepts::Cmplx,concepts::Real> fr2_c_r;
196
197 std::cout << sizeof(concepts::FrmE_Sum<concepts::Real, concepts::Real, concepts::Real, concepts::Real>);
198
199 std::cout << sizeof(concepts::FrmE_Product<concepts::Real, concepts::Real, concepts::Real, concepts::Real>);
200 std::cout << sizeof(concepts::FrmE_Product<concepts::Cmplx, concepts::Cmplx, concepts::Cmplx, concepts::Real>);
201 std::cout << sizeof(concepts::FrmE_Product<concepts::Point<double, 2>, concepts::Point<double, 2>, concepts::Real, concepts::Real>);
202 std::cout << sizeof(concepts::FrmE_Product<concepts::Point<concepts::Cmplx, 2>, concepts::Point<concepts::Cmplx, 2>, concepts::Cmplx, concepts::Real>);
203
204 std::cout << sizeof(concepts::PiecewiseFormulaBase<concepts::Real>);
205 std::cout << sizeof(concepts::PiecewiseFormulaBase<concepts::Cmplx>);
206 std::cout << sizeof(concepts::PiecewiseConstFormula<concepts::Real>);
207 std::cout << sizeof(concepts::PiecewiseConstFormula<concepts::Cmplx>);
208 std::cout << sizeof(concepts::PiecewiseConstFormula<concepts::Real2d>);
209 std::cout << sizeof(concepts::PiecewiseConstFormula<concepts::Cmplx2d>);
210 std::cout << sizeof(concepts::ParsedFormulaBase< 1u >);
211 std::cout << sizeof(concepts::ParsedFormulaBase< 2u >);
212
213 std::cout << sizeof(concepts::VectorElementFormulaBase<concepts::Real,2,concepts::Real>);
214 std::cout << sizeof(concepts::VectorElementFormulaBase<concepts::Cmplx,2,concepts::Real>);
215 std::cout << sizeof(concepts::VectorElementFormulaBase<concepts::Real,3,concepts::Real>);
216 std::cout << sizeof(concepts::VectorElementFormulaBase<concepts::Cmplx,3,concepts::Real>);
217 std::cout << sizeof(concepts::VectorElementFormula<concepts::Real,2,concepts::Real>);
218 std::cout << sizeof(concepts::VectorElementFormula<concepts::Cmplx,2,concepts::Real>);
219 std::cout << sizeof(concepts::VectorElementFormula<concepts::Real,3,concepts::Real>);
220 std::cout << sizeof(concepts::VectorElementFormula<concepts::Cmplx,3,concepts::Real>);
221
222 //geometry
223 std::cout << sizeof(concepts::Semantics<concepts::Boundary>);
224 std::cout << sizeof(concepts::Semantics<concepts::CellCondition>);
225 std::cout << sizeof(concepts::PointInCell<2>);
226 std::cout << sizeof(concepts::CellToCellMapping<2,2>);
227 std::cout << sizeof(concepts::CellDiameter);
228
229 //space
230 std::cout << sizeof(concepts::Space<concepts::Real>);
231 std::cout << sizeof(concepts::Subspace);
232 std::cout << sizeof(concepts::DummySpace<concepts::Real>);
233 std::cout << sizeof(concepts::SpaceOnCells<concepts::Real>);
234 std::cout << sizeof(concepts::SpaceOnCoarseCells<2,concepts::Real>);
235 std::cout << sizeof(concepts::Element<concepts::Real>);
236 std::cout << sizeof(concepts::Element<concepts::Cmplx>);
237 std::cout << sizeof(concepts::ElementWithCell<concepts::Real>);
238 std::cout << sizeof(concepts::ElementWithCell<concepts::Real> *);
239 std::cout << sizeof(concepts::ElementMatrix<concepts::Real>);
240 std::cout << sizeof(concepts::ElementMatrix<concepts::Cmplx>);
241 std::cout << sizeof(concepts::ElementPair<concepts::Real>);
242 std::cout << sizeof(concepts::ElementFunction<concepts::Real,concepts::Real>);
243 std::cout << sizeof(concepts::ElementFunction<concepts::Cmplx,concepts::Real>);
244 std::cout << sizeof(concepts::ElementFunction<concepts::Real,concepts::Cmplx>);
245 std::cout << sizeof(concepts::ElementFunction<concepts::Cmplx,concepts::Cmplx>);
246 std::cout << sizeof(concepts::CellPostprocess<concepts::Real>);
247 std::cout << sizeof(concepts::CellPostprocess<concepts::Cmplx>);
248 std::cout << sizeof(concepts::TColumn<concepts::Real>);
249 std::cout << sizeof(concepts::TColumn<bool>);
250 std::cout << sizeof(concepts::TMatrixBase<concepts::Real>);
251 std::cout << sizeof(concepts::TMatrixBase<concepts::Cmplx>);
252 std::cout << sizeof(concepts::TMatrix<concepts::Real>);
253 std::cout << sizeof(concepts::Dirichlet<concepts::Real>);
254 std::cout << sizeof(concepts::Dirichlet<concepts::Cmplx>);
255 std::cout << sizeof(concepts::AdaptiveAdjust<2>);
256 std::cout << sizeof(concepts::AdaptiveAdjustP<2>);
257 std::cout << sizeof(concepts::AdaptiveSpace<concepts::Real,concepts::AdaptiveAdjustP<2> >);
258 std::cout << sizeof(concepts::SpaceGraph<hp2D::hpAdaptiveSpaceH1>);
259 // free functions
260 concepts::Square square_mesh(1,1,1);
261 concepts::BoundaryConditions bc;
262 hp2D::hpAdaptiveSpaceH1 hpSpace(square_mesh, (uint)0, (uint)1, &bc);
263
264 // concepts::SpaceGraph<hp2D::hpAdaptiveSpaceH1> sgraph(hpSpace);
265 // sgraph.FillGraph(2);
266 // sgraph.SplitGraph(1);
267
268 concepts::ConstFormula<concepts::Real> const_formula(1.0);
269 concepts::ConstFormula<concepts::Real2d > const_formula_point(1.0);
270
271 //concepts::Formula<concepts::Real>* formula_ref = dynamic_cast<concepts::Formula<concepts::Real>* >(&const_formula);
272
273 concepts::ElementFormula<concepts::Real, concepts::Real>* element_formula = dynamic_cast< concepts::ElementFormula< concepts::Real, concepts::Real > * >(&const_formula);
274
275 concepts::ElementFormula<concepts::Real2d, concepts::Real>* element_formula_point = dynamic_cast< concepts::ElementFormula< concepts::Real2d, concepts::Real > * >(&const_formula_point);
276
277 //concepts::L2product(hpSpace, *(formula_ref));
278 concepts::L2product(hpSpace, *(element_formula));
279 concepts::L2product(hpSpace, *(element_formula_point));
280
281
282
283 // concepts::L2product(sgraph.getSubDomain(0), *(element_formula));
284 // concepts::L2product(sgraph.getSubDomain(0), *(element_formula_point));
285
286 // free operators
287 concepts::ElementFormulaContainer<concepts::Real, concepts::Real> el_form_cont(1.0);
288 el_form_cont + el_form_cont;
289 el_form_cont - el_form_cont;
290 concepts::operator+(el_form_cont, el_form_cont);
291 concepts::operator-(el_form_cont, el_form_cont);
292 //function
293 std::cout << sizeof(concepts::Function<concepts::Real>);
294 std::cout << sizeof(concepts::Function<std::complex<concepts::Real> >);
295 std::cout << sizeof(concepts::LinearForm<concepts::Real,concepts::Real>);
296 std::cout << sizeof(concepts::LinearForm<concepts::Cmplx,concepts::Real>);
297 std::cout << sizeof(concepts::Vector<concepts::Real>);
298 std::cout << sizeof(concepts::Vector<std::complex<concepts::Real> >);
299 std::cout << sizeof(concepts::Neumann<concepts::Real>);
300
301 std::cout << sizeof(concepts::Space<concepts::Real>);
302 std::cout << sizeof(concepts::Space<concepts::Cmplx>);
303 std::cout << sizeof(concepts::ElementFormulaVector<1,concepts::Real,concepts::Real,concepts::Real>);
304 std::cout << sizeof(concepts::ElementFormulaVector<2,concepts::Real,concepts::Real,concepts::Real>);
305 std::cout << sizeof(concepts::ElementFormulaVectorBase<concepts::Real,concepts::Real,concepts::Real,concepts::Real>);
306 std::cout << sizeof(concepts::ElementFormulaVectorBase<concepts::Real,
307 concepts::Real,
308 concepts::Point<concepts::Real,2>,
309 concepts::Real>);
310 //operator
311
312 std::cout << sizeof(concepts::Matrix<concepts::Real>);
313 std::cout << sizeof(concepts::Matrix<std::complex<concepts::Real> >);
314
315 std::cout << sizeof(concepts::SparseMatrix<concepts::Real>);
316 std::cout << sizeof(concepts::SparseMatrix<concepts::Cmplx>);
317
318 concepts::SparseMatrix<concepts::Real> SM_Real(3,3);
319 concepts::_HashedSMatrix_iterator<concepts::Real, const concepts::Real &, const concepts::Real *> SM_Real_cit = SM_Real.begin();
320 concepts::_HashedSMatrix_iterator<concepts::Real, concepts::Real &, concepts::Real *> SM_Real_it = SM_Real.begin();
321 std::cout << sizeof(SM_Real_cit) << sizeof(SM_Real_it) ;
322
323 concepts::SparseMatrix<concepts::Cmplx> SM_Cmplx(3,3);
324 concepts::_HashedSMatrix_iterator<concepts::Cmplx, const concepts::Cmplx &, const concepts::Cmplx *> SM_Cmplx_cit = SM_Cmplx.begin();
325 concepts::_HashedSMatrix_iterator<concepts::Cmplx, concepts::Cmplx &, concepts::Cmplx *> SM_Cmplx_it = SM_Cmplx.begin();
326 std::cout << sizeof(SM_Cmplx_cit) << sizeof(SM_Cmplx_it) ;
327
328 std::cout << sizeof(concepts::HashedSparseMatrix<concepts::Real>);
329 std::cout << sizeof(concepts::Operator<concepts::Real>);
330 std::cout << sizeof(concepts::VecOperator<concepts::Real>);
331 std::cout << sizeof(concepts::CG<concepts::Real>);
332 std::cout << sizeof(concepts::CG<concepts::Cmplx>);
333
334#ifdef HAS_MUMPS
335 std::cout << sizeof(concepts::Mumps<concepts::Real>);
336 std::cout << sizeof(concepts::Mumps<concepts::Cmplx>);
337#endif // HAS_MUMPS
338 std::cout << sizeof(concepts::DenseMatrix<concepts::Real>);
339 std::cout << sizeof(concepts::DenseMatrix<concepts::Cmplx>);
340
341 std::cout << sizeof(concepts::BilinearForm<concepts::Real,concepts::Real>);
342 // std::cout << sizeof(concepts::BilinearForm<concepts::Cmplx,concepts::Cmplx>);
343 // std::cout << sizeof(concepts::BilinearForm<concepts::Real,concepts::Cmplx>);
344 std::cout << sizeof(concepts::BilinearForm<concepts::Cmplx,concepts::Real>);
345
346 std::cout << sizeof(concepts::BilinearFormLiCo<concepts::Real,concepts::Real>);
347 std::cout << sizeof(concepts::BilinearFormLiCo<concepts::Cmplx,concepts::Real>);
348
349 std::cout << sizeof(concepts::BilinearFormContainer<concepts::Real,concepts::Real>);
350 std::cout << sizeof(concepts::BilinearFormContainer<concepts::Cmplx,concepts::Real>);
351
352 std::cout << sizeof(concepts::BilinearF_Sum<concepts::Real,concepts::Real,concepts::Real,concepts::Real>);
353 std::cout << sizeof(concepts::BilinearF_Sum<concepts::Cmplx,concepts::Cmplx,concepts::Cmplx,concepts::Real>);
354
355 std::cout << sizeof(concepts::BilinearF_W<concepts::Real,concepts::Real,concepts::Real,concepts::Real>);
356 std::cout << sizeof(concepts::BilinearF_W<concepts::Cmplx,concepts::Cmplx,concepts::Cmplx,concepts::Real>);
357
358 std::cout << sizeof(concepts::SolverFabric<concepts::Real>);
359
360#ifdef HAS_SuperLU
361 std::cout << sizeof(concepts::SuperLU<concepts::Real>);
362 std::cout << sizeof(concepts::SuperLU<concepts::Cmplx>);
363 std::cout << sizeof(concepts::SuperLUFabric<concepts::Real>);
364#endif // HAS_SuperLU
365
366 std::cout << sizeof(concepts::CGFabric<concepts::Real>);
367 std::cout << sizeof(concepts::DiagonalSolverFabric<concepts::Real>);
368 std::cout << sizeof(concepts::GMResFabric<concepts::Real>);
369
370 std::cout << sizeof(concepts::Estimator<concepts::Real>);
371 std::cout << sizeof(concepts::LocalEstimator<concepts::Real>);
372 std::cout << sizeof(concepts::ExplicitResidual<concepts::Real>);
373
374#endif // GEN_PYTHON
375}
376
377#endif
concepts::BilinearF_Sum
Definition bilinearForm.hh:211
concepts::BilinearF_W
Definition bilinearForm.hh:275
concepts::FrmE_Product
Definition frmE_product.hh:39
concepts::FrmE_Sum
Definition frmE_sum.hh:36
concepts::_HashedSMatrix_iterator
Definition hashedSMatrixIterator.hh:29
frmE_normalVector.hh
frmE_product.hh
frmE_sum.hh
concepts::operator-
ElementFormulaContainer< Real > operator-(const ElementFormulaContainer< Real > frm, const Real a)
Simple subtracting of a element formulas and a constant via –operator.
concepts::operator+
Frm_Sum< Real > operator+(const Formula< Real > &frm, const Real a)
concepts::Real
double Real
Definition typedefs.hh:39
concepts::isParallelRunning
bool isParallelRunning()
concepts::L2product
Real L2product(const ElementWithCell< G > &elm, const ElementFormula< F, G > &u, const ElementFormula< Real > *c=0, const Real t=0.0, IntegrationCell::intFormType form=IntegrationCell::ZERO)
Definition integral.hh:214
concepts::Cmplx
std::complex< Real > Cmplx
Type for a complex number. It also depends on the setting of Real.
Definition typedefs.hh:42
Generated on Wed Sep 13 2023 21:06:23 for Concepts by doxygen 1.9.8