From acbc5df64419f1c11011707358a211eb45be79af Mon Sep 17 00:00:00 2001
From: Mateusz Baran <mateuszbaran89@gmail.com>
Date: Mon, 18 Feb 2019 15:13:08 +0100
Subject: [PATCH] New benchmark for sqrt of small UpperTriangular matrices.

For tracking potential regressions to
https://github.com/JuliaLang/julia/pull/31100 .
---
 src/linalg/LinAlgBenchmarks.jl | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

diff --git a/src/linalg/LinAlgBenchmarks.jl b/src/linalg/LinAlgBenchmarks.jl
index ef9a6c3b..336afd7b 100644
--- a/src/linalg/LinAlgBenchmarks.jl
+++ b/src/linalg/LinAlgBenchmarks.jl
@@ -139,6 +139,20 @@ for s in SIZES
 
 end
 
+# Julia PR 31100
+begin
+    mf4 = randmat(4)
+    mc4 = randmat(Complex{Float64}, 4)
+
+    if VERSION >= v"0.7.0-DEV.1599"
+        g["sqrt", "LinearAlgebra.UpperTriangular{Float64} (#31100)", 4] = @benchmarkable sqrt($mf4)
+        g["sqrt", "LinearAlgebra.UpperTriangular{ComplexF64} (#31100)", 4] = @benchmarkable sqrt($mc4)
+    else
+        g["sqrt", "LinearAlgebra.UpperTriangular{Float64} (#31100)", 4] = @benchmarkable sqrtm($mf4)
+        g["sqrt", "LinearAlgebra.UpperTriangular{ComplexF64} (#31100)", 4] = @benchmarkable sqrtm($mc4)
+    end
+end
+
 for b in values(g)
     b.params.time_tolerance = 0.45
     b.params.samples = 100