Correct Residual Calc

cadet · Sep 10, 2024 · 1431ea7 · 1431ea7
1 parent da8d4bb
commit 1431ea7
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/CADETPythonSimulator/unit_operation.py b/CADETPythonSimulator/unit_operation.py
@@ -513,8 +513,8 @@ def compute_residual(
 
         """
         # Inlet DOFs are simply copied to the residual.
-
-        self.residuals['outlet']['c_poly'] = self.states['outlet']['c_poly']
+        t_poly = np.array([1, t, t**2, t**3])
+        self.residuals['outlet']['c'] = self.c_poly @ t_poly
         self.residuals['outlet']['viscosity'] = self.states['outlet']['viscosity']
 
 class Outlet(UnitOperationBase):
@@ -540,7 +540,8 @@ def compute_residual(
             Time at which to evaluate the residual.
 
         """
-        raise NotImplementedError()
+        self.residuals['inlet']['c'] = np.zeros(self.states['inlet']['c'].shape)
+        self.residuals['inlet']['viscosity'] = 0
 
 
 class Cstr(UnitOperationBase):