diff --git a/matrix/helper_functions.hpp b/matrix/helper_functions.hpp
index e13b1f297d..eee74ebeec 100644
--- a/matrix/helper_functions.hpp
+++ b/matrix/helper_functions.hpp
@@ -26,7 +26,7 @@ bool is_finite(Type x) {
  * @param x input possibly outside of the range
  * @param low lower limit of the allowed range
  * @param high upper limit of the allowed range
- * @return wrapped value inside the range
+ * @return wrapped value inside the range, or NAN if value is too far away from range.
  */
 template<typename Type>
 Type wrap(Type x, Type low, Type high) {