diff --git a/src/GT/GTeh_regularization.f90 b/src/GT/GTeh_regularization.f90
index e76886c..b651609 100644
--- a/src/GT/GTeh_regularization.f90
+++ b/src/GT/GTeh_regularization.f90
@@ -37,14 +37,14 @@ subroutine GTeh_regularization(nBas,nC,nO,nR,nS,e,Om,rhoL,rhoR)
 
       do i=nC+1,nO
         Dpim = e(p) - e(i) + Om(m)
-        kappa  = 1d0 - exp(-Dpim*Dpim*s)
+        kappa  = exp(-Dpim*Dpim*s)
         rhoL(i,p,m) = kappa*rhoL(i,p,m)
         rhoR(i,p,m) = kappa*rhoR(i,p,m)
       enddo
 
       do a=nO+1,nBas-nR
         Dpam = e(p) - e(a) - Om(m)
-        kappa  = 1d0 - exp(-Dpam*Dpam*s)
+        kappa  = exp(-Dpam*Dpam*s)
         rhoL(p,a,m) = kappa*rhoL(p,a,m)
         rhoR(p,a,m) = kappa*rhoR(p,a,m)
       enddo
diff --git a/src/GT/GTpp_regularization.f90 b/src/GT/GTpp_regularization.f90
index 7c2baa9..c4ab89b 100644
--- a/src/GT/GTpp_regularization.f90
+++ b/src/GT/GTpp_regularization.f90
@@ -40,7 +40,7 @@ subroutine GTpp_regularization(nBas,nC,nO,nR,nOO,nVV,e,Om1,rho1,Om2,rho2)
     do m=1,nVV
       do i=nC+1,nO
         Dpim = e(p) + e(i) - Om1(m)
-        kappa  = 1d0 - exp(-Dpim*Dpim*s)
+        kappa  = exp(-Dpim*Dpim*s)
         rho1(p,i,m) = kappa*rho1(p,i,m)
       end do
     end do
@@ -48,7 +48,7 @@ subroutine GTpp_regularization(nBas,nC,nO,nR,nOO,nVV,e,Om1,rho1,Om2,rho2)
     do m=1,nOO
       do a=nO+1,nBas-nR
         Dpam = e(p) + e(a) - Om2(m)
-        kappa  = 1d0 - exp(-Dpam*Dpam*s)
+        kappa  = exp(-Dpam*Dpam*s)
         rho2(p,a,m) = kappa*rho2(p,a,m)
       end do
     end do
diff --git a/src/GW/GW_regularization.f90 b/src/GW/GW_regularization.f90
index 3919db5..cff3ad9 100644
--- a/src/GW/GW_regularization.f90
+++ b/src/GW/GW_regularization.f90
@@ -36,13 +36,13 @@ subroutine GW_regularization(nBas,nC,nO,nR,nS,e,Om,rho)
 
       do i=nC+1,nO
         Dpim = e(p) - e(i) + Om(m)
-        kappa  = 1d0 - exp(-Dpim*Dpim*s)
+        kappa = exp(-Dpim*Dpim*s)
         rho(p,i,m) = kappa*rho(p,i,m)
       enddo
 
       do a=nO+1,nBas-nR
         Dpam = e(p) - e(a) - Om(m)
-        kappa  = 1d0 - exp(-Dpam*Dpam*s)
+        kappa = exp(-Dpam*Dpam*s)
         rho(p,a,m) = kappa*rho(p,a,m)
       enddo