# -*- coding: utf-8 -*-
# -------------------------------------------------------------
# Analysis from File
# -------------------------------------------------------------
path = u'theta-90.hdf5'
configuration = nlread(path, object_id='gID000')[0]

# -------------------------------------------------------------
# Spin Transfer Torque
# -------------------------------------------------------------
kpoint_grid = MonkhorstPackGrid(
    force_timereversal=False,
    )

spin_transfer_torque = SpinTransferTorque(
    configuration=configuration,
    energy=0*eV,
    kpoints=kpoint_grid,
    contributions=Left,
    energy_zero_parameter=AverageFermiLevel,
    infinitesimal=1e-06*eV,
    self_energy_calculator=RecursionSelfEnergy(),
    )
nlsave('theta-90.hdf5', spin_transfer_torque)
