On 11/06/2016 11:15 AM, Hamed Babaei wrote:
I am attaching the sequential and parallel codes as well as parameters.prm file.
It would be appreciated if you could compare tangent_matrix (the system_matrix
of the elasticity part of the problem) of parallel and sequential codes to
find out where I have ruined the S.P.D condition that the parallel code is no
longer solved by SolverCG+SSOR.
Hamed -- this goes beyond the level of support you can typically expect from
the volunteers on this mailing list. In your career, you will encounter
problems such as the ones you are having many times, and you will need to
develop the (mental and software) tools to deal with such situations yourself.
Use the opportunity to teach yourself about
- coming up with minimal testcases
- debugging in parallel
- printing output from different versions of the program and comparing
them.
In your case, I would try to reduce the problem as much as you can. For
example, is there a difference between the sequential code and the parallel
one only in 3d? Do you see the same problem in 2d? Does it happen in the first
iteration, and if so, can you just remove all of the code that only runs after
that and that clearly can't have anything to do with the issue (vastly
reducing the length of the program in which you suspect there to be a bug)?
What happens if you reduce the number of global refinements? Can you reduce it
to 1 (which means only 4 cells)? If so, if you have only 4 cells, the matrices
may be small enough to output and compare by hand. If you can't, write a small
script that takes the matrices and compares them internally.
These are the strategies successful programmers use all the time to find
problems. Use the opportunity to learn how to deal with such issues.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
