internal::band_matrix::lu_solve(std::vector<double,std::allocator<double>> const&,bool)