Thank you. I’ll keep that in mind. I use int32 because that is the int type in C++ and I tend to do @ccall so I keep the int types compatible. Yes I did try to pass a regular Julia int directly and it did not work correctly.
In order to do that I need to provide a mesh. The problematic code is highlighted in the first comment of this thread with a single threaded vs multi threaded approach, where only a mesh is missing to run.
This issue is related to my repo here: rkiefe/Femeko.jl: Not just a PDE toolbox. Adapt your ideas from a clean, modular code base with Femeko. - Codeberg.org
In the meshExample.jl you see an example of how to generate a mesh struct. You can then place the p1! Function call and test. I started my vacation today so I won’t be able to do that myself right now… Thank you.