/[escript]/branches/domexper/dudley/src/Assemble_PDE_System2_3D.c
ViewVC logotype

Diff of /branches/domexper/dudley/src/Assemble_PDE_System2_3D.c

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

revision 853 by gross, Wed Sep 20 05:56:36 2006 UTC revision 1811 by ksteube, Thu Sep 25 23:11:13 2008 UTC
# Line 1  Line 1 
1  /*  
2   ************************************************************  /*******************************************************
3   *          Copyright 2006 by ACcESS MNRF                   *  *
4   *                                                          *  * Copyright (c) 2003-2008 by University of Queensland
5   *              http://www.access.edu.au                    *  * Earth Systems Science Computational Center (ESSCC)
6   *       Primary Business: Queensland, Australia            *  * http://www.uq.edu.au/esscc
7   *  Licensed under the Open Software License version 3.0    *  *
8   *     http://www.opensource.org/licenses/osl-3.0.php       *  * 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    
# Line 33  Line 35 
35    
36  /**************************************************************/  /**************************************************************/
37    
 /*  Author: gross@access.edu.au */  
 /*  Version: $Id:$ */  
   
 /**************************************************************/  
   
38    
39  #include "Assemble.h"  #include "Assemble.h"
40  #include "Util.h"  #include "Util.h"
# Line 81  void  Finley_Assemble_PDE_System2_3D(Ass Line 78  void  Finley_Assemble_PDE_System2_3D(Ass
78                                                                                                                                                                                                                                                                                                                                                                                                            
79         if (!Finley_checkPtr(EM_S) && !Finley_checkPtr(EM_F) && !Finley_checkPtr(row_index) ) {         if (!Finley_checkPtr(EM_S) && !Finley_checkPtr(EM_F) && !Finley_checkPtr(row_index) ) {
80    
           #ifndef PASO_MPI  
81            for (color=elements->minColor;color<=elements->maxColor;color++) {            for (color=elements->minColor;color<=elements->maxColor;color++) {
82               /*  open loop over all elements: */               /*  open loop over all elements: */
83               #pragma omp for private(e) schedule(static)               #pragma omp for private(e) schedule(static)
84               for(e=0;e<elements->numElements;e++){               for(e=0;e<elements->numElements;e++){
85                  if (elements->Color[e]==color) {                  if (elements->Color[e]==color) {
           #else  
           {  
              for(e=0;e<elements->numElements;e++) {  
                 {  
           #endif  
86                     Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);                     Vol=&(p.row_jac->volume[INDEX2(0,e,p.numQuad)]);
87                     DSDX=&(p.row_jac->DSDX[INDEX4(0,0,0,e,p.row_NN,DIM,p.numQuad)]);                     DSDX=&(p.row_jac->DSDX[INDEX4(0,0,0,e,p.row_NN,DIM,p.numQuad)]);
88                     for (q=0;q<len_EM_S;++q) EM_S[q]=0;                     for (q=0;q<len_EM_S;++q) EM_S[q]=0;

Legend:
Removed from v.853  
changed lines
  Added in v.1811

  ViewVC Help
Powered by ViewVC 1.1.26