00001 #include <P0.h> 00002 00003 void P0:: compute_basis_functions() { 00004 // insert basis function 00005 Ns.insert(Ns.end(), GiNaC::numeric(1)); 00006 00007 // create and insert dof 00008 GiNaC::ex midpoint; 00009 for (int i=1; i<= p->no_vertices(); i++) { 00010 midpoint += p->vertex(i-1); 00011 } 00012 midpoint /= p->no_vertices(); 00013 dofs.insert(dofs.end(), midpoint); 00014 }