Commit 7b9fbe4e authored by Alvarez, Gonzalo's avatar Alvarez, Gonzalo
Browse files

expIpi: bug fix

parent 31deac61

src/Ainur/AinurComplex.hh

+1 −1
Original line number Diff line number Diff line
@@ -76,7 +76,7 @@ private: 


			if (isSqrtMinus1) {

				flag = true;

				real1 = atof(buffer.c_str());

				real1 = (buffer == "") ? 0 : atof(buffer.c_str());

				buffer = "";

				continue;

			}

src/Matrix.h

+2 −2
Original line number Diff line number Diff line
@@ -1053,9 +1053,9 @@ expIpi(Matrix<T>& A) 
			std::complex<RT> sum = 0;

			for (SizeType k = 0; k < n; ++k) {

				const RT arg = eigs[k]*M_PI;

				sum += PsimagLite::conj(A(k, i)) *

				sum += PsimagLite::conj(A(i, k)) *

				        std::complex<RT>(cos(arg), sin(arg)) *

				        A(k, j);

				        A(j, k);

			}



			result(i, j) = sum;