vcl_size_t internal_size1(matrix_base<NumericT, F> const & mat) { return mat.internal_size1(); }
vcl_size_t ld(matrix_base<NumericT> const & mat) { if (mat.row_major()) return mat.internal_size2(); return mat.internal_size1(); }