Dear developers,
I found that the hybridization matrix is rebuit everytime in remove_antisegment.
https://github.com/ALPSCore/CT-HYB-SEGMENT/blob/master/hybconfig.cpp#L69
Is it necessary? It seems that this line should be commented out as in other methods such as remove_segment.
The thing is that, at low temperature, the program often stops with an error message "in dgesv: info was not zero." This error is raised in invert function, which computes the inversion of large-size matrix, e.g., 100 x 100.
Best regards,
Junya Otsuki
Dear developers,
I found that the hybridization matrix is rebuit everytime in
remove_antisegment.https://github.com/ALPSCore/CT-HYB-SEGMENT/blob/master/hybconfig.cpp#L69
Is it necessary? It seems that this line should be commented out as in other methods such as
remove_segment.The thing is that, at low temperature, the program often stops with an error message "in dgesv: info was not zero." This error is raised in
invertfunction, which computes the inversion of large-size matrix, e.g., 100 x 100.Best regards,
Junya Otsuki