Bounded-operator realization of the resolvent #
This file contains only bounded operators and therefore avoids representing
1 - Δ, point masses, or
multiplication by a distribution as unbounded operators.
For bounded G and M, set K = G M. Invertibility of the Anderson
operator is represented by invertibility of 1 - K. On that event the
recentered resolvent satisfies
(1 - K)⁻¹ G - G = (1 - K)⁻¹ G M G.
The norm-small case is connected to the actual Neumann series in the complete operator algebra.
The bounded operator K = G ∘ M.
Equations
- Anderson4D.Kop G M = G * M
Instances For
The project definition of invertibility of the corresponding Anderson
operator: 1 - K is a unit of the bounded-operator algebra.
Equations
- Anderson4D.LopInvertible G M = IsUnit (1 - Anderson4D.Kop G M)
Instances For
The inverse Green operator on the invertibility event.
Equations
- Anderson4D.inverseGreen G M h = ↑(IsUnit.unit h)⁻¹ * G
Instances For
Matrix coefficient of the recentered inverse with the sign convention of paper (3.23), junk-totalized to zero off the invertibility event.
Equations
- Anderson4D.operatorModeCoeffH G M e α β = if h : Anderson4D.LopInvertible G M then inner ℂ (e (-α)) ((Anderson4D.inverseGreen G M h) (e β) - G (e β)) else 0
Instances For
The inverse unit is a two-sided inverse of 1 - K.
Algebraic core of the conditioned Hilbert--Schmidt difference argument.
Pointwise form of inverseGreen_sub_G.
A norm-small K gives the required invertibility event.
The canonical Neumann inverse of 1 - K when ‖K‖ < 1.
Equations
- Anderson4D.neumannInverse K hK = ↑(Units.oneSub K hK)⁻¹
Instances For
The Neumann inverse is the sum of the operator-valued geometric series.
The operator-valued Neumann series is summable.
Applying the Neumann inverse agrees with the pointwise sum of the iterated operators.
The canonical inverse used by the norm-small event is exactly the
inverse selected by its IsUnit witness.