[Users] ADMConstraints

[Erik Schnetter]

After testing for some time, I just applied corresponding changes to  
ADMAnalysis.

-erik

On May 31, 2007, at 07:56:00, Bela Szilagyi wrote:
> Should there be a similar patch applied to ADMAnalysis as well?
>
>
> On Thursday 31 May 2007 03:58:46 Erik Schnetter wrote:
>> While calculating the constraints with ADMConstraints in the presence
>> of mesh refinement works fine, there was a problem with calculating
>> norms of the constraints, which may require time interpolation.  In
>> order to do this correctly, it is necessary to keep several time
>> levels of the constraints, and to re-calculate the constraints after
>> changing the grid hierarchy.  If this is not done, then the values of
>> the constraints at each output grid point are correct, but norms of
>> the constraints may not be.
>>
>> I updated ADMConstraints so that it supports this behaviour.  In
>> order to remain backward compatible, this behaviour cannot be the
>> default.  Obtaining correct norms for the constraints requires
>> setting the parameters
>>
>> ADMConstraints::constraints_persist = yes
>> ADMConstraints::constraints_timelevels = 3
>>
>> where the number 3 depend on the time interpolation order; for linear
>> time interpolation, you only need 2 time levels.
>>
>> On the other hand, if you want to save space and time, and if you do
>> not need the norms of the constraints, then you should set the  
>> parameter
>>
>> ADMConstraints::constraints_prolongation_type = "none"
>>
>> -erik
>
>
>
>
>
> Bela Szilagyi
> ----------------------------------------------------
> Max-Planck-Institut für Gravitationsphysik
> Albert-Einstein-Institut
> Tel: +49 331 567 7189
> Fax: +49 331 567 7252
> ----------------------------------------------------
>
>
> _______________________________________________
> Users mailing list
> Users at cactuscode.org
> http://www.cactuscode.org/mailman/listinfo/users
>


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

My email is as private as my paper mail.  I therefore support encrypting
and signing email messages.  Get my PGP key from www.keyserver.net.



-------------- next part --------------
A non-text attachment was scrubbed...
Name: PGP.sig
Type: application/pgp-signature
Size: 194 bytes
Desc: This is a digitally signed message part
Url : http://www.cactuscode.org/pipermail/users/attachments/20080406/19340fb3/attachment.bin