From 70da3da1cc57259dfabc32693244996974f2cf0b Mon Sep 17 00:00:00 2001 From: Anthony Scemama Date: Thu, 10 Jun 2010 16:53:46 +0200 Subject: [PATCH] Bug in density --- src/density.irp.f | 4 ++++ src/test_1d.irp.f | 3 ++- 2 files changed, 6 insertions(+), 1 deletion(-) diff --git a/src/density.irp.f b/src/density.irp.f index 86f5669..9361f57 100644 --- a/src/density.irp.f +++ b/src/density.irp.f @@ -53,6 +53,8 @@ BEGIN_PROVIDER [ real, density_alpha_value_p ] exc(4) = exc(1)*exc(2) if (exc(4) /= 0) then exc(4) = exc(4)/abs(exc(4)) + else + exc(4) = 1 endif phase = dble(exc(4)) @@ -100,6 +102,8 @@ BEGIN_PROVIDER [ real, density_beta_value_p ] exc(4) = exc(1)*exc(2) if (exc(4) /= 0) then exc(4) = exc(4)/abs(exc(4)) + else + exc(4) = 1 endif phase = dble(exc(4)) diff --git a/src/test_1d.irp.f b/src/test_1d.irp.f index 3873c92..2fe848a 100644 --- a/src/test_1d.irp.f +++ b/src/test_1d.irp.f @@ -13,6 +13,7 @@ subroutine run do i=-40,40 point(3) = real(i)/20. TOUCH point - print *, point(3), eplf_value_p, eplf_up_up, eplf_up_dn +! print *, point(3), eplf_value_p, eplf_up_up, eplf_up_dn + print *, point(3), density_alpha_value_p enddo end