Improvement for Presyn and queue event transfer (#673)

* The issue was encountered while testing model of 250k cells
   from Ivan, see: neuronsimulator/nrn#1472
* Neuron queue transfer of PreSyn may be CoreNEURON InputPreSyn.
* Needed explicit transfer indication of PreSyn or InputPreSyn
* After direct mode initialize, clear the spike vectors.
  Those spikes have already been recorded in NEURON.

Author: nrnhines

coreneuron/apps/main1.cpp

@@ -342,6 +342,7 @@ void nrn_init_and_load_data(int argc,
         // to NEURON. Here there is some first time only
         // initialization and queue transfer.
         direct_mode_initialize();
+        clear_spike_vectors();  // PreSyn send already recorded by NEURON
         (*nrn2core_part2_clean_)();
     }
 

coreneuron/io/nrn2core_data_init.cpp

@@ -245,21 +245,28 @@ static void nrn2core_tqueue() {
                     } break;
 
                     case 4: {  // PreSyn
-                        int ps_index = ncte->intdata[idat++];
+                        int type = ncte->intdata[idat++];
+                        if (type == 0) {  // CoreNEURON PreSyn
+                            int ps_index = ncte->intdata[idat++];
 #if CORENRN_DEBUG_QUEUE
-                        printf("nrn2core_tqueue tid=%d i=%zd type=%d tdeliver=%g PreSyn %d\n",
-                               tid,
-                               i,
-                               ncte->type[i],
-                               ncte->td[i],
-                               ps_index);
+                            printf("nrn2core_tqueue tid=%d i=%zd type=%d tdeliver=%g PreSyn %d\n",
+                                   tid,
+                                   i,
+                                   ncte->type[i],
+                                   ncte->td[i],
+                                   ps_index);
 #endif
-                        PreSyn* ps = nt.presyns + ps_index;
-                        int gid = ps->output_index_;
-                        // Following assumes already sent to other machines.
-                        ps->output_index_ = -1;
-                        ps->send(ncte->td[i], net_cvode_instance, &nt);
-                        ps->output_index_ = gid;
+                            PreSyn* ps = nt.presyns + ps_index;
+                            int gid = ps->output_index_;
+                            // Following assumes already sent to other machines.
+                            ps->output_index_ = -1;
+                            ps->send(ncte->td[i], net_cvode_instance, &nt);
+                            ps->output_index_ = gid;
+                        } else {  // CoreNEURON InputPreSyn
+                            int gid = ncte->intdata[idat++];
+                            InputPreSyn* ps = gid2in[gid];
+                            ps->send(ncte->td[i], net_cvode_instance, &nt);
+                        }
                     } break;
 
                     case 6: {  // PlayRecordEvent