Thanks Matt. It works well for signal core. But is there any solution if I
need a MPI program?
On Tue, Sep 13, 2016 at 9:30 AM, Matthew Knepley <knep...@gmail.com> wrote:
> On Mon, Sep 12, 2016 at 8:24 PM, Ji Zhang <got...@gmail.com> wrote:
>> Dear all,
>> I'm using petsc4py and now face some problems.
>> I have a number of small petsc dense matrices mij, and I want to
>> construct them to a big matrix M like this:
>> [ m11 m12 m13 ]
>> M = | m21 m22 m23 | ,
>> [ m31 m32 m33 ]
>> How could I do it effectively?
>> Now I'm using the code below:
>> # get indexes of matrix mij
>> index1_begin, index1_end = getindex_i( )
>> index2_begin, index2_end = getindex_j( )
>> M[index1_begin:index1_end, index2_begin:index2_end] = mij[:, :]
>> which report such error messages:
>> petsc4py.PETSc.Error: error code 56
>>  MatGetValues() line 1818 in /home/zhangji/PycharmProjects/
>>  MatGetValues_MPIDense() line 154 in /home/zhangji/PycharmProjects/
> Make M a sequential dense matrix.
>>  No support for this operation for this object type
>>  Only local values currently supported
>> Zhang Ji
>> Beijing Computational Science Research Center
>> E-mail: got...@gmail.com
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener