Hi all, Using the print function of Sparsematrix class (Thanks to Daniel for letting me know that) I printed the elements of system_matrix for both sequential and parallel codes. I reduced the problem to only 54 DoFs. It seems that both system_matrices are symmetric and Identical except for some components which despite others, are very small numbers i.e. O(1e-14 or less). So it seems that problem comes from these very small components (probably as a numerical error) which are different for sequential and parallel codes and may destroy S.P.D condition.
I was wondering what is wrong in the parallel code's sytem_matrix that is no longer solved by SolverCG+SSOR. The print of the components( from row 37 to 54) are attached. It would be appreciated if you could take a look at it. Thanks, Hamed -- 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.
System_matrix elements.pdf
Description: Adobe PDF document
