Castro icon indicating copy to clipboard operation
Castro copied to clipboard

verify that MHD is tile safe

Open zingale opened this issue 3 years ago • 3 comments

PR #2039 turned tiling back on in the MHD solver. We previously disabled it (#289), but that was back when the solver was in Fortran and not ported to the GPU yet. It appears to work now, but we should double check that there are no threading issues.

The main place where an issue can come up is for nodal data when two threads might write to the same location at the same time. This is why we want to size temporary data to the tile box rather than to the full box, in general. We should look through the code to see if there are any potential issues.

zingale avatar Nov 13 '21 19:11 zingale

Note: since none of us really run on CPUs anymore, this has not been much of an issue, but it would be good to have full CPU support in this solver.

zingale avatar Nov 13 '21 19:11 zingale

If you point out the sections of code where the race condition might occur I could look at it.

joehellmers avatar Nov 15 '21 18:11 joehellmers

it's one of those "we'll know it when we see it" situations, so I am not sure I can point to anything, and in fact, all may be well.

zingale avatar Nov 16 '21 20:11 zingale