Access to the set of linear equations / Hines matrix

Anything that doesn't fit elsewhere.
Post Reply
roybens
Posts: 44
Joined: Fri Mar 14, 2008 7:57 am

Access to the set of linear equations / Hines matrix

Post by roybens » Wed May 09, 2012 1:56 pm

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

hines
Site Admin
Posts: 1580
Joined: Wed May 18, 2005 3:32 pm

Re: Access to the set of linear equations / Hines matrix

Post by hines » Thu May 10, 2012 7:34 am

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

roybens
Posts: 44
Joined: Fri Mar 14, 2008 7:57 am

Re: Access to the set of linear equations / Hines matrix

Post by roybens » Thu May 10, 2012 12:58 pm

Thanks a lot it is very helpfull

Post Reply