OdeLinearSystemSolver Class Reference

#include <OdeLinearSystemSolver.hpp>

Collaboration diagram for OdeLinearSystemSolver:

[legend]

List of all members.

Public Member Functions
	OdeLinearSystemSolver (unsigned systemSize, double timeStep)
	~OdeLinearSystemSolver ()
double	GetTimeStep ()
Mat &	rGetLhsMatrix ()
Vec &	rGetForceVector ()
void	SetInitialConditionVector (Vec initialConditionsVector)
Vec	SolveOneTimeStep ()
Private Attributes
double	mTimeStep
LinearSystem	mLinearSystem
Vec	mForceVector
Vec	mCurrentSolution

---

Detailed Description

Solve large systems of ODEs of the form

M dr/dt = f(t,r),

where M is a matrix, typically large and sparse, and r and f are vectors.

This differs from the OdeSolver classes in the ode component, which are for systems of the form dr/dt = f(t,r) and for small numbers of unknowns.

The solver uses forward euler to dicretise as M r^{n+1} = M r^{n} + dt f and solves this linear system.

The calling code is responsible with setting up M and f each timestep.

Definition at line 49 of file OdeLinearSystemSolver.hpp.

---

Constructor & Destructor Documentation

OdeLinearSystemSolver::OdeLinearSystemSolver	(	unsigned	systemSize,
		double	timeStep
	)

Constructor.

Parameters:

	systemSize	size of the ODE system
	timeStep	the time step used to integrate the ODE system

Definition at line 33 of file OdeLinearSystemSolver.cpp.

References PetscTools::CreateAndSetVec(), mCurrentSolution, mForceVector, and mTimeStep.

OdeLinearSystemSolver::~OdeLinearSystemSolver

(

)

Destructor.

Definition at line 44 of file OdeLinearSystemSolver.cpp.

References mCurrentSolution, and mForceVector.

---

Member Function Documentation

double OdeLinearSystemSolver::GetTimeStep

(

)

Get the timestep for the solver.

Definition at line 50 of file OdeLinearSystemSolver.cpp.

References mTimeStep.

Vec & OdeLinearSystemSolver::rGetForceVector

(

)

Get the force vector.

Definition at line 60 of file OdeLinearSystemSolver.cpp.

References mForceVector.

Referenced by NodeBasedCellPopulationWithBuskeUpdate< DIM >::UpdateNodeLocations().

Mat & OdeLinearSystemSolver::rGetLhsMatrix

(

)

Get the matrix.

Definition at line 55 of file OdeLinearSystemSolver.cpp.

References mLinearSystem, and LinearSystem::rGetLhsMatrix().

Referenced by NodeBasedCellPopulationWithBuskeUpdate< DIM >::UpdateNodeLocations().

void OdeLinearSystemSolver::SetInitialConditionVector

(

Vec

initialConditionsVector

)

Set the initial conditions.

Parameters:

initialConditionsVector

the initial condition

Definition at line 65 of file OdeLinearSystemSolver.cpp.

References mCurrentSolution.

Referenced by NodeBasedCellPopulationWithBuskeUpdate< DIM >::UpdateNodeLocations().

Vec OdeLinearSystemSolver::SolveOneTimeStep

(

)

Solve the ODE system over one time step.

Definition at line 70 of file OdeLinearSystemSolver.cpp.

References PetscVecTools::AddScaledVector(), mCurrentSolution, mForceVector, mLinearSystem, mTimeStep, LinearSystem::rGetLhsMatrix(), LinearSystem::rGetRhsVector(), and LinearSystem::Solve().

Referenced by NodeBasedCellPopulationWithBuskeUpdate< DIM >::UpdateNodeLocations().

---

Member Data Documentation

Vec OdeLinearSystemSolver::mCurrentSolution [private]

Solution at current timestep.

Definition at line 63 of file OdeLinearSystemSolver.hpp.

Referenced by OdeLinearSystemSolver(), SetInitialConditionVector(), SolveOneTimeStep(), and ~OdeLinearSystemSolver().

Vec OdeLinearSystemSolver::mForceVector [private]

Force vector (f in M dr/dt = f).

Definition at line 60 of file OdeLinearSystemSolver.hpp.

Referenced by OdeLinearSystemSolver(), rGetForceVector(), SolveOneTimeStep(), and ~OdeLinearSystemSolver().

LinearSystem OdeLinearSystemSolver::mLinearSystem [private]

The LHS matrix and the force vector. Solve() is called on this.

Definition at line 57 of file OdeLinearSystemSolver.hpp.

Referenced by rGetLhsMatrix(), and SolveOneTimeStep().

double OdeLinearSystemSolver::mTimeStep [private]

Timestep for solver.

Definition at line 54 of file OdeLinearSystemSolver.hpp.

Referenced by GetTimeStep(), OdeLinearSystemSolver(), and SolveOneTimeStep().

---

The documentation for this class was generated from the following files:

• cell_based/src/common/OdeLinearSystemSolver.hpp
• cell_based/src/common/OdeLinearSystemSolver.cpp