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