fn operator( local: &LocalMatrix, coupling: &DMatrix<Complex>, ) -> Result<DMatrix<Complex>>
I - T C for a finite coupling of matching size.
I - T C