about discretization limited

[xibaiw]

Hello,
I want to confirm if PeriLab has any limitations on mesh discretization. Is PeriLab only capable of simulating up to 100,000 points, or does the computable model size depend on the PC's performance? Is there has any estimation of the relationship between performance and model size?

[JTHesse]

Hi @xibaiw, there is no hard-coded limit to the number of nodes. Only the number itself as an Int64 can't be larger than 2^63-1.
But as you pointed out correctly the bottleneck here would be your hardware.
Take a look at Figure 6 of our initial paper, there you can see a few examples of performance, but keep in mind that this is an older PeriLab Version.
If you run into performance issues let us know.

[xibaiw]

Based on the above situation where computations work with low-density meshes but fail with refined meshes, I'd like to inquire whether PeriLab has any inherent limitations on the number of nodes it can handle.

[JTHesse]

It sounds to me like the RAM is not big enough for your example. However, julia should let you know why it is aborting, can you share the log?
As stated above, there are no limitations to the number of nodes.

[xibaiw]

I re-ran the simulation, and this time it executed successfully.However, I suspect insufficient memory may be the root cause. Thank you for your support Mr. Hesse(●'◡'●).

[xibaiw]

[1341] signal 15: Terminated
in expression starting at REPL[2]:1
Terminated
This is the information I got，after my program is terminated. It looks likes ther RAM is not big enough.

[JTHesse]

Alright, we are also working on reducing the memory allocations of PeriLab.

[xibaiw]

Additionally, when I tried to use MPI for parallel computing in Julia, Julia kept giving me an error that in the command $ mpiexecjl -n 2 julia -e "using PeriLab; PeriLab.main()" 1.yaml -v, the part 2 julia -e "using PeriLab; PeriLab.main()" 1.yaml -v was treated as extra tokens. I tried other commands but still couldn’t get it to work. How should MPI be properly used within PeriLab?

[JTHesse]

Try: mpiexecjl -n 2 julia -e 'using PeriLab; PeriLab.main("1.yaml")'