Calculate the determinant of any matrix

Boost.uBLAS can calculate the determinant of (at most) 2x2 matrices. I have a piece of code that can do this for any matrix size (which is http://www.richelbilderbeek.nl/CppUblasMatrixExample7.htm), that already has been tested.

If accepted, I volunteer to try to add this to Boost, but I'd need some help to up the quality even more.

Is this already been done? If no, would it be appreciated if I (create and) submit a Pull Request?

The issue asks whether Boost.uBLAS should support determinants for matrices larger than 2x2 and offers an externally linked implementation. Read the linked example and check the current uBLAS determinant support and project contribution guidance before proposing a direction. Done would mean a decided, accepted approach and—if pursued—a tested contribution integrated into uBLAS; the thread has had no comment for years.

Written by the indexing model from the issue text.