Hello,
Is there any way to access the Hines matrix in hoc?
That is, the almost tridiagonal set of linear equation, used in, e.g. the cnexp model.
Best regards,
Roy
Access to the set of linear equations / Hines matrix
Re: Access to the set of linear equations / Hines matrix
http://www.neuron.yale.edu/neuron/stati ... ml#fmatrix
One should call fcurrent() prior to fmatrix() to ensure the matrix is setup.
The matrix printed is for thread 0.
The 4 numbers printed on each line are
compartment index
jacobian element for parent in node equation (ClassicalNODEB)
jacobian element for node in parent equation (ClassicalNODEA)
diagonal element (NODED)
right hand side (NODERHS)
See nrn/src/nrnoc/fadvance.c
One should call fcurrent() prior to fmatrix() to ensure the matrix is setup.
The matrix printed is for thread 0.
The 4 numbers printed on each line are
compartment index
jacobian element for parent in node equation (ClassicalNODEB)
jacobian element for node in parent equation (ClassicalNODEA)
diagonal element (NODED)
right hand side (NODERHS)
See nrn/src/nrnoc/fadvance.c
Re: Access to the set of linear equations / Hines matrix
Thanks a lot it is very helpfull