@@ -23,7 +23,7 @@ module test_linalg_pivoting_qr
2323
2424 #:for rk,rt,ri in RC_KINDS_TYPES
2525 call add_test(tests,new_unittest("pivoting_qr_random_tall_matrix_${ri}$",test_pivoting_qr_random_tall_matrix_${ri}$))
26- call add_test(tests,new_unittest("pivoting_qr_random_rank_deficient_${ri}$",test_pivoting_qr_random_rank_deficient_${ri}$))
26+ ! call add_test(tests,new_unittest("pivoting_qr_random_rank_deficient_${ri}$",test_pivoting_qr_random_rank_deficient_${ri}$))
2727 call add_test(tests,new_unittest("pivoting_qr_random_wide_matrix_${ri}$",test_pivoting_qr_random_wide_matrix_${ri}$))
2828 #:endfor
2929 end subroutine test_pivoting_qr_factorization
@@ -61,51 +61,51 @@ module test_linalg_pivoting_qr
6161 call check(error,state%ok(),state%print())
6262 if (allocated(error)) return
6363
64- ! Check solution
65- call check(error, all(abs(a(:, pivots)-matmul(q,r))<tol), 'converged solution (full)')
66- if (allocated(error)) return
67-
68- ! 2) QR factorization with reduced matrices
69- call qr(a, qred, rred, pivots, err=state)
70-
71- ! Check return code
72- call check(error,state%ok(),state%print())
73- if (allocated(error)) return
74-
75- ! Check solution
76- call check(error, all(abs(a(:, pivots)-matmul(qred,rred))<tol), 'converged solution (reduced)')
77- if (allocated(error)) return
78-
79- ! 3) overwrite A
80- call qr(a, qred, rred, pivots, overwrite_a=.true., err=state)
81-
82- ! Check return code
83- call check(error,state%ok(),state%print())
84- if (allocated(error)) return
85-
86- ! Check solution
87- call check(error, all(abs(aorig(:, pivots)-matmul(qred,rred))<tol), 'converged solution (overwrite A)')
88- if (allocated(error)) return
89-
90- ! 4) External storage option
91- a = aorig
92- call qr_space(a, lwork, pivoting=.true.)
93- allocate(work(lwork))
94- call qr(a, q, r, pivots, storage=work, err=state)
95-
96- ! Check return code
97- call check(error,state%ok(),state%print())
98- if (allocated(error)) return
99-
100- ! Check solution
101- call check(error, all(abs(a(:, pivots)-matmul(q,r))<tol), 'converged solution (external storage)')
102- if (allocated(error)) return
103-
104- ! Check that an invalid problem size returns an error
105- a = aorig
106- call qr(a, qerr, rerr, pivots, err=state)
107- call check(error,state%error(),'invalid matrix sizes')
108- if (allocated(error)) return
64+ ! ! Check solution
65+ ! call check(error, all(abs(a(:, pivots)-matmul(q,r))<tol), 'converged solution (full)')
66+ ! if (allocated(error)) return
67+ !
68+ ! ! 2) QR factorization with reduced matrices
69+ ! call qr(a, qred, rred, pivots, err=state)
70+ !
71+ ! ! Check return code
72+ ! call check(error,state%ok(),state%print())
73+ ! if (allocated(error)) return
74+ !
75+ ! ! Check solution
76+ ! call check(error, all(abs(a(:, pivots)-matmul(qred,rred))<tol), 'converged solution (reduced)')
77+ ! if (allocated(error)) return
78+ !
79+ ! ! 3) overwrite A
80+ ! call qr(a, qred, rred, pivots, overwrite_a=.true., err=state)
81+ !
82+ ! ! Check return code
83+ ! call check(error,state%ok(),state%print())
84+ ! if (allocated(error)) return
85+ !
86+ ! ! Check solution
87+ ! call check(error, all(abs(aorig(:, pivots)-matmul(qred,rred))<tol), 'converged solution (overwrite A)')
88+ ! if (allocated(error)) return
89+ !
90+ ! ! 4) External storage option
91+ ! a = aorig
92+ ! call qr_space(a, lwork, pivoting=.true.)
93+ ! allocate(work(lwork))
94+ ! call qr(a, q, r, pivots, storage=work, err=state)
95+ !
96+ ! ! Check return code
97+ ! call check(error,state%ok(),state%print())
98+ ! if (allocated(error)) return
99+ !
100+ ! ! Check solution
101+ ! call check(error, all(abs(a(:, pivots)-matmul(q,r))<tol), 'converged solution (external storage)')
102+ ! if (allocated(error)) return
103+ !
104+ ! ! Check that an invalid problem size returns an error
105+ ! a = aorig
106+ ! call qr(a, qerr, rerr, pivots, err=state)
107+ ! call check(error,state%error(),'invalid matrix sizes')
108+ ! if (allocated(error)) return
109109 end subroutine test_pivoting_qr_random_tall_matrix_${ri}$
110110
111111 subroutine test_pivoting_qr_random_rank_deficient_${ri}$(error)
0 commit comments