cuIBM/NavierStokes_2generateQT_8inl_source.html

 #include <solvers/kernels/generateQT.h>


 template <>
 void NavierStokesSolver<device_memory>::generateQT()
 {
     logger.startTimer("generateQT");

     int nx = domInfo->nx,
         ny = domInfo->ny;

     int numU  = (nx-1)*ny;
     int numUV = numU + nx*(ny-1);
     int numP  = numU + ny;

     cooH QTHost(numP, numUV, 4*numP-2*(nx+ny));

     int *QTRows = thrust::raw_pointer_cast(&(QTHost.row_indices[0])),
         *QTCols = thrust::raw_pointer_cast(&(QTHost.column_indices[0]));
     real *QTVals = thrust::raw_pointer_cast(&(QTHost.values[0]));

     kernels::generateQT(QTRows, QTCols, QTVals, nx, ny);
     QT = QTHost;
     cusp::transpose(QT, Q);

     logger.stopTimer("generateQT");
 } // generateQT
real
double real
Is a float or a double depending on the machine precision.
Definition: types.h:116

kernels::generateQT
void generateQT(int *QTRows, int *QTCols, real *QTVals, int nx, int ny)
Generates the divergence matrix (on the host).
Definition: generateQT.cu:110

cooH
cusp::coo_matrix< int, real, host_memory > cooH
COO matrix stored on the host.
Definition: types.h:126

generateQT.h
Declaration of the kernels to generate gradient matrix and interpolation matrix.

NavierStokesSolver::generateQT
virtual void generateQT()