you may set maxfct ( maximum #of factors with identical sparsity structure) accordingly your case, make symbolic factorizations only once and then do numerical factorization maxfct times. we have a couple of examples show these cases.
Could you please direct me to an example to achieve this?
My problem involve solving many linear systems (to obtain solution to non linear system)
For example, I would like to solve
A1 x1 = b1
A2 x2 = b2
where A2 depends on x1, but the sparsity of A1 and A2 are the same (however I do not have A2 until I solve for x1)