COBS -- Constrained B-splines (Sparse matrix based)