Simplifying everything. MPI code now does not care about actual units. Works only in dimensionless ones.
Basic imag time propagation working.
Cleaned up main loop a bit. Moved unitaries to their respective functions.
Added code to calculate expectation values.