PetriMetaSens: dyn.py

import sys

import diana
try:
    dmain=diana.GetDianaMain()
    # get Diana model manger
    mmanger=dmain.GetModelManager()
    # create model
    model=mmanger.CreateModel(diana.CAPE_CONTINUOUS, "Example.so")
    model.Initialize()
    eso  = model.GetActiveESO()
    epar = eso.GetParameters()
    evar = eso.GetStateVariables()

    # create sensitivity model based on base model
    #sensmodel = model.GetSensitivityModel(['p1', 'p2', 'p3'], [])
    sensmodel = model.GetSensitivityModel(['p1'], [])
    sensmodel.Initialize()

    # create solver
    sfactory = dmain.GetSolverFactory()
    solver = sfactory.CreateSolver(diana.CAPE_DAE, sensmodel, "petri_meta_sens")
    solver.Initialize()
    spar = solver.GetParameters()
    spar['VerboseLevel'].SetValue(6)
    spar['Tend'].SetValue(40)

    # create report
    ri=dmain.CreateReportingInterface("basic")
    ri.SetComponentName("data")
    solver.SetReportingInterface(ri)
    ri.Add(spar['T'])
    ri.Add(evar['x1'])
    ri.Add(evar['x2'])
    ri.Add(epar['state_l'])
    ri.Add(epar['from_l_to_v'])
    ri.Add(epar['from_v_to_l'])

    solver.Solve()
    spar['T'].GetValue()

    z = sensmodel.GetSensESO()
    print "=" * 20
    print z.GetAllSensitivities()
    
    # store results into directory
    import os
    if not os.path.isdir('results'): os.mkdir('results');
    ri.WriteData("results/dyn")
    ri.WriteDataMatlab("results/dyn.m")

except diana.ECapeUser, exc:
    print '! Exception:', exc