[petsc-users] Galerkin multigrid coarsening

2014-01-30 Thread Boris Kaus
Hi, While doing multigrid it is possible to explicitly set restriction and prolongation/interpolation operators in PETSC using PCMGSetRestriction and PCMGSetInterpolation. In many cases, the restriction (R) and prolongation (P) operators are simply the transpose of each other (R=P’). Yet,

Re: [petsc-users] Galerkin multigrid coarsening

2014-01-30 Thread Jed Brown
Boris Kaus k...@uni-mainz.de writes: While doing multigrid it is possible to explicitly set restriction and prolongation/interpolation operators in PETSC using PCMGSetRestriction and PCMGSetInterpolation. In many cases, the restriction (R) and prolongation (P) operators are simply the

Re: [petsc-users] Configuration of Hybrid MPI-OpenMP

2014-01-30 Thread Danyang Su
I made a second check on initialization of PETSc, and found that the initialization does not take effect. The codes are as follows. call PetscInitialize(Petsc_Null_Character,ierrcode) call MPI_Comm_rank(Petsc_Comm_World,rank,ierrcode) call

Re: [petsc-users] Configuration of Hybrid MPI-OpenMP

2014-01-30 Thread Jed Brown
Danyang Su danyang...@gmail.com writes: I made a second check on initialization of PETSc, and found that the initialization does not take effect. The codes are as follows. call PetscInitialize(Petsc_Null_Character,ierrcode) call

Re: [petsc-users] Galerkin multigrid coarsening

2014-01-30 Thread Kaus, Boris
Thanks Jed, I'll give it a try and let you know. Might take a bit though. Boris On Jan 30, 2014, at 17:45, Jed Brown j...@jedbrown.org wrote: Boris Kaus k...@uni-mainz.de writes: While doing multigrid it is possible to explicitly set restriction and prolongation/interpolation

[petsc-users] Input arguments of DMPlexCreateFromDAG

2014-01-30 Thread Cedric Doucet
Hello, I tried to use DMPlexCreateFromDAG function to create a DM structure from a hybrid mesh. To understand how it works, I look at ex5.c file. Unfortunately, there are some things that I do not understand. 1. Why does coneSize lists faces of cells first, then faces of vertices and finally

Re: [petsc-users] Input arguments of DMPlexCreateFromDAG

2014-01-30 Thread Matthew Knepley
On Thu, Jan 30, 2014 at 12:19 PM, Cedric Doucet cedric.dou...@inria.frwrote: Hello, I tried to use DMPlexCreateFromDAG function to create a DM structure from a hybrid mesh. To understand how it works, I look at ex5.c file. Unfortunately, there are some things that I do not understand. 1.

Re: [petsc-users] Configuration of Hybrid MPI-OpenMP

2014-01-30 Thread Danyang Su
On 30/01/2014 9:30 AM, Jed Brown wrote: Danyang Su danyang...@gmail.com writes: I made a second check on initialization of PETSc, and found that the initialization does not take effect. The codes are as follows. call PetscInitialize(Petsc_Null_Character,ierrcode) call

Re: [petsc-users] Galerkin multigrid coarsening

2014-01-30 Thread Torquil Macdonald Sørensen
On 30/01/14 17:45, Jed Brown wrote: Boris Kaus k...@uni-mainz.de writes: While doing multigrid it is possible to explicitly set restriction and prolongation/interpolation operators in PETSC using PCMGSetRestriction and PCMGSetInterpolation. In many cases, the restriction (R) and

Re: [petsc-users] Galerkin multigrid coarsening

2014-01-30 Thread Jed Brown
Torquil Macdonald Sørensen torq...@gmail.com writes: This question was similar to my recent unanswered one: http://lists.mcs.anl.gov/pipermail/petsc-users/2014-January/020382.html That's why I Cc'd you. I just tried jed/pcmg-galerkin-rap, and it seems to work fineI verified that the coarse

Re: [petsc-users] Galerkin multigrid coarsening

2014-01-30 Thread Torquil Macdonald Sørensen
On 30/01/14 22:05, Jed Brown wrote: Torquil Macdonald Sørensen torq...@gmail.com writes: This question was similar to my recent unanswered one: http://lists.mcs.anl.gov/pipermail/petsc-users/2014-January/020382.html That's why I Cc'd you. Thanks, I didn't notice that. I am actually very

