initial commit

parity symmetry added in TPM and SPM, bar function works.
added Q condition, tested it and works
primal dual works for PQ
G up condition works.
bar function of PHM to SPM works
G down works
program works with PQG

but still an error in the spin expectation value function TPM::spin()
spin fixed
added the DPM class

started with the bookkeeping, tomorrow K > 0!
did the bookkeeping
get_phase_order function of DPM class is done

this does the bookkeeping of the parity states
T1 up condition works
S = 0 part of the TPM::bar(DPM) function works

S = 1 still contains bugs
the S=1 part of the TPM::bar(DPM) function works too now
T down condition works!
the full program works with the T1 condition!
Added PPHM class, did the new bookkeeping, dimensions work out
found a bug in de constructor, fixed it

wrong dimension
T up condition works!
the TPM::bar(const PPHM &) function works