<div dir="ltr"><div><div><div>Dear PETSc users, <br><br></div>I am a bit confused  about blocksizes in MPI matrices. <br><br></div>I define a MPI matrix A with a certain blocksize, let say bs = 3. <br></div><div>As far as I understand, the sparse pattern of the matrix A is made of square blocks of size 3x3, regardless of the actual values of the terms of  A.  <br></div><div>Am I right ? <br><br></div><div>The matrix A is distributed among several processors: does the local number of rows of A on each processor have to be a multiple of bs   ? <br><div><br></div></div><div>I want to use a multigrid preconditionner to  solve the linear system of matrix A : A x = b<br>Is it mandatory to define the right hand side b as a vector with the same blocksize bs=3 ? </div><div><br></div><div>Thank you very much of your help, <br></div><div>Best regards, <br></div><div>Natacha <br></div></div>