[petsc-users] ublas sparse matrix bindings?

[Jed Brown]

On Tue, 29 Jun 2010 12:43:36 -0500, Matthew Knepley <knepley at gmail.com> wrote:
> 3) What you really want is to make MatScalar to be float. Then use
>     doubles for the residual calculation.

I have no idea what Luke's usage is, but MatScalar != PetscScalar is not
currently supported.

  http://lists.mcs.anl.gov/pipermail/petsc-dev/2009-June/001417.html

I'm intrigued by it because my matrices are only for preconditioning
purposes, but haven't caught the bug to try fixing it since the above
thread.

Jed