/[escript]/temp_trunk_copy/finley/src/Assemble.h
ViewVC logotype

Diff of /temp_trunk_copy/finley/src/Assemble.h

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

trunk/esys2/finley/src/finleyC/Assemble.h revision 82 by jgs, Tue Oct 26 06:53:54 2004 UTC trunk/finley/src/Assemble.h revision 1312 by ksteube, Mon Sep 24 06:18:44 2007 UTC
# Line 1  Line 1 
1    
2  /* $Id$ */  /* $Id$ */
3    
4  #ifndef INC_FINLEY_ASSEMBLE  /*******************************************************
5  #define INC_FINLEY_ASSEMBLE   *
6     *           Copyright 2003-2007 by ACceSS MNRF
7     *       Copyright 2007 by University of Queensland
8     *
9     *                http://esscc.uq.edu.au
10     *        Primary Business: Queensland, Australia
11     *  Licensed under the Open Software License version 3.0
12     *     http://www.opensource.org/licenses/osl-3.0.php
13     *
14     *******************************************************/
15    
16  /**************************************************************/  /**************************************************************/
17    
# Line 9  Line 19 
19    
20  /**************************************************************/  /**************************************************************/
21    
22  /*   Copyrights by ACcESS Australia, 2003,2004 */  #ifndef INC_FINLEY_ASSEMBLE
23  /*   author: gross@access.edu.au */  #define INC_FINLEY_ASSEMBLE
 /*   Version: $Id$ */  
