P1_symbolic.cc

Go to the documentation of this file.

//
// P1 approximation
//
#include "basis_symbolic.h"
using namespace rheolef;
using namespace std;
 
class P1_symbolic : public basis_symbolic_nodal
{
public:
    P1_symbolic ();
};
P1_symbolic::P1_symbolic ()
: basis_symbolic_nodal("P1",1)
{
  basis_symbolic_nodal::_family_name = "P";
  basis_symbolic_nodal::set_degree_parameter();
  basis_symbolic_nodal::set_continuous_feature();
  on('p') << node(0) << poly (1) << end;
  on('e') << node (0)
          << node (1)
          << poly (1)
          << poly (x)
          << end;
  on('t') << node (0, 0)
          << node (1, 0)
          << node (0, 1)
          << poly (1)
          << poly (x)
          << poly (y)
          << end;
  on('q') << node (-1, -1)
          << node ( 1, -1)
          << node ( 1,  1)
          << node (-1,  1)
          << poly (1)
          << poly (x)
          << poly (y)
          << poly (x*y)
          << end;
  on('T') << node (0, 0, 0)
          << node (1, 0, 0)
          << node (0, 1, 0)
          << node (0, 0, 1)
          << poly (1)
          << poly (x)
          << poly (y)
          << poly (z)
          << end;
  on('P') << node(0, 0,-1)
          << node(1, 0,-1)
          << node(0, 1,-1)
          << node(0, 0, 1)
          << node(1, 0, 1)
          << node(0, 1, 1)
          << poly (1)
          << poly (x)
          << poly (y)
          << poly (z)
          << poly (x*z)
          << poly (y*z)
          << end;
  on('H') << node(-1,-1,-1)
          << node( 1,-1,-1)
          << node( 1, 1,-1)
          << node(-1, 1,-1)
          << node(-1,-1, 1)
          << node( 1,-1, 1)
          << node( 1, 1, 1)
          << node(-1, 1, 1)
          << poly (1)
          << poly (x)
          << poly (y)
          << poly (z)
          << poly (x*y)
          << poly (y*z)
          << poly (z*x)
          << poly (x*y*z)
          << end;
}
int main (int argc, char **argv) {
        P1_symbolic P1;
        P1.put_cxx_main (argc,argv);
}