On Fri, 1 Jan 2010 13:32:08 -0600, Matthew Knepley <knepley at gmail.com> wrote: > I would rather see then hierarchy change.
Maybe, but many things in the DM interface depend on Mat and Vec. Where would DMGetInterpolation and DMGetMatrix end up? > I think it is natural for the operator (Mat) to depend on the space > (DM) on which it is discretized. Agreed. > Just because shortsightedness in the past has confined us to really > simple spaces (R^N) does not mean we can't change that. KSP operates in a finite dimensional linear space, i.e. isometric to R^n, so at least up to that level, we're not really exploiting the more general interpretation. Jed