[ALPS-users] construct_matrix in update.h in the Hybridization Expansion

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, -- Kuang-Shing Chen