sistema lineal más rápido para matrices cuadradas pequeñas (10x10)

Estoy muy interesado en optimizar al máximo la resolución de sistemas lineales para matrices pequeñas (10x10), a veces llamadas matrices pequeñas . ¿Hay una solución lista para esto? La matriz se puede suponer no singular. Este solucionador se ejecutará más de 1 000 000 de veces en microsegundos...