14 std::string vCellModel;
16 std::vector<StimulusSharedPtr> vStimulus;
24 vSession = LibUtilities::SessionReader::CreateInstance(argc, argv);
35 vSession->LoadSolverInfo(
"CELLMODEL", vCellModel,
"");
36 ASSERTL0(vCellModel !=
"",
"Cell Model not specified.");
39 vCellModel, vSession, vExp);
43 vStimulus = Stimulus::LoadStimuli(vSession, vExp);
48 vDeltaT = vSession->GetParameter(
"TimeStep");
50 nSteps = vSession->GetParameter(
"NumSteps");
53 vSession->GetFunction(
"InitialConditions",
"u");
54 vSol[0][0] = e->Evaluate(0.0, 0.0, 0.0, 0.0);
57 for (
unsigned int i = 0; i < vCell->GetNumCellVariables(); ++i)
59 cout <<
" " << vCell->GetCellVarName(i);
64 for (
unsigned int i = 0; i < nSteps; ++i)
67 vCell->TimeIntegrate(vSol, vWsp, vTime);
70 for (
unsigned int i = 0; i < vStimulus.size(); ++i)
72 vStimulus[i]->Update(vWsp, vTime);
76 Vmath::Svtvp(1, vDeltaT, vWsp[0], 1, vSol[0], 1, vSol[0], 1);
82 cout << vTime <<
" " << vSol[0][0];
83 for (
unsigned int j = 0; j < vCell->GetNumCellVariables(); ++j)
85 cout <<
" " << vCell->GetCellSolution(j)[0];
90 for (
unsigned int i = 0; i < vCell->GetNumCellVariables(); ++i)
92 cout <<
"# " << vCell->GetCellVarName(i) <<
" "
93 << vCell->GetCellSolution(i)[0] << endl;
98 cerr <<
"An error occured" << endl;
#define ASSERTL0(condition, msg)
tBaseSharedPtr CreateInstance(tKey idKey BOOST_PP_COMMA_IF(MAX_PARAM) BOOST_PP_ENUM_BINARY_PARAMS(MAX_PARAM, tParam, x))
Create an instance of the class referred to by idKey.
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
void Svtvp(int n, const T alpha, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
svtvp (scalar times vector plus vector): z = alpha*x + y
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
boost::shared_ptr< CellModel > CellModelSharedPtr
A shared pointer to an EquationSystem object.
boost::shared_ptr< ExpList > ExpListSharedPtr
Shared pointer to an ExpList object.
boost::shared_ptr< Equation > EquationSharedPtr
CellModelFactory & GetCellModelFactory()
boost::shared_ptr< PointGeom > PointGeomSharedPtr