[petsc-users] tensor-product mesh

2017-11-25 Thread Mohammad Hassan Baghaei
Hi I am going to solve the PDEs on a domain of a circle, which I need to use the polar coordinates. As dm is supposed for Cartesian coordinate system. Is it possible to define a new coordinate system in that or I need to use dmplex for that. Thanks for your great time answering me.

[petsc-users] storing many petsc objects in a single file

2017-11-25 Thread Storm Weiner
Hey there, For simulations, its useful to store the history as a series of state vectors. For simulations with many time-steps it can get annoying to store each state vector as a separate file. It would be useful if there were some way to manage a database of petsc vectors. To save the current

Re: [petsc-users] crash with PCASM in parallel

2017-11-25 Thread Daniel Stone
The PFLOTRAN wiki seems to want us to install a very old version of PETSc: http://documentation.pflotran.org/user_guide/how_to/installation/linux.html#linux-install The relevant part is: " Install PETSc 3.1. Clone petsc and check out the supported version: git clone https://bitbucket.org/petsc/pe

Re: [petsc-users] crash with PCASM in parallel

2017-11-25 Thread Smith, Barry F.
I cannot find that routine ~/Src/petsc ((v3.6.4)) arch-basic $ git grep MatDestroy_MPIBAIJ_MatGetSubmatrices ~/Src/petsc ((v3.6.4)) arch-basic $ git checkout v3.7.5 Previous HEAD position was 401b1b531b... Increase patchlevel to 3.6.4 HEAD is now at b827f1350a... Increase patchlevel to 3.7.5 ~/

Re: [petsc-users] crash with PCASM in parallel

2017-11-25 Thread Daniel Stone
Thanks for the quick response. I tried Valgrind. Apart from a couple of other warnings in other parts of my code, now fixed, it shows the same stack I described: ==22498== Invalid read of size 4 ==22498==at 0x55A5BFF: MatDestroy_MPIBAIJ_MatGetSubmatrices (baijov.c:609) ==22498==by 0x538A20