Compatibility with multithreading

[mjyshin]

I wanted to look at a range of parameters on the same system, so I tried using Threads.@threads on the for loop, and I got the following error:

    nested task error: Allocations in assembly loop: cells: 208, bfaces: 0
    See the documentation of `check_allocs!` for more information

It runs just fine if I run it in series, but I thought I would try to speed up the process any way I could. Is this something VoronoiFVM cannot handle (yet), or did I fumble somewhere? Thank you!

[j-fu]

Hi, thanks for asking. Parallelization is in the plans, but this will need some work:

• grid partitioning + partition graph coloring
• multithreaded sparse matrix assembly
• parallel iterative solvers
  The later point possibly can be handeled relatively independently, there is a master thesis running on this and we will see if this functionality could be added first.

In the moment the multithreaded loop would end up in write conflicts.

[mjyshin]

Thank you for your reply, that sounds great, and I will patiently await the updates!

Two unrelated questions I had were 1) if there was a way to do unstructured/random Delaunay triangulation, even on a domain like (0,1) x (0,1), and 2) if you had any plans for an adaptive (i.e., time dependent) mesh at any point in the future?

[j-fu]

Hi,

as for 1) Have a look at SimplexGridFactory.jl - this package has been made for this purpose.
For 2): this has lesser priority than other things (and would need some input from math regarding error estimators etc.)