00001 #include <FE.h> 00002 00003 int StandardFE:: nbf() { 00004 return Ns.size(); 00005 } 00006 00007 void StandardFE:: compute_basis_functions() { 00008 // FIXME: give some valuable tips, can for instance 00009 // say that "standard Lagrange elements are implemented in LagrangeFE 00010 cout <<"StandardFE compute_basis_functions not implemented"<<endl; 00011 } 00012 00013 Polygon& StandardFE:: getPolygon() { 00014 // FIXME this is not nice. 00015 // should check that p is ok!! 00016 return *p; 00017 } 00018 00019 void StandardFE:: set(Polygon& p_) { 00020 Ns.clear(); 00021 dofs.clear(); 00022 p = &p_; 00023 } 00024 00025 00026 void StandardFE:: set(int order_) { 00027 Ns.clear(); 00028 dofs.clear(); 00029 order = order_; 00030 } 00031 00032 00033 GiNaC::ex StandardFE:: dof(int i) { 00034 //FIXME could be point or dirac on point 00035 return dofs[i-1]; 00036 } 00037 00038 00039 GiNaC::ex StandardFE::N(int i) { 00040 return Ns[i-1]; 00041 } 00042 00043 00044