Example t8_step5_element_data.jl fails when writing VTK

[ranocha]

I am using Julia v1.9.3 on macOS (ARM architecture). I cloned T8code.jl, activated the examples folder as project with T8code.jl and MPI.jl as dependencies:

(examples) pkg> st
Status `~/.julia/dev/T8code/examples/Project.toml`
  [da04e1cc] MPI v0.20.15
  [d0cc0030] T8code v0.4.1 `..`

The examples step0 to step4 run without problems. However, julia --project=. t8_step4_partition_balance_ghost.jl crashes with the following message

[t8]  [step5] 
[t8]  [step5] Hello, this is the step5 example of t8code.
[t8]  [step5] In this example we will store data on our elements and exchange the data of ghost elements.
[t8]  [step5] 
[t8]  [step5] 
[t8]  [step5] Creating an adapted forest as in step3.
[t8]  [step5] 
[t8] Constructed uniform forest with 8192 global elements.
[t8] Into t8_forest_adapt from 8192 total elements
[t8] Done t8_forest_adapt with 3089 total elements
[t8] Enter  forest partition.
[t8] Done forest partition.
[t8] Into t8_forest_balance with 3089 global elements.
[t8]  Computed maximum occurring level:	4
[t8]  Into t8_forest_ghost with 3089 local elements.
[t8]  Done t8_forest_ghost with 3089 local elements and 0 ghost elements.
[t8]  Profiling: 0
[t8]  Into t8_forest_adapt from 3089 total elements
[t8]  Done t8_forest_adapt with 3089 total elements
[t8] Done t8_forest_balance with 3089 global elements.
[t8]  [step5] Wrote forest to vtu files: *
[t8]  [step5] Computed level and volume data for local elements.
[t8]  [step5] Element 0 has level 2 and volume 3.255208e-04.
[t8]  [step5] Exchanged ghost data.
[libsc 0] Caught signal SEGV
[libsc 0] Abort: Obtained 21 stack frames
[libsc 0] Stack 0: 0   libsc.2.dylib                       0x0000000116895d44 sc_abort_handler + 100
[libsc 0] Stack 1: 1   libsc.2.dylib                       0x0000000116895e54 sc_abort + 20
[libsc 0] Stack 2: 2   libsc.2.dylib                       0x0000000116897994 sc_signal_handler + 100
[libsc 0] Stack 3: 3   libsystem_platform.dylib            0x0000000182d56a24 _sigtramp + 56
[libsc 0] Stack 4: 4   libt8.2.dylib                       0x0000000116f53d9c t8_forest_vtk_write_file + 2796
[libsc 0] Stack 5: 5   ???                                 0x000000011696c11c 0x0 + 4673945884
[libsc 0] Stack 6: 6   ???                                 0x00000001169743d0 0x0 + 4673979344
[libsc 0] Stack 7: 7   libjulia-internal.1.9.dylib         0x000000010531a4ec do_call + 188
[libsc 0] Stack 8: 8   libjulia-internal.1.9.dylib         0x0000000105318d10 eval_body + 1476
[libsc 0] Stack 9: 9   libjulia-internal.1.9.dylib         0x000000010531935c jl_interpret_toplevel_thunk + 260
[libsc 0] Stack 10: 10  libjulia-internal.1.9.dylib         0x000000010533059c jl_toplevel_eval_flex + 4620
[libsc 0] Stack 11: 11  libjulia-internal.1.9.dylib         0x00000001053304c0 jl_toplevel_eval_flex + 4400
[libsc 0] Stack 12: 12  libjulia-internal.1.9.dylib         0x00000001053312e4 ijl_toplevel_eval_in + 156
[libsc 0] Stack 13: 13  sys.dylib                           0x000000013c40d5f4 japi1_include_string_36282.clone_3 + 520
[libsc 0] Stack 14: 14  libjulia-internal.1.9.dylib         0x00000001053015c4 ijl_apply_generic + 1732
[libsc 0] Stack 15: 15  sys.dylib                           0x000000013b0d0b48 japi1__include_47977 + 824
[libsc 0] Stack 16: 16  sys.dylib                           0x000000013af6c00c julia_exec_options_50544 + 21780
[libsc 0] Stack 17: 17  libjulia-internal.1.9.dylib         0x0000000105358644 true_main + 192
[libsc 0] Stack 18: 18  libjulia-internal.1.9.dylib         0x0000000105358538 jl_repl_entrypoint + 180
[libsc 0] Stack 19: 19  julia                               0x0000000104d57f6c main + 12
[libsc 0] Stack 20: 20  dyld                                0x00000001829cff28 start + 2236
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
[unset]: PMIU_write error; fd=-1 buf=:cmd=abort exitcode=1 message=application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
:
system msg for write_line failure : Bad file descriptor

[jmark]

Noted. Looks like a problem on MacOS platforms?

[ranocha]

Could very well be the case. I haven't checked it on other systems. The previous examples using pre-defined VTK output functions work well, though.