From 616d3ef9d96759352ccc46ce6821dd69e42d6bac Mon Sep 17 00:00:00 2001
From: James Tocknell <aragilar@gmail.com>
Date: Tue, 21 Jan 2025 15:34:21 +1100
Subject: [PATCH] =?UTF-8?q?Remove=20the=20solver=20exiting=20when=20v=5F?=
 =?UTF-8?q?=CE=B8=20<=200?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

This may be stopping the solver at higher θ when it doesn't need to.
---
 src/disc_solver/solve/solution.py | 5 -----
 1 file changed, 5 deletions(-)

diff --git a/src/disc_solver/solve/solution.py b/src/disc_solver/solve/solution.py
index 250e5677..98d29b62 100644
--- a/src/disc_solver/solve/solution.py
+++ b/src/disc_solver/solve/solution.py
@@ -137,11 +137,6 @@ def rhs_equation(x, params, derivs):
                 # pylint: disable=unsubscriptable-object
                 problems[θ].append("negative density")
             return 1
-        if not derivs_post_solution and v_θ < 0:
-            if store_internal:
-                # pylint: disable=unsubscriptable-object
-                problems[θ].append("negative velocity")
-            return -1
 
         B_mag = sqrt(B_r**2 + B_φ**2 + B_θ**2)
         b_r, b_φ, b_θ = B_r/B_mag, B_φ/B_mag, B_θ/B_mag