> The problem is that you are already using a Trilinos vector, whereas the > code expects you to have a vector with a begin() method (like > dealii::Vector<double>). But then it is already easier: All you need to > do is to replace the two lines > solver.SetRHS(&b); > solver.SetLHS(&x); > by > solver.SetRHS(&right_hand_side.trilinos_vector()); > solver.SetLHS(&newton_update.trilinos_vector()); > > Note how you can access the underlying Epetra_Vector from > TrilinosWrappers::Vector. Does this help?
One remark: You need of course to delete the code where x and b are created, too. Best, Martin _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
