<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Timo,<div><br></div><div>I just learned a bunch of svn commands from Eric, that allowed me to check for other files</div><div>that were different, and I found two that I had forgotten. And, I found the problem.</div><div><br></div><div>First, I forgot that I had changed the temperature initial conditions, this is why when</div><div>you ran the model it finished in two time-steps.</div><div><br></div><div><div>Here's the change I made:</div><div><div>--- source/initial_conditions/box.cc<span class="Apple-tab-span" style="white-space: pre; "> </span>(revision 1118)</div><div>+++ source/initial_conditions/box.cc<span class="Apple-tab-span" style="white-space: pre; "> </span>(working copy)</div><div>@@ -44,7 +44,8 @@</div><div> double perturbation = 1;</div><div> for (unsigned int d=0; d<dim; ++d)</div><div> perturbation *= std::sin(numbers::PI*position[d]/geometry->get_extents()[d]);</div><div>- return 1 + perturbation/10;</div><div>+ // MIB return 1 + perturbation/10;</div><div>+<span class="Apple-tab-span" style="white-space: pre; "> </span>return 1400*(1 + perturbation/10);</div></div></div><div><br></div><div>(Note, I had to "hard-wire" this 1400 value, because as I understand it, given how Aspect is set-up right now,</div><div>there's no easy way to use information about the boundary conditions in defining the initial conditions. It</div><div>would certainly be helpful to make this easier in the future.).</div><div><br></div><div>Second, in Stokes solver I had changed the accuracy from 1e-12 to 1e-15, (this was leftover from looking</div><div>for ways to decrease the heat flux on the zero heat-flux side-walls, but didn't have any effect).</div><div><br></div><div>- When I change this back to 1e-12, the model runs fine on 4 processors. Its interesting that</div><div>using this smaller accuracy is only a problem on 4 processors, does that make any sense to you?</div><div><div><br></div><div>Thanks for responding to my e-mail…and sorry to have sent you on a wild goose-chase.</div><div>Magali</div><div><br></div><div><br></div><div><div>On Aug 30, 2012, at 1:53 PM, Timo Heister wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear Magali,<br><br><blockquote type="cite">Exception on processing:<br></blockquote><blockquote type="cite">Iterative method reported convergence failure in step 7551 with residual 3.15628e+07<br></blockquote><blockquote type="cite">Aborting!<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">(note it was only on step 59, but is stating it has an error in step 7551 ??)<br></blockquote><br>This means that the linear solver did not converge after 7551 iterations.<br><br>There can be several reasons for this happening. Does this happen in<br>the same time step every time you run it with the same number of<br>cores? Can you increase the number of graphical outputs and look at a<br>solution directly before this crash?<br><br><blockquote type="cite">I've attached the run output file (where the error is), the parameter file,<br></blockquote><blockquote type="cite">and the c and h files for the temperature boundary conditions for the box convection<br></blockquote><blockquote type="cite">(same as earlier this summer).<br></blockquote><br>I tried it here and it runs fine for me with 4 processors. Do you have<br>any other changes by any chance? The code finished in 2 time steps for<br>me:<br>*** Timestep 1: t=3.35488e+09 years<br>*** Timestep 2: t=1.1674e+10 years<br><br>I increased the final time and it worked for a couple hundred time steps.<br><br>Best,<br>Timo<br><br>-- <br>Timo Heister<br><a href="http://www.math.tamu.edu/~heister/">http://www.math.tamu.edu/~heister/</a><br>_______________________________________________<br>Aspect-devel mailing list<br>Aspect-devel@geodynamics.org<br>http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel<br></div></blockquote></div><br><div apple-content-edited="true">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div><div><div><div><div><div><div>--------------------------------------------------</div><div>Associate Professor & Chancellor Fellow, U.C. Davis </div><div>Chair, Geology Graduate Program</div><div>Department of Geology & KeckCAVES</div><div>2129 Earth & Physical Sciences Bldg</div><div>Davis, CA 95616</div><div><br></div><div><div>E-mail: <a href="mailto:mibillen@ucdavis.edu">mibillen@ucdavis.edu</a></div><div>Phone: no extension, please e-mail</div></div><div>--------------------------------------------------</div></div></div></div></div></div><br></div></div></div></span></div></span></div></span></div></span></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br></div></body></html>