Check for negative discriminant in time marcher point quadratic.

scikit-fmm · Jul 16, 2017 · 5a1e1db · 5a1e1db
1 parent 2354c65
commit 5a1e1db
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/skfmm/travel_time_marcher.cpp b/skfmm/travel_time_marcher.cpp
@@ -3,6 +3,7 @@
 #include "travel_time_marcher.h"
 #include "math.h"
 #include "heap.h"
+#include <stdexcept>
 
 void travelTimeMarcher::initalizeFrozen()
 {
@@ -82,9 +83,7 @@ double travelTimeMarcher::solveQuadratic(int i, const double &a,
   }
   else
   {
-    return 0;
+    throw std::runtime_error("Negative discriminant in time marcher quadratic.");
   }
   return r0;
 }
-
-