Skip to content

Commit

Permalink
Rep. %re & %im with real & aimag for nvfortran bug
Browse files Browse the repository at this point in the history
  • Loading branch information
gha3mi committed Jan 3, 2024
1 parent 6a59886 commit 2cb5aee
Showing 1 changed file with 12 additions and 12 deletions.
24 changes: 12 additions & 12 deletions src/forunittest.f90
Original file line number Diff line number Diff line change
Expand Up @@ -335,16 +335,16 @@ subroutine unit_test_cx1(this, res, expected, tol, msg)
this%msg = 'forunittest'
end if

if (norm2(expected%re)<tiny(0.0_rk)) then
rel_err_re = norm2(res%re-expected%re)
if (norm2(real(expected,rk))<tiny(0.0_rk)) then
rel_err_re = norm2(real(res,rk)-real(expected,rk))
else
rel_err_re = norm2(res%re-expected%re)/norm2(expected%re)
rel_err_re = norm2(real(res,rk)-real(expected,rk))/norm2(real(expected,rk))
end if

if (norm2(expected%im)<tiny(0.0_rk)) then
rel_err_im = norm2(res%im-expected%im)
if (norm2(aimag(expected))<tiny(0.0_rk)) then
rel_err_im = norm2(aimag(res)-aimag(expected))
else
rel_err_im = norm2(res%im-expected%im)/norm2(expected%im)
rel_err_im = norm2(aimag(res)-aimag(expected))/norm2(aimag(expected))
end if

if (present(tol)) then
Expand Down Expand Up @@ -374,16 +374,16 @@ subroutine unit_test_cx2(this, res, expected, tol, msg)
this%msg = 'forunittest'
end if

if (norm2(expected%re)<tiny(0.0_rk)) then
rel_err_re = norm2(res%re-expected%re)
if (norm2(real(expected,rk))<tiny(0.0_rk)) then
rel_err_re = norm2(real(res,rk)-real(expected,rk))
else
rel_err_re = norm2(res%re-expected%re)/norm2(expected%re)
rel_err_re = norm2(real(res,rk)-real(expected,rk))/norm2(real(expected,rk))
end if

if (norm2(expected%im)<tiny(0.0_rk)) then
rel_err_im = norm2(res%im-expected%im)
if (norm2(aimag(expected))<tiny(0.0_rk)) then
rel_err_im = norm2(aimag(res)-aimag(expected))
else
rel_err_im = norm2(res%im-expected%im)/norm2(expected%im)
rel_err_im = norm2(aimag(res)-aimag(expected))/norm2(aimag(expected))
end if

if (present(tol)) then
Expand Down

0 comments on commit 2cb5aee

Please sign in to comment.