Madderom’s Fortran subroutine
DILUCG
[
107
]
, which implements the method of
[
94
]
, has been used by a number of
numerical relativists for both this and other purposes.