i am a newbie in using the MKL.
If I build up a stiffness matrix on a FEM problem the stiffness matrix is singular. To solve a problem - no matter which one - you have to make the siffness matrix nonsingular/regular by inserting a boundary condition on the equation A x = b.
Does somebody of you know if there is a pre-written *.f or *.f90 routine for inserting a boundary condition in a CSR format of the BLAS/PARDISO package?
Do I have to write it by my own? :-)