From f4f3ac9b0ee5fbd9ec79ea5caa4b88516663992a Mon Sep 17 00:00:00 2001
From: Even Solbraa <41290109+EvenSol@users.noreply.github.com>
Date: Sun, 26 Nov 2023 17:22:21 +0000
Subject: [PATCH] update

---
 .../flashOps/TPmultiflash.java                   | 16 ++++++++++------
 1 file changed, 10 insertions(+), 6 deletions(-)

diff --git a/src/main/java/neqsim/thermodynamicOperations/flashOps/TPmultiflash.java b/src/main/java/neqsim/thermodynamicOperations/flashOps/TPmultiflash.java
index 800841f27c..b703d029a0 100644
--- a/src/main/java/neqsim/thermodynamicOperations/flashOps/TPmultiflash.java
+++ b/src/main/java/neqsim/thermodynamicOperations/flashOps/TPmultiflash.java
@@ -371,12 +371,14 @@ public void stabilityAnalysis() {
       int iter = 0;
       double errOld = 1.0e100;
       boolean useaccsubst = true;
+      int maxsucssubiter = 50;
+      int maxiter = 100;
       do {
         errOld = err;
         iter++;
         err = 0;
 
-        if (iter <= 150 || !system.isImplementedCompositionDeriativesofFugacity()) {
+        if (iter <= maxsucssubiter || !system.isImplementedCompositionDeriativesofFugacity()) {
           if (iter % 7 == 0 && useaccsubst) {
             double vec1 = 0.0;
 
@@ -513,7 +515,7 @@ public void stabilityAnalysis() {
             clonedSystem.get(0).getPhase(1).getComponents()[i].setx(1e-50);
           }
         }
-      } while ((Math.abs(err) > 1e-9 || err > errOld) && iter < 600);
+      } while ((Math.abs(err) > 1e-9 || err > errOld) && iter < maxiter);
       // logger.info("err: " + err + " ITER " + iter);
       double xTrivialCheck0 = 0.0;
       double xTrivialCheck1 = 0.0;
@@ -530,7 +532,7 @@ public void stabilityAnalysis() {
         xTrivialCheck0 += Math.abs(x[j][i] - system.getPhase(0).getComponent(i).getx());
         xTrivialCheck1 += Math.abs(x[j][i] - system.getPhase(1).getComponent(i).getx());
       }
-      if (iter >= 599) {
+      if (iter >= maxiter) {
         // logger.info("iter > maxiter multiphase stability ");
         // logger.info("error " + Math.abs(err));
         // logger.info("tm: " + tm[j]);
@@ -716,12 +718,14 @@ public void stabilityAnalysis3() {
       int iter = 0;
       double errOld = 1.0e100;
       boolean useaccsubst = true;
+      int maxsucssubiter = 50;
+      int maxiter = 100;
       do {
         errOld = err;
         iter++;
         err = 0;
 
-        if (iter <= 150 || !system.isImplementedCompositionDeriativesofFugacity()) {
+        if (iter <= maxsucssubiter || !system.isImplementedCompositionDeriativesofFugacity()) {
           if (iter % 7 == 0 && useaccsubst) {
             double vec1 = 0.0;
 
@@ -858,7 +862,7 @@ public void stabilityAnalysis3() {
             clonedSystem.get(0).getPhase(1).getComponents()[i].setx(1e-50);
           }
         }
-      } while ((Math.abs(err) > 1e-9 || err > errOld) && iter < 600);
+      } while ((Math.abs(err) > 1e-9 || err > errOld) && iter < maxiter);
       // logger.info("err: " + err + " ITER " + iter);
       double xTrivialCheck0 = 0.0;
       double xTrivialCheck1 = 0.0;
@@ -875,7 +879,7 @@ public void stabilityAnalysis3() {
         xTrivialCheck0 += Math.abs(x[j][i] - system.getPhase(0).getComponent(i).getx());
         xTrivialCheck1 += Math.abs(x[j][i] - system.getPhase(1).getComponent(i).getx());
       }
-      if (iter >= 599) {
+      if (iter >= maxiter - 1) {
         // logger.info("iter > maxiter multiphase stability ");
         // logger.info("error " + Math.abs(err));
         // logger.info("tm: " + tm[j]);