def METHOD_NAME(self, path, n_pts=100, method="cubic"):
 ''' 
 :param path: a list of points defining a path in the BZ [array(x,y), ...]
 :param n_pts: the number of points on the path
 :param method: interpolation method (cubic|linear|nearest)
 :return: L,Ltot
  with L: list of momenta on path
  and Ltot: Ltot[0][0] : list of eps(k)(0,0) for each k
 '''
 x,y,z,zmin,zmax = make_plottable(self, method=method)
 #print z
 #where_are_NaNs = np.isnan(z)
 #z[where_are_NaNs] = -10
 z=np.nan_to_num(z)
 sp_real = interpolate.RectBivariateSpline(x, y, z[0,0,:,:].real, kx=2, ky=2, s=0)
 sp_imag = interpolate.RectBivariateSpline(x, y, z[0,0,:,:].imag, kx=2, ky=2, s=0)
 L, high_sym=generate_points_on_path(path,n_pts)
 Lz_on_path = [sp_real(x0,y0)[0][0]+1j*sp_imag(x0,y0)[0][0] for x0,y0 in L]
 return L, np.array(Lz_on_path), high_sym