/[escript]/branches/arrayview_from_1695_trunk/finley/src/Assemble_AverageElementData.c
ViewVC logotype

Diff of /branches/arrayview_from_1695_trunk/finley/src/Assemble_AverageElementData.c

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

revision 1116 by gross, Mon Apr 23 03:36:34 2007 UTC revision 1312 by ksteube, Mon Sep 24 06:18:44 2007 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       */  /* $Id$ */
3    
4    /*******************************************************
5     *
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    
18  /*   author: gross@access.edu.au */  /*    assemblage routines: copies data between elements       */
 /*   Version: $Id$ */  
19    
20  /**************************************************************/  /**************************************************************/
21    
# Line 33  void Finley_Assemble_AverageElementData( Line 32  void Finley_Assemble_AverageElementData(
32      double *in_array,*out_array, vol, volinv, *wq;      double *in_array,*out_array, vol, volinv, *wq;
33      register double rtmp;      register double rtmp;
34      dim_t numComps=getDataPointSize(out);      dim_t numComps=getDataPointSize(out);
35        size_t numComps_size;
36    
37      Finley_resetError();      Finley_resetError();
38      if( elements == NULL )      if( elements == NULL )
# Line 84  void Finley_Assemble_AverageElementData( Line 84  void Finley_Assemble_AverageElementData(
84                   }                   }
85               }               }
86           } else {           } else {
87                 numComps_size=numComps*sizeof(double);
88               # pragma omp parallel for private(q,n,out_array,in_array) schedule(static)               # pragma omp parallel for private(q,n,out_array,in_array) schedule(static)
89               for (n=0;n<numElements;n++) {               for (n=0;n<numElements;n++) {
90                   in_array=getSampleData(in,n);                   in_array=getSampleData(in,n);
91                   out_array=getSampleData(out,n);                   out_array=getSampleData(out,n);
92                   for (q=0;q<numQuad_out;q++) Finley_copyDouble(numComps,in_array,out_array+q*numComps);                   for (q=0;q<numQuad_out;q++) memcpy(out_array+q*numComps,in_array,numComps_size);
93               }               }
94           }           }
95      }      }

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

  ViewVC Help
Powered by ViewVC 1.1.26