BBS: pickle is now a std::vector<char>

src/nrniv/nrnpy.h

@@ -1,5 +1,6 @@
 #pragma once
 #include <string_view>
+#include <vector>
 /**
  * Declarations of global symbols in NEURON that have to be populated with python-version-specific
  * values when dynamic Python is enabled. These are set by the nrnpython_reg_real function, and
@@ -24,7 +25,7 @@ namespace neuron::python {
 struct impl_ptrs {
     Object* (*callable_with_args)(Object*, int narg){};
     double (*call_func)(Object*, int, int*){};
-    char* (*call_picklef)(char*, std::size_t size, int narg, std::size_t* retsize){};
+    std::vector<char> (*call_picklef)(const std::vector<char>&, int narg){};
     void (*call_python_with_section)(Object*, Section*){};
     void (*cmdtool)(Object*, int type, double x, double y, int kd){};
     int (*guigetstr)(Object*, char**){};
@@ -43,9 +44,9 @@ struct impl_ptrs {
     Object* (*mpi_alltoall_type)(int, int){};
     double (*object_to_double)(Object*){};
     void* (*opaque_obj2pyobj)(Object*){};
-    Object* (*pickle2po)(char*, std::size_t size){};
+    Object* (*pickle2po)(const std::vector<char>&){};
     Object* (*po2ho)(PyObject*){};
-    char* (*po2pickle)(Object*, std::size_t* size){};
+    std::vector<char> (*po2pickle)(Object*){};
     double (*praxis_efun)(Object* pycallable, Object* hvec){};
     int (*pysame)(Object* o1, Object* o2){};
     void (*py2n_component)(Object*, Symbol*, int, int){};

src/nrnpython/nrnpy_p2h.cpp

@@ -611,26 +611,26 @@
 static PyObject* dumps;
 
 static void setpickle() {
-    PyObject* pickle;
-    if (!dumps) {
-        pickle = PyImport_ImportModule("pickle");
-        if (pickle) {
-            Py_INCREF(pickle);
-            dumps = PyObject_GetAttrString(pickle, "dumps");
-            loads = PyObject_GetAttrString(pickle, "loads");
-            if (dumps) {
-                Py_INCREF(dumps);
-                Py_INCREF(loads);
-            }
-        }
-        if (!dumps || !loads) {
-            hoc_execerror("Neither Python cPickle nor pickle are available", 0);
+    if (dumps) {
+        return;
+    }
+    PyObject* pickle = PyImport_ImportModule("pickle");
+    if (pickle) {
+        Py_INCREF(pickle);
+        dumps = PyObject_GetAttrString(pickle, "dumps");
+        loads = PyObject_GetAttrString(pickle, "loads");
+        if (dumps) {
+            Py_INCREF(dumps);
+            Py_INCREF(loads);
         }
     }
+    if (!dumps || !loads) {
+        hoc_execerror("Neither Python cPickle nor pickle are available", 0);
+    }
 }
 
 // note that *size includes the null terminating character if it exists
-static char* pickle(PyObject* p, size_t* size) {
+static std::vector<char> pickle(PyObject* p) {
     PyObject* arg = PyTuple_Pack(1, p);
     PyObject* r = nrnpy_pyCallObject(dumps, arg);
     Py_XDECREF(arg);
@@ -639,29 +639,25 @@
     }
     assert(r);
     assert(PyBytes_Check(r));
-    *size = PyBytes_Size(r);
-    char* buf1 = PyBytes_AsString(r);
-    char* buf = new char[*size];
-    for (size_t i = 0; i < *size; ++i) {
-        buf[i] = buf1[i];
-    }
+    std::size_t size = PyBytes_Size(r);
+    char* buf = PyBytes_AsString(r);
+    std::vector<char> ret(buf, buf + size);
     Py_XDECREF(r);
-    return buf;
+    return ret;
 }
 
-static char* po2pickle(Object* ho, std::size_t* size) {
+static std::vector<char> po2pickle(Object* ho) {
     setpickle();
     if (ho && ho->ctemplate->sym == nrnpy_pyobj_sym_) {
         PyObject* po = nrnpy_hoc2pyobject(ho);
-        char* buf = pickle(po, size);
-        return buf;
+        return pickle(po);
     } else {
-        return 0;
+        return {};
     }
 }
 
-static PyObject* unpickle(char* s, size_t size) {
-    PyObject* ps = PyBytes_FromStringAndSize(s, size);
+static PyObject* unpickle(const char* s, std::size_t len) {
+    PyObject* ps = PyBytes_FromStringAndSize(s, len);
     PyObject* arg = PyTuple_Pack(1, ps);
     PyObject* po = nrnpy_pyCallObject(loads, arg);
     assert(po);
@@ -670,9 +666,13 @@
     return po;
 }
 
-static Object* pickle2po(char* s, std::size_t size) {
+static PyObject* unpickle(const std::vector<char>& s) {
+    unpickle(s.data(), s.size());
+}
+
+static Object* pickle2po(const std::vector<char>& s) {
     setpickle();
-    PyObject* po = unpickle(s, size);
+    PyObject* po = unpickle(s);
     Object* ho = nrnpy_pyobject_in_obj(po);
     Py_XDECREF(po);
     return ho;
@@ -739,7 +739,7 @@
     return NULL;
 }
 
-char* call_picklef(char* fname, std::size_t size, int narg, std::size_t* retsize) {
+std::vector<char> call_picklef(const std::vector<char>& fname, int narg) {
     // fname is a pickled callable, narg is the number of args on the
     // hoc stack with types double, char*, hoc Vector, and PythonObject
     // callable return must be pickleable.
@@ -748,7 +748,7 @@
     PyObject* callable;
 
     setpickle();
-    PyObject* ps = PyBytes_FromStringAndSize(fname, size);
+    PyObject* ps = PyBytes_FromStringAndSize(fname.data(), fname.size());
     args = PyTuple_Pack(1, ps);
     callable = nrnpy_pyCallObject(loads, args);
     assert(callable);
@@ -777,7 +777,7 @@
             PyErr_Print();
         }
     }
-    char* rs = pickle(result, retsize);
+    auto rs = pickle(result);
     Py_XDECREF(result);
     return rs;
 }
@@ -811,17 +811,15 @@
 #if NRNMPI
 static PyObject* py_allgather(PyObject* psrc) {
     int np = nrnmpi_numprocs;
-    size_t sz;
-    char* sbuf = pickle(psrc, &sz);
+    auto sbuf = pickle(psrc);
     // what are the counts from each rank
     int* rcnt = new int[np];
-    rcnt[nrnmpi_myid] = int(sz);
+    rcnt[nrnmpi_myid] = sbuf.size();
     nrnmpi_int_allgather_inplace(rcnt, 1);
     int* rdispl = mk_displ(rcnt);
     char* rbuf = new char[rdispl[np]];
 
-    nrnmpi_char_allgatherv(sbuf, rbuf, rcnt, rdispl);
-    delete[] sbuf;
+    nrnmpi_char_allgatherv(sbuf.data(), rbuf, rcnt, rdispl);
 
     PyObject* pdest = char2pylist(rbuf, np, rcnt, rdispl);
     delete[] rbuf;
@@ -832,10 +830,9 @@
 
 static PyObject* py_gather(PyObject* psrc, int root) {
     int np = nrnmpi_numprocs;
-    size_t sz;
-    char* sbuf = pickle(psrc, &sz);
+    auto sbuf = pickle(psrc);
     // what are the counts from each rank
-    int scnt = int(sz);
+    int scnt = sbuf.size();
     int* rcnt = NULL;
     if (root == nrnmpi_myid) {
         rcnt = new int[np];
@@ -848,8 +845,7 @@
         rbuf = new char[rdispl[np]];
     }
 
-    nrnmpi_char_gatherv(sbuf, scnt, rbuf, rcnt, rdispl, root);
-    delete[] sbuf;
+    nrnmpi_char_gatherv(sbuf.data(), scnt, rbuf, rcnt, rdispl, root);
 
     PyObject* pdest = Py_None;
     if (root == nrnmpi_myid) {
@@ -865,25 +861,23 @@
 
 static PyObject* py_broadcast(PyObject* psrc, int root) {
     // Note: root returns reffed psrc.
-    char* buf = NULL;
+    std::vector<char> buf{};
     int cnt = 0;
     if (root == nrnmpi_myid) {
-        size_t sz;
-        buf = pickle(psrc, &sz);
-        cnt = int(sz);
+        buf = pickle(psrc);
+        cnt = buf.size();
     }
     nrnmpi_int_broadcast(&cnt, 1, root);
     if (root != nrnmpi_myid) {
-        buf = new char[cnt];
+        buf.resize(cnt);
     }
-    nrnmpi_char_broadcast(buf, cnt, root);
+    nrnmpi_char_broadcast(buf.data(), cnt, root);
     PyObject* pdest = psrc;
     if (root != nrnmpi_myid) {
-        pdest = unpickle(buf, size_t(cnt));
+        pdest = unpickle(buf);
     } else {
         Py_INCREF(pdest);
     }
-    delete[] buf;
     return pdest;
 }
 #endif
@@ -1003,25 +997,23 @@
                     Py_DECREF(p);
                     continue;
                 }
-                size_t sz;
-                char* b = pickle(p, &sz);
+                auto b = pickle(p);
                 if (size >= 0) {
-                    if (curpos + sz >= bufsz) {
-                        bufsz = bufsz * 2 + sz;
+                    if (curpos + b.size() >= bufsz) {
+                        bufsz = bufsz * 2 + b.size();
                         char* s2 = new char[bufsz];
                         for (size_t i = 0; i < curpos; ++i) {
                             s2[i] = s[i];
                         }
                         delete[] s;
                         s = s2;
                     }
-                    for (size_t j = 0; j < sz; ++j) {
+                    for (size_t j = 0; j < b.size(); ++j) {
                         s[curpos + j] = b[j];
                     }
                 }
-                curpos += sz;
-                scnt[i] = sz;
-                delete[] b;
+                curpos += b.size();
+                scnt[i] = b.size();
                 Py_DECREF(p);
             }
             Py_DECREF(iterator);
@@ -1048,7 +1040,6 @@
             // exchange
             sdispl = mk_displ(scnt);
             rdispl = mk_displ(rcnt);
-            char* r = 0;
             if (size < 0) {
                 pdest = PyTuple_New(2);
                 PyTuple_SetItem(pdest, 0, Py_BuildValue("l", (long) sdispl[np]));
@@ -1058,7 +1049,7 @@
                 delete[] rcnt;
                 delete[] rdispl;
             } else {
-                r = new char[rdispl[np] + 1];  // force > 0 for all None case
+                char* r = new char[rdispl[np] + 1];  // force > 0 for all None case
                 nrnmpi_char_alltoallv(s, scnt, sdispl, r, rcnt, rdispl);
                 delete[] s;
                 delete[] scnt;
@@ -1076,13 +1067,13 @@
             // destination counts
             rcnt = new int[1];
             nrnmpi_int_scatter(scnt, rcnt, 1, root);
-            r = new char[rcnt[0] + 1];  // rcnt[0] can be 0
+            std::vector<char> r(rcnt[0] + 1);  // rcnt[0] can be 0
 
             // exchange
             if (nrnmpi_myid == root) {
                 sdispl = mk_displ(scnt);
             }
-            nrnmpi_char_scatterv(s, scnt, sdispl, r, rcnt[0], root);
+            nrnmpi_char_scatterv(s, scnt, sdispl, r.data(), rcnt[0], root);
             if (s)
                 delete[] s;
             if (scnt)
@@ -1091,13 +1082,12 @@
                 delete[] sdispl;
 
             if (rcnt[0]) {
-                pdest = unpickle(r, size_t(rcnt[0]));
+                pdest = unpickle(r);
             } else {
                 pdest = Py_None;
                 Py_INCREF(pdest);
             }
 
-            delete[] r;
             delete[] rcnt;
             assert(rdispl == NULL);
         }