EOS.C

Go to the documentation of this file.

#include "EOS.H"

namespace eos {

int eos::ComputePressureBuffer(const pcpp::IndexIntervalType &regionInterval,
                               const std::vector<size_t>     &bufferSizes){
  
  numDim = bufferSizes.size();

  yStart = 0;
  yEnd = 0;
  zStart = 0;
  zEnd = 0;

  xStart = regionInterval[0].first;
  xEnd   = regionInterval[0].second;
  xSize  = bufferSizes[0];
  numPointsBuffer = xSize;

  if(numDim > 1){ 
    yStart = regionInterval[1].first;
    yEnd   = regionInterval[1].second;
    ySize  = bufferSizes[1];
    numPointsBuffer *= ySize;
  }   

  if(numDim > 2){ 
    zStart = regionInterval[2].first;
    zEnd   = regionInterval[2].second;
    zSize  = bufferSizes[2];
    numPointsBuffer *= zSize;
  }  

  size_t nPlane = xSize*ySize;
  size_t nPoints = xSize*ySize*zSize;

  for(size_t iZ = zStart;iZ <= zEnd;iZ++){
    size_t zIndex = iZ*nPlane;
    for(size_t iY = yStart;iY <= yEnd;iY++){
      size_t yzIndex = iY*xSize+zIndex;
      for(size_t iX = xStart;iX <= xEnd;iX++){
        size_t xyzIndex = yzIndex + iX;

        //ComputePressurePtr(xyzIndex);
        ComputePressure(xyzIndex);
      }
    }
  }
  return(0);
}

int eos::ComputeTemperatureBuffer(const pcpp::IndexIntervalType &regionInterval,
                                  const std::vector<size_t>     &bufferSizes) {

  numDim = bufferSizes.size();

  yStart = 0;
  yEnd = 0;
  zStart = 0;
  zEnd = 0;

  xStart = regionInterval[0].first;
  xEnd   = regionInterval[0].second;
  xSize  = bufferSizes[0];
  numPointsBuffer = xSize;

  if(numDim > 1){ 
    yStart = regionInterval[1].first;
    yEnd   = regionInterval[1].second;
    ySize  = bufferSizes[1];
    numPointsBuffer *= ySize;
  }   

  if(numDim > 2){ 
    zStart = regionInterval[2].first;
    zEnd   = regionInterval[2].second;
    zSize  = bufferSizes[2];
    numPointsBuffer *= zSize;
  }  
  
  size_t nPlane = xSize*ySize;
  size_t nPoints = xSize*ySize*zSize;

  for(size_t iZ = zStart;iZ <= zEnd;iZ++){
    size_t zIndex = iZ*nPlane;
    for(size_t iY = yStart;iY <= yEnd;iY++){
      size_t yzIndex = iY*xSize+zIndex;
      for(size_t iX = xStart;iX <= xEnd;iX++){
        size_t xyzIndex = yzIndex + iX;
    
        //ComputeTemperaturePtr(xyzIndex);
        ComputeTemperature(xyzIndex);
      }
    }
  }
  return(0);
}

} // namespace eos