`mul!(C,A,B)` hits method from DoubleFloats even after loading this package

[baggepinnen]

Hello and thank you for this package :)

I just tried it out, and noticed that

using DoubleBLAS
LinearAlgebra.mul!(C, A, B);

hits a method defined in DoubleFloats.jl https://github.com/JuliaMath/DoubleFloats.jl/blob/master/src/math/linearalgebra/matmul.jl#L81

How would I go about replacing that method with optimized ones from this package? Calling https://github.com/RalphAS/DoubleBLAS.jl/blob/master/src/gemm.jl#L6 directly is not really feasible if I'd like to keep implementations unchanged.

[RalphAS]

For the time being, the elaborated version mul!(C, A, B, true, false) should call into this package instead. I think you could put a 3-arg version of mul! into your own project that overrides the one in DoubleFloats to call the 5-arg one.

Secondly, I would not call the methods here "optimized"! They are a proof of principle until the experts have time and interest to do better.