Added small explanation about the mech_type enum

Author: iomaganaris

coreneuron/nrniv/nrniv_decl.h

@@ -37,6 +37,8 @@ THE POSSIBILITY OF SUCH DAMAGE.
 namespace coreneuron {
 
 /// Mechanism type to be used from stdindex2ptr and nrn_dblpntr2nrncore (in Neuron)
+/// Values of the mechanism types should be negative numbers to avoid any conflict with
+/// mechanism types of Memb_list(>0) or time(0) passed from Neuron
 enum mech_type {voltage = -1, i_membrane_ = -2};
 
 extern int cvode_active_;

coreneuron/nrniv/tqueue.h

@@ -60,7 +60,7 @@ namespace coreneuron {
 #define STRCMP(a, b) (a - b)
 
 class TQItem;
-#define SPBLK TQItem
+#define SPBLK SPC
 #define leftlink left_
 #define rightlink right_
 #define uplink parent_

coreneuron/nrnoc/fast_imem.cpp

@@ -73,9 +73,25 @@ void nrn_calc_fast_imem(NrnThread* _nt) {
 
     double* pd = _nt->nrn_fast_imem->nrn_sav_d;
     double* prhs = _nt->nrn_fast_imem->nrn_sav_rhs;
+    FILE *fp_rhs, *fp_d;
+    char rhs_filename[20], d_filename[20];
+    sprintf(rhs_filename,"rhs.CORENEURON.%d",nrnmpi_myid);
+    sprintf(d_filename,"d.CORENEURON.%d",nrnmpi_myid);
+    fp_rhs = fopen(rhs_filename, "a");
+    fp_d = fopen(d_filename, "a");
+    fprintf(fp_rhs, "\n%.8e time\n", _nt->_t);
+    fprintf(fp_d, "\n%.8e time\n", _nt->_t);
     for (i = i1; i < i3 ; ++i) {
         prhs[i] = (pd[i]*vec_rhs[i] + prhs[i])*vec_area[i]*0.01;
+        fprintf(fp_rhs, "%.8e, ", prhs[i]);
     }
+    fprintf(fp_rhs, "\n");
+    for (i = i1; i < i3 ; ++i) {
+        fprintf(fp_d, "%.8e, ", pd[i]);
+    }
+    fprintf(fp_d, "\n");
+    fclose(fp_rhs);
+    fclose(fp_d);
 }
 
 }

coreneuron/utils/reports/nrnreport.cpp

@@ -88,6 +88,7 @@ class ReportEvent : public DiscreteEvent {
 #pragma omp critical
         {
             // each thread needs to know its own step
+            printf("Real time: %lf\n",t);
             records_nrec(step, gids_to_report.size(), &gids_to_report[0], report_path);
             send(t + dt, nc, nt);
             step++;