MatGetLocalSize() and MatDenseGetLocalMatrix()

[Yujie]

Dear Matthew:

Two processors are used. The matrix dimension is 105*108;

MatGetOwnershipRange()

proc1: 0->56

proc2: 56->105

Mat->camp.rstart; Mat->cmap.rend

proc1: 0->54

proc2: 54->108

MatGetLocalSize()

proc1: row 56 col 54

porc2: row 49 col 54

MatDenseGetLocalMatrix()

proc1: 56*108

proc2: 49*108

thanks.

Yujie

On Tue, Jan 6, 2009 at 11:33 AM, Matthew Knepley <knepley at gmail.com> wrote:

> On Tue, Jan 6, 2009 at 11:52 AM, Yujie <recrusader at gmail.com> wrote:
>
>> Dear PETSc developers:
>>
>> I am trying to use MatGetArray() and MatSetValues() to combine several
>> MPIDense matrices into one matrix. At the beginning, I use
>>
>> MatGetOwnershipRange() and Mat->camp.rstart; Mat->cmap.rend (2.3.3-p8
>> version) to get the start, end row and column. I can calculate the local
>> rows and columns.
>>
>> I also use MatGetLocalSize() to confirm the accuracy. However, I always
>> find some data loses in the combined matrix.
>>
>> And then, I try to use MatDenseGetLocalMatrix() to get the lcoal matrix
>> and output it. I find column information by MatGetLocalSize() is not
>> consistent with by MatDenseGetLocalMatrix(), is it bug? could you give me
>> some advice? thanks a lot.
>>
> What information?
>
>   Matt
>
>
>> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090106/2e8f8230/attachment.htm>