def get_relative_trajectory(self):
if self._relative_trajectory is None:
self._relative_trajectory = relativize_trajectory(self)
return self._relative_trajectory
calculation.plot_trajectory(atoms=[0], coordinates=[2]) calculation.plot_velocity(atoms=[0], coordinates=[2]) #print(structure.get_cell()) #exit() calculation.write_trajectory_distribution([0, 0, 1], 'distribution.out') calculation.plot_trajectory_distribution([1, 0, 0]) calculation.plot_trajectory_distribution([0, 1, 0]) calculation.plot_trajectory_distribution([0, 0, 1]) exit() import dynaphopy.analysis.coordinates as coortest relative_trajectory = coortest.relativize_trajectory(trajectory) coortest.trajectory_projection(relative_trajectory, trajectory, np.array([1,0,0])) coortest.trajectory_projection(relative_trajectory, trajectory, np.array([0,1,0])) coortest.trajectory_projection(relative_trajectory, trajectory, np.array([0,0,1])) exit() #calculation.plot_vc(atoms=[0,1]) #calculation.plot_vq(modes=[0,1,2,3,4]) #print(structure.get_number_of_atoms()) #calculation.print_phonon_dispersion_spectrum() #calculation.get_phonon_dispersion_spectra() #calculation.set_band_ranges([[[0.2,0.0,0.2],[0.5,0.5,0.5]], [[0.5, 0.5, 0.5], [0.2, 0.0, 0.2]]])
