Improving info to screen output of phononH0 class.

phoebe-team · Mar 23, 2020 · 920f772 · 920f772
1 parent 432291c
commit 920f772
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 4 deletions.
diff --git a/src/harmonic/phononH0.cpp b/src/harmonic/phononH0.cpp
@@ -532,8 +532,6 @@ void PhononH0::dyndiag(Eigen::Tensor<std::complex<double>,4>& dyn,
 
 	Eigen::VectorXd w2 = eigensolver.eigenvalues();
 
-	std::cout << w2.transpose() << "\n";
-
 	for ( int i=0; i<numBands; i++ ) {
 		if ( w2(i) < 0 ) {
 			energies(i) = -sqrt(-w2(i));
@@ -742,6 +740,8 @@ void PhononH0::setAcousticSumRule(const std::string sumRule) {
 		Error e("invalid Acoustic Sum Rule", 1);
 	}
 
+	std::cout << "Start imposing " << sumRule << " acoustic sum rule.\n";
+
 	if ( sumRule == "simple" ) {
 
 		// Simple Acoustic Sum Rule on effective charges
@@ -1281,6 +1281,7 @@ void PhononH0::setAcousticSumRule(const std::string sumRule) {
 		forceConstants = frc_new;
 
 	}
+	std::cout << "Finished imposing " << sumRule << " acoustic sum rule.\n";
 }
 
 

diff --git a/src/parser/qe_input_parser.cpp b/src/parser/qe_input_parser.cpp
@@ -545,7 +545,6 @@ void QEParser::parsePhHarmonic(std::string fileName) {
 
 	//	Now, let's try to diagonalize some points, and start debugging at q=0
 
-
 	PhononH0 dynamicalMatrix(directUnitCell,
 			reciprocalUnitCell,
 			alat,
@@ -567,7 +566,6 @@ void QEParser::parsePhHarmonic(std::string fileName) {
 
 	dynamicalMatrix.diagonalize(q, omega, z);
 
-	std::cout << "Finito\n";
 	std::cout << omega.transpose() * ryToCmm1 << std::endl;
 
 	return;