修订版 | c5ae46b810735e050407e4cfed0eeb081aaab857 (tree) |
---|---|
时间 | 2013-01-22 05:48:45 |
作者 | Katsuhiko Nishimra <ktns.87@gmai...> |
Commiter | Katsuhiko Nishimra |
tmp output diis/damp
@@ -889,6 +889,7 @@ void Cndo2::DoDIIS(double** orbitalElectronPopulation, | ||
889 | 889 | } |
890 | 890 | |
891 | 891 | if(diisNumErrorVect <= step && diisEndError<eMax && eMax<diisStartError){ |
892 | + cout << "Using DIIS" << endl; | |
892 | 893 | MolDS_wrappers::Lapack::GetInstance()->Dsysv(diisErrorProducts, |
893 | 894 | diisErrorCoefficients, |
894 | 895 | diisNumErrorVect+1); |
@@ -911,6 +912,7 @@ void Cndo2::DoDamp(double rmsDensity, | ||
911 | 912 | double dampingThresh = Parameters::GetInstance()->GetDampingThreshSCF(); |
912 | 913 | double dampingWeight = Parameters::GetInstance()->GetDampingWeightSCF(); |
913 | 914 | if(0.0 < dampingWeight && dampingThresh < rmsDensity){ |
915 | + cout << "Using Damp" << endl; | |
914 | 916 | stringstream ompErrors; |
915 | 917 | #pragma omp parallel for schedule(auto) |
916 | 918 | for(int j=0; j<molecule.GetTotalNumberAOs(); j++){ |