Raise error when dividing vector by zero (#2605)

neuronsimulator · Nov 9, 2023 · 9700e91 · 9700e91
1 parent c1ac797
commit 9700e91
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/ivoc/ivocvect.cpp b/src/ivoc/ivocvect.cpp
@@ -2369,7 +2369,11 @@ static Object** v_mul(void* v1) {
 static Object** v_div(void* v1) {
     Vect* x = (Vect*) v1;
     if (hoc_argtype(1) == NUMBER) {
-        std::for_each(x->begin(), x->end(), [](double& d) { d /= *getarg(1); });
+        if (*getarg(1) == 0.0) {
+            hoc_execerror("Vector", "Division by zero");
+        } else {
+            std::for_each(x->begin(), x->end(), [](double& d) { d /= *getarg(1); });
+        }
     }
     if (hoc_is_object_arg(1)) {
         Vect* y = vector_arg(1);