Skip to content

Commit 406e170

Browse files
committed
typo: set complex NaN
1 parent 40062c4 commit 406e170

File tree

1 file changed

+8
-1
lines changed

1 file changed

+8
-1
lines changed

src/stdlib_linalg_inverse.fypp

Lines changed: 8 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -163,7 +163,14 @@ submodule (stdlib_linalg) stdlib_linalg_inverse
163163
! Return NaN on issues
164164
if (err%error()) then
165165
if (allocated(inva)) deallocate(inva)
166-
allocate(inva(size(a,1,kind=ilp),size(a,2,kind=ilp)),source=ieee_value(1.0_${rk}$,ieee_quiet_nan))
166+
allocate(inva(size(a,1,kind=ilp),size(a,2,kind=ilp)))
167+
168+
#:if rt.startswith('complex')
169+
inva = ieee_value(1.0_${rk}$,ieee_quiet_nan)
170+
#:else
171+
inva = cmplx(ieee_value(1.0_${rk}$,ieee_quiet_nan), &
172+
ieee_value(1.0_${rk}$,ieee_quiet_nan), kind=${rk}$)
173+
#:endif
167174
endif
168175

169176
end function stdlib_linalg_inverse_${ri}$_operator

0 commit comments

Comments
 (0)