Rheolef  7.2
an efficient C++ finite element environment
P1qd_symbolic.cc
Go to the documentation of this file.
1 //
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"
27 using namespace rheolef;
28 using namespace std;
29 using namespace GiNaC;
30 
31 class P1qd_symbolic : public basis_symbolic_nodal
32 {
33 public:
34  P1qd_symbolic ();
35 };
36 P1qd_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 }
56 int main (int argc, char **argv) {
57  P1qd_symbolic P1qd;
58  P1qd.put_cxx_main (argc,argv);
59 }
int main(int argc, char **argv)
This file is part of Rheolef.