public function penta_periodic_columns(tdsops) result(z)
Solves A_np*z_k = e_{p_k} for p = 1, 2, n-1, n with the non-periodic
pentadiagonal LU factors, as thread 1 of der_penta_periodic does in
the CUDA backend.
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.
Where possible, edges connecting nodes are
given different colours to make them easier to distinguish in
large graphs.