Is MKL/gesv an appropriate replacement for IMSL/lin_sol_gen? I believe both do an LU decomposition to get the solution. If I merely replace CALL LIN_SOL_GEN(A,B,X) with CALL GESVX(A,B,X), how close should I reasonably expect the solutions to match for double precision problems with, say 100x100, coefficient matrix.
↧