Исправление
dikiy,
(текущая версия)
:
def f(x,y,z):
return np.matmul(B(x,y,z),np.matmul(D,B(x,y,z))) # твоя функция. Насколько я понимаю, матрицы B функциональные
xx,yy,zz = np.meshgrid(x,y,z) #здесь подразумевается, что x,y,z уже содержат нужные наборты опорных точек значения
f_meshed=f(xx,yy,zz)
np.einsum('i,j,k,ijk',a,a,a,f_meshed)
Исходная версия
dikiy,
:
def f(x,y,z):
return np.matmul(B(x,y,z),np.matmul(D,B(x,y,z))) # твоя функция. Насколько я понимаю, матрицы B функциональные
xx,yy,zz = np.meshgrid(x,y,z) #здесь подразумевается, что x,y,z уже содержат нужные наборты опорных точек значения
f_meshed=f(xx,yy,zz)
np.einsum('i,j,k,ijk',a,a,a,f)