Hi,
I'm using the Pardiso sparse solver (i.e. with the ia, ja and aa arrays storage structure)for a finite element code I have written. The problem I solve with this ishighly non-linear so the stiffness matrix changes on almost every iteration of the algorithm.
The stiffness matrix (K)is established by assembly of smallerelement stiffness matrices (Ke)to appropriate positions in aa. TheKe matrices areusually fully dense. The way I do this now is by setting up an index array for each element at the start of the algorithm, which I then use to assemble rows of Ke to K (stored in aa) at each iteration.
I suspect that this indexing might be slow, but any other method I have thought of includes a lot of index searching and I suspect that to be even slower (at least if I write it myself).
Does anyone know or have any experience with this, and if it is possible to speed this up?
Cheers!