d3q Code

Minor changes, code still need polishing

Phonon dos subroutines moved here

I'm not sure this code work, putting it here so it does not get lost

An mpi barrier woould cause the code to stop if print_percent_wall is called inside a loop when the number of q-points is not a multiple of the number of CPUs

File was missing, neutrons subroutines moved here

QHA example, reference and testing still required

relatively decent QHA implementation added, documentation still lacking. Fun to use in conjunction with d3_db.x

The EOS from PW/tools/ev.f0 converted to library and interfaced with the Powell minimization (grabbed from the internet and made into a library as well). The original very simple minimization still exists

Dynbubble code polished a bit and removed the use of Bubble_sort in favour of quicksort

Oops