diff --git a/glm/detail/func_geometric.inl b/glm/detail/func_geometric.inl
index b3590152..9cde28fe 100644
--- a/glm/detail/func_geometric.inl
+++ b/glm/detail/func_geometric.inl
@@ -114,11 +114,10 @@ namespace detail
 	{
 		GLM_FUNC_QUALIFIER static vec<L, T, Q> call(vec<L, T, Q> const& I, vec<L, T, Q> const& N, T eta)
 		{
-			assert(eta >= static_cast<T>(-1) && eta <= static_cast<T>(1));
-
 			T const dotValue(dot(N, I));
 			T const k(static_cast<T>(1) - eta * eta * (static_cast<T>(1) - dotValue * dotValue));
-			vec<L, T, Q> const Result = (eta * I - (eta * dotValue + std::sqrt(k)) * N) * static_cast<T>(k >= static_cast<T>(0));
+			vec<L, T, Q> const Result =
+                (k >= static_cast<T>(0)) ? (eta * I - (eta * dotValue + std::sqrt(k)) * N) : vec<L, T, Q>(0);
 			return Result;
 		}
 	};
diff --git a/glm/geometric.hpp b/glm/geometric.hpp
index 28d76199..e0380830 100644
--- a/glm/geometric.hpp
+++ b/glm/geometric.hpp
@@ -99,8 +99,6 @@ namespace glm
 	/// and the ratio of indices of refraction eta,
 	/// return the refraction vector.
 	///
-	/// @param eta Indice of refraction. Must be a value between -1 and 1 inclusively.
-	///
 	/// @tparam L An integer between 1 and 4 included that qualify the dimension of the vector.
 	/// @tparam T Floating-point scalar types.
 	///