@@ -26,7 +26,7 @@ module test_linalg_svd
2626 #:endfor
2727
2828 #:for ck,ct,ci in CMPLX_KINDS_TYPES
29- #:if not ck in [ "xdp","qp"]
29+ #:if ck!= "xdp"
3030 tests = [tests,new_unittest("test_complex_svd_${ci}$",test_complex_svd_${ci}$)]
3131 #:endif
3232 #:endfor
@@ -180,18 +180,20 @@ module test_linalg_svd
180180 type(error_type), allocatable, intent(out) :: error
181181
182182 !> Reference solution
183- real(${ck}$), parameter :: tol = sqrt(epsilon(0.0_${ck}$))
184- real(${ck}$), parameter :: one = 1.0_${ck}$
185- real(${ck}$), parameter :: zero = 0.0_${ck}$
186- real(${ck}$), parameter :: sqrt2 = sqrt(2.0_${ck}$)
187- real(${ck}$), parameter :: rsqrt2 = one/sqrt2
188- ${ct}$, parameter :: cone = (1.0_${ck}$,0.0_${ck}$)
189- ${ct}$, parameter :: cimg = (0.0_${ck}$,1.0_${ck}$)
190- ${ct}$, parameter :: czero = (0.0_${ck}$,0.0_${ck}$)
183+ real(${ck}$), parameter :: tol = sqrt(epsilon(0.0_${ck}$))
184+ real(${ck}$), parameter :: one = 1.0_${ck}$
185+ real(${ck}$), parameter :: zero = 0.0_${ck}$
186+ real(${ck}$), parameter :: sqrt2 = sqrt(2.0_${ck}$)
187+ real(${ck}$), parameter :: rsqrt2 = one/sqrt2
188+ ${ct}$, parameter :: csqrt2 = (rsqrt2,zero)
189+ ${ct}$, parameter :: isqrt2 = (zero,rsqrt2)
190+ ${ct}$, parameter :: cone = (1.0_${ck}$,0.0_${ck}$)
191+ ${ct}$, parameter :: cimg = (0.0_${ck}$,1.0_${ck}$)
192+ ${ct}$, parameter :: czero = (0.0_${ck}$,0.0_${ck}$)
191193
192194 real(${ck}$), parameter :: s_sol(2) = [sqrt2,sqrt2]
193195 ${ct}$, parameter :: A_mat(2,2) = reshape([cone,cimg,cimg,cone],[2,2])
194- ${ct}$, parameter :: u_sol(2,2) = reshape(rsqrt2*[cone,cimg,cimg,cone ],[2,2])
196+ ${ct}$, parameter :: u_sol(2,2) = reshape([csqrt2,isqrt2,isqrt2,csqrt2 ],[2,2])
195197 ${ct}$, parameter :: vt_sol(2,2) = reshape([cone,czero,czero,cone],[2,2])
196198
197199 !> Local variables
0 commit comments