nedelec_ex1.cpp

Go to the documentation of this file.
00001 #include <SyFi.h>
00002 
00003 void print_out(FE& fe) { 
00004   for (int i=0; i< fe.nbf(); i++) { 
00005     cout <<"fe.N("<<i<<")        =       "<<fe.N(i)<<endl; 
00006   }
00007 }
00008 
00009 int main() {
00010   ReferenceTriangle triangle; 
00011   Nedelec fe; 
00012   fe.set(triangle); 
00013   fe.set(0); 
00014   fe.compute_basis_functions(); 
00015   cout <<"-------- 0 order ---------- 2D "<<endl; 
00016   print_out(fe); 
00017 
00018   fe.set(1); 
00019   fe.compute_basis_functions(); 
00020   cout <<"-------- 1 order ---------- 2D"<<endl; 
00021   print_out(fe); 
00022 
00023   fe.set(2); 
00024   fe.compute_basis_functions(); 
00025   cout <<"-------- 2 order ---------- 2D"<<endl; 
00026   print_out(fe); 
00027 
00028   fe.set(3); 
00029   fe.compute_basis_functions(); 
00030   cout <<"-------- 2 order ---------- 2D"<<endl; 
00031   print_out(fe); 
00032 
00033 
00034 
00035 
00036 
00037   nsd = 3; 
00038   ReferenceTetrahedron tetrahedon; 
00039   fe.set(tetrahedon); 
00040   fe.set(0); 
00041   fe.compute_basis_functions(); 
00042   cout <<"-------- 0 order ---------- 3D "<<endl; 
00043   print_out(fe); 
00044 
00045   fe.set(1); 
00046   fe.compute_basis_functions(); 
00047   cout <<"-------- 1 order ---------- 3D "<<endl; 
00048   print_out(fe); 
00049 
00050   fe.set(2); 
00051   fe.compute_basis_functions(); 
00052   cout <<"-------- 2 order ---------- 3D "<<endl; 
00053   print_out(fe); 
00054 
00055   fe.set(3); 
00056   fe.compute_basis_functions(); 
00057   cout <<"-------- 2 order ---------- 3D "<<endl; 
00058   print_out(fe); 
00059 
00060 
00061 
00062 
00063 }
00064 
00065 

Generated on Tue Jun 13 13:18:39 2006 for SyFi by  doxygen 1.4.4