Open vdrouint opened 1 year ago
This is partially solved - a routine has been built, but the find_zero function from Roots.jl does not like the intervals given to it. It is quite sensitive to the noise in the green's functions themselves. Tricky to simplify it all. Work in progress.
At each iteration of the process, one needs to find the solution to n_B(lambda) = q0 - this needs to be done numerically through some process. This is possibly pre-embedded in the Lehman rep package, or can be formulated directly.