LinearAlgebra/addmul/mul!(::Matrix{Float64}, ::LinearAlgebra.UpperTriangular{Float64, S} where S<:AbstractMatrix{Float64}, ::LinearAlgebra.Diagonal{Float64, V} where V<:AbstractVector{Float64}, α, β)/α = 1.0, β = -1.4383664597569412/adjoint and transpose/fa = identity, fb = transpose @test ≈(collect(returned_mat), exp_val, rtol = rtol, atol = atol)