Can you send this code to petsc-maint at mcs.anl.gov? I tried in C and this runs fine.
Matt On Thu, Jan 29, 2009 at 9:39 AM, Panyasantisuk Jarunan <Jarunan.Panyasantisuk at eleves.ec-nantes.fr> wrote: > Hi, > > Is anybody using MatCreateMPIAIJWithSplitArrays()? > I got the error below, eventhough, the matrix is built but not completed. > > [0]PETSC ERROR: --------------- Error Message --------------- > [0]PETSC ERROR: Argument out of range! > [0]PETSC ERROR: Column entry number 1 (actual colum 0) in row 1 is not > sorted! > > I think there is a problem with column indices but I thought that it has > been sorted, please see below. > > I am using this command > > call MatCreateMPIAIJWithSplitArrays(PETSC_COMM_WORLD,2,2, > $ PETSC_DETERMINE,PETSC_DETERMINE,pointer,column,v,opointer, > $ ocolumn,ov,D,ierr) > > for processor 0: > pointer = [0 2 4] ! row pointer into column column = [0 1 0 1] > ! local column index > v = [5.49 -2.74 -2.74 5.49] ! diagonal value > > off-diagonal processor 0: opointer = [0 1 3] ocolumn = [2 3] > ov = [-2.74 -2.74] > > The matrix: > row 0: (0, 5.49395) (1, -2.74697) (2, -2.74697) > row 1: (0, 0) (1, 5.49395) (3, -2.74697) > > I am using fortran with petsc-3.0.0. > > > Thank you, > Jarunan > > > > -- > Jarunan PANYASANTISUK > MSc. in Computational Mechanics > Erasmus Mundus Master Program > Ecole Centrale de Nantes > 1, rue de la no?, 44321 NANTES, FRANCE > > > > > > > -- 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
