<div dir="ltr"><div><div>Dear all, <br><br>I'm using petsc4py and now face some problems.<br>I have a number of small petsc dense matrices mij, and I want to construct them to a big matrix M like this:<br><br>         [  m11  m12  m13  ]<br>    M =  |  m21  m22  m23  |   ,<br>         [  m31  m32  m33  ]<br>How could I do it effectively?<br><br>Now I'm using the code below:<br><br>    # get indexes of matrix mij<br>    index1_begin, index1_end = getindex_i( )<br>    index2_begin, index2_end = getindex_j( )<br>    M[index1_begin:index1_end, index2_begin:index2_end] = mij[:, :]<br>which report such error messages:<br><br>    petsc4py.PETSc.Error: error code 56<br>    [0] MatGetValues() line 1818 in /home/zhangji/PycharmProjects/petsc-petsc-31a1859eaff6/src/mat/interface/matrix.c<br>    [0] MatGetValues_MPIDense() line 154 in /home/zhangji/PycharmProjects/petsc-petsc-31a1859eaff6/src/mat/impls/dense/mpi/mpidense.c<br>    [0] No support for this operation for this object type<br>    [0] Only local values currently supported<br><br></div>Thanks. <br><br></div><br><div>2016-09-13</div><span class="gmail-im"><span><span><div>Best, </div><div>Regards,</div><div>Zhang Ji </div><div>Beijing Computational Science Research Center </div><div>E-mail: <a target="_blank" href="mailto:gotofd@gmail.com">gotofd@gmail.com</a><br><br><br></div></span></span></span></div>