profiling PETSc code

[Matthew Knepley]

On 8/1/06, Matt Funk <mafunk at nmsu.edu> wrote:
> Hi,
>
> i don't think it is the mallocs since it says things like:
> [0] MatAssemblyEnd_SeqAIJMatrix size: 2912 X 2912; storage space: 0
> unneeded,2912 used
> [0] MatAssemblyEnd_SeqAIJNumber of mallocs during MatSetValues() is 0

This is only on one processor.

> However, i do get errors. They look like:
> [0]PETSC ERROR: StageLogRegister() line 95 in src/sys/plog/stageLog.c
> [0]PETSC ERROR: Invalid pointer!
> [0]PETSC ERROR: Null Pointer: Parameter # 3!
> [0]PETSC ERROR: PetscLogStageRegister() line 375 in src/sys/plog/plog.c

You gave an invalid pointer to the call. You should have

  int stage;
  PetscLogStageRegister(&stage, "MyStage");

> which happens during the call PETSCInitialize(...);
>
> After that i get an error like:
> [0] PetscCommDuplicateDuplicating a communicator 91 164 max tags = 1073741823
> [0] PetscCommDuplicateUsing internal PETSc communicator 91 164
> [0] PetscCommDuplicateUsing internal PETSc communicator 91 164
> [0] PetscCommDuplicateUsing internal PETSc communicator 91 164
> [0]PETSC ERROR: MatGetVecs() line 6283 in src/mat/interface/matrix.c
> [0]PETSC ERROR: Null argument, when expecting valid pointer!
> [0]PETSC ERROR: Null Object: Parameter # 1!
> [0]PETSC ERROR: KSPGetVecs() line 555 in src/ksp/ksp/interface/iterativ.c
> [0]PETSC ERROR: KSPDefaultGetWork() line 597 in
> src/ksp/ksp/interface/iterativ.c
> [0]PETSC ERROR: KSPSetUp_CG() line 75 in src/ksp/ksp/impls/cg/cg.c
> [0]PETSC ERROR: KSPSetUp() line 198 in src/ksp/ksp/interface/itfunc.c
>
> so i suppose that is a problem. I am just not sure what it means.
> any ideas?

It looks like you have not called KSPSetOperators().

   Matt

-- 
"Failure has a thousand explanations. Success doesn't need one" -- Sir
Alec Guiness