changed some abs() to fabs() just to be sure; however, that shouldn't be necessary in C++, see: https://stackoverflow.com/questions/3118165/when-do-i-use-fabs-and-when-is-it-sufficient-to-use-stdabs and https://en.cppreference.com/w/cpp/numeric/math/fabs