Added PRIVATE/PUBLIC attribute explicitly.

Author: h-m0r1

sparse_ir.f90

@@ -4,6 +4,8 @@ MODULE sparse_ir
   !
   IMPLICIT NONE
   !
+  PRIVATE
+  !
   INTEGER, PARAMETER :: DP = KIND(0d0)
   REAL(KIND = DP), PARAMETER :: one = 1.0D0
   REAL(KIND = DP), PARAMETER :: zero = 0.0D0
@@ -57,6 +59,13 @@ MODULE sparse_ir
     MODULE PROCEDURE evaluate_matsubara_b_from_dlr_zz, evaluate_matsubara_b_from_dlr_dz
   END INTERFACE evaluate_matsubara_b_from_dlr
   !
+  PUBLIC :: DecomposedMatrix_z, DecomposedMatrix_d, IR
+  PUBLIC :: init_ir, set_beta, finalize_ir
+  PUBLIC :: evaluate_tau, evaluate_matsubara_f, evaluate_matsubara_b
+  PUBLIC :: fit_tau, fit_matsubara_f, fit_matsubara_b
+  PUBLIC :: to_dlr
+  PUBLIC :: evaluate_tau_from_dlr, evaluate_matsubara_f_from_dlr, evaluate_matsubara_b_from_dlr
+  !
   !-----------------------------------------------------------------------
   TYPE DecomposedMatrix_z
   !-----------------------------------------------------------------------
@@ -393,6 +402,7 @@ SUBROUTINE finalize_ir(obj)
   !
   TYPE(IR) :: obj
   !! contains all the IR-basis objects
+  !
   IF (ALLOCATED(obj%x)) DEALLOCATE(obj%x)
   IF (ALLOCATED(obj%tau)) DEALLOCATE(obj%tau)
   IF (ALLOCATED(obj%s)) DEALLOCATE(obj%s)
@@ -430,7 +440,7 @@ SUBROUTINE finalize_dmat_z(dmat)
   IF (ALLOCATED(dmat%a_odd)) DEALLOCATE(dmat%a_odd)
   IF (ALLOCATED(dmat%a_even)) DEALLOCATE(dmat%a_even)
   IF (ALLOCATED(dmat%inv_s)) DEALLOCATE(dmat%inv_s)
-  IF (ALLOCATED(dmat%inv_s)) DEALLOCATE(dmat%inv_s_dl)
+  IF (ALLOCATED(dmat%inv_s_dl)) DEALLOCATE(dmat%inv_s_dl)
   IF (ALLOCATED(dmat%ut)) DEALLOCATE(dmat%ut)
   IF (ALLOCATED(dmat%v)) DEALLOCATE(dmat%v)
   IF (ALLOCATED(dmat%ut_real)) DEALLOCATE(dmat%ut_real)
@@ -455,7 +465,7 @@ SUBROUTINE finalize_dmat_d(dmat)
   IF (ALLOCATED(dmat%a)) DEALLOCATE(dmat%a)
   IF (ALLOCATED(dmat%a_real)) DEALLOCATE(dmat%a_real)
   IF (ALLOCATED(dmat%inv_s)) DEALLOCATE(dmat%inv_s)
-  IF (ALLOCATED(dmat%inv_s)) DEALLOCATE(dmat%inv_s_dl)
+  IF (ALLOCATED(dmat%inv_s_dl)) DEALLOCATE(dmat%inv_s_dl)
   IF (ALLOCATED(dmat%ut)) DEALLOCATE(dmat%ut)
   IF (ALLOCATED(dmat%v)) DEALLOCATE(dmat%v)
   IF (ALLOCATED(dmat%ut_real)) DEALLOCATE(dmat%ut_real)
@@ -1271,6 +1281,7 @@ SUBROUTINE evaluate_matsubara_f_zz(obj, arr, res)
   ENDIF
   DEALLOCATE(res_r)
   DEALLOCATE(res_i)
+  DEALLOCATE(arr_half)
   !
   !-----------------------------------------------------------------------
   END SUBROUTINE evaluate_matsubara_f_zz
