ILUCG3: Subprograms for the solution of a linear asymmetric matrix equation arising from A 7, 15, 19 or 27 point 3d discretization - CPC, Mendeley (Mendeley, 1983)