M x S x T (M) умножение с BLAS

Возможный дубликат:
Эффективный способ вычисления матрицы продукт AXA'?

У меня есть две матрицы двойной точности:

• M, который имеет размер l ⨯ m и не имеет определенного свойства;
• S, который имеет размер m ⨯ m и является симметричным.

Я хотел бы вычислить:

R = M ⨯ S ⨯ t(M)

где t(M) - транспонирование M.

Просматривая список подпрограмм (C) BLAS, я не нашел ничего, что непосредственно касалось бы этого случая. Какова наиболее эффективная комбинация подпрограмм, которая позволила бы мне вычислить R?