/[escript]/trunk/dudley/src/Assemble_CopyElementData.cpp
ViewVC logotype

Diff of /trunk/dudley/src/Assemble_CopyElementData.cpp

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

trunk/finley/src/Assemble_CopyElementData.c revision 1028 by gross, Wed Mar 14 00:15:24 2007 UTC branches/domexper/dudley/src/Assemble_CopyElementData.c revision 3199 by jfenwick, Thu Sep 23 06:08:56 2010 UTC
# Line 1  Line 1 
 /*  
  ************************************************************  
  *          Copyright 2006 by ACcESS MNRF                   *  
  *                                                          *  
  *              http://www.access.edu.au                    *  
  *       Primary Business: Queensland, Australia            *  
  *  Licensed under the Open Software License version 3.0    *  
  *     http://www.opensource.org/licenses/osl-3.0.php       *  
  *                                                          *  
  ************************************************************  
 */  
 /**************************************************************/  
1    
2  /*    assemblage routines: copies data between elements       */  /*******************************************************
3    *
4    * Copyright (c) 2003-2010 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    
14    
15  /**************************************************************/  /**************************************************************/
16    
17  /*   author: gross@access.edu.au */  /*    assemblage routines: copies data between elements       */
 /*   Version: $Id$ */  
18    
19  /**************************************************************/  /**************************************************************/
20    
# Line 26  Line 24 
24  #include <omp.h>  #include <omp.h>
25  #endif  #endif
26  /******************************************************************************************************/  /******************************************************************************************************/
27    #include "ShapeTable.h"
28    
29    void Dudley_Assemble_CopyElementData(Dudley_ElementFile* elements,escriptDataC* out,escriptDataC* in) {
 void Finley_Assemble_CopyElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in) {  
30      dim_t n,q, numElements, numQuad;      dim_t n,q, numElements, numQuad;
31        __const double *in_array;
32        double *out_array;
33      dim_t numComps=getDataPointSize(out);      dim_t numComps=getDataPointSize(out);
34      double *in_array,*out_array;      size_t len_size;
35      Finley_resetError();  
36        Dudley_resetError();
37        if( elements == NULL )
38        {
39           return;
40        }
41    
     if (elements==NULL) return;  
42      numElements=elements->numElements;      numElements=elements->numElements;
43      numQuad=elements->ReferenceElement->numQuadNodes;      if (Dudley_Assemble_reducedIntegrationOrder(in)) {
44    //       numQuad=elements->referenceElementSet->referenceElementReducedQuadrature->BasisFunctions->numQuadNodes;
45        numQuad=QuadNums[elements->numDim][0];
46        } else {
47    //       numQuad=elements->referenceElementSet->referenceElement->BasisFunctions->numQuadNodes;
48        numQuad=QuadNums[elements->numDim][1];
49    
50        }
51    
52      /* check out and in */      /* check out and in */
53      if (numComps!=getDataPointSize(in)) {      if (numComps!=getDataPointSize(in)) {
54         Finley_setError(TYPE_ERROR,"__FILE__: number of components of input and output Data do not match.");         Dudley_setError(TYPE_ERROR,"Dudley_Assemble_CopyElementData: number of components of input and output Data do not match.");
55      } else if (!numSamplesEqual(in,numQuad,numElements)) {      } else if (!numSamplesEqual(in,numQuad,numElements)) {
56         Finley_setError(TYPE_ERROR,"__FILE__: illegal number of samples of input Data object");         Dudley_setError(TYPE_ERROR,"Dudley_Assemble_CopyElementData: illegal number of samples of input Data object");
57      } else if (!numSamplesEqual(out,numQuad,numElements)) {      } else if (!numSamplesEqual(out,numQuad,numElements)) {
58         Finley_setError(TYPE_ERROR,"__FILE__: illegal number of samples of output Data object");         Dudley_setError(TYPE_ERROR,"Dudley_Assemble_CopyElementData: illegal number of samples of output Data object");
59      } else if (!isExpanded(out)) {      } else if (!isExpanded(out)) {
60         Finley_setError(TYPE_ERROR,"__FILE__: expanded Data object is expected for output data.");         Dudley_setError(TYPE_ERROR,"Dudley_Assemble_CopyElementData: expanded Data object is expected for output data.");
61      }      }
62    
63      /* now we can start */      /* now we can start */
64    
65      if (Finley_noError()) {      if (Dudley_noError()) {
66           if (isExpanded(in)) {           if (isExpanded(in)) {
67               # pragma omp parallel for private(n) schedule(static)               len_size=numComps*numQuad*sizeof(double);
68               for (n=0;n<numElements;n++)           requireWrite(out);
69                   Finley_copyDouble(numComps*numQuad,getSampleData(in,n),getSampleData(out,n));           #pragma omp parallel private(n)
70             {
71                   # pragma omp for schedule(static)
72                   for (n=0;n<numElements;n++)
73                     memcpy(getSampleDataRW(out,n),getSampleDataRO(in,n), len_size);
74             }
75           } else {           } else {
76               # pragma omp parallel for private(q,n,out_array,in_array) schedule(static)               len_size=numComps*sizeof(double);
77               for (n=0;n<numElements;n++) {           requireWrite(out);
78                   in_array=getSampleData(in,n);           #pragma omp parallel private(q,n,out_array,in_array)
79                   out_array=getSampleData(out,n);           {
80                   for (q=0;q<numQuad;q++) Finley_copyDouble(numComps,in_array,out_array+q*numComps);                 # pragma omp for schedule(static)
81               }                 for (n=0;n<numElements;n++) {
82                     in_array=getSampleDataRO(in,n);
83                     out_array=getSampleDataRW(out,n);
84                     for (q=0;q<numQuad;q++) memcpy(out_array+q*numComps,in_array,len_size);
85                   }
86             }
87           }           }
88      }      }
89      return;      return;
90  }  }
 /*  
  * $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/02 05:29:11  gross  
  * bug in finley/src/Assemble_CopyElementData fixed  
  *  
  * 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.2  2004/07/21 05:00:54  gross  
  * name changes in DataC  
  *  
  * Revision 1.1  2004/07/02 04:21:13  gross  
  * Finley C code has been included  
  *  
  *  
  */  

Legend:
Removed from v.1028  
changed lines
  Added in v.3199

  ViewVC Help
Powered by ViewVC 1.1.26