#include <cppad/cppad.hpp>
#include <boost/numeric/ublas/lu.hpp>

typedef CppAD::AD<double> my_dbl;

int main() {

    boost::numeric::ublas::matrix<my_dbl> a(3, 3);
    boost::numeric::ublas::permutation_matrix<int> pert(3);

    const int singular = lu_factorize(a, pert);

    return 0;
}