/[escript]/branches/3.4.1/ripley/src/Rectangle.cpp
ViewVC logotype

Diff of /branches/3.4.1/ripley/src/Rectangle.cpp

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

revision 4587 by jfenwick, Wed Dec 11 06:17:09 2013 UTC revision 4588 by sshaw, Wed Dec 11 23:29:29 2013 UTC
# Line 4047  escript::Data Rectangle::randomFill(long Line 4047  escript::Data Rectangle::randomFill(long
4047          throw RipleyException("Radius of gaussian filter must be less than half the width/height of a rank");          throw RipleyException("Radius of gaussian filter must be less than half the width/height of a rank");
4048      }      }
4049            
4050      size_t inset=2*radius+1;    
     size_t Eheight=ext[1]-2*inset;  // how high the E (shared) region is  
     size_t Swidth=ext[0]-2*inset;  
       
4051      double* src=new double[ext[0]*ext[1]];      double* src=new double[ext[0]*ext[1]];
4052      esysUtils::randomFillArray(seed, src, ext[0]*ext[1]);        esysUtils::randomFillArray(seed, src, ext[0]*ext[1]);  
4053            
4054        
4055  #ifdef ESYS_MPI      #ifdef ESYS_MPI    
4056        size_t inset=2*radius+1;
4057        size_t Eheight=ext[1]-2*inset;  // how high the E (shared) region is
4058        size_t Swidth=ext[0]-2*inset;
4059    
4060      double* SWin=new double[inset*inset];  memset(SWin, 0, inset*inset*sizeof(double));      double* SWin=new double[inset*inset];  memset(SWin, 0, inset*inset*sizeof(double));
4061      double* SEin=new double[inset*inset];  memset(SEin, 0, inset*inset*sizeof(double));      double* SEin=new double[inset*inset];  memset(SEin, 0, inset*inset*sizeof(double));

Legend:
Removed from v.4587  
changed lines
  Added in v.4588

  ViewVC Help
Powered by ViewVC 1.1.26