Multigrid Preconditioner Crashes

More
5 years 1 month ago #1952 by dfoiles
Hello,

I am using adaptive refinement with the multigrid preconditioner and I have found that after each iteration, the program will stop without giving any indication of an error. I save the mesh after each iteration, so I can rerun the program and improve the mesh manually one iteration at a time, but I would like to know how to fix this issue. When running the program, I can see that the bilinear and linear forms are being assembled, but the program crashes when it gets to building the preconditioner. I have not had this problem with any other preconditioners. If could you help me figure this out, I would be appreciative.

Thanks
More
5 years 1 month ago #1955 by christopher
Hi,
can you send a minimal example reproducing that behaviour? In the example from the itutorials the multigrid preconditioner is working.
ngsolve.org/docu/latest/i-tutorials/unit...vity/adaptivity.html
Best
Christopher
More
5 years 1 month ago - 5 years 1 month ago #1960 by dfoiles
I have attached an example. I find that it will successfully complete one iteration, but will unexpectedly stop as soon as assembly is completed during the second iteration.

Thanks

File Attachment:

File Name: LSPR_minim...10-15.py
File Size:4 KB
Last edit: 5 years 1 month ago by dfoiles.
More
5 years 1 month ago #1971 by christopher
I see, thanks for reporting this, there seems to be some needed functionality for complex geometric multigrid missing for the HCurl space. We will investigate this. For the time being you could use the bddcc preconditioner, it needs more iterations but every iteration is faster, so for me it's about the same speed as the multigrid.
I'll let you know when the multigrid is fixed.
Best
Christopher
Time to create page: 0.104 seconds