polh.cpp

Go to the documentation of this file.

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

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

int main() {

    archive ar; 
    int order = 3; 

    initSyFi(1);
    cout <<"third order homogenous polynomial in 1D"<<endl; 
    ex polh = homogenous_pol(order, nsd, "a");  
    cout <<"polh "<<polh<<endl; 
    ar.archive_ex(polh, "polh1_3"); 


    initSyFi(2);
    cout <<"third order homogenous polynomial in 2D"<<endl; 
    polh = homogenous_pol(order, nsd, "a");  
    cout <<"polh "<<polh<<endl; 
    ar.archive_ex(polh, "polh2_3"); 

    initSyFi(3);
    cout <<"third order homogenous polynomial in 3D"<<endl; 
    polh = homogenous_pol(order, nsd, "a");  
    cout <<"polh "<<polh<<endl; 
    ar.archive_ex(polh, "polh3_3"); 

    ofstream vfile("polh.gar.v"); 
    vfile << ar; vfile.close(); 
    if(!compare_archives("polh.gar.v", "polh.gar.r")) { 
            cerr << "Failure!" << endl;
            return -1;
    }

    return 0; 


}