When running the program in single threaded mode I get the error message "Parameter 5 was incorrect on entry to DORMQR". The message is not present if I either run the program with multiple threads or if I compile it with reference LAPACK. My compiltervar settings are intel64 ilp64 and the program is compiled with
ifort testgsvd.f90 -lmkl_rt
testgsvd.f90:
program testgsvd
implicit none
integer*8 :: k, l, alpha(10), beta(10), iwork(20), info
double precision :: A(10,20), B(10,20), U, V, Q, work(80)
open(1, file="a.txt", action = "read")
read(1, *) A
open(2, file="b.txt", action = "read")
read(2, *) B
close(1)
close(2)
call dggsvd('N', 'N', 'N', 10_8, 20_8, 10_8, k, l, A, 10_8, B, 10_8, alpha, beta, U, 1_8, V, 1_8, Q, 1_8, work, iwork, info)
end program testgsvd
The matrices A and B are attached.
