On 07/06/2017 09:02 AM, Weixiong Zheng wrote:
Both are generated with subdivided hyper rectangle with same reps and diagonals.

And no adaptive refinement?

In that case, the solution *should* be the same, assuming your linear solver tolerance is small enough. If it isn't, there is probably a bug somewhere. But it should be possible to debug this on a relatively coarse mesh, since you have such a nice invariant (the symmetry) that you can check on a very coarse mesh as well as in your matrix and rhs vectors. For example, you could look graphically at the rhs vector to see whether the quantity it represents is symmetric or not.

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.

Reply via email to