00001 #include <SyFi.h> 00002 00003 void usage(FE& fe) { 00004 ex Ni; 00005 ex gradNi; 00006 ex dofi; 00007 for (int i=1; i<= fe.nbf() ; i++) { 00008 Ni = fe.N(i); 00009 gradNi = grad(Ni); 00010 dofi = fe.dof(i); 00011 cout <<"The basis function, N("<<i<<")="<<Ni<<endl; 00012 cout <<"The gradient of N("<<i<<")="<<gradNi<<endl; 00013 cout <<"The corresponding dof, L("<<i<<")="<<dofi<<endl; 00014 } 00015 } 00016 00017 00018 00019 int main(){ 00020 int order = 2; 00021 Triangle triangle(lst(0,0), lst(1,0), lst(0,1)); 00022 LagrangeFE fe; 00023 fe.set(order); 00024 fe.set(triangle); 00025 fe.compute_basis_functions(); 00026 00027 usage(fe); 00028 00029 } 00030 00031