Ignored LLPO in the removed nodes; using 0.0 for floating point error

causing/graph.py

@@ -238,6 +238,8 @@ def recalc_graphs(graphs, model, xdat) -> Iterable[networkx.DiGraph]:
     for i, approx_graph in enumerate(graphs):
         individual_xdat = xdat[:, i : i + 1]
         removed_nodes = set(model.graph.nodes) - set(approx_graph.nodes)
+        if "LLPO" in removed_nodes:
+            removed_nodes.remove("LLPO")
 
         # Calc effects on shrunken model
         individual_model = model.shrink(removed_nodes)

causing/model.py

@@ -1,4 +1,5 @@
 from __future__ import annotations
+
 from dataclasses import dataclass, field
 from typing import Iterable, Callable
 from functools import cached_property
@@ -82,8 +83,22 @@ def compute(
                     eq_inputs[:, fixed_from_ind] = fixed_vals
 
                 try:
+                    # print(f"Comuting variable: {self.yvars[i]}")
+                    # yhat[i] = np.array(
+                    #     [eq(*eq_in, *parameters.values()) for eq_in in eq_inputs],
+                    #     dtype=np.float64,
+                    # )
+                    computed_yvars = []
+                    for eq_in in eq_inputs:
+                        try:
+                            computed_yvars.append(eq(*eq_in, *parameters.values()))
+                        except FloatingPointError:
+                            # Floating Point Error for self.yvars[i]
+                            # Adding 0.0 to overcome this.
+                            computed_yvars.append(0.0)
+
                     yhat[i] = np.array(
-                        [eq(*eq_in, *parameters.values()) for eq_in in eq_inputs],
+                        computed_yvars,
                         dtype=np.float64,
                     )
                 except Exception as e: