<div dir="ltr">Hi Klaus,<div><br></div><div>  That's exactly what was happening! I've implemented your ChildOrder suggestion (switching HYPRE and Multigrid ranking), and now the potential is correct, and I am not having the non-convergence issues anymore. </div><div><br></div><div>  Also, thank you for the comments on the boundary conditions. I've always been a bit uncertain how these can be paired together, and your email cleared that up. </div><div><br></div><div>Take care,</div><div>  -Mark</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 3 March 2016 at 19:45, Klaus Weide <span dir="ltr"><<a href="mailto:klaus@flash.uchicago.edu" target="_blank">klaus@flash.uchicago.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Mark,<br>
<br>
On going back to one of your older messages, I noticed something<br>
suspicious in your log file excerpt:<br>
<div><div class="h5"><br>
<br>
On Mon, 15 Feb 2016, Mark Richardson wrote:<br>
> I'm having one other problem, which is maybe related. Originally I was<br>
> running this with a self-made Eos (Tillotson for planetary bodies, that<br>
> includes radiative transfer). I was getting the following errors<br>
> immediately after initialisation:<br>
><br>
>  [ 02-15-2016  09:57:48.810 ] [gr_hypreInit]: HYPRE_RELEASE_NAME=hypre,<br>
> HYPRE_RELEASE_VERSION=2.10.0b, HYPRE_RELEASE_DATE=2015/01/22<br>
><br>
>  [ 02-15-2016  09:57:48.812 ] memory: /proc vsize    (MB):     2006.58<br>
> (min)       2008.35 (max)       2007.17 (avg)<br>
>  [ 02-15-2016  09:57:48.813 ] memory: /proc rss      (MB):     1284.28<br>
> (min)       1286.47 (max)       1285.21 (avg)<br>
>  [ 02-15-2016  09:57:48.815 ] memory: /proc vsize    (MB):     2006.58<br>
> (min)       2008.35 (max)       2007.17 (avg)<br>
>  [ 02-15-2016  09:57:48.816 ] memory: /proc rss      (MB):     1284.37<br>
> (min)       1286.56 (max)       1285.29 (avg)<br>
><br>
>  [ 02-15-2016  09:57:58.820 ] [gr_hypreSolve]: Nonconv./failure: ierr=<br>
>       0<br>
>  [ 02-15-2016  09:57:58.820 ] [gr_hypreSolve]: Nonconv./failure: component=<br>
>           0<br>
>  [ 02-15-2016  09:57:58.820 ] [gr_hypreSolve]: Nonconv./failure: converged=<br>
> F<br>
>  [ 02-15-2016  09:57:58.820 ] [gr_hypreSolve]: Nonconv./failure: |initial<br>
> RHS|^2=  1.2520E+35<br>
>  [ 02-15-2016  09:57:58.820 ] [gr_hypreSolve]: Nonconv./failure:<br>
> num_iterations=         500<br>
>  [ 02-15-2016  09:57:58.820 ] [gr_hypreSolve]: Nonconv./failure:<br>
> final_res_norm=   88143026.2676523<br>
<br>
</div></div>The fact that these messages appear this early - before writing of the<br>
initial checkpoint file, before the message "Entering evolution loop" -<br>
indicates that they come from some call from Driver_initFlash, not<br>
Driver_evolveFlash. Unless you have made significant changes, there should<br>
be no calls to Diffuse_solveScalar from Driver_initFlash (directly or<br>
indirectly). But there IS an initial call to<br>
Gravity_potentialListOfBlocks. So I strongly suspect that your setup is<br>
using GridSolvers/HYPRE/Grid_solvePoisson.F90 as the Poisson solver<br>
implementation, rather than the intended (I assume)<br>
GridSolvers/Multigrid/Grid_solvePoisson.F90.<br>
<br>
Q1. Am I correct to assume that you are intending to use Multigrid as<br>
Poisson solver for Gravity?<br>
<br>
Q2. Can you confirm that Grid_solvePoisson.F90 in your object directory<br>
was pointing to GridSolvers/HYPRE/Grid_solvePoisson.F90, not<br>
GridSolvers/Multigrid/Grid_solvePoisson.F90 as it should?<br>
<br>
Q3. Is the situation still the same in your most recent configuration?<br>
<br>
If the answers are Y/Y/Y, the solution may be simple - try the following:<br>
in source/Grid/GridSolvers/Config, change the line that begins with<br>
<br>
 CHILDORDER<br>
<br>
to<br>
<br>
 CHILDORDER Pfft HYPRE Multigrid<br>
<br>
The important part here is changing the relative order of Multigrid and<br>
HYPRE, so Multigrid comes later.<br>
(I have omitted various alternative implementations from the directive<br>
that are not included in the released code anyway.)<br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
Klaus<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><br>Mark Richardson<br><div>Beecroft Postdoctoral Fellow, Department of Physics<div>Denys Wilkinson Building - 555F</div><div>University of Oxford<div><a href="mailto:Mark.Richardson@physics.ox.ac.uk" target="_blank">Mark.Richardson@physics.ox.ac.uk</a><br><div>+44 1865 280793</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div>