On 6/17/23 13:58, Abbas Ballout wrote:
I plotted the stress in Paraview along a line for both MPI runs 1 and 3. The
results matched.
That's good news!
TheMatrixTools::apply_boundary_values(boundary_values, system_matrix, tmp,
system_rhs, false)
function appears to be the culprit .
With MPI 1
rhs.l2() before calling apply_boundary_values 3032.65
rhs.l2() after calling apply_boundary_values 1.86145e+10
with MPI 3
rhs.l2() before calling apply_boundary_values 3032.65
rhs.l2() after calling apply_boundary_values 3.67161e+10
Is it because apply_boundary_values isdoing something special to ghost nodes?
Perhaps. I don't know. But you can find out if you output the elements of
system_rhs that correspond to boundary nodes or constrained nodes or in fact
any other node. It shouldn't be an issue with ghost nodes because the rhs
vector has no ghost entries.
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.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/e608b161-302e-c3fc-374c-6454d74843d0%40colostate.edu.
