[BUG] - The sequence of adding the multiple solvers in the deformed mesh will change the results

NengLu commented 2 weeks ago

Describe the bug

What is the bug? The sequence of adding the multiple solvers in the deformed mesh will change the results. Screenshots mesh._equation_systems_register: SNES_Poisson SNES_Projection SNES_Stokes mesh._equation_systems_register: SNES_Poisson SNES_Stokes SNES_Projection
What version code? latest dev branch, python-3.11
What system are you running the code on? linux/amd64, ubuntu:22.04
Steps to reproduce? Ex_2DRelaxation_FreeSurf_y_stokes_projection.txt Ex_2DRelaxation_FreeSurf_y_projection_stokes.txt

NengLu commented 2 weeks ago

Short scripts for reproducing the bug: Test_solver_Proj_Stokes.txt Test_solver_Stokes_Proj.txt

lmoresi commented 5 days ago

I did find a bug that might have been causing problems with the DM updates. It might be worth trying with the code from my recent PR (branch is joss-setup). I can't really tell if this is fixing anything because the examples are very complicated and I don't really see what to do to reproduce the error / how to know if it is fixed.

NengLu commented 4 days ago

Thanks for that. I tried the latest version code, but the issue remains.

underworldcode / underworld3

[BUG] - The sequence of adding the multiple solvers in the deformed mesh will change the results #205