get several submatrices from several matrices and combine submatrices into another matrix

Matthew Knepley knepley at gmail.com
Wed Dec 31 11:29:16 CST 2008


Yes.

   Matt

On Wed, Dec 31, 2008 at 11:25 AM, Yujie <recrusader at gmail.com> wrote:

> Dear Matthew:
>
> thank you very much. Regarding my problem about MatAXPY() for MPIDense, I
> have checked SEQDense has this function. To my knowledge, in each node,
> MPIDense's format is SEQDense, right? if it is, I should ba able to write a
> MPIDense-based function using SEQDense? thanks  a lot.
>
> Regards,
> Yujie
>
>
> On Wed, Dec 31, 2008 at 8:56 AM, Matthew Knepley <knepley at gmail.com>wrote:
>
>> Yes, you can have 0 length IS arguments.
>>
>>   Matt
>>
>>
>> On Wed, Dec 31, 2008 at 10:21 AM, Yujie <recrusader at gmail.com> wrote:
>>
>>> Hi, PETSc Developers
>>>
>>> In parallel mode, I have 3 MPIDense matrices, A1, A2, A3. Now, I need to
>>> exact submatrices B1, B2, B3 in parallel mode from A1, A2, A3 respectively.
>>> I know I should use MatGetSubMatrix(). My problem is in
>>>
>>> MatGetSubMatrix(Mat mat,IS isrow,IS iscol,PetscInt csize,MatReuse cll,Mat *newmat)
>>>
>>> , the parameter "isrow" means "rows this processor should obtain". if Bi
>>> in some nodes (cpus) of the cluster (Ai have rows in them) don't have rows,
>>> is it work?
>>>
>>> After getting Bi, I will create MPIDense-based C and combine them into C
>>> like
>>>
>>>      B1
>>> C= B2
>>>      B3
>>> MatGetRow();
>>> MatSetValues() should work, right?
>>>
>>> Could you give me some comments about these operations? thanks a lot.
>>> Happy new year!
>>>
>>> Regards,
>>> Yujie
>>>
>>
>>
>>
>> --
>> 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
>>
>
>


-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20081231/ab23b000/attachment.htm>


More information about the petsc-users mailing list