From ded7bed6fff8cad4eeddb0bfe896d681f1495062 Mon Sep 17 00:00:00 2001 From: Ardavan Oskooi Date: Sun, 2 Oct 2022 21:51:22 -0700 Subject: [PATCH] fix output in fields::get_eigenmode for MPB mode wavevector in Cartesian coordinates converted from non-orthogonal reciprocal lattice --- src/mpb.cpp | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/mpb.cpp b/src/mpb.cpp index f2e39d84d..f10b7e828 100644 --- a/src/mpb.cpp +++ b/src/mpb.cpp @@ -572,8 +572,10 @@ void *fields::get_eigenmode(double frequency, direction d, const volume where, c if (verbosity > 0) master_printf("MPB solved for frequency_%d(%g,%g,%g) " "= %g after %d iters\n", - band_num, G[0][0] * k[0], G[1][1] * k[1], G[2][2] * k[2], - sqrt(eigvals[band_num - 1]), num_iters); + band_num, G[0][0] * k[0] + G[1][0] * k[1] + G[2][0] * k[2], + G[0][1] * k[0] + G[1][1] * k[1] + G[2][1] * k[2], + G[0][2] * k[0] + G[1][2] * k[1] + G[2][2] * k[2], sqrt(eigvals[band_num - 1]), + num_iters); // copy desired single eigenvector into scratch arrays evectmatrix_resize(&W[0], 1, 0);