假设我有以下网格omega、phiy和phix。

import numpy as np

phi=np.linspace(-0.75,0.75,num=4)
omega=np.linspace(-100,100,num=5)
omega,phiy,phix=np.meshgrid(omega,phi,phi,sparse=True,indexing='ij')

rand=np.random.rand(5,4)

然后我想得到下面的C

def f(j,k):
   return np.argmin(np.abs(phi[j]**2+np.sin(phi[k])-phi))

C=np.empty((5,4,4))
for i in range(5):
  for j in range(4):
    for k in range(4):
       C[i,j,k]=rand[i,f(j,k)]

有没有更紧凑更好的方法来得到它，使用一些meshgrids的工具？