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.
Probably best to ask in the MKL forum.
Thanks. I will do that.
For more complete information about compiler optimizations, see our Optimization Notice.