Thanks. I’ll take a look at figuring out what is going on. We keep shrinking the domain beyond anything reasonable.
Barry [0] SNESSolve_NEWTONTR(): fnorm=1826.65, gnorm=1826.65, ynorm=5.56247e-74 [0] SNESSolve_NEWTONTR(): gpred=1826.65, rho=0, delta=1.66874e-74 [0] SNESSolve_NEWTONTR(): Trying again in smaller region [0] SNESSolve_NEWTONTR(): Scaling direction by 7.94052e-81 On Apr 25, 2014, at 8:53 AM, Norihiro Watanabe <[email protected]> wrote: > attached is the output written until I killed the program. Please grep the > file with "Time step: 7" to go to log messages where the problem happens. > > Thanks, > Nori > > > On Fri, Apr 25, 2014 at 3:09 PM, Barry Smith <[email protected]> wrote: > > Run with -info and send ALL the output. -info triggers the printing of all > the messages with from the calls to PetscInfo() thus we’ll be able to see > what is happening to the values. In theory it should always eventually get > out but there code be either a bug in our code or an error in your function > evaluation that prevents it from ever getting out. > > Barry > > On Apr 25, 2014, at 7:59 AM, Norihiro Watanabe <[email protected]> wrote: > > > I mean "it's keep running but not print anything". I guess the program is > > running in while(1) loop after linear solve in SNESSolve_NEWTONTR() in > > src/snes/impls/tr/tr.c. But I'm not sure if the breaking condition (if (rho > > > neP->sigma) break; ) can be always satisfied at the end. > > > > > > > > On Fri, Apr 25, 2014 at 2:49 PM, Barry Smith <[email protected]> wrote: > > > > On Apr 25, 2014, at 7:31 AM, Norihiro Watanabe <[email protected]> wrote: > > > > > Hi, > > > > > > In my simulation, nonlinear solve with the trust regtion method got > > > stagnent after linear solve (see output below). > > > > What do you mean, get stagnant? Does the code just hang, that is keep > > running but not print anything. > > > > > Is it possible that the method goes to inifite loop? > > > > This is not suppose to be possible. > > > > > Is there any parameter to avoid this situation? > > > > You need to determine what it is “hanging” on. Try running with > > -start_in_debugger and when it “hangs” hit control C and type where to > > determine where it is. > > > > Barry > > > > > > > > 0 SNES Function norm 1.828728087153e+03 > > > 0 KSP Residual norm 91.2735 > > > Linear solve converged due to CONVERGED_ITS iterations 1 > > > Linear solve converged due to CONVERGED_RTOL iterations 3 > > > 1 KSP Residual norm 3.42223 > > > Linear solve converged due to CONVERGED_STEP_LENGTH iterations 1 > > > > > > > > > Thank you in advance, > > > Nori > > > > > > > > > > -- > > Norihiro Watanabe > > > > > -- > Norihiro Watanabe > <log.tar.gz>
