testm

ipqa-research · Oct 22, 2024 · 69c3219 · 69c3219
1 parent 663c5a3
commit 69c3219
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 70 deletions.
diff --git a/src/equilibria/boundaries/phase_envelopes_pt.f90 b/src/equilibria/boundaries/phase_envelopes_pt.f90
@@ -238,10 +238,10 @@ subroutine update_spec(X, ns, S, dS, dXdS, step_iters)
             ] &
             )
 
-         do while(maxval(abs(dXdS(:nc)*dS)) > 0.07)
+         do while(maxval(abs(dXdS(:nc)*dS)) > 0.1 * maxval(abs(X(:nc))))
             dS = 0.7*dS
          end do
-
+         
          call save_point(X, step_iters)
          call detect_critical(X, dXdS, ns, S, dS)
       end subroutine update_spec
@@ -316,9 +316,8 @@ subroutine detect_critical(X, dXdS, ns, S, dS)
 
          inner = 0
          do while (&
-                  maxval(abs(X(:nc))) < 0.1 &
-            .and. inner < 500&
-            .and. abs(Vz - Vy) < 0.1)
+                  maxval(abs(X(:nc))) < 0.06 &
+            .and. inner < 5000)
             ! If near a critical point, jump over it
             inner = inner + 1
 

diff --git a/test/test_psrk.f90 b/test/test_psrk.f90