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

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

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

trunk/finley/src/finley/Assemble.h revision 201 by jgs, Wed Nov 23 04:10:21 2005 UTC trunk/finley/src/Assemble.h revision 2548 by jfenwick, Mon Jul 20 06:20:06 2009 UTC
# Line 1  Line 1 
 /*  
  ******************************************************************************  
  *                                                                            *  
  *       COPYRIGHT  ACcESS 2003,2004,2005 -  All Rights Reserved              *  
  *                                                                            *  
  * This software is the property of ACcESS. No part of this code              *  
  * may be copied in any form or by any means without the expressed written    *  
  * consent of ACcESS.  Copying, use or modification of this software          *  
  * by any unauthorised person is illegal unless that person has a software    *  
  * license agreement with ACcESS.                                             *  
  *                                                                            *  
  ******************************************************************************  
 */  
1    
2  /**************************************************************/  /*******************************************************
3    *
4    * Copyright (c) 2003-2009 by University of Queensland
5    * Earth Systems Science Computational Center (ESSCC)
6    * http://www.uq.edu.au/esscc
7    *
8    * Primary Business: Queensland, Australia
9    * Licensed under the Open Software License version 3.0
10    * http://www.opensource.org/licenses/osl-3.0.php
11    *
12    *******************************************************/
13    
 /*    assemblage routines: header file */  
14    
15  /**************************************************************/  /**************************************************************/
16    
17  /*  Copyrights by ACcESS Australia 2003,2004,2005 */  /*    assemblage routines: header file */
 /*  Author: gross@access.edu.au */  
 /*  Version: $Id$ */  
