tetrahedron_ex3.cpp

Go to the documentation of this file.

#include <SyFi.h>

using namespace GiNaC; 
using namespace SyFi; 
using namespace std; 

int main() {


    initSyFi(3); 

    //example that check scaling

    numeric h(1,100); //1.0/100 
    numeric a(1,2); 
    ex p0 = lst(a,a,a);
    ex p1 = lst(a+h,a,a);
    ex p2 = lst(a,a+h,a);
    ex p3 = lst(a,a,a+h);

    Tetrahedron tetrahedron(p0,p1,p2,p3);


    ex f = 1;
    ex intf = tetrahedron.integrate(f); 
    cout <<"intf "<<intf<<endl; 
    EQUAL_OR_DIE(intf, "1/6000000"); 

    return 0; 
}