> > [0]PETSC ERROR: No support for this operation for this object type!
> > [0]PETSC ERROR: No support for omega != 1.0; use -mat_no_inode!
>
> I don't know why PETSc changed this, but it appears to me from the error
> message that you can only use a relaxation factor of 1.0 with SOR now. You
> set it to 1.2 above. What happens if you leave it at 1.0?
This is indeed the case as Wolfgang says (I also do not know why).
This changed in petsc-3.0.0.p0 and affects anything after that.
Feel free to "safely" comment out:
// PETScWrappers::PreconditionSOR::AdditionalData addldata(1.2);
since the deal.II interface defaults to this value.
On the other hand, if you really want to use a damping
factor and petsc-3.0.0-px, if you figure out how to do it please let me
know and we can try to update the deal.II wrappers!
It is in my interests to do things like this.... ;-)
Hint: Check for differences between petsc versions 3.0.0 and 2.3.3
in the subroutine first:
MatRelax () around line 3000 of in src/mat/interface/matrix.c
(the other errors follow from this). I guess there must be an alternative
way to invoke damping in petsc-3.0.0? Can you find it?
Best,
Toby
-----
Toby D. Young
Philosopher-Physicist
Adiunkt (Assistant Professor)
Polish Academy of Sciences
Warszawa, Polska
www: http://www.ippt.gov.pl/~tyoung
skype: stenografia
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
