On 2/25/22 19:44, Ali Seddiq wrote:

Thank you for your reply. In order to reconstruct the following final discretized form as
(M+A)u + A'p=F+G,
 a 2*2 matrix is subdivided into blocks in a way that (block)components (0,0) and (1,1) are M+A and A' respectively and components  (0,1) and (1,0) set to zero. Solution vector is declared as a BlockVector but contains/subdivided into two scalar components u and p.  Right hand side is also a BlockVector containing F_u and G_p. I hope I could have addressed your point here. Would you please advise me if there are any discrepancies in my structuring?

I have no idea. The error message says
  You are providing a matrix whose subdivision into blocks in either row or
  column direction does not use the same block sizes as the solution vector or
  right hand side vectors, respectively.
You need to find out why that is so. Some debug output right before you call MatrixTools::apply_boundary_values() will be able to tell you whether the error message is correct or not. If the error message is correct, you ought to find out why the subdivision is how it is.

Best
 W.


--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           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 [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/cd49fe8d-95c2-5c10-1162-9f25ff78d914%40colostate.edu.

Reply via email to