Dear Gerry,
Check please arguments of the function call. I think at least one of them is incorrect.
According to the gerfs interface, arguments should be:
call gerfs( a, af, ipiv, b, x [,trans] [,ferr] [,berr] [,info] )
Where:
a Holds the matrix A of size (n, n).
af Holds the matrix AF of size (n, n).
ipiv Holds the vector of length (n).
b Holds the matrix B of size (n, nrhs).
x Holds the matrix X of size (n, nrhs).
ferr Holds the vector of length (nrhs).
berr Holds the vector of length (nrhs).
trans Must be 'N', 'C', or 'T'. The default value is 'N'.
If nrhs=1, then b,x could be vectors of length n; ferr,berr variables.
Simple test for gerfs with wrong arguments gives the same error:
program test
use mkl95_lapack
use mkl95_precision
implicit none
call gerfs()
pause
end !test
fortcom: Error: source/gbsv.f90, line 5: There is no matching specific subroutine for this generic subroutine call. [GERFS]
call gerfs()
-----^
-Vladimir
