more debugging

CarlosP24 · May 13, 2024 · f036dbc · f036dbc
1 parent 027dfd2
commit f036dbc
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/utilities/LP_usadel.jl b/src/utilities/LP_usadel.jl
@@ -18,7 +18,7 @@ function LP_lobe(n, ξd, R, d)
     RLP = R + d/2
     pre = 1.76 * π * exp(π/4) / 4
     frac1 = 1 / (1 + (d / (2 * RLP))^2)
-    frac2 = 1 / (1 + (2 * RLP / d)^2)
+    frac2 = d^2 / (d^2 + 4 * RLP^2)
     frac3 = RLP^2 / ξd^2
     root = pre * frac3 * frac1 + n^2 * frac2 * (frac2 - 4/3)
     ΦLPa = n * frac1 - sqrt(root)