check for sigma added

lanl · Mar 7, 2024 · fcb26d0 · fcb26d0
1 parent 0a2f59c
commit fcb26d0
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/analysis/analysis.cpp b/src/analysis/analysis.cpp
@@ -6,6 +6,9 @@ std::shared_ptr<StateDescriptor> Initialize(ParameterInput *pin) {
   auto analysis_pkg = std::make_shared<StateDescriptor>("analysis");
   Params &params = analysis_pkg->AllParams();
   Real sigma = pin->GetOrAddReal("analysis", "sigma", 2e-3);
+  if (sigma < 0){
+    PARTHENON_THROW("sigma must be greater than 0");
+  }
   params.Add("sigma", sigma);
 
   return analysis_pkg;