История изменений

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)

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)