raviartthomas_ex2.cpp

Go to the documentation of this file.
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 

Generated on Tue Apr 18 12:00:26 2006 for SyFi by  doxygen 1.4.4