Skip to content

Commit

Permalink
Fix some more loops
Browse files Browse the repository at this point in the history
(cherry picked from commit c74cefc)
  • Loading branch information
njansson committed May 30, 2024
1 parent 4d03281 commit e91b172
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions src/fluid/bcknd/cpu/pnpn_res_cpu.f90
Original file line number Diff line number Diff line change
Expand Up @@ -92,7 +92,7 @@ subroutine pnpn_prs_res_cpu_compute(p, p_res, u, v, w, u_e, v_e, w_e, f_x, &

call Ax%compute(p_res%x,p%x,c_Xh,p%msh,p%Xh)

do i = 1, n
do concurrent (i = 1:n)
p_res%x(i,1,1,1) = (-p_res%x(i,1,1,1)) &
+ wa1%x(i,1,1,1) + wa2%x(i,1,1,1) + wa3%x(i,1,1,1)
end do
Expand Down Expand Up @@ -161,7 +161,7 @@ subroutine pnpn_vel_res_cpu_compute(Ax, u, v, w, u_res, v_res, w_res, &

call opgrad(ta1%x, ta2%x, ta3%x, p%x, c_Xh)

do i = 1, n
do concurrent (i = 1:n)
u_res%x(i,1,1,1) = (-u_res%x(i,1,1,1)) - ta1%x(i,1,1,1) + f_x%x(i,1,1,1)
v_res%x(i,1,1,1) = (-v_res%x(i,1,1,1)) - ta2%x(i,1,1,1) + f_y%x(i,1,1,1)
w_res%x(i,1,1,1) = (-w_res%x(i,1,1,1)) - ta3%x(i,1,1,1) + f_z%x(i,1,1,1)
Expand Down

0 comments on commit e91b172

Please sign in to comment.