On Sat, Nov 10, 2012 at 4:30 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote: > On Sat, Nov 10, 2012 at 3:25 PM, Matthew Knepley <knepley at gmail.com> wrote: >> >> Everything is just a point, so I get back mixed dimension things all the >> time. I >> usually separate by stratum. > > > Okay, but what is the compelling use case for DMComplexGetStratumIS() > returning mixed-dimension that cannot be reasonably accomplished by getting > dimensional ISs and union/concatenating them?
I am not sure I understand. Labels just define point sets. If you want to segregate point sets by dimension, that is fine. Nothing in the interface is sensitive to this, nor should it be. >> >> You don't need coordinates for interpolation? > > > Are we talking about the same interpolation? If I have cell-to-vertex > connectivity, I can create the faces without coordinates, yes. Ah, we are not. Mesh interpolation. The way I have done interpolation, you have to know the mesh dimension, DMComplexGetDimension(), and then it assumes that height 0 stuff is cells of that dimension. Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
