title: determinant calculation with SIMD
body: <p>Does there exist an approach for calculating the determinant of matrices with low dimensions (about 4), that works well with SIMD (neon, SSE, SSE2)? I am using a hand-expansion formula, which does not work so well. I am using  SSE all the way to SSE3 and neon, both under linux. The matrix elements are all floats.</p>