[petsc-dev] Mat/DM dependency

[Jed Brown]

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