@@ -49,15 +49,24 @@ module test_linalg_inverse
49
49
if (allocated(error)) return
50
50
51
51
call check(error,all(abs(a-inva)<epsilon(0.0_${rk}$)),'inverse_${ri}$_eye (function): data converged')
52
- if (allocated(error)) return
53
-
54
- !> Inverse subroutine
55
- call invert(a,err=state)
52
+ if (allocated(error)) return
53
+
54
+ !> Inverse subroutine: split
55
+ call invert(a,inva, err=state)
56
56
57
57
call check(error,state%ok(),'inverse_${ri}$_eye (subroutine): '//state%print())
58
58
if (allocated(error)) return
59
59
60
60
call check(error,all(abs(a-inva)<epsilon(0.0_${rk}$)),'inverse_${ri}$_eye (subroutine): data converged')
61
+ if (allocated(error)) return
62
+
63
+ !> Inverse subroutine in-place
64
+ call invert(a,err=state)
65
+
66
+ call check(error,state%ok(),'inverse_${ri}$_eye (in-place): '//state%print())
67
+ if (allocated(error)) return
68
+
69
+ call check(error,all(abs(a-inva)<epsilon(0.0_${rk}$)),'inverse_${ri}$_eye (in-place): data converged')
61
70
if (allocated(error)) return
62
71
63
72
end subroutine test_${ri}$_eye_inverse
0 commit comments