cudaktmethod.cpp

Go to the documentation of this file.

#include "cudaktmethod.h"
#include "numericmethods.h"
extern "C" {
void integrate_hyperbolic_code(double m_CFL, double t, double tEnd, double *m_vSol, double *m_vX0, double *m_vX1, double *m_vY0, double *m_vY1, double *m_VZ0, double *m_VZ1, double *m_Vx, double *m_Vy, double *m_Vz);
}


CUDAKTMethod::CUDAKTMethod(OrthoMesh &mesh,Function3D &fInitU,const VecDouble &cPor,Function3D &fPrescribedU,double v1,double v2,double theta,double rs1,double rs2,double CFL)
  :m_mesh(mesh)
{

  //Set the initial values for the saturation;
  m_vSol.reinit(mesh.numCells(),1);
  mesh.setCentralValuesFromFunction(fInitU,m_vSol);
  m_vFluxes.reinit(mesh.numFaces());
  m_fValues.reinit(mesh.numFaces());
  VecIndex vIX0,vIX1,vIY0,vIY1,vIZ0,vIZ1;
  m_CFL = CFL;
  
  mesh.getFacesInSlab(mesh.getP()[0],vIX0,0);
  mesh.getFacesInSlab(mesh.getQ()[0],vIX1,0);
  mesh.getFacesInSlab(mesh.getP()[1],vIY0,1);
  mesh.getFacesInSlab(mesh.getQ()[1],vIY1,1);
  mesh.getFacesInSlab(mesh.getP()[2],vIZ0,2);
  mesh.getFacesInSlab(mesh.getQ()[2],vIZ1,2);
  
  NumericMethods::fillVector(m_fValues,INFINITY);
  mesh.setFacesValuesFromFunction(fPrescribedU,m_fValues,0,true);


  NumericMethods::getSubVector(m_fValues,vIX0,m_vX0);
  NumericMethods::getSubVector(m_fValues,vIX1,m_vX1);
  NumericMethods::getSubVector(m_fValues,vIY0,m_vY0);
  NumericMethods::getSubVector(m_fValues,vIY1,m_vY1);
  NumericMethods::getSubVector(m_fValues,vIZ0,m_vZ0);
  NumericMethods::getSubVector(m_fValues,vIZ1,m_vZ1);

  //
  //mesh.setFacesValuesFromFunction(fPrescribedU,m_fValues,0,true);
  m_nX = mesh.numElemX();
  m_nY = mesh.numElemY();
  m_nZ = mesh.numElemZ();

  m_Q = mesh.getQ();

  m_rs1=rs1;
  m_rs2=rs2;
  m_theta=theta;
  m_v1 = v1;
  m_v2=v2;
  m_por = cPor(0);
}


void CUDAKTMethod::updateVelocities(DynamicBase &dynMod)
{
  //Get normal velocities
  dynMod.getNormalVelocityAtFaces(m_vNV);

  //split the velocity vector field at faces in three parts
  //one corresponding to each dimension
  m_vX = &(m_vNV(0));
  m_vY = m_vX+m_mesh.getFacesInYOffset();
  m_vZ = m_vX+m_mesh.getFacesInZOffset();
    
  //Store the velocitieson GPU
}

ArrayOfVecDouble&  CUDAKTMethod::getSolutionAtCells()
{
  return m_vSol;

}


 void CUDAKTMethod::printOutput() 
{

}


 void CUDAKTMethod::updateDataForDynamicModule() 
{

}


 void CUDAKTMethod::iterateN(unsigned n,double dt) 
{
  throw new Exception("Cuda not working for now\n");


  //m_vSol has the initial value parameter
  //m_vX0, m_vX1 LEFT AND RIGHT FACE
  //m_vZ0, m_vZ1 BOTTOM AND UP  FACE
  //m_vY0, m_vY1 FRONT AND BACK FACE
  //m_vX, m_vY, m_vZ normal components of the velocities at the faces

  //To avoid linking problems this line of code was commented.
  //integrate_hyperbolic_code(m_CFL, t, tEnd, &(m_vSol(0)), &(m_vX0(0)), &(m_vX1(0)), &(m_vY0(0)), &(m_vY1(0)), &(m_vZ0(0)), &(m_vZ1(0)), m_vX, m_vY, m_vZ);
}




 double CUDAKTMethod::getDt(double t,double tEnd) 
{
  return NAN;
  throw new Exception("Cuda not working for now\n");

}

---