distribute_mask<TITLE> <META NAME="description" CONTENT="distribute_mask"> <META NAME="keywords" CONTENT="S2HATdocs"> <META NAME="resource-type" CONTENT="document"> <META NAME="distribution" CONTENT="global"> <META NAME="Generator" CONTENT="LaTeX2HTML v2002-2-1"> <META HTTP-EQUIV="Content-Style-Type" CONTENT="text/css"> <LINK REL="STYLESHEET" HREF="S2HATdocs.css"> </HEAD> <BODY >  <BR> <B> Next:</B> <A NAME="tex2html231" HREF="./CdistrWghts.html">distribute_w8ring</A> <B> Up:</B> <A NAME="tex2html229" HREF="./CgathScat.html">Gather/scatter routines</A> <B> Previous:</B> <A NAME="tex2html223" HREF="./CdistrPartMap.html">distribute_partialmap</A> <BR> <BR>  <HR> <table width=750 border=0 cellpadding=1 cellspacing=0> <tbody> <tr><td> <H3><A NAME="SECTION00015300000000000000"> distribute_map</A> </H3> subroutine <tt>distribute_mask( pixeltype pixelization, int nmasks, int masknum, int ncomp, int first_ring, int last_ring, int mask_size, int *local_mask, int *mask, int myid, int numprocs, int root, MPI_Comm comm) </tt><BR> <BR> <U>Purpose:</U><br><br> The routine distributes a full sky mask <tt><b>mask</b></tt> made of <tt><b>ncomp</b></tt> components stored in the memory of the proc <tt><b>root</b></tt> over all processors of the communicator <tt><b>comm</b></tt>. The result is stored in a vector <tt><b>local_mask</b></tt> as a <tt><b>masknum</b></tt> mask. To distribute <tt><b>nmasks</b></tt> the routines needs to be called <tt><b>nmasks</b></tt> times, with <tt><b>masknum</b></tt> changing from <tt>0</tt> to <tt>nmask-1</tt>. In this way the need for storing all the <tt><b>nmasks</b></tt> full sky masks in the memory of a single proc is avoided. <BR> The mask is assumed to be a vector of 0 and 1 integers delineating the observed patch of the sky. For the float-valued map-like objects (windows, apodization) use <tt><b>distribute_map</b></tt>. <BR> To be called by all procs of the communicator <tt><b>comm</b></tt>. Blocking. <BR> <BR> <BR><U>Parameters:</U> <BR> <BR> input: <DL> <DD><tt><b>pixelization(*)--</b> [pixeltype]</tt>, an S2HAT structure containing the info about the pixelization; </DD> <BR> <DD><tt><b>masknum(*)--</b> [int]</tt>, a current number of a mask which is being now distributed; </DD> <BR> <DD><tt><b>nmasks(*)--</b></tt> [int], a total number of masks to be distributed and later processed; </DD> <br> <DD><tt><b>ncomp(*)--</b> [int]</tt>, a number of components of the masks to be distributed <tt>(=1,2,3)</tt>; </DD> <BR> <DD><tt><b>first_ring, last_ring--</b> [int]</tt>, a number of a first and last ring which define a distibuted map (both included) to be stored on a given proc; </DD> <BR> <DD><tt><b>mask_size--</b> [int]</tt>, a number of pixels corresponding to a set of rings defined by <tt>first_ring</tt> and <tt>last_ring</tt> above; </DD> <BR> <DD><tt><b>mask--</b> [int*]</tt>, a vector storing to a 2dim array of dimensions <tt>(0:npixall-1, 1:ncomp)</tt> in the column-wise order. It needs to be defined and allocated only on the proc <tt>root</tt> and it will be distributed over all procs of the communicator <tt>comm</tt> as a result of the routine. </DD> <BR> <DD><tt><b>myid--</b> [int]</tt>, proc id <tt>(0,...,numprocs-1)</tt>; </DD> <BR> <DD><tt><b>numprocs(*)--</b> [int]</tt>, total number of procs used; </DD> <BR> <DD><tt><b>root(*)--</b> [int]</tt>, defines a proc <tt>root</tt> on which the full mask <tt>mask</tt> is stored; </DD> <BR> <DD><tt><b>comm(*)--</b> [MPI_Comm]</tt>, an mpi communicator for the numprocs procs. </DD> </DL> <font size = 2pt> (*)The value must be the same on all procs of the communicator <tt>comm</tt>. </font> <P> output: <DL> <DD><tt><b>local_mask--</b> [int*]</tt>, a vector storing a 3dim array of a size <tt>(0:mask_size-1,1:ncomp,1:nmasks)</tt> in the column-wise order to contain the distributed mask. </DD> </DL> <U>Comments:</U> <BR><BR><DD> <LI>On the proc <tt>root</tt> a full sky mask is always assumed, even if only parts of it (as defined by <tt>first_ring</tt> and <tt>last_ring</tt>) may actually be distributed. </DD><br> <DD><LI>On all procs the sub-mask, <tt>local_mask,</tt> is made of complete rings of both northern and southern hemispheres with numbers falling between <tt>first_ring</tt> and the <tt>last_ring</tt>, inclusive. </DD> <BR> </tr></td> </tbody> </table> <HR> <BR> <BR> <B> Next:</B> <A NAME="tex2html231" HREF="./CdistrWghts.html">distribute_w8ring</A> <B> Up:</B> <A NAME="tex2html229" HREF="./CgathScat.html">Gather/scatter routines</A> <B> Previous:</B> <A NAME="tex2html223" HREF="./CdistrPartMap.html">distribute_partialmap</A> <ADDRESS> <BR> radek stompor 2009-10-15 </ADDRESS> </BODY> </HTML>