Faster matrix sqrt for upper triangular matrices.

Replacing Val(::Bool) in an argument list with an explicit if makes the 
call type stable. The issue was discussed in JuliaLang#31007.

Author: mateuszbaran

stdlib/LinearAlgebra/src/triangular.jl

@@ -2392,7 +2392,13 @@ function sqrt(A::UpperTriangular)
             end
         end
     end
-    sqrt(A,Val(realmatrix))
+    # Writing an explicit if instead of using Val(realmatrix) below
+    # makes the calls to sqrt(::UpperTriangular,::Val) type stable.
+    if realmatrix
+        return sqrt(A,Val(true))
+    else
+        return sqrt(A,Val(false))
+    end
 end
 function sqrt(A::UpperTriangular{T},::Val{realmatrix}) where {T,realmatrix}
     B = A.data