Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
P1qd_symbolic.cc
Go to the documentation of this file.
1
21//
22// P1qd approximation: linear polynoms in quadrangles or hexahedra
23// used for Stokes : incompressible velocity-pressure element P2-P1q
24// P1qd is globally discontinuous (see P1qd_numbering)
25// see Brezzi-Fortin 1991, page 266.
26#include "basis_symbolic.h"
27using namespace rheolef;
28using namespace std;
29using namespace GiNaC;
30
31class P1qd_symbolic : public basis_symbolic_nodal
32{
33public:
34 P1qd_symbolic ();
35};
36P1qd_symbolic::P1qd_symbolic ()
37: basis_symbolic_nodal("P1qd",1)
38{
39 on('q') << node ( 0.5, 0.5)
40 << node ( 0.5, -0.5)
41 << node (-0.5, -0.5)
42 << poly (1)
43 << poly (x)
44 << poly (y)
45 << end;
46 on('H') << node ( 0.5, 0.5, 0.5)
47 << node ( 0.5, 0.5, -0.5)
48 << node ( 0.5, -0.5, -0.5)
49 << node (-0.5, -0.5, -0.5)
50 << poly (1)
51 << poly (x)
52 << poly (y)
53 << poly (z)
54 << end;
55}
56int main (int argc, char **argv) {
57 P1qd_symbolic P1qd;
58 P1qd.put_cxx_main (argc,argv);
59}
int main()
Definition field2bb.cc:58
This file is part of Rheolef.
STL namespace.