00001 #include <SyFi.h> 00002 00003 using namespace GiNaC; 00004 00005 int main () { 00006 int order = 3; 00007 00008 ReferenceTriangle triangle("t"); 00009 RaviartThomas vfe; 00010 vfe.set(triangle); 00011 vfe.set(order); 00012 vfe.compute_basis_functions(); 00013 00014 DiscontinuousLagrangeFE pfe; 00015 pfe.set(triangle); 00016 pfe.set(order); 00017 pfe.compute_basis_functions(); 00018 00019 for (int i=0; i< vfe.nbf(); i++) 00020 cout <<"vfe.N("<<i<<")="<<vfe.N(i)<<endl; 00021 00022 for (int i=0; i< pfe.nbf(); i++) 00023 cout <<"pfe.N("<<i<<")="<<pfe.N(i)<<endl; 00024 00025 return 0; 00026 } 00027 00028 00029 00030 00031 00032 00033