Средства и методы разработки

Я обратных матриц не считаю (т.к. они у меня не всегда полного ранга), а потому использую сингулярное разложение (на собственные значения и вектора). Для матрицы 500х500 (типичный мой размер :)) это занимает около 70 мкс (колеблется в интервале 45-100 мкс в зависимости от содержимого). Т.е. вместо обратной матрицы пользуюсь псевдообратной, исключая вклады с собственными значениями, меньшими, чем 1E-10.