From 027dfd234112c9d5e6d2367f6976fb6310e7ac98 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Carlos=20Pay=C3=A1?= <carlos.paya@csic.es>
Date: Mon, 13 May 2024 15:05:36 +0200
Subject: [PATCH] debug destructive

---
 src/utilities/LP_usadel.jl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/utilities/LP_usadel.jl b/src/utilities/LP_usadel.jl
index de284e5..ae8f4ea 100644
--- a/src/utilities/LP_usadel.jl
+++ b/src/utilities/LP_usadel.jl
@@ -20,7 +20,7 @@ function LP_lobe(n, ξd, R, d)
     frac1 = 1 / (1 + (d / (2 * RLP))^2)
     frac2 = 1 / (1 + (2 * RLP / d)^2)
     frac3 = RLP^2 / ξd^2
-    root = pre * frac3 * frac1 - n^2 * frac2 * (1/3 + frac1)
+    root = pre * frac3 * frac1 + n^2 * frac2 * (frac2 - 4/3)
     ΦLPa = n * frac1 - sqrt(root)
     ΦLPb = n * frac1 + sqrt(root)
     return ΦLPa, ΦLPb