[Users] Evolving grid arrays using MoL

[Ian Hinder]

On 23 Sep 2013, at 18:06, Peter Diener <diener at cct.lsu.edu> wrote:

> Hi,
> 
> No, you're right. The examples I mentioned are using either
> uni-grid or multi-patch without AMR. I didn't think of the complication of adding AMR to the mix.

So unless anything has changed, the latest information I have is from Erik's email from September 2007:

> On 28 Sep 2007, at 18:10, Erik Schnetter <schnetter at cct.lsu.edu> wrote:
> 
>> I am now evolving geodesics with the following work-around:
>> 
>> I use MoL.  Each time the RHS is calculated, I check whether the current refinement level is the finest level.  If so, the RHS is calculated, otherwise the RHS is set to zero.
>> 
>> The current refinement level is accessed by
>> 
>> interface.ccl:
>> USES INCLUDE: carpet.hh
>> 
>> #include <carpet.hh>
>> 
>> if (Carpet::reflevel == Carpet::reflevels - 1) // do calculation
>> 
>> This works only from C++.  The RHS calculating routine has to be scheduled in level mode.



However, I get stuck well before this stage.  I try to register the array variable with MoL, but it complains that it is not an array.  

-- 
Ian Hinder
http://numrel.aei.mpg.de/people/hinder

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 203 bytes
Desc: Message signed with OpenPGP using GPGMail
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20130924/0767109a/attachment.bin