Hi,<br><br>I look through the source code of hybridization in alps/applications/dmft/qmc and have a question:<br><br>The template function construct_matrix in update.h is defined but NEVER used in the code.<br><br>I use grep to search for that:<br>
grep  construct_matrix * -R<br clear="all"><br>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?<br>
<br>Thank you,<br>-- <br>Kuang-Shing Chen<br><br>