SCOREC / pumi-pic

support libraries for unstructured mesh particle in cell simulations on GPUs and CPUs
BSD 3-Clause "New" or "Revised" License
36 stars 14 forks source link

gyro scs failure #34

Closed cwsmith closed 4 years ago

cwsmith commented 4 years ago

mesh: 2M triangles planes: 6 processes-per-group: 6 picParts: 16

221 Particle not in element! ptcl 0 elem 39722 => 39722 orig 0.000000000000000 0.000000000000000 bcc -1109.295 1373.866 -263.570
221 Particle not in element! ptcl 0 elem 39720 => 39720 orig 0.000000000000000 0.000000000000000 bcc 1306.852 -1153.194 -152.659
221 Particle not in element! ptcl 0 elem 39717 => 39717 orig 0.000000000000000 0.000000000000000 bcc -1256.801 1312.209 -54.408
221 Particle not in element! ptcl 0 elem 39716 => 39716 orig 0.000000000000000 0.000000000000000 bcc -219.453 1307.392 -1086.939
221 Particle not in element! ptcl 0 elem 39714 => 39714 orig 0.000000000000000 0.000000000000000 bcc -1104.396 1378.750 -273.353

and other errors

cwsmith commented 4 years ago

Element centroid computation is returning 0.0 with a 24 picpart run on 288 gpus with 6 planes and 2 processes/picpart. The vtu files produced by xgcp mesh render(...) are invalid for picparts 22 and 23 (the last two). If only the coordinate information were incorrect in these vtu files they would still render. This may indicate that the mesh topology is the problem; I'm guessing the vtu files only contain element to vertex downward adjacencies.

update: The centroid check was broken.