Add missing partial derivative block in jacobian of heteroazeotrope_t (#320)

g-bauer · web-flow · commit 051888cbd51d · 2025-11-18T16:24:31.000+01:00
diff --git a/crates/feos-core/src/phase_equilibria/phase_diagram_binary.rs b/crates/feos-core/src/phase_equilibria/phase_diagram_binary.rs
@@ -442,10 +442,10 @@ impl<E: Residual> PhaseEquilibrium<E, 3> {
 
             // calculate Jacobian
             let jacobian = stack![
-                    dmu_drho_l1, 0, 0;
-                    0, dmu_drho_l2, -dmu_drho_v;
-                    dp_drho_l1, 0, -&dp_drho_v;
-                    0, dp_drho_l2, -dp_drho_v
+                dmu_drho_l1, 0          , -&dmu_drho_v;
+                0          , dmu_drho_l2, -dmu_drho_v;
+                dp_drho_l1 , 0          , -&dp_drho_v;
+                0          , dp_drho_l2 , -dp_drho_v
             ];
 
             // calculate Newton step
