adding drop tolerance to ilut inverse construction
separating richardson extrapolation relaxation parameter from neumann series relaxation parameter