preeq_wrms and posteq_wrms values are zero with FSA

[plakrisenko]

The values should be equal to "weighted root-mean-square of the rhs when steady state was reached".

I get values equal to "0.0" in case
solver.setSensitivityOrder(amici.SensitivityOrder_first)
amici_solver.setSensitivityMethod(SensitivityMethod.forward)
amici_solver.setNewtonMaxSteps(0)
amici_model.setSteadyStateSensitivityMode(amici.SteadyStateSensitivityMode.integrationOnly).
However, in case I use SensitivityMethod.adjoint I get values different from "0.0", e.g. "0.91".

• AMICI version: 0.16.1

I suppose, it's related to the line here
https://github.com/AMICI-dev/AMICI/blob/master/src/steadystateproblem.cpp#L711

I think, in the end the output of getWrmsFSA is saved to preeq_wrms (or posteq_wrms), not output of getWrms. Not sure why it's always "0.0", though. I tried with different tolerances, fsa tolerances.

[dweindl]

I think, in the end the output of getWrmsFSA is saved to preeq_wrms (or posteq_wrms), not output of getWrms.

Correct.

Agreed that this doesn't match the current documentation at

AMICI/include/amici/rdata.h

Lines 332 to 336 in 1a99308

	/**
	* weighted root-mean-square of the rhs when steadystate
	* was reached (preequilibration)
	*/
	realtype preeq_wrms = NAN;

or

AMICI/include/amici/rdata.h

Lines 343 to 347 in 1a99308

	/**
	* weighted root-mean-square of the rhs when steadystate
	* was reached (postequilibration)
	*/
	realtype posteq_wrms = NAN;