/[escript]/trunk/finley/src/NodeFile.h
ViewVC logotype

Annotation of /trunk/finley/src/NodeFile.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 751 - (hide annotations)
Mon Jun 26 01:46:34 2006 UTC (13 years, 5 months ago) by bcumming
File MIME type: text/plain
File size: 5259 byte(s)
Changes relating to the MPI version of escript
The standard OpenMP version of escript is unchanged

- updated data types (Finley_Mesh, Finley_NodeFile, etc) to store meshes
  over multiple MPI processes.
- added CommBuffer code in Paso for communication of Data associated
  with distributed meshes
- updates in Finley and Escript to support distributed data and operations
  on distributed data (such as interpolation).
- construction of RHS in MPI, so that simple explicit schemes (such as
  /docs/examples/wave.py without IO and the Locator) can run in MPI.
- updated mesh generation for first order line, rectangle and brick
  meshes and second order line meshes in MPI.        
- small changes to trunk/SConstruct and trunk/scons/ess_options.py to
  build the MPI version, these changes are turned off by default.

1 jgs 150 /*
2 elspeth 626 ************************************************************
3     * Copyright 2006 by ACcESS MNRF *
4     * *
5     * http://www.access.edu.au *
6     * Primary Business: Queensland, Australia *
7     * Licensed under the Open Software License version 3.0 *
8     * http://www.opensource.org/licenses/osl-3.0.php *
9     * *
10     ************************************************************
11 jgs 150 */
12 jgs 82
13 jgs 150 /* Version: $Id$ */
14    
15 jgs 82 #ifndef INC_FINLEY_NODEFILE
16     #define INC_FINLEY_NODEFILE
17    
18     #define MAX_numDim 3
19    
20 jgs 150 #include "Finley.h"
21 robwdcock 682 #include "escript/DataC.h"
22 jgs 82
23 bcumming 730 #ifdef PASO_MPI
24     #include "Distribution.h"
25 bcumming 751 #include "./paso/CommBuffer.h"
26 bcumming 730 #endif
27    
28 jgs 82 struct Finley_NodeFile {
29 bcumming 730 #ifdef PASO_MPI
30     Paso_MPIInfo *MPIInfo; /* MPI information */
31     Finley_NodeDistribution *degreeOfFreedomDistribution; /* information about the distribution of degrees of freedom
32     on this subdomain and over other subdomains */
33 bcumming 751 Finley_NodeDistribution *reducedDegreeOfFreedomDistribution;
34     Paso_CommBuffer *CommBuffer;
35     Paso_CommBuffer *reducedCommBuffer;
36 bcumming 730 #endif
37 jgs 123 dim_t numNodes; /* number of nodes */
38     dim_t numDim; /* spatial dimension */
39     index_t *Id; /* Id[i] is the id number of node i. this number is not really
40 jgs 82 used but useful when the nodes are
41     relabled. see also Finley_resolveNodeIds. */
42 jgs 123 /* in the entire code the term 'node id' refers to i but nor to Id[i] */
43     /* if not explicitly stated otherwise. */
44     index_t *Tag; /* Tag[i] is the tag of node i. */
45     double *Coordinates; /* Coordinates[INDEX2(k,i,numDim)] is the k-th coordinate of the */
46     /* node i. */
47     index_t* degreeOfFreedom; /* degreeOfFreedom[i] is the degree of freedom assigned to node i */
48     /* this index is used to consider periodic boundary conditions by assigning */
49     /* the same degreeOfFreedom to the same node */
50 jgs 82 /* the following data are set by Finley_NodeFile_? */
51 jgs 123 dim_t numDegreesOfFreedom; /* number of degrees of freedom in the mesh (<=numNodes)*/
52     /* notice that numDegreesOfFreedom=max(degreeOfFreedom[:]) */
53     index_t* reducedDegreeOfFreedom; /* reducedDegreeOfFreedom[i] is the degree of freedom in the reduced version of the mesh */
54     /* assigned to node i. reducedDegreeOfFreedom[i]<0 indicates that the node is not appearing */
55     /* as a degree of freedom in the reduced version of the mesh. */
56     dim_t reducedNumDegreesOfFreedom; /* number of degrees of freedom in the reduced version of the mesh */
57     /* notice that reducedNumDegreesOfFreedom=max(reducedDegreeOfFreedom[:]) */
58 jgs 82 /* this information is used by interfaces to postprocessing tools that can deal with linear elements only: */
59 jgs 123 dim_t reducedNumNodes; /* returnes the number of nodes in the mesh when looking at the corners */
60     /* of elements only */
61     index_t *toReduced; /* toReduced[i] is the node id in the reduced mesh. if toReduced[i]<0 it means that the node does not appear in the reduced mesh */
62 jgs 82 };
63    
64 bcumming 730
65    
66 jgs 82 typedef struct Finley_NodeFile Finley_NodeFile;
67    
68 bcumming 730 #ifdef PASO_MPI
69     Finley_NodeFile* Finley_NodeFile_alloc(dim_t, Paso_MPIInfo *MPIInfo);
70     void Finley_NodeFile_allocTable(Finley_NodeFile*,dim_t);
71     #else
72     Finley_NodeFile* Finley_NodeFile_alloc(dim_t numDim);
73     void Finley_NodeFile_allocTable(Finley_NodeFile*,dim_t);
74     #endif
75    
76     /* Finley_NodeFile * Finley_NodeFile_alloc(dim_t); */
77 jgs 82 void Finley_NodeFile_dealloc(Finley_NodeFile*);
78 jgs 123 void Finley_NodeFile_setIdRange(index_t*,index_t*,Finley_NodeFile*);
79     void Finley_NodeFile_scatter(index_t*,Finley_NodeFile*,Finley_NodeFile*);
80     void Finley_NodeFile_gather(index_t*,Finley_NodeFile*,Finley_NodeFile*);
81 jgs 82 void Finley_NodeFile_setCoordinates(Finley_NodeFile*,escriptDataC*);
82 jgs 123 void Finley_NodeFile_copyTable(dim_t,Finley_NodeFile*,dim_t,dim_t,Finley_NodeFile*);
83 bcumming 730
84 jgs 82 void Finley_NodeFile_deallocTable(Finley_NodeFile*);
85    
86     #endif
87    
88     /*
89     * $Log$
90 jgs 150 * Revision 1.3 2005/09/15 03:44:23 jgs
91     * Merge of development branch dev-02 back to main trunk on 2005-09-15
92     *
93     * Revision 1.2.2.1 2005/09/07 06:26:20 gross
94     * the solver from finley are put into the standalone package paso now
95     *
96 jgs 123 * Revision 1.2 2005/07/08 04:07:55 jgs
97     * Merge of development branch back to main trunk on 2005-07-08
98 jgs 82 *
99 jgs 123 * Revision 1.1.1.1.2.1 2005/06/29 02:34:54 gross
100     * some changes towards 64 integers in finley
101     *
102     * Revision 1.1.1.1 2004/10/26 06:53:57 jgs
103     * initial import of project esys2
104     *
105 jgs 82 * Revision 1.2 2004/07/02 04:21:13 gross
106     * Finley C code has been included
107     *
108     * Revision 1.1.1.1 2004/06/24 04:00:40 johng
109     * Initial version of eys using boost-python.
110     *
111     *
112     */

Properties

Name Value
svn:eol-style native
svn:keywords Author Date Id Revision

  ViewVC Help
Powered by ViewVC 1.1.26