add test: init NaNs, check correct return

Author: perazz

test/linalg/test_linalg_qr.fypp

@@ -53,22 +53,20 @@ module test_linalg_qr
         #:endif
         aorig = a
 
-        ! 1) QR factorization with full matrices
-        q = 1.0_${rk}$
+        ! 1) QR factorization with full matrices. Input NaNs to be sure Q and R are OK on return
+        q = 0.0_${rk}$
+        q = 1.0_${rk}$/q
+        r = 0.0_${rk}$
+        r = 1.0_${rk}$/r       
         call qr(a,q,r,err=state)
 
         ! Check return code
         call check(error,state%ok(),state%print())
         if (allocated(error)) return        
 
         ! Check solution
-        if (.not.all(abs(a-matmul(q,r))<tol)) then 
-            state = linalg_state_type('qr',LINALG_VALUE_ERROR,'converged solution (full) max error: ',&
-                                      maxval(abs(a-matmul(q,r))))
-            call check(error, state%ok(), state%print())
-            if (allocated(error)) return     
-        endif
-             
+        call check(error, all(abs(a-matmul(q,r))<tol), 'converged solution (fulle)')
+        if (allocated(error)) return        
 
         ! 2) QR factorization with reduced matrices
         call qr(a,qred,rred,err=state)