update herk for consistency

Author: araujoms

src/matmul.jl

@@ -782,8 +782,7 @@ Base.@constprop :aggressive function syrk_wrapper!(C::StridedMatrix{T}, tA::Abst
         else
             generic_syrk!(C, A, false, tA_uc == 'N', alpha, beta)
         end
-        copytri!(C, 'U')
-        return C
+        return copytri!(C, 'U')
     end
     return gemm_wrapper!(C, tA, tAt, A, A, α, β)
 end
@@ -829,10 +828,11 @@ Base.@constprop :aggressive function herk_wrapper!(C::Union{StridedMatrix{T}, St
                 stride(A, 1) == stride(C, 1) == 1 &&
                 _fullstride2(A) && _fullstride2(C)) &&
                 max(nA, mA) ≥ 4
-            return copytri!(BLAS.herk!('U', tA, alpha, A, beta, C), 'U', true)
+            BLAS.herk!('U', tA, alpha, A, beta, C)
         else
-            return copytri!(generic_syrk!(C, A, true, tA_uc == 'N', alpha, beta), 'U', true)
+            generic_syrk!(C, A, true, tA_uc == 'N', alpha, beta)
         end
+        return copytri!(C, 'U', true)
     end
     return gemm_wrapper!(C, tA, tAt, A, A, α, β)
 end