modification of MKL source code

[xilin1]

hello,
i would like to perform the operation A*D*B, where D is diagonal, and A,B are appropriate dimensions using BLAS. there are no routines specific to diagonal matrices, so i was able to make a one-line modification to the DGEMM code available on netlib to implement this. it runs much slower than the BLAS in MKL (of course). is the source code available to make such a modification?

[TimP]

If you compiled netlib DGEMM with a vectorizing compiler, it ought to perform well on problems which don't need cache blocking. By removing the if tests which skip loops involving 0. operations, you would give ifort -O3 more freedom to optimize.