18    
19  /**************************************************************/  /**************************************************************/
20    
# Line 33  Line 27 
27  #include "Finley.h"  #include "Finley.h"
28  #include "ElementFile.h"  #include "ElementFile.h"
29  #include "NodeFile.h"  #include "NodeFile.h"
30  #include "escript/Data/DataC.h"  #include "escript/DataC.h"
31  #include "paso/SystemMatrix.h"  #include "paso/SystemMatrix.h"
32    
33  struct Assemble_Parameters {  struct Assemble_Parameters {
34     dim_t numQuad;     dim_t numQuad;
35     dim_t numDim;     dim_t numDim;
    dim_t numElementDim;  
   
36     dim_t NN;     dim_t NN;
    dim_t NS;  
    Finley_RefElement* referenceElement;  
37    
38     dim_t numEqu;     dim_t numEqu;
39     index_t* label_row;     index_t* row_DOF;
40     Finley_RefElement* referenceElement_row;     dim_t row_DOF_UpperBound;
41       Finley_ElementFile_Jacobeans* row_jac;
42     index_t* row_node;     index_t* row_node;
43     dim_t NN_row;     dim_t row_NN;
44     dim_t NS_row;     dim_t row_NS;
45    
46     dim_t numComp;     dim_t numComp;
47     index_t * label_col;     index_t * col_DOF;
48     Finley_RefElement* referenceElement_col;     dim_t col_DOF_UpperBound;
49       Finley_ElementFile_Jacobeans* col_jac;
50     index_t* col_node;     index_t* col_node;
51     dim_t NN_col;     dim_t col_NN;
52     dim_t NS_col;     dim_t col_NS;
53    
54       /* added by Ben Cumming for MPI version */
55    
56     index_t id[MAX_numNodes]; /* used to hold a reordering vector, referenced by row_node and col_node */     index_t id[MAX_numNodes]; /* used to hold a reordering vector, referenced by row_node and col_node */
57  };  };
# Line 65  struct Assemble_Parameters { Line 59  struct Assemble_Parameters {
59  typedef struct Assemble_Parameters Assemble_Parameters;  typedef struct Assemble_Parameters Assemble_Parameters;
60    
61    
62  typedef void (Finley_Assemble_handelShapeMissMatch) (dim_t, dim_t,dim_t, double*,dim_t, dim_t);  #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) )
63    
64  void Finley_Assemble_PDE(Finley_NodeFile*,Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,  void Finley_Assemble_PDE(Finley_NodeFile*,Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
65                                      escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*) ;                                      escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*) ;
66  void Finley_Assemble_PDE_RHS(Finley_NodeFile*,Finley_ElementFile*,escriptDataC*,escriptDataC*,escriptDataC*) ;  void Finley_Assemble_LumpedSystem(Finley_NodeFile*,Finley_ElementFile*, escriptDataC*, escriptDataC*);
67  void Finley_Assemble_RobinCondition(Finley_NodeFile*,Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,  
68                                      escriptDataC*,escriptDataC*,Finley_Assemble_handelShapeMissMatch) ;  void Assemble_getAssembleParameters(Finley_NodeFile*,Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,bool_t, Assemble_Parameters*);
69  void Finley_Assemble_RobinCondition_RHS(Finley_NodeFile*,Finley_ElementFile*,escriptDataC*,escriptDataC*,Finley_Assemble_handelShapeMissMatch);  void  Finley_Assemble_PDE_System2_3D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
70  /* void Finley_Assemble_Points(Finley_Mesh*,Paso_SystemMatrix*,escriptDataC*,escriptDataC*,escriptDataC*) ;*/                                       escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
71    void  Finley_Assemble_PDE_System2_2D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
72                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
73    void  Finley_Assemble_PDE_System2_1D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
74                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
75    void  Finley_Assemble_PDE_System2_C(Assemble_Parameters p, Finley_ElementFile*, Paso_SystemMatrix*, escriptDataC*, escriptDataC*, escriptDataC*);
76    void  Finley_Assemble_PDE_Single2_3D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
77                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
78    void  Finley_Assemble_PDE_Single2_2D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
79                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
80    void  Finley_Assemble_PDE_Single2_1D(Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
81                                         escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
82    void  Finley_Assemble_PDE_Single2_C(Assemble_Parameters p, Finley_ElementFile*, Paso_SystemMatrix*, escriptDataC*, escriptDataC*, escriptDataC*);
83    
84    
85  void Finley_Assemble_NodeCoordinates(Finley_NodeFile*,escriptDataC*);  void Finley_Assemble_NodeCoordinates(Finley_NodeFile*,escriptDataC*);
86  void Finley_Assemble_setNormal(Finley_NodeFile*, Finley_ElementFile*, escriptDataC*);  void Finley_Assemble_setNormal(Finley_NodeFile*, Finley_ElementFile*, escriptDataC*);
87  void Finley_Assemble_interpolate(Finley_NodeFile*,Finley_ElementFile*,escriptDataC*, escriptDataC*);  void Finley_Assemble_interpolate(Finley_NodeFile*,Finley_ElementFile*,escriptDataC*, escriptDataC*);
# Line 82  void Finley_Assemble_integrate(Finley_No Line 90  void Finley_Assemble_integrate(Finley_No
90  void Finley_Assemble_getSize(Finley_NodeFile*,Finley_ElementFile*, escriptDataC*);  void Finley_Assemble_getSize(Finley_NodeFile*,Finley_ElementFile*, escriptDataC*);
91  void Finley_Assemble_CopyNodalData(Finley_NodeFile* nodes,escriptDataC* out,escriptDataC* in);  void Finley_Assemble_CopyNodalData(Finley_NodeFile* nodes,escriptDataC* out,escriptDataC* in);
92  void Finley_Assemble_CopyElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in);  void Finley_Assemble_CopyElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in);
93  void Finley_Assemble_PDEMatrix_System2(dim_t,dim_t,dim_t,dim_t,dim_t,double*,double*, double*,dim_t, double*, double*,dim_t, double*,dim_t,double*,dim_t,double*,dim_t);  void Finley_Assemble_AverageElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in);
 void Finley_Assemble_PDEMatrix_Single2(dim_t,dim_t,dim_t,double*,double*, double*,dim_t, double*, double*,dim_t, double*,dim_t,double*,dim_t,double*,dim_t);  
 void Finley_Assemble_RHSMatrix_System(dim_t,dim_t,dim_t,dim_t,double*,double*,double*,dim_t, double*,double*,dim_t,double*,dim_t);  
 void Finley_Assemble_RHSMatrix_Single(dim_t,dim_t,dim_t,double*,double*,double*,dim_t, double*,double*,dim_t,double*,dim_t);  
   
   
 void Assemble_getAssembleParameters(Finley_NodeFile*,Finley_ElementFile*,Paso_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;  
94  void Finley_Assemble_addToSystemMatrix(Paso_SystemMatrix*,dim_t,index_t*, dim_t,dim_t,index_t*,dim_t, double*);  void Finley_Assemble_addToSystemMatrix(Paso_SystemMatrix*,dim_t,index_t*, dim_t,dim_t,index_t*,dim_t, double*);
95    void Assemble_jacobeans_1D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
96    void Assemble_jacobeans_2D(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
97    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*);
98    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*);
99    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*);
100    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*);
101    void Assemble_jacobeans_3D(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_M2D_E3D(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_C(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_E2D(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_C(double*, dim_t,double*, dim_t, dim_t, dim_t, index_t* , double*, dim_t, double*, double*, double*, index_t*);
106    void Finley_Assemble_LumpedSystem(Finley_NodeFile* nodes,Finley_ElementFile* elements, escriptDataC* lumpedMat, escriptDataC* D);
107    
108    
109  #endif /* #ifndef INC_FINLEY_ASSEMBLE */  #endif /* #ifndef INC_FINLEY_ASSEMBLE */
110    
 /*  
  * $Log$  
  * Revision 1.4  2005/09/15 03:44:21  jgs  
  * Merge of development branch dev-02 back to main trunk on 2005-09-15  
  *  
  * Revision 1.3  2005/08/12 01:45:42  jgs  
  * erge of development branch dev-02 back to main trunk on 2005-08-12  
  *  
  * Revision 1.2.2.2  2005/09/07 06:26:17  gross  
  * the solver from finley are put into the standalone package paso now  
  *  
  * Revision 1.2.2.1  2005/08/04 22:41:11  gross  
  * some extra routines for finley that might speed-up RHS assembling in some cases (not actived right now)  
  *  
  * Revision 1.2  2005/07/08 04:07:45  jgs  
  * Merge of development branch back to main trunk on 2005-07-08  
  *  
  * Revision 1.1.1.1.2.1  2005/06/29 02:34:46  gross  
  * some changes towards 64 integers in finley  
  *  
  * Revision 1.1.1.1  2004/10/26 06:53:56  jgs  
  * initial import of project esys2  
  *  
  * Revision 1.1  2004/07/02 04:21:13  gross  
  * Finley C code has been included  
  *  
  *  
  */  

Legend:
Removed from v.201  
changed lines
  Added in v.2548

  ViewVC Help
Powered by ViewVC 1.1.26