lorenzo-mechbau
commented
5 years ago Branch renamed to PR_develop_backsubstitute
Closed lorenzo-mechbau closed 5 years ago
lorenzo-mechbau
commented
5 years ago Branch renamed to PR_develop_backsubstitute
A change has been implemented in the file equations_set_routines.F90 to correctly perform backsubstitute (in EQUATIONS_SET_BACKSUBSTITUTE) in the case of sparse matrices (compressed row storage). Backsubstitute performs a matrix-vector multiplication e.g. when a flux is computed based on the solution of a Laplace problem. The previous implementation was correct for full matrices (block storage), but led to wrong results in sequential/parallel simulations with sparse equation matrices.
The fix involves global-to-local maps which should be avoided/replaced in the next future to avoid replicated data.