diff --git a/src/equations/compressible_euler_3d.jl b/src/equations/compressible_euler_3d.jl
index d6de349fd91..2910765bee4 100644
--- a/src/equations/compressible_euler_3d.jl
+++ b/src/equations/compressible_euler_3d.jl
@@ -1331,7 +1331,7 @@ function flux_hllc(u_ll, u_rr, normal_direction::AbstractVector,
 
     norm_ = norm(normal_direction)
     norm_sq = norm_ * norm_
-    inv_norm_sq = 1.0 / norm_sq
+    inv_norm_sq = inv(norm_sq)
 
     c_ll = sqrt(equations.gamma * p_ll / rho_ll) * norm_
     c_rr = sqrt(equations.gamma * p_rr / rho_rr) * norm_