diff options
Diffstat (limited to 'src/Dissipation_6_5_min_err_coeff.F90')
-rw-r--r-- | src/Dissipation_6_5_min_err_coeff.F90 | 32 |
1 files changed, 13 insertions, 19 deletions
diff --git a/src/Dissipation_6_5_min_err_coeff.F90 b/src/Dissipation_6_5_min_err_coeff.F90 index 623efad..4fe6120 100644 --- a/src/Dissipation_6_5_min_err_coeff.F90 +++ b/src/Dissipation_6_5_min_err_coeff.F90 @@ -69,24 +69,20 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & call set_dmatrix_half ( d(1:ir,1:ir), 1 ) - tmp(1:ir,1:ir) = -transpose ( & - matmul ( h(1:ir,1:ir), & - matmul ( & - transpose ( d(1:ir,1:ir) ), & - matmul ( b(1:ir,1:ir), & - d(1:ir,1:ir) ) ) ) ) + tmp(1:ir,1:ir) = matmul ( b(1:ir,1:ir), d(1:ir,1:ir) ) + tmp(1:ir,1:ir) = matmul ( transpose ( d(1:ir,1:ir) ), tmp(1:ir,1:ir) ) + tmp(1:ir,1:ir) = matmul ( h(1:ir,1:ir), tmp(1:ir,1:ir) ) + tmp(1:ir,1:ir) = -transpose ( tmp(1:ir,1:ir) ) atmp(1:ir,1:center) = tmp(1:ir,1:center) il = center + 1 - max(gsize(1),3) d = zero call set_dmatrix_half ( d(il:ni,il:ni), 2 ) - tmp(il:ni,il:ni) = -transpose ( & - matmul ( h(il:ni,il:ni), & - matmul ( & - transpose ( d(il:ni,il:ni) ), & - matmul ( b(il:ni,il:ni), & - d(il:ni,il:ni) ) ) ) ) + tmp(il:ni,il:ni) = matmul ( b(il:ni,il:ni), d(il:ni,il:ni) ) + tmp(il:ni,il:ni) = matmul ( transpose ( d(il:ni,il:ni) ), tmp(il:ni,il:ni) ) + tmp(il:ni,il:ni) = matmul ( h(il:ni,il:ni), tmp(il:ni,il:ni) ) + tmp(il:ni,il:ni) = -transpose ( tmp(il:ni,il:ni) ) atmp(il:ni,center+1:ni) = tmp(il:ni,center+1:ni) @@ -281,12 +277,10 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & call set_dmatrix_half ( d(1:jr,1:jr), 1 ) - tmp(1:jr,1:jr) = -transpose ( & - matmul ( h(1:jr,1:jr), & - matmul ( & - transpose ( d(1:jr,1:jr) ), & - matmul ( b(1:jr,1:jr), & - d(1:jr,1:jr) ) ) ) ) + tmp(1:jr,1:jr) = matmul ( b(1:jr,1:jr), d(1:jr,1:jr) ) + tmp(1:jr,1:jr) = matmul ( transpose ( d(1:jr,1:jr) ), tmp(1:jr,1:jr) ) + tmp(1:jr,1:jr) = matmul ( h(1:jr,1:jr), tmp(1:jr,1:jr) ) + tmp(1:jr,1:jr) = -transpose ( tmp(1:jr,1:jr) ) atmp(1:jr,1:center) = tmp(1:jr,1:center) @@ -744,7 +738,7 @@ contains implicit none CCTK_REAL, dimension(:,:), intent(out) :: d - CCTK_INT, intent(in) :: part + integer, intent(in) :: part CCTK_INT :: n CCTK_REAL, dimension(4), save :: ac = (/ -1.0_wp, 3.0_wp, -3.0_wp, 1.0_wp /) CCTK_INT :: i |