GiNaCDE/test/test1.cpp at master · mithun218/GiNaCDE · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79

/** @file test1.cpp
 *
 *   Program to test GiNaCDE library. This solve following odes
        i.   linear ode - the damped harmonic oscillator.
        ii.  nonlinear ode - modified Painlev-Ince equation */


#include "GiNaCDE.h"

 int main()
 {
    const ex u=reader("u"),x=reader("x"),a=reader("a"),b=reader("b"),
             c=reader("c"),w=reader("w"),A_0=reader("A_0"),A_1=reader("A_1"),A_2=reader("A_2");
    ex ode,res;
    size_t solu_num1,solu_num2;
    stringstream diffStr,algSoluStr,diffSoluStr;
    string str;

    depend(u, {x});

    twcPhase = lst{lst{},lst{}};

    ode = Diff(u,x,2) + Diff(u,x,1) + w*w*u; // the damped harmonic oscillator
    output = ginac;// Outputs are saved in ginac format;
    filename = "damped_FIM.txt";
    desolve(ode,{u},FIM,true);
    /* Checking all solutions*/
    diffStr.str("");
    diffStr<<ode;
    solu_num1 = (solutionClt).size();
    for(size_t i=0;i<solu_num1;i++)
    {
        algSoluStr.str("");
        algSoluStr<<solutionClt[i][0];
        solu_num2 = nops(solutionClt[i]);
        for(size_t j=1;j<solu_num2;j++)
        {
            diffSoluStr.str("");
            diffSoluStr<<solutionClt[i][j];
            res = checkSolu(diffStr.str(),diffSoluStr.str(),algSoluStr.str());
            if(res!=_ex0)
                return -1;
        }

    }


    ode = Diff(u,x,2) + a*u*Diff(u,x,1) + b*u*u*u; // modified Painlev-Ince equation
    output = ginac;
    filename = "Painlev_FIMextravar.txt";
    paraInDiffSolve = lst{a,b};
    desolve(ode,{u},FIM,true);
    /* Checking all solutions*/
    diffStr.str("");
    diffStr<<ode;
    solu_num1 = (solutionClt).size();
    for(size_t i=0;i<solu_num1;i++)
    {
        algSoluStr.str("");
        algSoluStr<<solutionClt[i][0];
        solu_num2 = nops(solutionClt[i]);
        for(size_t j=1;j<solu_num2;j++)
        {
            diffSoluStr.str("");
            diffSoluStr<<solutionClt[i][j];
            res = checkSolu(diffStr.str(),diffSoluStr.str(),algSoluStr.str());
            if(res!=_ex0)
                return -1;
        }

    }

    return 0;

 }