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

revision 1311 by ksteube, Fri Aug 3 00:27:20 2007 UTC revision 1312 by ksteube, Mon Sep 24 06:18:44 2007 UTC
# Line 1  Line 1
/*
************************************************************
*          Copyright 2006,2007 by ACcENULLNULL MNRF              *
*                                                          *
*              http://www.access.edu.au                    *
*       Primary Business: Queensland, Australia            *
*                                                          *
************************************************************
*/
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
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 119  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
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 168  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
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 219  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
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 272  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