Hi,
I look through the source code of hybridization in alps/applications/dmft/qmc and have a question:
The template function construct_matrix in update.h is defined but NEVER used in the code.
I use grep to search for that: grep construct_matrix * -R
The result is only the definition in update.h. Now the question is: How is the matrix M (inverse of hybridization function) constructed? M is only initialized in the constructor of HybridizationRun but it is never assigned to any value after that. How is the code working without defining M?
Thank you,