/[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 trunk/dudley/src/Assemble_CopyElementData.c revision 3981 by jfenwick, Fri Sep 21 02:47:54 2012 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-2012 by University of Queensland
5    * http://www.uq.edu.au
6    *
7    * Primary Business: Queensland, Australia
8    * Licensed under the Open Software License version 3.0
9    * http://www.opensource.org/licenses/osl-3.0.php
10    *
11    * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12    * Development since 2012 by School of Earth Sciences
13    *
14    *****************************************************************************/
15    
16  /**************************************************************/  /************************************************************************************/
17    
18  /*   author: gross@access.edu.au */  /*    assemblage routines: copies data between elements       */
 /*   Version: $Id$ */  
19    
20  /**************************************************************/  /************************************************************************************/
21    
22  #include "Assemble.h"  #include "Assemble.h"
23  #include "Util.h"  #include "Util.h"
24  #ifdef _OPENMP  #ifdef _OPENMP
25  #include <omp.h>  #include <omp.h>
26  #endif  #endif
27  /******************************************************************************************************/  /****************************************************************************************************************************/
28    #include "ShapeTable.h"
29    
30    void Dudley_Assemble_CopyElementData(Dudley_ElementFile * elements, escriptDataC * out, escriptDataC * in)
31    {
32        dim_t n, q, numElements, numQuad;
33        __const double *in_array;
34        double *out_array;
35        dim_t numComps = getDataPointSize(out);
36        size_t len_size;
37    
38        Dudley_resetError();
39        if (elements == NULL)
40        {
41        return;
42        }
43    
44  void Finley_Assemble_CopyElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in) {      numElements = elements->numElements;
45      dim_t n,q, numElements, numQuad;      if (Dudley_Assemble_reducedIntegrationOrder(in))
46      dim_t numComps=getDataPointSize(out);      {
47      double *in_array,*out_array;      numQuad = QuadNums[elements->numDim][0];
48      Finley_resetError();      }
49        else
50      if (elements==NULL) return;      {
51      numElements=elements->numElements;      numQuad = QuadNums[elements->numDim][1];
52      numQuad=elements->ReferenceElement->numQuadNodes;  
53        }
54    
55      /* check out and in */      /* check out and in */
56      if (numComps!=getDataPointSize(in)) {      if (numComps != getDataPointSize(in))
57         Finley_setError(TYPE_ERROR,"__FILE__: number of components of input and output Data do not match.");      {
58      } else if (!numSamplesEqual(in,numQuad,numElements)) {      Dudley_setError(TYPE_ERROR,
59         Finley_setError(TYPE_ERROR,"__FILE__: illegal number of samples of input Data object");              "Dudley_Assemble_CopyElementData: number of components of input and output Data do not match.");
60      } else if (!numSamplesEqual(out,numQuad,numElements)) {      }
61         Finley_setError(TYPE_ERROR,"__FILE__: illegal number of samples of output Data object");      else if (!numSamplesEqual(in, numQuad, numElements))
62      } else if (!isExpanded(out)) {      {
63         Finley_setError(TYPE_ERROR,"__FILE__: expanded Data object is expected for output data.");      Dudley_setError(TYPE_ERROR, "Dudley_Assemble_CopyElementData: illegal number of samples of input Data object");
64        }
65        else if (!numSamplesEqual(out, numQuad, numElements))
66        {
67        Dudley_setError(TYPE_ERROR, "Dudley_Assemble_CopyElementData: illegal number of samples of output Data object");
68        }
69        else if (!isExpanded(out))
70        {
71        Dudley_setError(TYPE_ERROR,
72                "Dudley_Assemble_CopyElementData: expanded Data object is expected for output data.");
73      }      }
74    
75      /* now we can start */      /* now we can start */
76    
77      if (Finley_noError()) {      if (Dudley_noError())
78           if (isExpanded(in)) {      {
79               # pragma omp parallel for private(n) schedule(static)      if (isExpanded(in))
80               for (n=0;n<numElements;n++)      {
81                   Finley_copyDouble(numComps*numQuad,getSampleData(in,n),getSampleData(out,n));          len_size = numComps * numQuad * sizeof(double);
82           } else {          requireWrite(out);
83               # pragma omp parallel for private(q,n,out_array,in_array) schedule(static)  #pragma omp parallel private(n)
84               for (n=0;n<numElements;n++) {          {
85                   in_array=getSampleData(in,n);  # pragma omp for schedule(static)
86                   out_array=getSampleData(out,n);          for (n = 0; n < numElements; n++)
87                   for (q=0;q<numQuad;q++) Finley_copyDouble(numComps,in_array,out_array+q*numComps);              memcpy(getSampleDataRW(out, n), getSampleDataRO(in, n), len_size);
88               }          }
89           }      }
90        else
91        {
92            len_size = numComps * sizeof(double);
93            requireWrite(out);
94    #pragma omp parallel private(q,n,out_array,in_array)
95            {
96    # pragma omp for schedule(static)
97            for (n = 0; n < numElements; n++)
98            {
99                in_array = getSampleDataRO(in, n);
100                out_array = getSampleDataRW(out, n);
101                for (q = 0; q < numQuad; q++)
102                memcpy(out_array + q * numComps, in_array, len_size);
103            }
104            }
105        }
106      }      }
107      return;      return;
108  }  }
 /*  
  * $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.3981

  ViewVC Help
Powered by ViewVC 1.1.26