fix a bug

numericalEFT · Jun 16, 2023 · 495633d · 495633d · kunyuan · Jun 16, 2023
1 parent e1235c2
commit 495633d
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "ElectronGas"
 uuid = "8c98e03e-8755-48b1-a5ef-1e2e7812a7c9"
 authors = ["Xiansheng Cai, Pengcheng Hou, Kun Chen"]
-version = "0.2.4"
+version = "0.2.5"
 
 [deps]
 CompositeGrids = "b5136c89-beeb-4521-9139-60d2cac8be56"

diff --git a/src/polarization.jl b/src/polarization.jl
@@ -90,7 +90,8 @@ function finitetemp_kgrid_ladder(μ::Float64, m::Float64, q::Float64, kF::Float6
     mink = (q < 1e-16 / minterval) ? minterval * kF : minterval * min(q, kF)
     mixX = mink / (1 + mink)
     if q >= sqrt(8 * μ * m)
-        kgrid = CompositeGrid.LogDensedGrid(:gauss, [0.0, 1.0], [0.0, kF], scaleN, mink, gaussN)
+        x = kF / (1 + kF)
+        kgrid = CompositeGrid.LogDensedGrid(:gauss, [0.0, 1.0], [0.0, x], scaleN, mink, gaussN)
     else
         k1 = abs(q - sqrt(8 * μ * m - q^2)) / 2
         k2 = (q + sqrt(8 * μ * m - q^2)) / 2