Rheolef  7.1
an efficient C++ finite element environment
bubble_symbolic.cc
Go to the documentation of this file.
1 //
22 // bubble approximation
23 //
24 #include "basis_symbolic.h"
25 using namespace rheolef;
26 using namespace std;
27 using namespace GiNaC;
28 
29 class bubble_symbolic : public basis_symbolic_nodal
30 {
31 public:
32  bubble_symbolic ();
33 };
34 bubble_symbolic::bubble_symbolic ()
35 : basis_symbolic_nodal("bubble",4)
36 {
37  on('p') << node(0)
38  << poly (1)
39  << end;
40  on('e') << node (0.5)
41  << poly (4*x*(1-x))
42  << end;
43  on('t') << node (ex(1)/3, ex(1)/3)
44  << poly (27*x*y*(1-x-y))
45  << end;
46  on('T') << node (1./4, 1./4, 1./4)
47  << poly (256*x*y*z*(1-x-y-z))
48  << end;
49 }
50 int main (int argc, char **argv) {
51  bubble_symbolic bubble;
52  bubble.put_cxx_main (argc,argv);
53 }
int main(int argc, char **argv)
This file is part of Rheolef.