trunk/ascxx/idatest.cpp

#include "simulation.h"
#include "library.h"
#include "solver.h"
#include "solverreporter.h"
#include "integrator.h"
#include "solverparameters.h"
#include "integratorreporter.h"

#include <iostream>
using namespace std;

int main(void){

    Library L;
    L.load("johnpye/idadenx.a4c");

    Type T = L.findType("idadenx");
    Simulation M = T.getSimulation("sim");

//  M.setSolver(Solver("QRSlv"));
    
    M.build();

    Integrator I(M);

    I.setEngine("IDA");
    SolverParameters P = I.getParameters();
    P.getParameter("safeeval").setBoolValue(false);
    P.getParameter("linsolver").setStrValue("DENSE");
    I.setParameters(P);

    IntegratorReporterConsole R(&I);
    I.setReporter(&R);

    I.setLogTimesteps(UnitsM("s"), 0.4, 4e10, 11);

    I.analyse();
    I.solve();
}
1	#include "simulation.h"
2	#include "library.h"
3	#include "solver.h"
4	#include "solverreporter.h"
5	#include "integrator.h"
6	#include "solverparameters.h"
7	#include "integratorreporter.h"
8
9	#include <iostream>
10	using namespace std;
11
12	int main(void){
13
14	Library L;
15	L.load("johnpye/idadenx.a4c");
16
17	Type T = L.findType("idadenx");
18	Simulation M = T.getSimulation("sim");
19
20	// M.setSolver(Solver("QRSlv"));
21
22	M.build();
23
24	Integrator I(M);
25
26	I.setEngine("IDA");
27	SolverParameters P = I.getParameters();
28	P.getParameter("safeeval").setBoolValue(false);
29	P.getParameter("linsolver").setStrValue("DENSE");
30	I.setParameters(P);
31
32	IntegratorReporterConsole R(&I);
33	I.setReporter(&R);
34
35	I.setLogTimesteps(UnitsM("s"), 0.4, 4e10, 11);
36
37	I.analyse();
38	I.solve();
39	}