Class NewtonSolverDense

• Defined in File newton_solver.h

Inheritance Relationships

Base Type

• public amici::NewtonSolver (Class NewtonSolver)

Class Documentation

class NewtonSolverDense : public amici::NewtonSolver

The NewtonSolverDense provides access to the dense linear solver for the Newton method.

Public Functions

explicit NewtonSolverDense(Model const &model)

constructor for sparse solver

Parameters:

model – model instance that provides problem dimensions

NewtonSolverDense(NewtonSolverDense const&) = delete

NewtonSolverDense &operator=(NewtonSolverDense const &other) = delete

~NewtonSolverDense() override

virtual void solveLinearSystem(AmiVector &rhs) override

Solves the linear system for the Newton step.

Parameters:

rhs – containing the RHS of the linear system, will be overwritten by solution to the linear system

virtual void prepareLinearSystem(Model &model, SimulationState const &state) override

Writes the Jacobian for the Newton iteration and passes it to the linear solver.

Parameters:

• model – pointer to the model instance

• state – current simulation state

virtual void prepareLinearSystemB(Model &model, SimulationState const &state) override

Writes the Jacobian (JB) for the Newton iteration and passes it to the linear solver

Parameters:

• model – pointer to the model instance

• state – current simulation state

virtual void reinitialize() override

Reinitialize the linear solver.

virtual bool is_singular(Model &model, SimulationState const &state) const override

Checks whether linear system is singular.

Parameters:

• model – pointer to the model instance

• state – current simulation state

Returns:

boolean indicating whether the linear system is singular (condition number < 1/machine precision)