parallel/pico/SubDomain.cpp

// -*- C++ -*-
//
//  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
//                                   Fehmi Cirak
//                        California Institute of Technology
//                           (C) 2004 All Rights Reserved
//
//  <LicenseText>
//
//  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
// EDITED BY MATT GALLOWAY FOR USE WITH DAPTA

#include "DAPTA.hpp"

#include <limits>
#include <cassert>
#include <cstdlib>
#include <algorithm>

#include <mpi.h>

namespace DAPTA { // Define namespace DAPTA

namespace pico { // Define namespace pico

template <typename V>
SubDomain<V>::SubDomain(const V& p, const V& q, double maxElemSize) {
   for (int i=0; i<V::dimension; ++i) {
      _low.coords[i] = std::min(p.coords[i], q.coords[i]);
      _upp.coords[i] = std::max(p.coords[i], q.coords[i]);
   }
   
   _maxElemSize = maxElemSize;
   
   // initialize Buffers
   _sendMsgSize = 0;
   _sBuf=NULL;
   _rBuf=NULL;
}

template <typename V>
SubDomain<V>::~SubDomain(){
   if (_sBuf) delete _sBuf;
   if (_rBuf) delete _rBuf;
}

template <typename V>
bool SubDomain<V>::checkOverlap (const SubDomain * const b) const {
   VertexType lp = b->lowerPoint();
   VertexType up = b->upperPoint();
   
   if ( lp.coords[0] > _upp.coords[0] ||
       lp.coords[1] > _upp.coords[1] ||
       lp.coords[2] > _upp.coords[2] ||
       up.coords[0] < _low.coords[0] ||
       up.coords[1] < _low.coords[1] ||
       up.coords[2] < _low.coords[2] ) return false;
   
   return true;    
}

template <typename V>
bool SubDomain<V>::isIn (const V& point) const {
   for (int i=0; i<V::dimension; ++i) 
      if (point.coords[i] < _low.coords[i] || _upp.coords[i] < point.coords[i]) 
         return false;
   
   return true;    
}

template <typename V>
void SubDomain<V>::registerSendBuffer(const char * const start, const size_t&  size) {
   if (!_sBuf) _sBuf = new CBuffer();
   _sendMsgSize = size;
   _sBuf->copyToBuffer(start, size);
   
   return;
}

template <typename V>
void SubDomain<V>::allocateReceiveBuffer(const size_t& size) {    
   if (!_rBuf) _rBuf = new CBuffer();
   assert(_rBuf != NULL);
   assert(size>0);
   _rBuf->allocateBuffer(size);
   
   return;
}
   
} // namespace pico

} // namespace DAPTA

---