Ошибка ODEINT: ValueError: установка элемента массива с последовательностью [duplicate]

Z=np.array([1.0,1.0,1.0,1.0])  

def func(TempLake,Z):
    A=TempLake
    B=Z
    return A*B
Nlayers=Z.size
N=3
TempLake=np.zeros((N+1,Nlayers))
kOUT=np.vectorize(func)(TempLake,Z)

Это тоже работает, вместо того, чтобы зацикливать, просто векторизовать, но читать ниже заметки из скудной документации: https://docs.scipy.org/doc/numpy/reference/generated/numpy.vectorize. html

Функция векторизации предоставляется в первую очередь для удобства, а не для производительности. Реализация по существу является циклом for.

Если otypes не указывается, то вызов функции с первым аргументом будет использоваться для определения количества выходов. Результаты этого вызова будут кэшироваться, если кеш имеет значение True, чтобы предотвратить вызов функции дважды. Однако для реализации кеша должна быть завершена оригинальная функция, которая замедлит последующие вызовы, поэтому сделайте это только в том случае, если ваша функция будет дорогой.

Я считаю, что массивы python просто принимают значения. Поэтому преобразуем его в список:

kOUT = np.zeros(N+1)
kOUT = kOUT.tolist()

KOUT [i] - это один элемент списка. Но вы назначаете список этому элементу. ваш func генерирует список.

Чтобы поместить последовательность или другой массив numpy в массив numpy, просто измените эту строку:

kOUT=np.zeros(N+1)

на:

kOUT=np.asarray([None]*(N+1))

Или:

kOUT=np.zeros((N+1), object)

Жаль, что оба ответа определяют проблему, но не дали заключения для ее решения. Давайте посмотрим на код.

Z = np.array([1.0, 1.0, 1.0, 1.0])  

def func(TempLake, Z):
    A = TempLake
    B = Z
    return A * B
Nlayers = Z.size
N = 3
TempLake = np.zeros((N+1, Nlayers))
kOUT = np.zeros(N + 1)

for i in xrange(N):
    # store the i-th result of
    # function "func" in i-th item in kOUT
    kOUT[i] = func(TempLake[i], Z)

Ошибка показывает, что вы устанавливаете i-й элемент kOUT (dtype: int) в массив, каждый элемент в kOUT является всего лишь объектом int, не может указывать на другой тип данных, вы должны изменить оператор, чтобы изменить тип данных kOUT. Например:

Измените следующую инструкцию:

kOUT = np.zeros(N + 1)

в:

kOUT = np.zeros(N + 1, dtype=object)

или:

kOUT = np.zeros((N + 1, N + 1))

Весь код:

import numpy as np
Z = np.array([1.0, 1.0, 1.0, 1.0])

def func(TempLake, Z):
    A = TempLake
    B = Z
    return A * B

Nlayers = Z.size
N = 3
TempLake = np.zeros((N + 1, Nlayers))

kOUT = np.zeros(N + 1, dtype=object)
for i in xrange(N):
    kOUT[i] = func(TempLake[i], Z)

Надеюсь, он вам поможет.