24    
25  /**************************************************************/  /**************************************************************/
26    
27  #include "ReferenceElements.h"  #include "ReferenceElements.h"
28  #include "System.h"  #include "Finley.h"
29  #include "ElementFile.h"  #include "ElementFile.h"
30  #include "NodeFile.h"  #include "NodeFile.h"
31  #include "escript/Data/DataC.h"  #include "escript/DataC.h"
32    #include "paso/SystemMatrix.h"
33    
34  struct Assemble_Parameters {  struct Assemble_Parameters {
35     int numQuad;     dim_t numQuad;
36     int numDim;     dim_t numDim;
37     int numElementDim;     dim_t NN;
38    
39     int NN;     dim_t numEqu;
40     int NS;     index_t* row_DOF;
41     Finley_RefElement* referenceElement;     dim_t row_DOF_UpperBound;
42       Finley_ElementFile_Jacobeans* row_jac;
43     int numEqu;     index_t* row_node;
44     maybelong* label_row;     dim_t row_NN;
45     Finley_RefElement* referenceElement_row;     dim_t row_NS;
46     maybelong* row_node;  
47     int NN_row;     dim_t numComp;
48     int NS_row;     index_t * col_DOF;
49       dim_t col_DOF_UpperBound;
50     int numComp;     Finley_ElementFile_Jacobeans* col_jac;
51     maybelong * label_col;     index_t* col_node;
52     Finley_RefElement* referenceElement_col;     dim_t col_NN;
53     maybelong* col_node;     dim_t col_NS;
    int NN_col;  
    int NS_col;  
54    
55     maybelong id[MAX_numNodes]; /* used to hold a reordering vector, referenced by row_node and col_node */     /* added by Ben Cumming for MPI version */
56    
57       index_t id[MAX_numNodes]; /* used to hold a reordering vector, referenced by row_node and col_node */
58  };  };
59    
60  typedef struct Assemble_Parameters Assemble_Parameters;  typedef struct Assemble_Parameters Assemble_Parameters;
61    
62    
63  typedef void (Finley_Assemble_handelShapeMissMatch) (int, int,int, double*,int, int);  #define Finley_Assemble_reducedIntegrationOrder(__in__) ( (getFunctionSpaceType(__in__) == FINLEY_REDUCED_ELEMENTS) || (getFunctionSpaceType(__in__) == FINLEY_REDUCED_FACE_ELEMENTS) || (getFunctionSpaceType(__in__) == FINLEY_REDUCED_CONTACT_ELEMENTS_1) || (getFunctionSpaceType(__in__) == FINLEY_REDUCED_CONTACT_ELEMENTS_2) )
64    
65  void Finley_Assemble_PDE(Finley_NodeFile*,Finley_ElementFile*,Finley_SystemMatrix*,escriptDataC*,  void Finley_Assemble_PDE(Finley_NodeFile*,Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
66                                      escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*) ;                                      escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*) ;
67  void Finley_Assemble_RobinCondition(Finley_NodeFile*,Finley_ElementFile*,Finley_SystemMatrix*,escriptDataC*,  void Finley_Assemble_LumpedSystem(Finley_NodeFile*,Finley_ElementFile*, escriptDataC*, escriptDataC*);
68                                      escriptDataC*,escriptDataC*,Finley_Assemble_handelShapeMissMatch) ;  
69  /* void Finley_Assemble_Points(Finley_Mesh*,Finley_SystemMatrix*,escriptDataC*,escriptDataC*,escriptDataC*) ;*/  void Assemble_getAssembleParameters(Finley_NodeFile*,Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,bool_t, Assemble_Parameters*);
70    void  Finley_Assemble_PDE_System2_3D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
71                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
72    void  Finley_Assemble_PDE_System2_2D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
73                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
74    void  Finley_Assemble_PDE_System2_1D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
75                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
76    void  Finley_Assemble_PDE_System2_C(Assemble_Parameters p, Finley_ElementFile*, Paso_SystemMatrix*, escriptDataC*, escriptDataC*, escriptDataC*);
77    void  Finley_Assemble_PDE_Single2_3D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
78                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
79    void  Finley_Assemble_PDE_Single2_2D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
80                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
81    void  Finley_Assemble_PDE_Single2_1D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
82                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
83    void  Finley_Assemble_PDE_Single2_C(Assemble_Parameters p, Finley_ElementFile*, Paso_SystemMatrix*, escriptDataC*, escriptDataC*, escriptDataC*);
84    
85    
86  void Finley_Assemble_NodeCoordinates(Finley_NodeFile*,escriptDataC*);  void Finley_Assemble_NodeCoordinates(Finley_NodeFile*,escriptDataC*);
87  void Finley_Assemble_setNormal(Finley_NodeFile*, Finley_ElementFile*, escriptDataC*);  void Finley_Assemble_setNormal(Finley_NodeFile*, Finley_ElementFile*, escriptDataC*);
88  void Finley_Assemble_interpolate(Finley_NodeFile*,Finley_ElementFile*,escriptDataC*, escriptDataC*);  void Finley_Assemble_interpolate(Finley_NodeFile*,Finley_ElementFile*,escriptDataC*, escriptDataC*);
# Line 64  void Finley_Assemble_integrate(Finley_No Line 91  void Finley_Assemble_integrate(Finley_No
91  void Finley_Assemble_getSize(Finley_NodeFile*,Finley_ElementFile*, escriptDataC*);  void Finley_Assemble_getSize(Finley_NodeFile*,Finley_ElementFile*, escriptDataC*);
92  void Finley_Assemble_CopyNodalData(Finley_NodeFile* nodes,escriptDataC* out,escriptDataC* in);  void Finley_Assemble_CopyNodalData(Finley_NodeFile* nodes,escriptDataC* out,escriptDataC* in);
93  void Finley_Assemble_CopyElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in);  void Finley_Assemble_CopyElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in);
94  void Finley_Assemble_PDEMatrix_System2(int,int,int,int,int,double*,double*, double*,int, double*, double*,int, double*,int,double*,int,double*,int);  void Finley_Assemble_AverageElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in);
95  void Finley_Assemble_PDEMatrix_Single2(int,int,int,double*,double*, double*,int, double*, double*,int, double*,int,double*,int,double*,int);  void Finley_Assemble_addToSystemMatrix(Paso_SystemMatrix*,dim_t,index_t*, dim_t,dim_t,index_t*,dim_t, double*);
96  void Finley_Assemble_RHSMatrix_System(int,int,int,int,double*,double*,double*,int, double*,double*,int,double*,int);  void Assemble_jacobeans_1D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
97  void Finley_Assemble_RHSMatrix_Single(int,int,int,double*,double*,double*,int, double*,double*,int,double*,int);  void Assemble_jacobeans_2D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
98    void Assemble_jacobeans_2D_M1D_E2D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
99    void Assemble_jacobeans_2D_M1D_E2D_C(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
100    void Assemble_jacobeans_2D_M1D_E1D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
101    void Assemble_jacobeans_2D_M1D_E1D_C(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
102    void Assemble_jacobeans_3D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
103    void Assemble_jacobeans_3D_M2D_E3D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
104    void Assemble_jacobeans_3D_M2D_E3D_C(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
105    void Assemble_jacobeans_3D_M2D_E2D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
106    void Assemble_jacobeans_3D_M2D_E2D_C(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
107    void Finley_Assemble_LumpedSystem(Finley_NodeFile* nodes,Finley_ElementFile* elements, escriptDataC* lumpedMat, escriptDataC* D);
108    
109    
 void Assemble_getAssembleParameters(Finley_NodeFile*,Finley_ElementFile*,Finley_SystemMatrix*,escriptDataC*,Assemble_Parameters*);  
 Finley_Assemble_handelShapeMissMatch Finley_Assemble_handelShapeMissMatch_Step_out;  
 Finley_Assemble_handelShapeMissMatch Finley_Assemble_handelShapeMissMatch_Step_in;  
 Finley_Assemble_handelShapeMissMatch Finley_Assemble_handelShapeMissMatch_Mean_out;  
 Finley_Assemble_handelShapeMissMatch Finley_Assemble_handelShapeMissMatch_Mean_in;  
   
110  #endif /* #ifndef INC_FINLEY_ASSEMBLE */  #endif /* #ifndef INC_FINLEY_ASSEMBLE */
111    
 /*  
  * $Log$  
  * Revision 1.1  2004/10/26 06:53:56  jgs  
  * Initial revision  
  *  
  * Revision 1.1  2004/07/02 04:21:13  gross  
  * Finley C code has been included  
  *  
  *  
  */  

Legend:
Removed from v.82  
changed lines
  Added in v.1312

  ViewVC Help
Powered by ViewVC 1.1.26