[petsc-users] Fwd: Sorted CSR Matrix and Multigrid PC.

[Barry Smith]

  Sorry, this was due to our bug, the fortran function for PCGAMGSetType() was wrong. I have fixed this in the maint and master branch of PETSc in the git repository. But you can simply remove the call to PCGAMGSetType() from your code since what you are setting is the default type.

  BTW: there are other problems with your code after that call that you will have to work through.

  Barry

> On Sep 7, 2016, at 8:46 PM, Manuel Valera <mvalera at mail.sdsu.edu> wrote:
> 
> 
> ---------- Forwarded message ----------
> From: Manuel Valera <mvalera at mail.sdsu.edu>
> Date: Wed, Sep 7, 2016 at 6:40 PM
> Subject: Re: [petsc-users] Sorted CSR Matrix and Multigrid PC.
> To: 
> Cc: PETSc users list <petsc-users at mcs.anl.gov>
> 
> 
> Hello,
> 
> I was able to sort the data but the PCGAMG does not seem  to be working. 
> 
> I reconfigured everything from scratch as suggested and updated to the latest PETSc version, same results,
> 
> I get the following error:
> 
> [0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
> [...]
> [0]PETSC ERROR: ---------------------  Stack Frames ------------------------------------
> [0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,
> [0]PETSC ERROR:       INSTEAD the line number of the start of the function
> [0]PETSC ERROR:       is given.
> [0]PETSC ERROR: [0] PetscStrcmp line 524 /home/valera/sergcemv4/bitbucket/serucoamv4/petsc-3.7.3/src/sys/utils/str.c
> [0]PETSC ERROR: [0] PetscFunctionListFind_Private line 352 /home/valera/sergcemv4/bitbucket/serucoamv4/petsc-3.7.3/src/sys/dll/reg.c
> [0]PETSC ERROR: [0] PCGAMGSetType_GAMG line 1157 /home/valera/sergcemv4/bitbucket/serucoamv4/petsc-3.7.3/src/ksp/pc/impls/gamg/gamg.c
> [0]PETSC ERROR: [0] PCGAMGSetType line 1102 /home/valera/sergcemv4/bitbucket/serucoamv4/petsc-3.7.3/src/ksp/pc/impls/gamg/gamg.c
> [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
> ....
> 
> You can find the working program attached, sorry for the large files and messy makefile,
> 
> Many thanks,
> 
> Manuel Valera.
> 
> 
> ​
>  MGpetscSolver.tar.gz
> ​
>