This introduces the following options to the KKRhost code:

* NSIMPLEMIXFIRST= integer: number of simple mixing steps that are done before Broyden or Anderson mixing starts
* SPINMIXQBOUND= real: threshold of the rms of the nonco angles (in degrees) that is used to fix the directions if the angles change less than this value
* <use_boyden_spinmix>= boolean: activate broyden mixing for noncollinear angles
* <write_angles_alliter>= boolean: activate writeout of the nonco angles for all iterations in a new file

The following parameter control the new broyden mixing for the nonco angles:

* SPINMIXALPHA= real: mixing parameter use in the Broyden spinmixing
* SPINMIXNSIMPLE= integer: number of simple mixing steps before Broyden spin mixing starts
* SPINMIXMEMLEN= integer: number of iterations that are remembered in the Broyden mixing

- NSIMPLEMIXFIRST>0 activates simple mixing for the first NSIMPLEMIXFIRST if IMIX>3
- SPINMIXQBOUND is the spin mixing threshold after which all angles are assumed to be converged and will be fixed

Bug fix: corrected calculation of spherical angles from the spin moment vectors for the cases theta=0 and pi

Changes involving read_angles:
-- wunfiles: sets theta, phi and fixdir in t_params, new override with fix_nonco_angles, broadcasts values
-- main1a, main1b, main1c, writekkrflex: replaced call to read_angles with assignment from t_params

Changes involving fix_nonco_angles:
-- wunfiles: used to override fixdir
-- rhovalnew: now uses fixdir to decide if the local frame has to be corrected, new angles always computed
-- main1c: decides on nonco_angle_out.dat, now always write information on new angles

Take kkrparser_functions.py from masci-tools, commit 2df8019

Otherwise from 50x50x50 k mesh on the NofKs number exceeds the field width

force at least two iterations for LDA+U (otherwise initial step will stop if starting potential is converged)

- fix initialization of gfint
- implement QBOUND_LDAU to control how well LDAU potential should be converged

replaces cutoff compared to qbound

Add EFSET option for voronoi which sets Fermi level to disired value and shifts core states accordingly