[ET Trac] [Einstein Toolkit] #1222: Reduction weight for periodic domains.

Einstein Toolkit trac-noreply at einsteintoolkit.org
Thu Jan 17 02:55:17 CST 2013


#1222: Reduction weight for periodic domains.
-------------------------+--------------------------------------------------
  Reporter:  bentivegna  |       Owner:  eschnett           
      Type:  defect      |      Status:  new                
  Priority:  minor       |   Milestone:                     
 Component:  Carpet      |     Version:  development version
Resolution:              |    Keywords:                     
-------------------------+--------------------------------------------------

Comment (by bentivegna):

 If I understand your argument correctly, I do not agree. Let's say we have
 a periodic domain where we identify x=+1 with x=-1; we cover this with two
 grids, one with spacing dx=0.2 and one with dx/2, and with a lower
 shiftout of 1 and an upper shiftout of 0. With three boundary points, this
 will look like this:

 {{{
            -1                                      +1
             |                                       |
             |                                       |
             |                                       |
 B   B   B   I   I   I   I   I   I   I   I   I   I   B   B   B
       b b b i i i i i i i i i i i i i i i i i i i i b b b
 }}}

 The coordinate of the last internal point will indeed be different on the
 two levels (0.8 on the coarse and 0.9 on the fine), but things like the
 volume of the two levels (assigning weight of 1 to the interior points and
 0 otherwise) is the same:

 C: #I*dx=10*0.2=2
 F: #i*dx/2=20*0.1=2

 Notice that the volume of the two levels is the same with the current
 weight assignment too  -- the patched code produces identical results as
 the current one, when the boundaries have been correctly set. It only
 makes a difference when the boundaries are junk, in which case it prevents
 this junk from polluting the reduction operations.

-- 
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1222#comment:4>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit


More information about the Trac mailing list