/[MITgcm]/manual/s_phys_pkgs/text/exch2.tex

Diff of /manual/s_phys_pkgs/text/exch2.tex

Parent Directory | Revision Log | View Revision Graph Revision Graph | View Patch Patch

-revision 1.14 by afe,
Wed Mar 17 21:44:02 2004 UTC
+revision 1.15 by afe,
Thu Mar 18 14:56:25 2004 UTC
 Line 281 
 generated by \file{driver.m}.\\
  The scalar parameters \varlink{exch2\_domain\_nxt}{exch2_domain_nxt}
  and \varlink{exch2\_domain\_nyt}{exch2_domain_nyt} express the number
  of tiles in the $x$ and $y$ global indices.  For example, the default
- setup of six tiles (Fig. \ref{fig:6tile}) has \code{exch2\_domain\_nxt=6} and
+ setup of six tiles (Fig. \ref{fig:6tile}) has
- \code{exch2\_domain\_nyt=1}.  A topology of twenty-four square tiles,
+ \code{exch2\_domain\_nxt=6} and \code{exch2\_domain\_nyt=1}.  A
- four per subdomain (as in figure \ref{fig:24tile}), will have
+ topology of twenty-four square tiles, four per subdomain (as in figure
- \code{exch2\_domain\_nxt=12} and \code{exch2\_domain\_nyt=2}.  Note
+ \ref{fig:24tile}), will have \code{exch2\_domain\_nxt=12} and
- that these parameters express the tile layout to allow global data
+ \code{exch2\_domain\_nyt=2}.  Note that these parameters express the
- files that are tile-layout-neutral and have no bearing on the internal
+ tile layout to allow global data files that are tile-layout-neutral
- storage of the arrays.  The tiles are internally stored in a range
+ and have no bearing on the internal storage of the arrays.  The tiles
- from \code{(1:\varlink{bi}{bi})} the $x$ axis, and $y$ axis variable
+ are internally stored in a range from \code{(1:\varlink{bi}{bi})} the
- \varlink{bj}{bj} is generally ignored within the package. \\
+ $x$ axis, and the $y$ axis variable \varlink{bj}{bj} is generally
+ ignored within the package. \\
  \subsubsection{Arrays Indexed to Tile Number}
- The following arrays are of length \code{NTILES}, are indexed to the
+ The following arrays are of length \code{NTILES}and are indexed to the
- tile number, and the indices are omitted in their descriptions. \\
+ tile number, which is indicated in the diagrams with the notation
+ \textsf{t}$n$.  The indices are omitted in the descriptions. \\
  The arrays \varlink{exch2\_tnx}{exch2_tnx} and
  \varlink{exch2\_tny}{exch2_tny} express the $x$ and $y$ dimensions of
-Line 324 
 the subdomain of each tile, in a range \
+Line 326 
 the subdomain of each tile, in a range \
  standard cube topology and indicated by \textbf{\textsf{f}}$n$ in
  figures \ref{fig:12tile} and
  \ref{fig:24tile}. \varlink{exch2\_nNeighbours}{exch2_nNeighbours}
- contains a count of how many neighboring tiles each tile has, and is
+ contains a count the  neighboring tiles each tile has, and is
  used for setting bounds for looping over neighboring tiles.
  \varlink{exch2\_tProc}{exch2_tProc} holds the process rank of each
  tile, and is used in interprocess communication.  \\