@@ -1355,6 +1366,7 @@ SUBROUTINE evaluate_matsubara_f_dz(obj, arr, res)
   res = CMPLX(res_r, res_i, KIND = DP)
   DEALLOCATE(res_r)
   DEALLOCATE(res_i)
+  DEALLOCATE(arr_half)
   !
   !-----------------------------------------------------------------------
   END SUBROUTINE evaluate_matsubara_f_dz
@@ -1472,6 +1484,7 @@ SUBROUTINE evaluate_matsubara_b_zz(obj, arr, res)
   ENDIF
   DEALLOCATE(res_r)
   DEALLOCATE(res_i)
+  DEALLOCATE(arr_half)
   !
   !-----------------------------------------------------------------------
   END SUBROUTINE evaluate_matsubara_b_zz
@@ -1556,6 +1569,7 @@ SUBROUTINE evaluate_matsubara_b_dz(obj, arr, res)
   res = CMPLX(res_r, res_i, KIND = DP)
   DEALLOCATE(res_r)
   DEALLOCATE(res_i)
+  DEALLOCATE(arr_half)
   !
   !-----------------------------------------------------------------------
   END SUBROUTINE evaluate_matsubara_b_dz
@@ -1896,7 +1910,7 @@ SUBROUTINE fit_tau_zd(obj, arr, res)
   DEALLOCATE(ut_arr, arr_tmp)
   !
   !-----------------------------------------------------------------------
-  END SUBROUTINE
+  END SUBROUTINE fit_tau_zd
   !-----------------------------------------------------------------------
   !
   !-----------------------------------------------------------------------
@@ -1989,7 +2003,7 @@ SUBROUTINE fit_tau_dd(obj, arr, res)
   DEALLOCATE(ut_arr)
   !
   !-----------------------------------------------------------------------
-  END SUBROUTINE
+  END SUBROUTINE fit_tau_dd
   !-----------------------------------------------------------------------
   !
   !-----------------------------------------------------------------------

sparse_ir_io.f90

@@ -6,6 +6,12 @@ MODULE sparse_ir_io
   !
   IMPLICIT NONE
   !
+  PRIVATE
+  !
+  INTEGER, PARAMETER :: DP = KIND(0d0)
+  !
+  PUBLIC :: read_ir
+  !
   CONTAINS
   !
   !-----------------------------------------------------------------------
@@ -46,7 +52,6 @@ FUNCTION read_ir(unit, beta, positive_only) RESULT(obj)
   END FUNCTION read_ir
   !-----------------------------------------------------------------------
   !
-  !
   !-----------------------------------------------------------------------
   FUNCTION read_v1(unit, beta, positive_only) RESULT(obj)
   !-----------------------------------------------------------------------
@@ -207,7 +212,7 @@ FUNCTION read_v1(unit, beta, positive_only) RESULT(obj)
     CALL init_ir(obj, beta, lambda, eps, s, tau, freq_f, freq_b, u, uhat_f, uhat_b, omega, v, dlr, 1d-16, positive_only)
   ENDIF
   !
-  DEALLOCATE(u, uhat_f, uhat_b, v, dlr)
+  DEALLOCATE(s, tau, u, freq_f, uhat_f, freq_b, uhat_b, omega, v, dlr)
   !
   !-----------------------------------------------------------------------
   END FUNCTION read_v1

test_sparse_ir.f90

@@ -7,6 +7,14 @@ PROGRAM test
   !
   IMPLICIT NONE
   !
+  INTEGER, PARAMETER :: DP = KIND(0d0)
+  REAL(KIND = DP), PARAMETER :: one = 1.0D0
+  REAL(KIND = DP), PARAMETER :: zero = 0.0D0
+  REAL(KIND = DP), PARAMETER :: pi = 4.D0*ATAN(1.D0)
+  COMPLEX(KIND = DP), PARAMETER :: cone = (1.0D0, 0.0D0)
+  COMPLEX(KIND = DP), PARAMETER :: ci = (0.0D0, 1.0D0)
+  COMPLEX(KIND = DP), PARAMETER :: czero = (0.0D0, 0.0D0)
+  !
   ! use preset
   CALL test_fermion(.true., .false., .false., .false.)
   CALL test_boson  (.true., .false., .false., .false.)