import numpy as np import sys from phonopy.vasp import read_vasp from enthought.mayavi import mlab def line_plot( m, n, pt ): mlab.plot3d( [pt[m][0], pt[n][0]], [pt[m][1], pt[n][1]], [pt[m][2], pt[n][2]], tube_radius=0.025, colormap='Spect…