diff --git a/src/geometry/GeometricalFunctions.cpp b/src/geometry/GeometricalFunctions.cpp
index e904312d97b483058ee6f8551142898b2c6ae3aa..3ae39d9ee041d1ccef4ca2bfed7e65ce52e8c389 100644
--- a/src/geometry/GeometricalFunctions.cpp
+++ b/src/geometry/GeometricalFunctions.cpp
@@ -64,6 +64,7 @@ namespace geometry {
 void getClosestLineBoxPoints( const Vector3<real_t>& p1, const Vector3<real_t>& p2, const Vector3<real_t>& c, const Matrix3<real_t>& R,
                               const Vector3<real_t>& side, Vector3<real_t>& lret, Vector3<real_t>& bret )
 {
+   WALBERLA_ABORT("This function produces incorrect results! See https://i10git.cs.fau.de/walberla/walberla/-/issues/25");
    using namespace walberla::math;
    //----- Note: All computations will be performed in the box-relative coordinate-system -----