> PETScWrappers::PreconditionSOR::AdditionalData addldata(1.2); > PETScWrappers::PreconditionSOR preconditioner(flux_matrix, addldata); > > bicgstab.solve (flux_matrix, solution_H0, > flux_rhs, preconditioner); > > This never caused problems before. Now I get a PETSc error: > > [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? Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
