I need to either invert or solve for g iteratively in the most efficient way possible. Can anyone suggesta routine that they know of themselves or know of in the NetLib repository (preferably in fortran).