LinearAlgebra/addmul/mul!(::Matrix{Float64}, ::Matrix{Float64}, ::LinearAlgebra.Diagonal{Float64, V} where V<:AbstractVector{Float64}, α, β)/α = 0.04046995577353876, β = 0.0/adjoint and transpose/fa = transpose, fb = transpose @test ≈(collect(returned_mat), exp_val, rtol = rtol, atol = atol)