/[escript]/trunk/ripley/src/Rectangle.cpp
ViewVC logotype

Diff of /trunk/ripley/src/Rectangle.cpp

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

revision 4696 by jfenwick, Wed Feb 19 07:29:50 2014 UTC revision 4705 by jfenwick, Fri Feb 21 02:36:15 2014 UTC
# Line 1883  namespace Line 1883  namespace
1883              for (int x=-r;x<=r;++x)              for (int x=-r;x<=r;++x)
1884              {                      {        
1885                  arr[(x+r)+(y+r)*(r*2+1)]=common*exp(-(x*x+y*y)/(2*sigma*sigma));                  arr[(x+r)+(y+r)*(r*2+1)]=common*exp(-(x*x+y*y)/(2*sigma*sigma));
 // cout << (x+y*(r*2+1)) << " " << arr[(x+r)+(y+r)*(r*2+1)] << endl;  
1886                  total+=arr[(x+r)+(y+r)*(r*2+1)];                  total+=arr[(x+r)+(y+r)*(r*2+1)];
1887              }              }
1888          }          }
1889          double invtotal=1/total;          double invtotal=1/total;
 //cout << "Inv total is "        << invtotal << endl;  
1890          for (size_t p=0;p<(radius*2+1)*(radius*2+1);++p)          for (size_t p=0;p<(radius*2+1)*(radius*2+1);++p)
1891          {          {
1892              arr[p]*=invtotal;              arr[p]*=invtotal;
 //cout << p << "->" << arr[p] << endl;        
1893          }          }
1894          return arr;          return arr;
1895      }      }
# Line 2067  escript::Data Rectangle::randomFillWorke Line 2064  escript::Data Rectangle::randomFillWorke
2064  #endif  #endif
2065                    
2066      esysUtils::patternFillArray2D(ext[0], ext[1], src, 4, basex, basey, numvals);      esysUtils::patternFillArray2D(ext[0], ext[1], src, 4, basex, basey, numvals);
2067  */      */
2068    
2069            
     
2070  #ifdef ESYS_MPI    #ifdef ESYS_MPI  
2071            
2072      BlockGrid2 grid(m_NX[0]-1, m_NX[1]-1);      BlockGrid2 grid(m_NX[0]-1, m_NX[1]-1);
# Line 2092  escript::Data Rectangle::randomFillWorke Line 2089  escript::Data Rectangle::randomFillWorke
2089      grid.generateOutNeighbours(X, Y, outcoms);      grid.generateOutNeighbours(X, Y, outcoms);
2090            
2091      block.copyAllToBuffer(src);        block.copyAllToBuffer(src);  
   
 // for (int i=0;i<9;++i)  
 // {  
 //     if (i!=4)  
 //     {  
 //       for (int j=0;j<block.getBuffSize(i);++j)  
 //       {  
 //    block.getOutBuffer(i)[j]=100+i;  
 //       }  
 //     }  
 // }  
       
2092            
2093      int comserr=0;          int comserr=0;    
2094      for (size_t i=0;i<incoms.size();++i)      for (size_t i=0;i<incoms.size();++i)
# Line 2117  escript::Data Rectangle::randomFillWorke Line 2102  escript::Data Rectangle::randomFillWorke
2102      for (size_t i=0;i<outcoms.size();++i)      for (size_t i=0;i<outcoms.size();++i)
2103      {      {
2104          message& m=outcoms[i];          message& m=outcoms[i];
 // cout << "Sending " <<    (int)m.srcbuffid << " with tag " << m.tag << endl;  
2105          comserr|=MPI_Isend(block.getOutBuffer(m.srcbuffid), block.getBuffSize(m.srcbuffid) , MPI_DOUBLE, m.destID, m.tag, m_mpiInfo->comm, reqs+(rused++));          comserr|=MPI_Isend(block.getOutBuffer(m.srcbuffid), block.getBuffSize(m.srcbuffid) , MPI_DOUBLE, m.destID, m.tag, m_mpiInfo->comm, reqs+(rused++));
 // for (int i=0;i<block.getBuffSize(m.srcbuffid);++i)  
 // {  
 //     cout << block.getOutBuffer(m.srcbuffid)[i] << " ";  
 // }  
 // cout << endl;      
2106      }          }    
2107            
2108      if (!comserr)      if (!comserr)
2109      {          {    
2110          comserr=MPI_Waitall(rused, reqs, stats);          comserr=MPI_Waitall(rused, reqs, stats);    
       
       
 //     for (size_t i=0;i<incoms.size();++i)  
 //     {  
 //         message& m=incoms[i];        
 // cout << "Gettinging " <<     (int)m.destbuffid << " with tag " << m.tag << endl;  
 // for (int i=0;i<block.getBuffSize(m.destbuffid);++i)  
 // {  
 //     cout << block.getInBuffer(m.destbuffid)[i] << " ";  
 // }  
 // cout << endl;  
 //    
 //     }      
       
2111      }      }
2112    
2113      if (comserr)      if (comserr)

Legend:
Removed from v.4696  
changed lines
  Added in v.4705

  ViewVC Help
Powered by ViewVC 1.1.26