From 2a9e482e44cd3b49bdd292ed43e443a89c2612a5 Mon Sep 17 00:00:00 2001
From: Jan Petykiewicz <jan@mpxd.net>
Date: Tue, 23 May 2023 12:51:12 -0700
Subject: [PATCH] Z is y0 transposed

---
 meanas/fdfd/bloch.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/meanas/fdfd/bloch.py b/meanas/fdfd/bloch.py
index 1dadea2..afada02 100644
--- a/meanas/fdfd/bloch.py
+++ b/meanas/fdfd/bloch.py
@@ -561,7 +561,7 @@ def eigsolve(
     if y0 is None:
         Z = numpy.random.rand(*y_shape) + 1j * numpy.random.rand(*y_shape)
     else:
-        Z = numpy.array(y0, copy=False)
+        Z = numpy.array(y0, copy=False).T
 
     while True:
         Z *= num_modes / norm(Z)