On Thu, Oct 31, 2013 at 1:21 PM, Jose David Bermeol <[email protected]>wrote:

> Hi small question. In the method  MatCreateSeqDense(MPI_Comm comm,PetscInt
> m,PetscInt n,PetscScalar *data,Mat *A) I'm giving the local number of rows
> and columns, so I should pass the total number of local columns or I should
> pass columns/number_of_MPI process.
>

If its sequential, you pass the total number of columns. If its
MatCreateMPIDense(), pass the number of local columns. This
is used to create a vector with the same layout as y = A^T x.

  Matt


> Thanks
>



-- 
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

Reply via email to