strange behaviour with PetscViewerBinary on MATIS

[Andreas Grassl]

Jed Brown schrieb:
> You can put this in MatView_IS if you really need it, but I doubt it
> will actually be useful.  Unfortunately, you cannot change the domain
> decomposition with Neumann preconditioners, hence they will have limited
> use for solving a system with a saved matrix.  Why do you want to save
> the matrix, it's vastly slower and less useful than a function which
> assembles that matrix?

I assemble the Matrix by reading out from a data structure produced by a
proprietary program and just used this easy approach to compare the solvers on
different machines, where this program is not installed.

Since the implementation of the NN-preconditioner is suboptimal at all, I will
not waste much time on this issues and my post at the list was lead mostly by my
curiosity.

thanks for the explanation

cheers,

ando

-- 
 /"\                               Grassl Andreas
 \ /    ASCII Ribbon Campaign      Uni Innsbruck Institut f. Mathematik
  X      against HTML email        Technikerstr. 13 Zi 709
 / \                               +43 (0)512 507 6091