Add the missing lead_dim argument for a Lapack routine.

parent e0511a41
Pipeline #146437162 passed with stage
in 0 seconds
......@@ -732,8 +732,8 @@ subroutine dlinsyssolve(n, nrhs, a, b, x)
SAFE_ALLOCATE(c(1:n))
SAFE_ALLOCATE(af(1:n, 1:n))
call X(gesvx) ("N", "N", n, nrhs, a(1, 1), n, af(1, 1), n, ipiv(1), equed, r(1), c(1), b(1, 1), n, x(1, 1), n, &
rcond, ferr(1), berr(1), work(1), iwork(1), info)
call X(gesvx) ("N", "N", n, nrhs, a(1, 1), lead_dim(a), af(1, 1), n, ipiv(1), equed, r(1), c(1), &
b(1, 1), lead_dim(b), x(1, 1), lead_dim(x), rcond, ferr(1), berr(1), work(1), iwork(1), info)
if(info /= 0) then
write(message(1), '(3a, i5)') 'In dlinsyssolve, LAPACK ', TOSTRING(X(gesvx)), ' returned info = ', info
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment