I am trying to using pardiso (sparse direct solver) as my subdomain (substructure) solver in my domain decomposition application.
The first solve with dense right hand side (rhs) coming from source terms and boundary conditions.
The following sloves using the factorized sparse matrices with sparse right hand side only assotiated with boundary conditions.
With pardiso, the first solve I set IPARM(31) = 0 and PERM(i)=0 with phase = 13.
The following solves, I set IPARM(31) = 1, and only some PERM(i) = 1.
It seems like the entries b( m ) affected my solution even if PERM(m) =0 when phase = 33.
Here are my questions.
1. Is PERM usefull for function pardiso( ) when phase = 33 ?
2. Is there a better way to acheive my goal which is
the first solve (phase = 13) is with dense right hand side, but the following solves (phase = 33) using the factorized matrices with sparse right hand side ?