OdeLinearSystemSolver Class Reference

#include <OdeLinearSystemSolver.hpp>

Collaboration diagram for OdeLinearSystemSolver:

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 56 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 38 of file OdeLinearSystemSolver.cpp.

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

OdeLinearSystemSolver::~OdeLinearSystemSolver

(

)

Destructor.

Definition at line 49 of file OdeLinearSystemSolver.cpp.

References PetscTools::Destroy(), mCurrentSolution, and mForceVector.

Member Function Documentation

double OdeLinearSystemSolver::GetTimeStep

(

)

Returns

the timestep for the solver.

Definition at line 55 of file OdeLinearSystemSolver.cpp.

References mTimeStep.

Vec & OdeLinearSystemSolver::rGetForceVector

(

)

Returns

the force vector.

Definition at line 65 of file OdeLinearSystemSolver.cpp.

References mForceVector.

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

Mat & OdeLinearSystemSolver::rGetLhsMatrix

(

)

Returns

the matrix.

Definition at line 60 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 70 of file OdeLinearSystemSolver.cpp.

References mCurrentSolution.

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

Vec OdeLinearSystemSolver::SolveOneTimeStep

(

)

Solve the ODE system over one time step.

Returns

solution vector

Definition at line 75 of file OdeLinearSystemSolver.cpp.

References PetscVecTools::AddScaledVector(), PetscTools::Destroy(), 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 70 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 67 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 64 of file OdeLinearSystemSolver.hpp.

Referenced by rGetLhsMatrix(), and SolveOneTimeStep().

double OdeLinearSystemSolver::mTimeStep

private

Timestep for solver.

Definition at line 61 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