-Line 334 
 The arrays \varlink{exch2\_isWedge}{exch
+Line 336 
 The arrays \varlink{exch2\_isWedge}{exch
  \varlink{exch2\_isEedge}{exch2_isEedge},
  \varlink{exch2\_isSedge}{exch2_isSedge}, and
  \varlink{exch2\_isNedge}{exch2_isNedge} are set to \code{1} if the
- indexed tile lies on the edge of a subdomain, \code{0} if not.  The
+ indexed tile lies on the respective edge of a subdomain, \code{0} if
- values are used within the topology generator to determine the
+ not.  The values are used within the topology generator to determine
- orientation of neighboring tiles, and to indicate whether a tile lies
+ the orientation of neighboring tiles, and to indicate whether a tile
- on the corner of a subdomain.  The latter case requires special
+ lies on the corner of a subdomain.  The latter case requires special
  exchange and numerical handling for the singularities at the eight
  corners of the cube. \\
-Line 350 
 orientations between the the tiles. \\
+Line 352 
 orientations between the the tiles. \\
  The array \code{exch2\_neighbourId(a,T)} holds the tile number
  \code{Tn} for each of the tile number \code{T}'s neighboring tiles
- \code{a}.  The neighbor tiles are indexed \code{(1:MAX\_NEIGHBOURS)}
+ \code{a}.  The neighbor tiles are indexed
- in the order right to left on the north then south edges, and then top
+ \code{(1:exch2\_NNeighbours(T))} in the order right to left on the
- to bottom on the east and west edges.  Maybe throw in a fig here, eh?
+ north then south edges, and then top to bottom on the east and west
- \\
+ edges.  Maybe throw in a fig here, eh?  \\
- \sloppy
+ \sloppy The \code{exch2\_opposingSend\_record(a,T)} array holds the
- The \code{exch2\_opposingSend\_record(a,T)} array holds the index
+ index \code{b} of the element in \texttt{exch2\_neighbourId(b,Tn)}
- \code{b} in \texttt{exch2\_neighbourId(b,Tn)} that holds the tile
+ that holds the tile number \code{T}, given
- number \code{T}.  In other words,
+ \code{Tn=exch2\_neighborId(a,T)}.  In other words,
  \begin{verbatim}
     exch2_neighbourId( exch2_opposingSend_record(a,T),
                        exch2_neighbourId(a,T) ) = T
-Line 366 
 number \code{T}.  In other words,
+Line 368 
 number \code{T}.  In other words,
  This provides a back-reference from the neighbor tiles. \\
  The arrays \varlink{exch2\_pi}{exch2_pi} and
- \varlink{exch2\_pj}{exch2_pj} specify the transformations of variables
+ \varlink{exch2\_pj}{exch2_pj} specify the transformations of indices
  in exchanges between the neighboring tiles.  These transformations are
- necessary in exchanges between subdomains because a physical vector
+ necessary in exchanges between subdomains because the array index in
- component in one direction may map to one in a different direction in
+ one dimension may map to the other index in an adjacent subdomain, and
- an adjacent subdomain, and may be have its indexing reversed. This
+ may be have its indexing reversed. This swapping arises from the
- swapping arises from the ``folding'' of two-dimensional arrays into a
+ ``folding'' of two-dimensional arrays into a three-dimensional cube.
- three-dimensional cube.
  The dimensions of \code{exch2\_pi(t,N,T)} and \code{exch2\_pj(t,N,T)}
  are the neighbor ID \code{N} and the tile number \code{T} as explained
- above, plus a vector of length 2 containing transformation factors
+ above, plus a vector of length \code{2} containing transformation
- \code{t}.  The first element of the transformation vector indicates
+ factors \code{t}.  The first element of the transformation vector
- the factor \code{t} by which variables representing the same
+ holds the factor to multiply the index in the same axis, and the
- \emph{physical} vector component of a tile \code{T} will be multiplied
+ second element holds the the same for the orthogonal index.  To
- in exchanges with neighbor \code{N}, and the second element indicates
+ clarify, \code{exch2\_pi(1,N,T)} holds the mapping of the $x$ axis
- the transform to the physical vector in the other direction.  To
+ index of tile \code{T} to the $x$ axis of tile \code{T}'s neighbor
- clarify (hopefully), \code{exch2\_pi(1,N,T)} holds the transform of
+ \code{N}, and \code{exch2\_pi(2,N,T)} holds the mapping of \code{T}'s
- the $i$ component of a vector variable in tile \code{T} to the $i$
+ $x$ index to the neighbor \code{N}'s $y$ index. \\
- component of tile \code{T}'s neighbor \code{N}, and
- \code{exch2\_pi(2,N,T)} holds the transform of \code{T}'s $i$
- components to the neighbor \code{N}'s $j$ component. \\
- Under the current cube topology, one of the two elements of
+ One of the two elements of \code{exch2\_pi} or \code{exch2\_pj} for a
- \code{exch2\_pi} or \code{exch2\_pj} for a given tile \code{T} and
+ given tile \code{T} and neighbor \code{N} will be \code{0}, reflecting
- neighbor \code{N} will be \code{0}, reflecting the fact that the two
+ the fact that the two axes are orthogonal.  The other element will be
- vector components are orthogonal.  The other element will be \code{1}
+ \code{1} or \code{-1}, depending on whether the axes are indexed in
- or \code{-1}, depending on whether the components are indexed in the
+ the same or opposite directions.  For example, the transform vector of
- same or opposite directions.  For example, the transform vector of the
+ the arrays for all tile neighbors on the same subdomain will be
- arrays for all tile neighbors on the same subdomain will be
  \code{(1,0)}, since all tiles on the same subdomain are oriented
- identically.  A vector direction that corresponds to the orthogonal
+ identically.  An axis that corresponds to the orthogonal dimension
- dimension with the same index direction in a particular tile-neighbor
+ with the same index direction in a particular tile-neighbor
- orientation will have \code{(0,1)}, whereas those in the opposite
+ orientation will have \code{(0,1)}.  Those in the opposite index
- index direction will have \code{(0,-1)}. \\
+ direction will have \code{(0,-1)} in order to reverse the ordering. \\
  The arrays \varlink{exch2\_oi}{exch2_oi},
  \varlink{exch2\_oj}{exch2_oj}, \varlink{exch2\_oi\_f}{exch2_oi_f}, and

 Legend:



Removed from v.1.14
 


changed lines


 
Added in v.1.15
 Legend:



Removed from v.1.14
 


changed lines


 
Added in v.1.15
-Removed from v.1.14
+Added in v.1.15

	ViewVC Help
Powered by ViewVC 1.1.22