Merge pull request #59 from OrderN/hotfix-1.0.3pre

Hotfix for memory leak in Bessel transforms
OrderN · May 18, 2020 · 8dc94a8 · 8dc94a8
2 parents 0450523 + 50f1119
commit 8dc94a8
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/src/ol_bessel_transforms.f90 b/src/ol_bessel_transforms.f90
@@ -383,7 +383,7 @@ subroutine new_bessel_transform_evenFFTW(n,function_in,npts,npts2,rcut,delta_r,f
      use datatypes
      use numbers
      use fft_interface_module, ONLY: cosft_init_wrapper, cosft_exec_wrapper, &
-          sinft_init_wrapper, sinft_exec_wrapper
+          sinft_init_wrapper, sinft_exec_wrapper, cosft_dest_wrapper, sinft_dest_wrapper
 
      implicit none
 
@@ -483,6 +483,8 @@ subroutine new_bessel_transform_evenFFTW(n,function_in,npts,npts2,rcut,delta_r,f
            function_out(1) = function_out(1)+(r*r*delta_r*function_in(i))
         enddo
      end if
+     call cosft_dest_wrapper
+     call sinft_dest_wrapper
    end subroutine new_bessel_transform_evenFFTW
 !!***
 
@@ -511,7 +513,7 @@ subroutine new_bessel_transform_oddFFTW(n,function_in,npts,npts2,rcut,delta_r,fu
      use datatypes
      use numbers
      use fft_interface_module, ONLY: cosft_init_wrapper, cosft_exec_wrapper, &
-          sinft_init_wrapper, sinft_exec_wrapper
+          sinft_init_wrapper, sinft_exec_wrapper, cosft_dest_wrapper, sinft_dest_wrapper
 
      implicit none
 
@@ -602,6 +604,8 @@ subroutine new_bessel_transform_oddFFTW(n,function_in,npts,npts2,rcut,delta_r,fu
         end do
      end do
      function_out(1) = zero
+     call cosft_dest_wrapper
+     call sinft_dest_wrapper
    end subroutine new_bessel_transform_oddFFTW
 
  end module bessel_integrals