amici.conserved_quantities_rref

Find conserved quantities deterministically

Functions

nullspace_by_rref(mat)

Compute basis of the nullspace of mat based on the reduced row echelon form

pivots(mat)

Get indices of pivot columns in mat, assumed to be in reduced row echelon form

rref(mat[, round_ndigits])

Bring matrix mat to reduced row echelon form

amici.conserved_quantities_rref.nullspace_by_rref(mat)[source]

Compute basis of the nullspace of mat based on the reduced row echelon form

Return type

numpy.array

amici.conserved_quantities_rref.pivots(mat)[source]

Get indices of pivot columns in mat, assumed to be in reduced row echelon form

Return type

typing.List[int]

amici.conserved_quantities_rref.rref(mat, round_ndigits=None)[source]

Bring matrix mat to reduced row echelon form

see https://en.wikipedia.org/wiki/Row_echelon_form

Parameters
  • mat (numpy.array) – Numpy float matrix to operate on (will be copied)

  • round_ndigits (typing.Union[typing.Literal[False], int, None]) – Number of digits to round intermediary results to, or False to disable rounding completely. Helps to avoid numerical artifacts.

Return type

numpy.array

Returns

mat in rref form.