Switch to column-major layout

[gha3mi]

Currently uses row-major layout. Example:

Tgc(i,:) = kron( &
    basis_bspline(Xt(i,2), knot2, nc(2), degree(2)), &
    basis_bspline(Xt(i,1), knot1, nc(1), degree(1)) )

Column-major layout would look like:

Tgc(:,i) = kron( &
    basis_bspline(Xt(2,i), knot2, nc(2), degree(2)), &
    basis_bspline(Xt(1,i), knot1, nc(1), degree(1)) )

Switching requires consistent updates across the entire library. Consider for a future major release.