<html xmlns="http://www.w3.org/1999/xhtml" xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office"><head><!--[if gte mso 9]><xml><o:OfficeDocumentSettings><o:AllowPNG/><o:PixelsPerInch>96</o:PixelsPerInch></o:OfficeDocumentSettings></xml><![endif]--></head><body><div style="font-family:garamond, new york, times, serif;font-size:13px;"><div>Hello PetSc users,</div><div><br></div><div>1) I have problem with SNES convergence. I call SNESSolve in a time loop and use the inbuilt Jacobian feature. The code works fine for about 5380 time steps after which it breaks down. The solution till that point looks fine. I have used newtonls of type l2. (newtontr and others aren't working). Since I have used <span>inbuilt Jacobian feature and the code worked for about 5000 time steps I don't understand the reason for failure, is it an incorrect function evaluation? Attached are the outputs with -pc_type lu and ilu along with -snes_linesearch_type l2 -snes_converged_reason -snes_monitor -snes_view -ksp_converged_reason -ksp_monitor_true_residual</span></div><div><span><br></span></div><div><span>2) How to get hold of failure signal, like Nonlinear solve <span>DIVERGED_MAX_IT or <span>DIVERGED_LINEAR_SOLVE so that whenever it occurs I can use a reduced time step and see if the code converges.</span></span></span></div><div><span><span><span><br></span></span></span></div><div><span><span><span>Thank you,</span></span></span></div><div><span><span><span>Rahul.<br></span></span></span></div><div><span><span><span></span></span></span><br></div><div><a href="https://drive.google.com/file/d/12bAiJm0i2NEmv7NXld6LSuzUVpNrLwSA/view?usp=sharing" rel="nofollow" target="_blank" class="enhancr_card_5702848584">output_ilu.txt</a></div><div><br></div><div id="ydpd045cdb9enhancr_card_5702848584" class="ydpd045cdb9yahoo-link-enhancr-card ydpd045cdb9ymail-preserve-class ydpd045cdb9ymail-preserve-style" style="max-width:400px;font-family:"Helvetica Neue", "Segoe UI", Helvetica, Arial, sans-serif" data-url="https://drive.google.com/file/d/12bAiJm0i2NEmv7NXld6LSuzUVpNrLwSA/view?usp=sharing" data-type="YENHANCER" data-size="MEDIUM" contenteditable="false"><a href="https://drive.google.com/file/d/12bAiJm0i2NEmv7NXld6LSuzUVpNrLwSA/view?usp=sharing" style="text-decoration:none !important;color:#000 !important" class="ydpd045cdb9yahoo-enhancr-cardlink" rel="nofollow" target="_blank"><table class="ydpd045cdb9card-wrapper ydpd045cdb9yahoo-ignore-table" style="max-width:400px" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td width="400"><table class="ydpd045cdb9card ydpd045cdb9yahoo-ignore-table" style="max-width:400px;border-width:1px;border-style:solid;border-color:rgb(224, 228, 233);border-radius:2px" width="100%" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td class="ydpd045cdb9card-primary-image-cell" style="background-color: rgb(0, 0, 0); background-repeat: no-repeat; background-size: cover; position: relative; border-radius: 2px 2px 0px 0px; min-height: 175px;" valign="top" height="175" bgcolor="#000000" background="https://s.yimg.com/lo/api/res/1.2/rUpWl8RM2SABokN5iaWULg--~A/Zmk9ZmlsbDt3PTQwMDtoPTIwMDthcHBpZD1pZXh0cmFjdA--/https://lh5.googleusercontent.com/Z0ks0ak03LaePVJ3I5X1ejFVwqVcOrul6f2cvmuWXjdOn-PMJbR2mYhihLA=w1200-h630-p.cf.jpg"><!--[if gte mso 9]><v:rect fill="true" stroke="false" style="width:396px;height:175px;position:absolute;top:0;left:0;"><v:fill type="frame" color="#000000" src="https://s.yimg.com/lo/api/res/1.2/rUpWl8RM2SABokN5iaWULg--~A/Zmk9ZmlsbDt3PTQwMDtoPTIwMDthcHBpZD1pZXh0cmFjdA--/https://lh5.googleusercontent.com/Z0ks0ak03LaePVJ3I5X1ejFVwqVcOrul6f2cvmuWXjdOn-PMJbR2mYhihLA=w1200-h630-p.cf.jpg"/></v:rect><![endif]--><table class="ydpd045cdb9card-overlay-container-table ydpd045cdb9yahoo-ignore-table" style="width:100%" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td class="ydpd045cdb9card-overlay-cell" style="background-color: transparent; border-radius: 2px 2px 0px 0px; min-height: 175px;" valign="top" bgcolor="transparent" background="https://s.yimg.com/cv/ae/nq/storm/assets/enhancrV21/1/enhancr_gradient-400x175.png"><!--[if gte mso 9]><v:rect fill="true" stroke="false" style="width:396px;height:175px;position:absolute;top:-18px;left:0;"><v:fill type="pattern" color="#000000" src="https://s.yimg.com/cv/ae/nq/storm/assets/enhancrV21/1/enhancr_gradient-400x175.png"/><v:textbox inset="0,0,20px,0"><![endif]--><table class="ydpd045cdb9yahoo-ignore-table" style="width: 100%; min-height: 175px;" height="175" border="0"><tbody><tr><td class="ydpd045cdb9card-richInfo2" style="text-align:left;padding:15px 0 0 15px;vertical-align:top"></td><td class="ydpd045cdb9card-actions" style="text-align:right;padding:15px 15px 0 0;vertical-align:top"><div class="ydpd045cdb9card-share-container"></div></td></tr></tbody></table><!--[if gte mso 9]></v:textbox></v:rect><![endif]--></td></tr></tbody></table></td></tr><tr><td><table class="ydpd045cdb9card-info ydpd045cdb9yahoo-ignore-table" style="background-color: rgb(255, 255, 255); background-repeat: repeat; background-attachment: scroll; background-image: none; background-size: auto auto; position: relative; z-index: 2; width: 100%; max-width: 400px; border-radius: 0px 0px 2px 2px; border-top: 1px solid rgb(224, 228, 233);" cellspacing="0" cellpadding="0" border="0" align="center"><tbody><tr><td style="background-color:#ffffff;padding:16px 0 16px 12px;vertical-align:top;border-radius:0 0 0 2px"><img class="ydpd045cdb9card-object-1 ydpd045cdb9yahoo-ignore-inline-image ydpd045cdb9ymail-preserve-class" src="https://s.yimg.com/nq/storm/assets/enhancrV2/23/logos/google.png" style="min-width:36px;margin-top:3px" height="36"></td><td style="vertical-align:middle;padding:12px 24px 16px 12px;width:99%;font-family:"Helvetica Neue", "Segoe UI", Helvetica, Arial, sans-serif;border-radius:0 0 2px 0"><h2 class="ydpd045cdb9card-title" style="font-size:14px;line-height:19px;margin:0 0 6px 0;font-family:"Helvetica Neue", "Segoe UI", Helvetica, Arial, sans-serif;word-break:break-word;color:#26282a">output_ilu.txt</h2><p class="ydpd045cdb9card-description" style="font-size:12px;line-height:16px;margin:0;color:#979ba7;word-break:break-word"></p></td></tr></tbody></table></td></tr></tbody></table></td></tr></tbody></table></a></div><div><br></div><div><br></div><div><br></div><div><a href="https://drive.google.com/file/d/1CbQuLLIVtkACo6eedMmMQvX41TvjtoiQ/view?usp=sharing" rel="nofollow" target="_blank" class="enhancr_card_6437182933">output_lu.txt</a></div><div><br></div><div id="ydp352c65f3enhancr_card_6437182933" class="ydp352c65f3yahoo-link-enhancr-card ydp352c65f3ymail-preserve-class ydp352c65f3ymail-preserve-style" style="max-width:400px;font-family:"Helvetica Neue", "Segoe UI", Helvetica, Arial, sans-serif" data-url="https://drive.google.com/file/d/1CbQuLLIVtkACo6eedMmMQvX41TvjtoiQ/view?usp=sharing" data-type="YENHANCER" data-size="MEDIUM" contenteditable="false"><a href="https://drive.google.com/file/d/1CbQuLLIVtkACo6eedMmMQvX41TvjtoiQ/view?usp=sharing" style="text-decoration:none !important;color:#000 !important" class="ydp352c65f3yahoo-enhancr-cardlink" rel="nofollow" target="_blank"><table class="ydp352c65f3card-wrapper ydp352c65f3yahoo-ignore-table" style="max-width:400px" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td width="400"><table class="ydp352c65f3card ydp352c65f3yahoo-ignore-table" style="max-width:400px;border-width:1px;border-style:solid;border-color:rgb(224, 228, 233);border-radius:2px" width="100%" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td class="ydp352c65f3card-primary-image-cell" style="background-color: rgb(0, 0, 0); background-repeat: no-repeat; background-size: cover; position: relative; border-radius: 2px 2px 0px 0px; min-height: 175px;" valign="top" height="175" bgcolor="#000000" background="https://s.yimg.com/lo/api/res/1.2/HgnG8J__r3zsgWOMnL.MYg--~A/Zmk9ZmlsbDt3PTQwMDtoPTIwMDthcHBpZD1pZXh0cmFjdA--/https://lh4.googleusercontent.com/y7DZz1myxbsXhzL9zsNqDJ57VQlARxD3TCVrOA0TA19WFihCo31AfrlORHQ=w1200-h630-p.cf.jpg"><!--[if gte mso 9]><v:rect fill="true" stroke="false" style="width:396px;height:175px;position:absolute;top:0;left:0;"><v:fill type="frame" color="#000000" src="https://s.yimg.com/lo/api/res/1.2/HgnG8J__r3zsgWOMnL.MYg--~A/Zmk9ZmlsbDt3PTQwMDtoPTIwMDthcHBpZD1pZXh0cmFjdA--/https://lh4.googleusercontent.com/y7DZz1myxbsXhzL9zsNqDJ57VQlARxD3TCVrOA0TA19WFihCo31AfrlORHQ=w1200-h630-p.cf.jpg"/></v:rect><![endif]--><table class="ydp352c65f3card-overlay-container-table ydp352c65f3yahoo-ignore-table" style="width:100%" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td class="ydp352c65f3card-overlay-cell" style="background-color: transparent; border-radius: 2px 2px 0px 0px; min-height: 175px;" valign="top" bgcolor="transparent" background="https://s.yimg.com/cv/ae/nq/storm/assets/enhancrV21/1/enhancr_gradient-400x175.png"><!--[if gte mso 9]><v:rect fill="true" stroke="false" style="width:396px;height:175px;position:absolute;top:-18px;left:0;"><v:fill type="pattern" color="#000000" src="https://s.yimg.com/cv/ae/nq/storm/assets/enhancrV21/1/enhancr_gradient-400x175.png"/><v:textbox inset="0,0,20px,0"><![endif]--><table class="ydp352c65f3yahoo-ignore-table" style="width: 100%; min-height: 175px;" height="175" border="0"><tbody><tr><td class="ydp352c65f3card-richInfo2" style="text-align:left;padding:15px 0 0 15px;vertical-align:top"></td><td class="ydp352c65f3card-actions" style="text-align:right;padding:15px 15px 0 0;vertical-align:top"><div class="ydp352c65f3card-share-container"></div></td></tr></tbody></table><!--[if gte mso 9]></v:textbox></v:rect><![endif]--></td></tr></tbody></table></td></tr><tr><td><table class="ydp352c65f3card-info ydp352c65f3yahoo-ignore-table" style="background-color: rgb(255, 255, 255); background-repeat: repeat; background-attachment: scroll; background-image: none; background-size: auto auto; position: relative; z-index: 2; width: 100%; max-width: 400px; border-radius: 0px 0px 2px 2px; border-top: 1px solid rgb(224, 228, 233);" cellspacing="0" cellpadding="0" border="0" align="center"><tbody><tr><td style="background-color:#ffffff;padding:16px 0 16px 12px;vertical-align:top;border-radius:0 0 0 2px"><img class="ydp352c65f3card-object-1 ydp352c65f3yahoo-ignore-inline-image ydp352c65f3ymail-preserve-class" src="https://s.yimg.com/nq/storm/assets/enhancrV2/23/logos/google.png" style="min-width:36px;margin-top:3px" height="36"></td><td style="vertical-align:middle;padding:12px 24px 16px 12px;width:99%;font-family:"Helvetica Neue", "Segoe UI", Helvetica, Arial, sans-serif;border-radius:0 0 2px 0"><h2 class="ydp352c65f3card-title" style="font-size:14px;line-height:19px;margin:0 0 6px 0;font-family:"Helvetica Neue", "Segoe UI", Helvetica, Arial, sans-serif;word-break:break-word;color:#26282a">output_lu.txt</h2><p class="ydp352c65f3card-description" style="font-size:12px;line-height:16px;margin:0;color:#979ba7;word-break:break-word"></p></td></tr></tbody></table></td></tr></tbody></table></td></tr></tbody></table></a></div><div><br></div><div><br></div></div></body></html>