[petsc-users] saving log info

2014-01-30 Thread Mohammad Mirzadeh
Hi guys, I'd like to be able to post-process the values reported by -log_summary. Is there any programmatic way to extract the data and save it manually to a file? If not, is there any parser that go over stdout dump and extract the info? Mohammad

Re: [petsc-users] Prometheus vs GAMG for elasticity/plasticity problems

2014-01-30 Thread Mark Adams
Try running prometheus with -out_verbose 2 and gamg with -pc_gamg_verbose 2 and send me the output. Mark On Thu, Jan 30, 2014 at 4:40 AM, Thomas Gross tgr...@ilsb.tuwien.ac.atwrote: Please find enclosed the output for GAMG using -mg_levels_ksp_max_it 1: Prometheus: -ksp_type cg -pc_type

Re: [petsc-users] Prometheus vs GAMG for elasticity/plasticity problems

2014-01-30 Thread Mark Adams
And is this just a simple cube problem in elasticity as it looks like from your input file name? On Thu, Jan 30, 2014 at 4:40 AM, Thomas Gross tgr...@ilsb.tuwien.ac.atwrote: Please find enclosed the output for GAMG using -mg_levels_ksp_max_it 1: Prometheus: -ksp_type cg -pc_type prometheus

Re: [petsc-users] saving log info

2014-01-30 Thread Jed Brown
Mohammad Mirzadeh mirza...@gmail.com writes: Hi guys, I'd like to be able to post-process the values reported by -log_summary. Is there any programmatic way to extract the data and save it manually to a file? If not, is there any parser that go over stdout dump and extract the info? I

Re: [petsc-users] saving log info

2014-01-30 Thread Mohammad Mirzadeh
Thanks Jed. This looks interesting. On Thu, Jan 30, 2014 at 2:45 PM, Jed Brown j...@jedbrown.org wrote: Mohammad Mirzadeh mirza...@gmail.com writes: Hi guys, I'd like to be able to post-process the values reported by -log_summary. Is there any programmatic way to extract the data and

Re: [petsc-users] Prometheus vs GAMG for elasticity/plasticity problems

2014-01-30 Thread Mark Adams
We're seeing a bit lower performance in MatMult with GAMG, perhaps because we are not using block formats specialized for elasticity. The block info seems to be there. Mark, what else is different? * Prometheus seems to be coarsening slower and taking nearly twice the iterations. *

Re: [petsc-users] Prometheus vs GAMG for elasticity/plasticity problems

2014-01-30 Thread Jed Brown
Mark Adams mfad...@lbl.gov writes: We're seeing a bit lower performance in MatMult with GAMG, perhaps because we are not using block formats specialized for elasticity. The block info seems to be there. I meant that GAMG isn't using BAIJ, but Prometheus is. That affects kernel efficiency.

Re: [petsc-users] saving log info

2014-01-30 Thread Barry Smith
master now has a viewer that dumps the raw data in Python, then you can write python scripts to process it any way you want. -log_view ascii:filename:ascii_info_detail sample python script in bin/pythonscripts/petsclogformat.py I personally think that parsing the output of -log_summary

Re: [petsc-users] Prometheus vs GAMG for elasticity/plasticity problems

2014-01-30 Thread Mark Adams
I meant that GAMG isn't using BAIJ, but Prometheus is. That affects kernel efficiency. Ah yes, that must be it.

[petsc-users] How can PETSc configure with MSMPI

2014-01-30 Thread Danyang Su
Hi All, When configure petsc with msmpi ./configure --with-cc='win32fe cl' --with-fc='win32fe ifort' --with-cxx='win32fe cl' --download-f-blas-lapack --with-threadco mm --with-openmp --with-mpi-include=/cygdrive/c/Program Files/Microsoft HPC Pack 2008 R2/Inc --with-mpi-lib=/cygdrive/c/Prog

Re: [petsc-users] saving log info

2014-01-30 Thread Matthew Knepley
On Thu, Jan 30, 2014 at 5:56 PM, Barry Smith bsm...@mcs.anl.gov wrote: master now has a viewer that dumps the raw data in Python, then you can write python scripts to process it any way you want. -log_view ascii:filename:ascii_info_detail sample python script in