This class encapsulates a general timestep solver concept to solve linear PDE in time and space. More...

#include <timestepping.hh>

Inheritance diagram for timestepping::TimeStepping:

Collaboration diagram for timestepping::TimeStepping:

Public Member Functions
void	operator() (concepts::Vector< Real > &sol, uint n)
	Returns the solution of the n'th timestep.
uint	operator() (concepts::Vector< Real > &sol)
Real	time () const
	Return the time of the current solution.
	TimeStepping (TimeStepStrategy &timeStepStrategy)
	Constructor.
Protected Member Functions
virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream.
Private Attributes
uint	n_
	The number of the timestep of the actual solution.
TimeStepStrategy &	timeStepStrategy_
	Timestep strategy of the solver.

Detailed Description

This class encapsulates a general timestep solver concept to solve linear PDE in time and space.

Any linear partial differential equation of the form

$( D_k \partial_t^k + \dots + D_1 \partial_t + D_0 ) y(x,t) = f(x,t),$

where D_n, ..., D₀ are space operators that not depend on time. After the discretisation it can always be brought in the form

$( D + const.(\Delta t) D') Y_n = V(Y_{n-1}, Y_{n-2}, ...),$

where the LHS depends only on the timestep size and the RHS can be evalated from previous steps. Therefore the new solution Y_n is best obtained with an LU separation of the LHS.

Author:: Manuel Walser, 2002

Definition at line 39 of file timestepping.hh.

Constructor & Destructor Documentation

timestepping::TimeStepping::TimeStepping ( TimeStepStrategy & timeStepStrategy )

Constructor.

Parameters:

timeStepStrategy Timestep strategy for solving the PDE

Member Function Documentation

virtual std::ostream& timestepping::TimeStepping::info ( std::ostream & os ) const [protected, virtual]

Returns information in an output stream.

Reimplemented from concepts::OutputOperator.

uint timestepping::TimeStepping::operator() ( concepts::Vector< Real > & sol )

void timestepping::TimeStepping::operator()	(	concepts::Vector< Real > &	sol,
		uint	n
	)

Returns the solution of the n'th timestep.

With the parameter n the actual solution is return, otherwise the specified step in the futur. A step number n before the actual step number is not possible and throws an exception

Parameters: