ValueError решение обобщенной проблемы собственных значений с помощью scipy ⇐ Python

[Anonymous]

Привет, я пытаюсь решить обобщенную задачу собственных значений (GEVP) на матрице коррелятора.
Часть кода, которая вызывает у меня проблему:

Код: Выделить всё

    import scipy.linalg as la

[...]

corr_j = simple_corr(data_jack, mask)

#compute the eigenvalues
eig_j = np.zeros((nOps, L0, L0))

for t0 in range(2,L0//2-1):
for t1 in range(t0+1, L0//2):

#solve the generalized eigenvalue problem C(t) v = lambda(t) C(t0) v

C = corr_j[:,:,t1]
C0 = corr_j[:,:,t0]

print(f"C = {C.shape}")
print(f"C = {C}")
print(f"C0 = {C0.shape}")
print(f"C0 = {C0}")
eig_t = la.eig(C, b = C0)

Все вызовы печати предназначены для целей отладки. И чтобы показать вам, что матрицы имеют правильную форму, не являются пустыми или сингулярными.
Этот фрагмент кода находится внутри функции gevp_jack(data,mask), которая должна выполнять этот GEVP для множества различных образцы складного ножа.
Когда я запускаю код, я получаю следующий результат:

Код: Выделить всё

C = (7, 7)
C = [[ 1.59584972 -2.15633651  1.24114058  1.75321492 -1.81777217  1.81941323
-0.91749864]
[-2.15633651  2.9590494  -1.68747421 -2.40047725  2.50849557 -2.51730532
1.25051207]
[ 1.24114058 -1.68747421  0.98367127  1.37224169 -1.40936567  1.40846551
-0.72515541]
[ 1.75321492 -2.40047725  1.37224169  1.9320236  -1.99770539  2.00577185
-1.01332363]
[-1.81777217  2.50849557 -1.40936567 -1.99770539  2.17348089 -2.16582443
1.08016084]
[ 1.81941323 -2.51730532  1.40846551  2.00577185 -2.16582443  2.15743968
-1.09213161]
[-0.91749864  1.25051207 -0.72515541 -1.01332363  1.08016084 -1.09213161
0.51369853]]
C0 = (7, 7)
C0 = [[ 5.09266872 -6.49228588  4.14126736  5.37497589 -4.90238238  5.82926292
-3.01593039]
[-6.49228588  8.46178084 -5.22826411 -6.92532952  6.49513817 -7.44136205
3.78233514]
[ 4.14126736 -5.22826411  3.40637298  4.36963407 -3.89903505  4.7350309
-2.4841818 ]
[ 5.37497589 -6.92532952  4.36963407  5.68667154 -5.23611244  6.14753923
-3.16321806]
[-4.90238238  6.49513817 -3.89903505 -5.23611244  5.28974258 -5.64769213
2.79528734]
[ 5.82926292 -7.44136205  4.7350309   6.14753923 -5.64769213  6.65394343
-3.42858045]
[-3.01593039  3.78233514 -2.4841818  -3.16321806  2.79528734 -3.42858045
1.80604997]]
Unexpected exception formatting exception. Falling back to standard exception
Traceback (most recent call last):
File "/home/andrea/anaconda3/lib/python3.11/site-packages/IPython/core/interactiveshell.py", line 3526, in run_code
exec(code_obj, self.user_global_ns, self.user_ns)
File "/tmp/ipykernel_4321/229272726.py", line 2, in 
eig, err_eig = gu.gevp_jack(data, mask)
^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/andrea/Documents/[...]/gevp_utils.py", line 378, in gevp_jack
eig_t = np.sort(eig_t)
^^^^^^^^^^^^^^
File "", line 200, in sort
File "/home/andrea/anaconda3/lib/python3.11/site-packages/numpy/core/fromnumeric.py", line 1029, in sort
a = asanyarray(a).copy(order="K")
^^^^^^^^^^^^^
ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 2 dimensions.  The detected shape was (2, 7) + inhomogeneous part.

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "/home/andrea/anaconda3/lib/python3.11/site-packages/IPython/core/interactiveshell.py", line 2120, in showtraceback
stb = self.InteractiveTB.structured_traceback(
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/IPython/core/ultratb.py", line 1435, in structured_traceback
return FormattedTB.structured_traceback(
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/IPython/core/ultratb.py", line 1326, in structured_traceback
return VerboseTB.structured_traceback(
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/IPython/core/ultratb.py", line 1173, in structured_traceback
formatted_exception = self.format_exception_as_a_whole(etype, evalue, etb, number_of_lines_of_context,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/IPython/core/ultratb.py", line 1088, in format_exception_as_a_whole
frames.append(self.format_record(record))
^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/IPython/core/ultratb.py", line 970, in format_record
frame_info.lines, Colors, self.has_colors, lvals
^^^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/IPython/core/ultratb.py", line 792, in lines
return self._sd.lines
^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/stack_data/utils.py", line 145, in cached_property_wrapper
value = obj.__dict__[self.func.__name__] = self.func(obj)
^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/stack_data/core.py", line 698, in lines
pieces = self.included_pieces
^^^^^^^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/stack_data/utils.py", line 145, in cached_property_wrapper
value = obj.__dict__[self.func.__name__] = self.func(obj)
^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/stack_data/core.py", line 649, in included_pieces
pos = scope_pieces.index(self.executing_piece)
^^^^^^^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/stack_data/utils.py", line 145, in cached_property_wrapper
value = obj.__dict__[self.func.__name__] = self.func(obj)
^^^^^^^^^^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/stack_data/core.py", line 628, in executing_piece
return only(
^^^^^
File "/home/andrea/anaconda3/lib/python3.11/site-packages/executing/executing.py", line 164, in only
raise NotOneValueFound('Expected one value, found 0')
executing.executing.NotOneValueFound: Expected one value, found 0

Я не знаю, что делаю неправильно, кажется, это относится к массиву неоднородного размера, но весь мой массив представляет собой обычные невинные квадратные матрицы

Подробнее здесь: https://stackoverflow.com/questions/786 ... with-scipy