Add iaf_psc_exp_ps_lossless model by Krishnan et al

precise/CMakeLists.txt

@@ -23,6 +23,7 @@ set( precise_sources
     iaf_psc_alpha_canon.cpp iaf_psc_alpha_canon.h
     iaf_psc_alpha_presc.cpp iaf_psc_alpha_presc.h
     iaf_psc_exp_ps.cpp iaf_psc_exp_ps.h
+    iaf_psc_exp_ps_lossless.cpp iaf_psc_exp_ps_lossless.h
     poisson_generator_ps.cpp poisson_generator_ps.h
     parrot_neuron_ps.cpp parrot_neuron_ps.h
     precisemodule.cpp precisemodule.h

precise/iaf_psc_alpha_canon.h

@@ -102,9 +102,8 @@ The following parameters can be set in the status dictionary.
                         0-none, 1-linear, 2-quadratic, 3-cubic
 
 Remarks:
-If tau_m is very close to tau_syn_ex or tau_syn_in, the model
-will numerically behave as if tau_m is equal to tau_syn_ex or
-tau_syn_in, respectively, to avoid numerical instabilities.
+If tau_m is very close to tau_syn, the model will numerically behave as
+if tau_m is equal to tau_syn, to avoid numerical instabilities.
 For details, please see IAF_Neruons_Singularity.ipynb in
 the NEST source code (docs/model_details).
 

precise/iaf_psc_alpha_presc.h

@@ -81,9 +81,8 @@
   in addition to the on-grid spike times.
 
   Remarks:
-  If tau_m is very close to tau_syn_ex or tau_syn_in, the model
-  will numerically behave as if tau_m is equal to tau_syn_ex or
-  tau_syn_in, respectively, to avoid numerical instabilities.
+  If tau_m is very close to tau_syn, the model will numerically behave
+  as if tau_m is equal to tau_syn, to avoid numerical instabilities.
   For details, please see IAF_Neruons_Singularity.ipynb in
   the NEST source code (docs/model_details).