/[escript]/trunk/finley/src/Assemble_LumpedSystem.c
ViewVC logotype

Diff of /trunk/finley/src/Assemble_LumpedSystem.c

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

trunk/finley/src/Assemble_LumpedSystem.c revision 1220 by btully, Thu Aug 2 04:46:20 2007 UTC temp/finley/src/Assemble_LumpedSystem.c revision 1387 by trankine, Fri Jan 11 07:45:26 2008 UTC
# Line 1  Line 1 
 /*  
  ************************************************************  
  *          Copyright 2006,2007 by ACcENULLNULL MNRF              *  
  *                                                          *  
  *              http://www.access.edu.au                    *  
  *       Primary Business: Queensland, Australia            *  
  *  Licensed under the Open NULLoftware License version 3.0    *  
  *     http://www.opensource.org/licenses/osl-3.0.php       *  
  *                                                          *  
  ************************************************************  
 */  
1    
2    /* $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    
# Line 21  Line 23 
23    
24  /**************************************************************/  /**************************************************************/
25    
 /*  Author: gross@access.edu.au */  
 /*  Version: $Id$ */  
   
 /**************************************************************/  
   
26  #include "Assemble.h"  #include "Assemble.h"
27  #include "Util.h"  #include "Util.h"
28  #ifdef _OPENMP  #ifdef _OPENMP
# Line 33  Line 30 
30  #endif  #endif
31    
32    
33  #define NEW_LUMPING /* */  /* Disabled until the tests pass */
34    /* #define NEW_LUMPING */ /* */
35    
36  /**************************************************************/  /**************************************************************/
37    
# Line 118  void Finley_Assemble_LumpedSystem(Finley Line 116  void Finley_Assemble_LumpedSystem(Finley
116         if ( !Finley_checkPtr(EM_lumpedMat) && !Finley_checkPtr(row_index) ) {         if ( !Finley_checkPtr(EM_lumpedMat) && !Finley_checkPtr(row_index) ) {
117            if (p.numEqu == 1) {            if (p.numEqu == 1) {
118               if (expandedD) {               if (expandedD) {
                  #ifndef PASO_MPI  
119                   for (color=elements->minColor;color<=elements->maxColor;color++) {                   for (color=elements->minColor;color<=elements->maxColor;color++) {
120                      /*  open loop over all elements: */                      /*  open loop over all elements: */
121                      #pragma omp for private(e) schedule(static)                      #pragma omp for private(e) schedule(static)
122                      for(e=0;e<elements->numElements;e++){                      for(e=0;e<elements->numElements;e++){
123                         if (elements->Color[e]==color) {                         if (elements->Color[e]==color) {
                  #else  
                  {  
                     for(e=0;e<elements->numElements;e++) {  
                        {  
                  #endif  
124                            Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);                            Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);
125                            memset(EM_lumpedMat,0,len_EM_lumpedMat_size);                            memset(EM_lumpedMat,0,len_EM_lumpedMat_size);
126                            D_p=getSampleData(D,e);                            D_p=getSampleData(D,e);
# Line 167  void Finley_Assemble_LumpedSystem(Finley Line 159  void Finley_Assemble_LumpedSystem(Finley
159                      } /* end element loop */                      } /* end element loop */
160                    } /* end color loop */                    } /* end color loop */
161               } else  {               } else  {
                  #ifndef PASO_MPI  
162                   for (color=elements->minColor;color<=elements->maxColor;color++) {                   for (color=elements->minColor;color<=elements->maxColor;color++) {
163                   /*  open loop over all elements: */                   /*  open loop over all elements: */
164                   #pragma omp for private(e) schedule(static)                   #pragma omp for private(e) schedule(static)
165                   for(e=0;e<elements->numElements;e++){                   for(e=0;e<elements->numElements;e++){
166                      if (elements->Color[e]==color) {                      if (elements->Color[e]==color) {
                  #else  
                  {  
                     for(e=0;e<elements->numElements;e++) {  
                        {  
                  #endif  
167                            Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);                            Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);
168                            memset(EM_lumpedMat,0,len_EM_lumpedMat_size);                            memset(EM_lumpedMat,0,len_EM_lumpedMat_size);
169                            D_p=getSampleData(D,e);                            D_p=getSampleData(D,e);
# Line 218  void Finley_Assemble_LumpedSystem(Finley Line 204  void Finley_Assemble_LumpedSystem(Finley
204               }               }
205            } else {            } else {
206               if (expandedD) {               if (expandedD) {
                  #ifndef PASO_MPI  
207                   for (color=elements->minColor;color<=elements->maxColor;color++) {                   for (color=elements->minColor;color<=elements->maxColor;color++) {
208                      /*  open loop over all elements: */                      /*  open loop over all elements: */
209                      #pragma omp for private(e) schedule(static)                      #pragma omp for private(e) schedule(static)
210                      for(e=0;e<elements->numElements;e++){                      for(e=0;e<elements->numElements;e++){
211                         if (elements->Color[e]==color) {                         if (elements->Color[e]==color) {
                  #else  
                  {  
                     for(e=0;e<elements->numElements;e++) {  
                        {  
                  #endif  
212                            Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);                            Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);
213                            memset(EM_lumpedMat,0,len_EM_lumpedMat_size);                            memset(EM_lumpedMat,0,len_EM_lumpedMat_size);
214                            D_p=getSampleData(D,e);                            D_p=getSampleData(D,e);
# Line 271  void Finley_Assemble_LumpedSystem(Finley Line 251  void Finley_Assemble_LumpedSystem(Finley
251                      } /* end element loop */                      } /* end element loop */
252                  } /* end color loop */                  } /* end color loop */
253               } else {               } else {
254                   #ifndef PASO_MPI                   /*  open loop over all elements: */
255                   for (color=elements->minColor;color<=elements->maxColor;color++) {                   for (color=elements->minColor;color<=elements->maxColor;color++) {
                     /*  open loop over all elements: */  
256                      #pragma omp for private(e) schedule(static)                      #pragma omp for private(e) schedule(static)
257                      for(e=0;e<elements->numElements;e++){                      for(e=0;e<elements->numElements;e++){
258                         if (elements->Color[e]==color) {                         if (elements->Color[e]==color) {
                  #else  
                  {  
                     for(e=0;e<elements->numElements;e++) {  
                        {  
                  #endif  
259                            Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);                            Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);
260                            memset(EM_lumpedMat,0,len_EM_lumpedMat_size);                            memset(EM_lumpedMat,0,len_EM_lumpedMat_size);
261                            D_p=getSampleData(D,e);                            D_p=getSampleData(D,e);

Legend:
Removed from v.1220  
changed lines
  Added in v.1387

  ViewVC Help
Powered by ViewVC 1.1.26