24 if (displacement[0] == 0 && displacement[1] == 0 && displacement[2] == 0) {
37 real_t radius = std::sqrt(local_cartesian[0] * local_cartesian[0] +
38 local_cartesian[1] * local_cartesian[1] +
39 local_cartesian[2] * local_cartesian[2]);
40 return {radius, std::acos(local_cartesian[2] / radius),
41 std::atan2(local_cartesian[1], local_cartesian[0])};