Dg = colon(D,g)
DJ = colon(D,J)
R. Krone and A. Leykin, "Numerical algorithms for detecting embedded components.", arXiv:1405.7871
The object colon is a method function with options.