Re: [petsc-users] Weird behavior for log_summary

2015-06-08 Thread Michele Rosso
Is there any external software I should link at compile time to allow PETSc to provide log infos? Also, is there any define in petscconf.h that controls profiling? Michele On Mon, 2015-06-08 at 12:48 -0500, Barry Smith wrote: On Jun 8, 2015, at 12:16 PM, Michele Rosso mro...@uci.edu wrote:

Re: [petsc-users] Weird behavior for log_summary

2015-06-08 Thread Jed Brown
Michele Rosso mro...@uci.edu writes: Is there any external software I should link at compile time to allow PETSc to provide log infos? Also, is there any define in petscconf.h that controls profiling? There is PETSC_USE_LOG, but that should always be on, in which case everything is run-time

Re: [petsc-users] Guidance on GAMG preconditioning

2015-06-08 Thread Young, Matthew, Adam
Thanks for your helpful comments, Jed, Matt, and Mark. They've given me some things to think about, so I'll work on those then open a new email thread when I have further questions. --Matt -- Matthew Young Graduate Student Boston

Re: [petsc-users] Weird behavior for log_summary

2015-06-08 Thread Barry Smith
On Jun 8, 2015, at 12:53 PM, Michele Rosso mro...@uci.edu wrote: Is there any external software I should link at compile time to allow PETSc to provide log infos? No Also, is there any define in petscconf.h that controls profiling? Yes, but it would stop anything from being printed,

Re: [petsc-users] Weird behavior for log_summary

2015-06-08 Thread Michele Rosso
Jed, In the petscconf.h I have #ifndef PETSC_USE_LOG #define PETSC_USE_LOG 1 #endif so I guess that is not the problem. I run ex50: I attached the output. It does prints the summary. The I guess there is something wrong with my code. I call mpi_init before petsc_initialize and then I finalize

Re: [petsc-users] Guidance on GAMG preconditioning

2015-06-08 Thread Jed Brown
Young, Matthew, Adam m...@bu.edu writes: This is a problem from ionospheric plasma physics. The simulation treats ions via a particle-in-cell method and electrons as an inertialess fluid, the justification being that ionospheric ions are 10^4 times more massive than electrons. We further

Re: [petsc-users] Weird behavior for log_summary

2015-06-08 Thread Michele Rosso
Hi Barry, I run a small test case like you suggested: this results in no error, but the problem with log_summary still persists. Please find attached the output of log_summary. Thanks, Michele On Fri, 2015-06-05 at 21:34 -0500, Barry Smith wrote: [NID 04001] 2015-06-04 19:07:24 Apid 25022256:

Re: [petsc-users] DMs not transferred into PCCOMPOSITE?

2015-06-08 Thread Matthew Knepley
On Mon, Jun 8, 2015 at 4:31 AM, Lawrence Mitchell lawrence.mitch...@imperial.ac.uk wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi all, I have a multi-field system where the splits are defined by a DM. I can use a fieldsplit PC on its own with no trouble, but cannot do so inside

Re: [petsc-users] Tao iterations

2015-06-08 Thread Justin Chang
Hi Jason, One more question about BLMVM... if it only uses gradient information and does not require the definitition of a Hessian Matrix, can this method be applied to solve problems that are nonsymmetric by nature? (e.g., advection-diffusion equations). If I had wanted to solve the same

Re: [petsc-users] Guidance on GAMG preconditioning

2015-06-08 Thread Matthew Knepley
On Sun, Jun 7, 2015 at 9:08 PM, Justin Chang jychan...@gmail.com wrote: Matt (Knepley), I see what you're saying and it makes perfect sense. The point of my work isn't necessarily to compare CG/Jacobi with GAMG. Rather I am trying to compare both the numerical solution and the computational

[petsc-users] DMs not transferred into PCCOMPOSITE?

2015-06-08 Thread Lawrence Mitchell
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi all, I have a multi-field system where the splits are defined by a DM. I can use a fieldsplit PC on its own with no trouble, but cannot do so inside a composite PC. My matrices do not have a blocksize set (because the data layout doesn't match)

Re: [petsc-users] Tao iterations

2015-06-08 Thread Jed Brown
Justin Chang jychan...@gmail.com writes: Hi Jason, One more question about BLMVM... if it only uses gradient information and does not require the definitition of a Hessian Matrix, can this method be applied to solve problems that are nonsymmetric by nature? (e.g., advection-diffusion

Re: [petsc-users] petsc4py Build Problem

2015-06-08 Thread Lisandro Dalcin
On 8 June 2015 at 02:50, Mikhail Khodak mkho...@princeton.edu wrote: Hello, I am trying to build petsc4py-3.5.1 using Cygwin on 64-bit Windows 7. I have built PETSc 3.5.4 with shared and dynamic libraries using mpich2-1.2.1 and successfully ran the installation tests. I am using Python 2.7 and

Re: [petsc-users] Tao iterations

2015-06-08 Thread Justin Chang
Jed, Thank you for your response. I agree completely with all that you said. I just wonder what would happen if I attempted to use the TAO routines after forming the Jacobian J and residual r arising from the advection diffusion equation. In my current (linear) diffusion framework, i have the

Re: [petsc-users] Tao iterations

2015-06-08 Thread Jed Brown
Justin Chang jychan...@gmail.com writes: Jed, Thank you for your response. I agree completely with all that you said. I just wonder what would happen if I attempted to use the TAO routines after forming the Jacobian J and residual r arising from the advection diffusion equation. In my

Re: [petsc-users] Tao iterations

2015-06-08 Thread Justin Chang
Ah I see that makes sense. Thank you very much On Monday, June 8, 2015, Jed Brown j...@jedbrown.org wrote: Justin Chang jychan...@gmail.com javascript:; writes: Jed, Thank you for your response. I agree completely with all that you said. I just wonder what would happen if I attempted