diff --git a/src/celmath/geomutil.h b/src/celmath/geomutil.h
index f44e0fc1..2e0fa6a8 100644
--- a/src/celmath/geomutil.h
+++ b/src/celmath/geomutil.h
@@ -47,7 +47,7 @@ ZRotation(T radians)
  *  of the component of 'up' that is orthogonal to the z-axis.
  */
 template<class T> Eigen::Quaternion<T>
-LookAt(Eigen::Matrix<T, 3, 1> from, Eigen::Matrix<T, 3, 1> to, Eigen::Matrix<T, 3, 1> up)
+LookAt(const Eigen::Matrix<T, 3, 1>& from, const Eigen::Matrix<T, 3, 1>& to, const Eigen::Matrix<T, 3, 1>& up)
 {
     Eigen::Matrix<T, 3, 1> n = to - from;
     n.normalize();