Main Page | Namespace List | Class Hierarchy | Class List | Directories | File List | Namespace Members | Class Members | File Members | Search for

# tetrahedron_ex1.cpp

Go to the documentation of this file.

    00001 #include <SyFi.h>
    00002
    00003 using namespace GiNaC;
    00004
    00005 int main() {
    00006
    00007   ex p0 = lst(0.0,0.0,0.0);
    00008   ex p1 = lst(1.0,0.0,0.0);
    00009   ex p2 = lst(0.0,1.0,0.0);
    00010   ex p3 = lst(0.0,0.0,1.0);
    00011
    00012   Tetrahedron tetrahedron(p0,p1,p2,p3);
    00013
    00014   ex repr = tetrahedron.repr();
    00015   cout <<"t.repr "<<repr<<endl;
    00016   EQUAL_OR_DIE(repr, "{x==r,y==s,z==t,{r,0,1},{s,0,1-r},{t,0,1-s-r}}");
    00017
    00018   ex f = x*y*z;
    00019   ex intf = tetrahedron.integrate(f);
    00020   cout <<"intf "<<intf<<endl;
    00021   EQUAL_OR_DIE(intf, "1/720");
    00022
    00023   return 0;
    00024 }
    00025