ODE IV Control and Variable Scale

Kevin H. Hobbs hobbsk@ohiou.edu
Thu Mar 18 18:28:00 GMT 2010

I'm trying make sure that I am using the gsl_odeiv_control functions

I'm modelling neurons and my variables have different scales. Voltage is
usually between -100 and 50 mV. The membrane current gating variables
stay between 0 and 1. Internal calcium concentration is usually between
0 and a few hundred micro molar.

Each time I read the descriptions of the gsl_odeiv_control functions I
become more convinced that the parameters eps_rel, a_y, and a_dydt are
there to excuse some error in variables and derivatives that are far
from zero. Is that correct?

Since I have no reason to care more about values close to zero (unless I
do some complicated remapping of the variable values) should I set
eps_rel, a_y, and a_dydt to 0?

Since my variables have different ranges should I use
gsl_odeiv_control_scaled_new with an appropriate eps_abs and scale_abs
set with something like 1 for internal calcium concentration, 0.3 for
membrane voltage, and 0.001 for all of the gating variables?

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 252 bytes
Desc: OpenPGP digital signature
URL: <http://sourceware.org/pipermail/gsl-discuss/attachments/20100318/cf675036/attachment.sig>

More information about the Gsl-discuss mailing list