I am looking for a routine (or set of routines) in the Intel MKL that that can replace
dslucs (Incomplete LU BiConjugate Gradient Squared Ax=b Solver) in
Looking at the MKL documentation, I see that MKL has only CG and GMRES solvers, along with ILU preconditioners. It may be that the GMRES solver will be adequate for your purposes.
you can use in site for more information perhaps this help you. http://aam.mathematik.uni-freiburg.de/IAM/Research/projectskr/lin_solver/