restore all, avoid real128 complex parameter math

Author: perazz

test/linalg/test_linalg_svd.fypp

@@ -26,7 +26,7 @@ module test_linalg_svd
         #:endfor
 
         #:for ck,ct,ci in CMPLX_KINDS_TYPES
-        #:if not ck in ["xdp","qp"]
+        #:if ck!="xdp"
         tests = [tests,new_unittest("test_complex_svd_${ci}$",test_complex_svd_${ci}$)]        
         #:endif
         #:endfor
@@ -180,18 +180,20 @@ module test_linalg_svd
         type(error_type), allocatable, intent(out) :: error
 
         !> Reference solution
-        real(${ck}$), parameter :: tol     = sqrt(epsilon(0.0_${ck}$))
-        real(${ck}$), parameter :: one     = 1.0_${ck}$
-        real(${ck}$), parameter :: zero    = 0.0_${ck}$
-        real(${ck}$), parameter ::  sqrt2  = sqrt(2.0_${ck}$)
-        real(${ck}$), parameter :: rsqrt2  = one/sqrt2
-        ${ct}$, parameter :: cone  = (1.0_${ck}$,0.0_${ck}$)
-        ${ct}$, parameter :: cimg  = (0.0_${ck}$,1.0_${ck}$)
-        ${ct}$, parameter :: czero = (0.0_${ck}$,0.0_${ck}$)
+        real(${ck}$), parameter :: tol    = sqrt(epsilon(0.0_${ck}$))
+        real(${ck}$), parameter :: one    = 1.0_${ck}$
+        real(${ck}$), parameter :: zero   = 0.0_${ck}$
+        real(${ck}$), parameter :: sqrt2  = sqrt(2.0_${ck}$)
+        ${ct}$, parameter :: rsqrt2 = (1.0_${ck}$/sqrt(2.0_${ck}$),0.0_${ck}$)
+        ${ct}$, parameter :: isqrt2 = (0.0_${ck}$,1.0_${ck}$/sqrt(2.0_${ck}$))
+        ${ct}$, parameter :: cone   = (1.0_${ck}$,0.0_${ck}$)
+        ${ct}$, parameter :: cimg   = (0.0_${ck}$,1.0_${ck}$)
+        ${ct}$, parameter :: czero  = (0.0_${ck}$,0.0_${ck}$)
+        ${ct}$, parameter :: czero  = (0.0_${ck}$,0.0_${ck}$)
 
         real(${ck}$), parameter ::  s_sol(2)   = [sqrt2,sqrt2]
         ${ct}$, parameter ::  A_mat(2,2) = reshape([cone,cimg,cimg,cone],[2,2])
-        ${ct}$, parameter ::  u_sol(2,2) = reshape(rsqrt2*[cone,cimg,cimg,cone],[2,2])
+        ${ct}$, parameter ::  u_sol(2,2) = reshape([rsqrt2,isqrt2,isqrt2,rsqrt2],[2,2])
         ${ct}$, parameter :: vt_sol(2,2) = reshape([cone,czero,czero,cone],[2,2])
 
         !> Local variables