Segfault in markergen rhoto5d about ascot5 HOT 1 OPEN

Just noting that (in case if you're in a hurry) rhoto5d is not needed to generate alpha particles. You can sample markers directly from the fusion distribution and, if you need importance sampling for the wall loads, just pick ones that have rho > some_value. For practical purposes this is quite close to what rhoto5d would provide.

In any case this segfault should be fixed. I'm fairly certain this is stellarator specific as it has not been tested for stellarators. My guess is that it is this line where the segfault is coming from (though this is the only part I'm fairly sure should work for stellarators):

Could you test that by replicating that function call in a separate script? Just replace self._ascot with the Ascot object (a5). You could also try running valgrind to locate the segfault, but that might mean you'd have to use smaller distribution grid resolution for it to complete within a reasonable time.

Can you confirm that the magnetic field data covers the entire (R,z) grid you are using in the distribution?

from ascot5.