Trouble porting model from 1 to multiple mpi ranks about underworld2 HOT 3 OPEN

Hi Gilly,

Good to hear from you. Nice model!
The problem is to do with the multigrid solver and the indenter (internal BC). I'm working on a fix now. To get around the problem you can use the mumps solver. i.e.

Model.solver.set_inner_method("mumps")

For 2D models this works well and gets around the problem. For 3D solves it DOESN"T work well because mumps is a direct solve method that consumes far too much memory.
Let me know how it goes.

cheers,
Jules

from underworld2.

Thanks Jules! Sorry for the delayed response. I'll give it a go and let you know if I have any issue with this workaround.

Cheers,
Gilly

from underworld2.

Quick feedback on the workaround you pointed out:

• the model starts on multiple ranks with mumps (Model.solver.set_inner_method("mumps")) but it takes forever to converge... I tried adding some penalty too, but it doesn't do any good.
• following your original suggestion Jules I did play around a bit and found that using superludist for the inner solves (Model.solver.set_inner_method("superludist")) is doing a better job with this setup.

Thanks again for your help Jules!

from underworld2.