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 nsd = 3; 00011 ReferenceTetrahedron tetrahedon; 00012 Nedelec2Hdiv fe; 00013 fe.set(tetrahedon); 00014 fe.set(1); 00015 fe.compute_basis_functions(); 00016 cout <<"-------- 1 order ---------- 3D "<<endl; 00017 print_out(fe); 00018 00019 fe.set(2); 00020 fe.compute_basis_functions(); 00021 cout <<"-------- 2 order ---------- 3D "<<endl; 00022 print_out(fe); 00023 00024 fe.set(3); 00025 fe.compute_basis_functions(); 00026 cout <<"-------- 3 order ---------- 3D "<<endl; 00027 print_out(fe); 00028 00029 fe.set(4); 00030 fe.compute_basis_functions(); 00031 cout <<"-------- 4 order ---------- 3D "<<endl; 00032 print_out(fe); 00033 00034 00035 00036 00037 } 00038 00039 00040