MAINT: bsplines: use lapack banded solver instead of solve_banded

This is to reduce the memory footprint: solve_banded needs to make
a copy of the lhs matrix internally, but here we can directly construct
the collocation matrix in the LAPACK banded storage format.