delete laplacian function definition

Author: rouson

src/matcha/subdomain_s.F90

@@ -80,70 +80,6 @@
     my_dz = dz_
   end procedure
 
-  module procedure laplacian
-    integer i, j, k
-    real, allocatable :: halo_west(:,:), halo_east(:,:)
-
-    call_assert(allocated(rhs%s_))
-    call_assert(allocated(halo_x))
-
-    allocate(laplacian_rhs%s_(my_nx, ny, nz))
-    halo_west = merge(halo_x(west,:,:), rhs%s_(1,:,:), me/=1) ! conditionally use halo value
-    i = my_internal_west
-    call_assert_describe(i+1<=my_nx, "laplacian: westernmost subdomain too small")
-    ! Compute Laplacians throughout the low-x boundary subdomain using non-allocatable associations:
-    associate( laplacian_phi =>  laplacian_rhs%s_, inbox => halo_west, phi=>rhs%s_)
-#if HAVE_2018_LOCALITY_SPECIFIERS
-      do concurrent(j=2:ny-1, k=2:nz-1) & 
-        default(none) shared(laplacian_phi, inbox, phi, dx_, dy_, dz_, i) !Fortran 2018 loacality specifiers
-#else
-      do concurrent(j=2:ny-1, k=2:nz-1)
-#endif
-        laplacian_phi(i,j,k) = (inbox(j,k    ) - 2*phi(i,j,k) + phi(i+1,j  ,k  ))/dx_**2 + &
-                               (phi(i,j-1,k  ) - 2*phi(i,j,k) + phi(i  ,j+1,k  ))/dy_**2 + &
-                               (phi(i,j  ,k-1) - 2*phi(i,j,k) + phi(i  ,j  ,k+1))/dz_**2
-      end do
-    end associate
-    ! Compute Laplacians throughout non-boundary subdomains with non-allocatable associations:
-    associate(laplacian_phi => laplacian_rhs%s_, phi=>rhs%s_)
-#if HAVE_2018_LOCALITY_SPECIFIERS
-      do concurrent(i=my_internal_west+1:my_internal_east-1, j=2:ny-1, k=2:nz-1) &
-        default(none) shared(laplacian_phi, phi, dx_, dy_, dz_) ! Fortran 2018 locality specifiers
-#else
-      do concurrent(i=my_internal_west+1:my_internal_east-1, j=2:ny-1, k=2:nz-1)
-#endif
-        laplacian_phi(i,j,k) = (phi(i-1,j  ,k  ) - 2*phi(i,j,k) + phi(i+1,j  ,k  ))/dx_**2 + &
-                               (phi(i  ,j-1,k  ) - 2*phi(i,j,k) + phi(i  ,j+1,k  ))/dy_**2 + &
-                               (phi(i  ,j  ,k-1) - 2*phi(i,j,k) + phi(i  ,j  ,k+1))/dz_**2
-      end do
-    end associate
-
-    halo_east = merge(halo_x(east,:,:), rhs%s_(my_nx,:,:), me/=num_subdomains) !conditionally use halo value
-    i = my_internal_east
-    call_assert_describe(i-1>0, "laplacian: easternmost subdomain too small")
-    ! Compute Laplacians throughout the high-x boundary subdomain using non-allocatable associations:
-    associate(laplacian_phi =>  laplacian_rhs%s_, inbox => halo_east, phi=>rhs%s_)
-#if HAVE_2018_LOCALITY_SPECIFIERS
-      do concurrent(j=2:ny-1, k=2:nz-1) & ! compute Laplacian in low-x boundary subdomain
-        default(none) shared(laplacian_phi, inbox, phi, dx_, dy_, dz_, i) ! Fortran 2018 locality specifiers
-#else
-      do concurrent(j=2:ny-1, k=2:nz-1) ! compute Laplacian in low-x boundary subdomain
-#endif
-        laplacian_phi(i,j,k) = (phi(i-1,j  ,k  ) - 2*phi(i,j,k) + inbox(  j  ,k  ))/dx_**2 + &
-                               (phi(i  ,j-1,k  ) - 2*phi(i,j,k) + phi(i  ,j+1,k  ))/dy_**2 + &
-                               (phi(i  ,j  ,k-1) - 2*phi(i,j,k) + phi(i  ,j  ,k+1))/dz_**2
-      end do
-    end associate
-
-    laplacian_rhs%s_(:, 1,:) = 0. ! low-y boundary
-    laplacian_rhs%s_(:,ny,:) = 0. ! high-y boundary
-    laplacian_rhs%s_(:,:, 1) = 0. ! low-z boundary
-    laplacian_rhs%s_(:,:,nz) = 0. ! high-z boundary
-
-    if (me==1) laplacian_rhs%s_(1,:,:) = 0. ! low-x boundary
-    if (me==num_subdomains) laplacian_rhs%s_(my_nx,:,:) = 0. ! high-x boundary
-  end procedure
-
   module procedure multiply
     call_assert(allocated(rhs%s_))
     product%s_ =  lhs * rhs%s_