search

chaos-polymtl / lethe

Repository for the open-source lethe CFD/DEM/CFD-DEM project
https://chaos-polymtl.github.io/lethe/index.html
GNU Lesser General Public License v2.1
267 stars 59 forks source link

Fix NavierStokesOperatorBase::get_system_matrix() #1111

Closed peterrum closed 3 months ago

peterrum commented 3 months ago

At the moment, the matrix is cleared during reinit(), and only setup once although the parameters have changed.

blaisb commented 3 months ago

This will change the steady-state results significantly. @lpsaavedra I think some stuff will need to be rerun. I iwll merge after CI

lpsaavedra commented 3 months ago

I have added a small correction to avoid re-calculating the system matrix of the operator when we use AMG in the coarse-grid solver... after CI, we can merge this