I used the following code to generate a list of tdof's with a boundary attribute of 3 for each process ... inner_dofs[i], targetPos[0], targetPos[1], inner_gtdofs[i]); } The final output shows that ...

Currently the standard Fortran's print *, A prints a 2D array A as a 1D list of numbers. Rather, I would like stdlib to have a function print_array (we can discuss a better name) ... Pretty printing ...