Check that norm_val is positive before dividing

dknez · dknez · commit 36305d2903f2 · 2025-01-29T15:05:24.000-05:00
diff --git a/src/reduced_basis/rb_eim_construction.C b/src/reduced_basis/rb_eim_construction.C
@@ -684,7 +684,7 @@ void RBEIMConstruction::apply_normalization_to_solution_snapshots()
       // error tolerance in the Greedy training, and the L2-norm-scaled value of
       // _max_abs_value_in_training_set that we obtain here should be sufficient
       // for that purpose.
-      if (i == _max_abs_value_in_training_set_index)
+      if ((i == _max_abs_value_in_training_set_index) && (norm_val > 0.))
         _max_abs_value_in_training_set /= norm_val;
     }
 

Original file line number	Diff line number	Diff line change
`@@ -684,7 +684,7 @@ void RBEIMConstruction::apply_normalization_to_solution_snapshots()`
`684`	`684`	`// error tolerance in the Greedy training, and the L2-norm-scaled value of`
`685`	`685`	`// _max_abs_value_in_training_set that we obtain here should be sufficient`
`686`	`686`	`// for that purpose.`
`687`		`- if (i == _max_abs_value_in_training_set_index)`
	`687`	`+ if ((i == _max_abs_value_in_training_set_index) && (norm_val > 0.))`
`688`	`688`	`_max_abs_value_in_training_set /= norm_val;`
`689`	`689`	`}`
`690`	`690`