problem with patchDomain

[Julian C. Cummings]

Hi All,

This is probably most relevant for Scott and Stephen.
I almost have my new Particles SpatialLayout with
NewFields working, but there is a problem that occurs
with a cell-centered Field.  One of the things I need
to do is check whether a particle should be on the 
current patch by examining its position and comparing
it with a bounding box around the domain of this patch.
The bounding box that I create stretches one-half cell
beyond the first and last field points on this patch.

To get the domain of the current patch, I've been saying

layout().patchDomain(lid);

where layout() returns the FieldLayout and lid is the 
local id number of the current patch.  The problem is 
that the patch domain I get handed is always in terms 
of the vertex domain (i.e., the mesh domain) because 
that is all that the FieldLayout knows about.  This is
a problem if the actual Field is cell-centered (or 
something else besides vertex-centered).

Would it make sense to put these patchDomain() functions
into the Field interface, so that the Field could check
its centering and provide a properly adjusted patch
domain?  I guess an alternative for me is to try taking 
a view of the current patch of the Field and then ask
that patch-view for its domain.  Does that work, or do
you get a zero-based domain when you view a Field patch?

Julian C.


Dr. Julian C. Cummings
Staff Scientist, CACR/Caltech
(626) 395-2543
cummings at cacr.caltech.edu