GreenPolyProducts/cantera
5
0
Fork 0
mirror of https://github.com/Cantera/cantera.git synced 2025-02-25 18:55:29 -06:00
Code Issues Packages Projects Releases Wiki Activity

[Numerics] Make sparse Jacobian accessible from Python

Browse Source
This commit is contained in:
Ray Speth 2025-01-03 13:12:59 -05:00 committed by Ingmar Schoegl
parent 9fa3dda0a7
commit 74340f8e03
2 changed files with 9 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
interfaces/cython/cantera/jacobians.pxd
View File

@ -20,6 +20,7 @@ cdef extern from "cantera/numerics/EigenSparseJacobian.h" namespace "Cantera":
CxxEigenSparseJacobian() except +translate_exception
void printPreconditioner()
CxxSparseMatrix matrix() except +translate_exception
CxxSparseMatrix jacobian() except +translate_exception
cdef extern from "cantera/numerics/AdaptivePreconditioner.h" namespace "Cantera":
cdef cppclass CxxAdaptivePreconditioner "Cantera::AdaptivePreconditioner" \

8
interfaces/cython/cantera/jacobians.pyx
View File

@ -86,6 +86,14 @@ cdef class EigenSparseJacobian(SystemJacobian):
cdef CxxSparseMatrix smat = self.sparse_jac.matrix()
return get_from_sparse(smat, smat.rows(), smat.cols())
property jacobian:
"""
Property to retrieve the latest Jacobian.
"""
def __get__(self):
cdef CxxSparseMatrix smat = self.sparse_jac.jacobian()
return get_from_sparse(smat, smat.rows(), smat.cols())
cdef class EigenSparseDirectJacobian(EigenSparseJacobian):
"""