# -------------------------------------------------------------
# Load device configuration
# -------------------------------------------------------------
device_configuration = nlread('mgo_para.nc', DeviceConfiguration)[0]

# -------------------------------------------------------------
# Transmission Spectrum using a dense Monkhorst-Pack grid.
# -------------------------------------------------------------
transmission_spectrum = TransmissionSpectrum(
    configuration=device_configuration,
    energies=numpy.linspace(0,0,1)*eV,
    kpoints=MonkhorstPackGrid(500,500),
    )
nlsave('adaptive_grid.nc', transmission_spectrum)

# -------------------------------------------------------------
# AdaptiveGrid.
# -------------------------------------------------------------
adaptive_grid = AdaptiveGrid(
    tolerance=0.01,
    error_measure=Relative)

# -------------------------------------------------------------
# Transmission Spectrum using the AdaptiveGrid.
# -------------------------------------------------------------
transmission_spectrum = TransmissionSpectrum(
    configuration=device_configuration,
    energies=numpy.linspace(0,0,1)*eV,
    kpoints=adaptive_grid,
    )
nlsave('adaptive_grid.nc', transmission_spectrum)