Request GEMMT API

[edelsohn]

Intel MKL provides an additional GEMM variant, GEMMT, that updates only the upper or lower triangular part of the result matrix. This would be a great addition to OpenBLAS.

https://software.intel.com/en-us/node/590047

[brada4]

You mean - to wrap _symm with 2 transforms?

[edelsohn]

I believe that the difference is more than wrapping SYMM with two transforms, but maybe I am misunderstanding your question.

SYMM is described as
C := alpha_A_B + beta*C
https://software.intel.com/en-us/node/468488

GEMMT is described as
C := alpha_op(A)_op(B) + beta*C
https://software.intel.com/en-us/node/590047

with different limitations on A, B, and C. For GEMMT, the upper triangular or lower triangular part of C is overwritten by the respective part of the result.

[brada4]

I wrote DGEMMT in julia and in no-transform case it became plain DSYMM call...
It is really k--ka-s job they did with llvm

[zerothi]

This would be nice! :)

[xrq-phys]

Hi.

I'm curious if anyone's still tuned-in for this request.
To reillustrate the problem, GEMMT computes only a part of GEMM.
It's different from SYMM in my opinion as neither of the 2 operands A&B has symmetry restrictions.
Rather, it's more similar to SYR2K but with only A*transpose(B) term.

GEMMT is beyond the BLAS standard but I guess it's implementation could be very close to SRYK/HERK?

[martin-frbg]

Not implemented yet, but not forgotten either. I agree that it does look more like syr2k than symm.

[jgillis]

Fwiw, MUMPS solver ( http://mumps.enseeiht.fr/doc/userguide_5.4.1.pdf ) would benefit a lot: "We strongly recommend to use this ability if your BLAS library enables it"

[egaudry]

I'd like to second that : GEMMT feature's implementation would be great :).

[mohawk2]

ReLAPACK provides this: https://github.com/HPAC/ReLAPACK/blob/master/src/dgemmt.c

[egaudry]

Thanks for the pointer, I still believe this feature might be offered by the BLAS it-self (as with other BLAS implementations on the market).

[mohawk2]

Not the reference one: https://github.com/Reference-LAPACK/lapack/search?q=dgemmt shows no hits.

[egaudry]

agreed, however when looking at MKL and BLIS, you can see them supported as it enables extra performance on various workload (MUMPS Direct Sparse Solver for instance); it seems like a low hanging fruit.

Also, I understand why you gave the pointer to ReLAPACK, thanks for sharing it.

[martin-frbg]

Lots of low hanging fruit but never enough pickers (happens in real-world orchards too). ReLAPACK is included in OpenBLAS as a build-time option, but the gemmt there is not built by default (even in the original ReLAPACK source, see its config.h). I have not gotten around to checking if Peise's algorithm there actually works and is efficient.

[egaudry]

fair point @martin-frbg; I'll stop arguing as I cannot devote time to help :)