费米曲面图

Question

本质上，我试图做一个二维的费米曲面图，即k空间中某个平面的f（n，vec_k）=e_f的二维切割，用插值法。具体来说，我有一个纽比阵列：本征，有形状

在特征形状=（100100100,10），其中前三个索引在向量vec_k上，第三个索引是频带索引“n”。如何通过一个任意曲面特征值==e_f，其中e_f是标量数？在

Answer 1

这通常使用marching cube算法来解决。您应该在这里研究MLab的contour3d()函数：http://docs.enthought.com/mayavi/mayavi/auto/mlab_helper_functions.html#mayavi.mlab.contour3d。这正是你想要的。在

contour3d(Eigen[:,:,:, bandNo], contours=e_f)

应该为你做。也可以为所有标注栏指定多个曲面，如下所示：

^{pr2}$

你应该得到这样的图：http://docs.enthought.com/mayavi/mayavi/_images/enthought_mayavi_mlab_contour3d1.jpg