Re: [deal.II] Sparsity pattern discrepancy when applying constraints using Affine constraints

2023-06-29 Thread Wolfgang Bangerth



Vijay:


The constraint matrix that I get from the above code is:
0 = 0
1 = 0
2 = 0
3 = 0
4 = 0
5 = 0


OK, so these all set DoFs 0...6 to zero.


12 18: 1
13 19: 1
14 20: 1
15 21: 1
16 22: 1
17 23: 1


And these say x_12=x_18, x_13=x_19, etc.

Are these all of the constraints you have? The pictures you show add whole 
blocks of entries to the sparsity pattern, but the constraints you show above 
do not support this. Something does not match.


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/c534e801-4a56-7c55-bff0-5586f17a8f82%40colostate.edu.


Re: [deal.II] Is it normal for the results to vary this much when running in parallel.?

2023-06-29 Thread Abbas Ballout
This isn't a problem anymore.
At first I thought that what was happening with step 18 and with distort 
random were caused by the same issue.
I am applying BCs weakly so it's okay.  
Thanks for the help. 

On Wednesday, June 28, 2023 at 4:17:55 PM UTC+2 Wolfgang Bangerth wrote:

> On 6/24/23 09:13, Abbas Ballout wrote:
> > **
> > 
> > Bellow I have plotted the rhs before and after applying boundary 
> constraints 
> > when running with 1, 2 and 3 cores.
> > The rhs looks different at the boundary with MPI 3 as opposed to MPIs 1 
> and 2.
>
> Right. The question is whether it is *wrong*.
> Best
> W.
>
> -- 
> 
> Wolfgang Bangerth email: bang...@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/0655e9e3-7aaa-4454-8ac6-415b35f20564n%40googlegroups.com.