Pseudo code

C := alpha A × B + beta C or C := alpha B × A + beta C, where alpha and beta are scalars, A is a symmetric or hermitian matrix and B and C are m × n matrices.