Thanks Matt, it's clear now On 6/19/12 6:49 PM, Matthew Knepley wrote: > On Tue, Jun 19, 2012 at 10:06 AM, Anton Popov <popov at uni-mainz.de > <mailto:popov at uni-mainz.de>> wrote: > > Hi petsc team! > I've recently upgraded to 3.3, and currently experience problems. > One of them follows. > If I compile this simple code: > > > In 3.3, we now require that MatSetUp() is called before messing with > the storage scheme. > We should have been more specific in Changes. Moreover, this should be > a better error > message. Will fix, and see below. > > #include "petsc.h" > #undef __FUNCT__ > #define __FUNCT__ "main" > int main( int argc,char **argv ) > { > Mat A; > PetscErrorCode ierr; > PetscInt n_rows_Local=100; > PetscInt n_cols_Local=100; > ierr = PetscInitialize(&argc, &argv, PETSC_NULL, PETSC_NULL); > CHKERRQ(ierr); > ierr = MatCreate(PETSC_COMM_WORLD, &A); CHKERRQ(ierr); > ierr = MatSetSizes(A, n_rows_Local, n_cols_Local, > PETSC_DETERMINE, PETSC_DETERMINE); CHKERRQ(ierr); > ierr = MatSetType(A, MATAIJ); CHKERRQ(ierr); > > > ierr = MatSetUp(A);CHKERRQ(ierr); > > Matt > > ierr = MatSetOption(A, MAT_KEEP_NONZERO_PATTERN, PETSC_TRUE); > CHKERRQ(ierr); > ierr = PetscFinalize(); CHKERRQ(ierr); > PetscFunctionReturn(0); > } > > and run it in parallel (mpiexec -np 2 ./example) it gives me > errors like this: > > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [0]PETSC ERROR: Null argument, when expecting valid pointer! > [0]PETSC ERROR: Null Object: Parameter # 1! > [0]PETSC ERROR: > ------------------------------------------------------------------------ > > Commenting the line with "MatSetOption" removes the error, but > doesn't solve the problem, because I actually need to > KEEP_NONZERO_PATTERN > > Running the code sequentially, also vanishes the error. > > Please check whether it's a bug in petsc-3.3, or I should set this > option in a different way. > > Thanks, > > Anton Popov > > > > > -- > 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/20120620/451dc046/attachment-0001.html>
