undefined | Better HN

Better HN

Top Best Ask Show New Jobs

0 pointsKolenCh4y ago0 comments

Any reference to how Julia outperform Fortran in BLAS? I thought it is/was calling existing BLAS libraries?

0 pointsKolenCh4y ago0 comments

Any reference to how Julia outperform Fortran in BLAS? I thought it is/was calling existing BLAS libraries?

0 comments

4 comments · 1 top-level

adgjlsfhk14y ago· 3 in thread

Julia currently ships with fortran based Blas, but Octavian.jl is apure Julia matmul that is faster. (it's nowhere near finished though)

cbkeller4y ago

And the "how" behind Octavian.jl is basically LoopVectorization.jl [1], which helps make optimal use of your CPU's SIMD instructions.

Currently there can some nontrivial compilation latency with this approach, but since LV ultimately emits custom LLVM it's actually perfectly compatible with StaticCompiler.jl [2] following Mason's rewrite, so stay tuned on that front.

[1] https://github.com/JuliaSIMD/LoopVectorization.jl

[2] https://github.com/tshort/StaticCompiler.jl

KolenChOP4y ago

Thanks. But how LoopVectorization.jl is helping here, say comparing to C/Fortran optimized w.r.t. to the CPU? Is there somewhere in their doc mentioning this?

1 more reply

KolenChOP4y ago

Thanks! Very interesting.

Link here for others: https://octavian.julialinearalgebra.org/stable/

j / k navigate · click thread line to collapse

0 comments

4 comments · 1 top-level

adgjlsfhk14y ago· 3 in thread

Julia currently ships with fortran based Blas, but Octavian.jl is apure Julia matmul that is faster. (it's nowhere near finished though)

cbkeller4y ago

And the "how" behind Octavian.jl is basically LoopVectorization.jl [1], which helps make optimal use of your CPU's SIMD instructions.

[1] https://github.com/JuliaSIMD/LoopVectorization.jl

[2] https://github.com/tshort/StaticCompiler.jl

KolenChOP4y ago

Thanks. But how LoopVectorization.jl is helping here, say comparing to C/Fortran optimized w.r.t. to the CPU? Is there somewhere in their doc mentioning this?

1 more reply

KolenChOP4y ago

Thanks! Very interesting.

Link here for others: https://octavian.julialinearalgebra.org/stable/

j / k navigate · click thread line to collapse