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(const Model &model)

constructor for sparse solver

Parameters

model – model instance that provides problem dimensions

NewtonSolverDense(const NewtonSolverDense&) = delete

NewtonSolverDense &operator=(const NewtonSolverDense &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, const SimulationState &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, const SimulationState &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, const SimulationState &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)