[Users] Schedule change in Exact.

Roland Haas roland.haas at physics.gatech.edu
Tue May 4 22:23:14 CDT 2010


Hello Eloisa,

> But is this actually correct? Roland's example with "BSSNtoADM" being 
> naturally scheduled in PseudoEvolution makes it clear that there are at 
> least two levels of pseudo-evolution: one that translates the 
> gravitational/hydro variables between all the different formulations one 
> may need, and a later one where these variables (cast in whichever 
> formulation) are used to build derived objects. Also from Roland's 
> comments, I understand that the latter class must remain in 
> PseudoEvolution (rather than going in ANALYSIS which I would have found 
> more natural) because derivatives calculated in ANALYSIS will be wrong 
> at the refinement boundaries (did I get this point right?). 
Yes this is my understanding. There was some discussion about this wrt 
to WeylScal4 (between Tanja, Erik, Ian, Peter, myself, around 
2009-09-18). Essentially in ANALYSIS you get:

--8<--
What I wanted to say is that if the Psi4 grid function is set during 
ANALYSIS then the fine and coarse grid levels are not set in the proper 
order. To test this I wrote a tiny thorn with two functions that are 
both scheduled in ANALYSIS, one which sets a grid function to 
cctk_iteration+lvl/10, where lvl is the current refinement level and a 
second function to write one value from this grid function to stdout 
(see attachement). The output of those (for a six level system) looks 
like this:

#      iteration i*j*k lvl grid-function
debug: 0         0     0   0.000000
debug: 0         0     1   0.100000
debug: 0         0     2   0.200000
debug: 0         0     3   0.300000
debug: 0         0     4   0.400000
debug: 0         0     5   0.500000
debug: 0         0     6   0.600000
debug: 1         0     6   1.600000 <-
debug: 2         0     5   2.500000
debug: 2         0     6   2.600000


Those ones with iteration number 0 are just initial data and are fine. 
If prolongation was to be used in ANALYSIS and it works as in EVOL then 
in step 1 (marked with an <- above) on level 6 one would need values 
from step 2 of level 5 which at this point have not yet been calculated, 
resulting in poison (or old values) being prolongated into the slice.
--8<--
Essentially the time interpolation part of prolongation fails during 
ANALYSIS (so computing only every_coarse is almost fine [no restriction 
from fine to coarse happens]).

> Is there a 
> way to solve this other problem perhaps?
No idea.

Yours,
Roland

-- 
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://keys.gnupg.net.


More information about the Users mailing list