Как быстро производить матричные операции?

Доброго времени суток!

Встала проблема по максимально возможному ускорению таких операций:

1. перемножение, сложение матриц 2х2, 4х4; умножение такой матрицы на вектор (вычисления интеграла методом Симпсона). Числа комплексные, double.
2. произведение векторов. Один -- набор констант, другой -- вычисляется некоторой функцией. Числа вещественные, double.

На данный момент все матрицы и вектора хранятся в простых массивах. Пункт 1 реализую с помощью простейших и очевидных самописных макросов (по определению операций). Пункт 2 выполняю в цикле.
Компилирую с помощью gcc 4.4.3 при опции O2. Параметр -msse4.2 скорость не изменяет... Стоит вопрос, что можно сделать:

• Boost?
• STL?
• GSL?
• Явным образом переписать с использованием SSE?
• Иное? :)

К сожалению, не могу понять, используют ли упомянутые библиотеки SSE-инструкции, и вообще на сколько это будет оптимальнее, чем моя реализация. Экспериментов с использованием различных методов пока не проводил... Надеюсь на пинок консультирующего свойства от знающего человека :)

С уважением, V_oron.