from gpaw.transport.analysor import Transport_Plotter import numpy as np from pylab import * import sys plotter = Transport_Plotter() plotter.plot_setup() if len(sys.argv) <= 2: if len(sys.argv[1]) <= 2: nt = plotter.get_info('nty', int(sys.argv[1]), 0) else: tmp = sys.argv[1].split('-') sam = int(tmp[0]) ref = int(tmp[1]) nt = plotter.get_info('nty', sam, 0) - plotter.get_info('nty', ref, 0) else: nt = plotter.get_info('nty', int(sys.argv[1]), int(sys.argv[2])) matshow((nt[0])) colorbar() xlabel('Transport Direction') ylabel('Pseudo Density') show()
from gpaw.transport.analysor import Transport_Plotter from pylab import * plotter = Transport_Plotter('bias', 'Pt_H2_nsc.dat') plotter.set_ele_steps(None, 0) plotter.plot_ele_step_info('ham', [0], 0, 0)