barycenter_simplex.cpp

Go to the documentation of this file.

#include <SyFi.h>
#include <fstream>

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

int main() {

  /*
    initSyFi(3); 

    ex p0 = lst(0.0,0.0,0.0);
    ex p1 = lst(1.0,0.0,0.0);
    ex p2 = lst(0.0,1.0,0.0);
    ex p3 = lst(0.0,0.0,1.0);

    Simplex simplex(lst(p0,p1,p2,p3)); 

    ex coord = barycenter(simplex); 
    cout <<"coord "<<coord<<endl; 
    */

    initSyFi(4); 

    ex q0 = lst(0.0,0.0,0.0,0.0);
    ex q1 = lst(1.0,0.0,0.0,0.0);
    ex q2 = lst(0.0,1.0,0.0,0.0);
    ex q3 = lst(0.0,0.0,1.0,0.0);
    ex q4 = lst(0.0,0.0,0.0,1.0);

    Simplex simplex4(lst(q0,q1,q2,q3,q4)); 

    ex coord4 = barycenter(simplex4); 
    cout <<"coord "<<coord4<<endl; 

}