MatCreateMPIAIJWithSplitArrays

[Panyasantisuk Jarunan]

Oh, I could not use MatCreateMPIAIJWithArrays either but the mechanism 
below works.

      call MatCreate(PETSC_COMM_WORLD,D,ierr)
      call MatSetSizes(D,N,N,PETSC_DETERMINE,PETSC_DETERMINE,
     $     ierr)
      call MatSetType(D,MATMPIAIJ,ierr) ! to set type a parallel matrix
      call MatSetFromOptions(D,ierr)
      call MatMPIAIJSetPreallocationCSR(D,pointer,Column,v,ierr)

Where pointer is start-row indices a
      Column is local column indices
      v is value

Is there the different beteween the start-row indices in 
MatMPIAIJSetPreallocationCSR and row indices in 
MatCreateMPIAIJWithArrays ?



Regards,
Jarunan




Hello,

To define a matrix with arrays, I cannot use 
MatCreateMPIAIJWithSplitArrays in my program which is written in 
Fortran:

      call MatCreateMPIAIJWithSplitArrays(PETSC_COMM_WORLD,N,N,
     $     PETSC_DETERMINE,PETSC_DETERMINE,pointer,column,v,opointer,
     $     oColumn,ov,D,ierr)

The error is
F:246: undefined reference to `matcreatempiaijwithsplitarrays_'

I could use MatCreateMPIAIJWithArrays but the off diagonal values are 
missing with this command.

I would be appreciate for any advice. Thank you before hand.

Regards,
Jarunan




-- 
Jarunan PANYASANTISUK
MSc. in Computational Mechanics
Erasmus Mundus Master Program
Ecole Centrale de Nantes
1, rue de la noë, 44321 NANTES, FRANCE