Re: [Libmesh-users] Output Example 1 FEMSystem

2016-03-19 Thread Paul T. Bauman
Hi Adriano, Thanks for reporting this. The fix for this got put into the "primary" PETSc classes, but not this one (again, this isn't exercised very much). I'll put in a fix. (As an aside: Hopefully in the next few days, I'll have an interface to this in GRINS so we can exercise this in some test

Re: [Libmesh-users] Output Example 1 FEMSystem

2016-03-11 Thread Paul T. Bauman
Hi Adriano, On Fri, Mar 11, 2016 at 1:49 PM, Adriano Côrtes wrote: > Dear Roy, > > I discovered that the example 1 of FEMSystem (NS equation) uses > NewtonSolver, that is, _diff_solver in libMesh::TimeSolver points to > NewtonSolver. Where should I change the code to use PetscDiffSolver? Is > th

Re: [Libmesh-users] Output Example 1 FEMSystem

2016-03-11 Thread Adriano Côrtes
Dear Roy, I discovered that the example 1 of FEMSystem (NS equation) uses NewtonSolver, that is, _diff_solver in libMesh::TimeSolver points to NewtonSolver. Where should I change the code to use PetscDiffSolver? Is there an example where I can take a look? Thank you. Best, Adriano. 2016-03-10 14

Re: [Libmesh-users] Output Example 1 FEMSystem

2016-03-10 Thread Adriano Côrtes
Thank you Roy for the explanations. Best, Adriano. 2016-03-10 14:08 GMT-03:00 Roy Stogner : > > On Thu, 10 Mar 2016, Adriano Côrtes wrote: > > *1. Why only the timestep 0 printed the line bellow?* >> >> Nonlinear step: |du|/|u| = 0.99946, |du| = 194.712 >> > > That's only printed for interm

Re: [Libmesh-users] Output Example 1 FEMSystem

2016-03-10 Thread Adriano Côrtes
Thank you John for you reply and for the issue. Adriano. 2016-03-10 13:54 GMT-03:00 John Peterson : > > > On Thu, Mar 10, 2016 at 9:49 AM, John Peterson > wrote: > >> >> >> On Thu, Mar 10, 2016 at 9:33 AM, Adriano Côrtes >> wrote: >> >>> Dear developers, >>> >>> I am trying to understand the F

Re: [Libmesh-users] Output Example 1 FEMSystem

2016-03-10 Thread Roy Stogner
On Thu, 10 Mar 2016, Adriano Côrtes wrote: *1. Why only the timestep 0 printed the line bellow?* Nonlinear step: |du|/|u| = 0.99946, |du| = 194.712 That's only printed for intermediate nonlinear steps, IIRC. In later time steps the problem is linear enough to converge in one quasi-New

Re: [Libmesh-users] Output Example 1 FEMSystem

2016-03-10 Thread John Peterson
On Thu, Mar 10, 2016 at 9:49 AM, John Peterson wrote: > > > On Thu, Mar 10, 2016 at 9:33 AM, Adriano Côrtes > wrote: > >> Dear developers, >> >> I am trying to understand the FEMSystem framework, and the first thing I >> did was to run the example 1, but the output I'm obtaining from the >> defa

Re: [Libmesh-users] Output Example 1 FEMSystem

2016-03-10 Thread John Peterson
On Thu, Mar 10, 2016 at 9:33 AM, Adriano Côrtes wrote: > Dear developers, > > I am trying to understand the FEMSystem framework, and the first thing I > did was to run the example 1, but the output I'm obtaining from the default > run (using make and run.sh) is completely different from the one i