Numerically stable solution of coupled channel equations: The local reflection matrix