Edith, If I understand you correctly, you basically consider step-7 and replace the Neumann boundary conditions by inhomogeneous Dirichlet boundary conditions. The system matrix is symmetric and positive definite and hence CG should work.
Can you verify that the solution is correct if you use a direct solver (SparseDirectUMFPACK [1])? Best, Daniel [1] https://www.dealii.org/8.4.1/doxygen/deal.II/classSparseDirectUMFPACK.html -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
