Note: The output unit for <T> is hbar/2. ...this has to be changed in the future.
Note2: Switches have to be implemented to toggle the computation of these quantities.

Pre merge commit. Commented out the non-funtional source-free code. Numerical problems yet to be solved...

PBE and FFN now bring the BCC Iron testcase down to 1.06*10**-3 using straight mixing. => It seems that there are still some convergence Issues.

FFN converges now with GGA in 46 iterations. But the magnetic moment of BCC Iron is calculated to 2.01 which indicates that there are still error in the code.

Using GGA we achive Dist of 10**-3 in 150 iterations now. There still seem to be some error in the code which I'm trying to locate now.

Added 'Rotation' of den_work. Next=> Thoughts about grads of magnetic moments and how to implement this on radial grid. Also: Angles?

Fully fully noco not working (temporarily) with LDA anymore after this commit. Preparations to fix GGA which has never been working with fully fully noco.

Introduced rot of den before vgen and backwards rot before mixing. HIGHLY experimental. Not stable converging at all.