(To be removed) Column partition for sparse finite differences
color will be removed in a future release.
group = color(J,P)
group = color(J,P) returns
a partition of the column corresponding to a coloring of the column-intersection
GROUP(I) = J means column
All columns belonging to a color can be estimated in a single finite difference.