Custom functions in HH_traub.mod and IM_cortex.mod

[JustasB]

In the original version, the following functions are used in the mod files:

FUNCTION vtrap(x,y) {
    if (fabs(x/y) < 1e-6) {
        vtrap = y*(1 - x/y/2)
    }else{
        vtrap = x/(Exp(x/y)-1)
    }
}

FUNCTION Exp(x) {
    if (x < -100) {
        Exp = 0
    }else{
        Exp = exp(x)
    }
} 
FUNCTION exptable(x) { 
    TABLE  FROM -25 TO 25 WITH 10000

    if ((x > -25) && (x < 25)) {
        exptable = exp(x)
    } else {
        exptable = 0.
    }
}

They probably lead to small differences between LEMS -neuron and original versions.

[pgleeson]

Shouldn’t lead to any significant difference with the LEMS version, these are just used to optimise the exp calculation and catch infinities. You could reduce the 1e-6 even further and just use the standard exp expression instead of Exp/exptable in your version of the mod files to be sure these differences will be removed, but the models will run slightly slower.