DianaICSolver: Cstr.py

import diana 
try:
    dmain=diana.GetDianaMain()
    # get Diana model manger
    mmanger=dmain.GetModelManager();

    # create model
    model=mmanger.CreateModel(diana.CAPE_CONTINUOUS, "Cstr.so");
    model.Initialize();
    eso=model.GetActiveESO();
    esopar=eso.GetParameters();
    ri=dmain.CreateReportingInterface("basic");
    ri.SetComponentName("data");

    # create solver
    sfactory=dmain.GetSolverFactory();

    ## consistent initial conditions solver calling
    icsolver=sfactory.CreateSolver(diana.CAPE_IC, model, "icsolver.so");
    icsolver.Initialize();
    icpar=icsolver.GetParameters();
    icpar['VerboseLevel'].SetValue(0);
    icsolver.Solve();
    ##

    solver=sfactory.CreateSolver(diana.CAPE_DAE, model, "ida.so");
    solver.Initialize();
    solpar=solver.GetParameters();
    solver.SetReportingInterface(ri);

    ri.Add(solpar["T"]);
    ri.Add(eso.GetStateVariables().ItemByName("x1"));
    ri.Add(eso.GetStateVariables().ItemByName("x2"));

    solpar["VerboseLevel"].SetValue(1);
    solpar["Tend"].SetValue(100);
    solver.Solve();

    # create directory and store results
    import os
    if not os.path.isdir('results'): os.mkdir('results');
    ri.WriteData("results/Cstr");
    ri.WriteDataMatlab("results/Cstr.m");
except diana.ECapeUser, exc:
    print repr(exc)