corrections by Vadim

added corrections to sqrt and sq functions
passed vector by reference instead of by value in
ECHARM_crystal_interpolation