[Users] changing CarpetRegrid2 parameters for ETK Maxwell + mercurial Carpet?

[Erik Schnetter]

Bernard

I am looking at your report now.

needrecv is a bboxset, not an integer. Can you instead output it via

cout << "neecrecv=" << needrecv << "\n";

? This will output its content. Since the routine expects it (the set
of grid points) to be empty at this point, but it is not, it is
important to know what points are left over -- outer boundary, ghost
points, refinement boundary, ...

Your long output (with output_bboxes etc., all output turned on) is
fine. The more output the better!

-erik

On Thu, Feb 2, 2012 at 2:45 PM, Bernard Kelly <bernard.j.kelly at nasa.gov> wrote:
> Hi all (again).
>
> I mailed the ETK list back in November (!) because of a problem with our
> BHB evolutions using ETK "Maxwell"; said evolutions crashed after only a
> few M of evolution. Erik made some suggestions, I supplied more
> diagnostic information (output_bboxes=yes), Erik looked it over again,
> and suggested a further diagnostic message to nail down what the problem
> was.
>
> I thought I was doing this, when the problem magically disappeared, so I
> didn't get back to the list.
>
> Unfortunately, it turned out that what had happened was I accidentally
> went back to ETK "Curie" (and hence my older, non-crashing problem).
> ["Accidentally" here means I was led astray by some symbolic links I
> created to point to our own Hahndol thorns within the Cactus
> arrangements directory.] Now that I've discovered my own stupidity, I
> still need ysome help. I've implemented Erik's suggested diagnostic
> message in CarpetLib. Here's my new code, inserted at
> line 948 of the repository CarpetLib/src/dh.cc:
>
> --------------------------------------------------------------
> 948,954d947
> <          /* [2011-11-21] diagnostic code suggested by Erik */
> <          if (not(needrecv.empty())) {
> <            cout << "rl=" << rl << " c=" << c << " needrecv=" << needrecv
> << endl;
> <            CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, "rl=%d,
> c=%d, needrecv=%d\n",rl,c,needrecv);
> <          }
> <          /* end of diagnostic code */
> --------------------------------------------------------------
>
>
> I've run both with and without ouput_bboxes, in a relatively
> low-resolution run. I've again put the screen output files on webdrive
> again, along with the source parameter file. Here's the URL to access it:
>
> https://webdrive.gsfc.nasa.gov/longauth/600/bernard.j.kelly/QSDKYC3
>
> As before, you'll need the temporary username (beanykelly) and password
> (helpmercurial).
>
> Thanks in advance -- again -- for your help,
>
> Bernard
>
> On 11/17/11 12:33 PM, Kelly, Bernard J. (GSFC-660.0)[UNIVERSITY OF
> MARYLAND BALTIMORE COUNTY] wrote:
>> Hi all.
>>
>> I followed Erik's suggestion with the extra diagnostic information. The
>> resulting screen output is about 200 MB in size, but zips down to 20MB.
>>
>> I've placed it on Goddard's "webdrive" for anyone to download. Here's the
>> link:
>>
>> https://webdrive.gsfc.nasa.gov/shortauth/bernard.j.kelly/jvXqUYC
>>
>> You'll need the temporary username (beanykelly) and password
>> (helpmercurial) to access it.
>>
>> All suggestions (especially the blindingly obvious) are appreciated.
>>
>> Thanks, Bernard
>>
>>
>> ---------------------------------------------------------------------------
>
> --
> -----------------------------------------------
> Bernard J. Kelly
>
> NASA Goddard Space Flight Center, Code 663
> 8800 Greenbelt Road
> Greenbelt, MD 20771, U.S.A.
>
> phone: +1 (301) 286-7243
> -----------------------------------------------
> _______________________________________________
> Users mailing list
> Users at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users



-- 
Erik Schnetter <schnetter at cct.lsu.edu>   http://www.cct.lsu.edu/~eschnett/