tetrahedron_ex3.cpp

Go to the documentation of this file.
00001 #include <SyFi.h>
00002 
00003 using namespace GiNaC; 
00004 
00005 int main() {
00006   
00007   //example that check scaling
00008 
00009   numeric h(1,100); //1.0/100 
00010   numeric a(1,2); 
00011   ex p0 = lst(a,a,a);
00012   ex p1 = lst(a+h,a,a);
00013   ex p2 = lst(a,a+h,a);
00014   ex p3 = lst(a,a,a+h);
00015   
00016   Tetrahedron tetrahedron(p0,p1,p2,p3);
00017   
00018  
00019   ex f = 1;
00020   ex intf = tetrahedron.integrate(f); 
00021   cout <<"intf "<<intf<<endl; 
00022   EQUAL_OR_DIE(intf, "1/6000000"); 
00023  
00024   return 0; 
00025 }
00026 

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