Files
2016-04-27 01:07:35 -04:00

53 lines
1.8 KiB
Plaintext

//******************************************************************************
//
// File: Test06.cu
// Package: edu.rit.gpu.test
// Unit: Test06 kernel function
//
// This C/CUDA source file is copyright (C) 2014 by Alan Kaminsky. All rights
// reserved. For further information, contact the author, Alan Kaminsky, at
// ark@cs.rit.edu.
//
// This C/CUDA source file is part of the Parallel Java 2 Library ("PJ2"). PJ2
// is free software; you can redistribute it and/or modify it under the terms of
// the GNU General Public License as published by the Free Software Foundation;
// either version 3 of the License, or (at your option) any later version.
//
// PJ2 is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
// A PARTICULAR PURPOSE. See the GNU General Public License for more details.
//
// A copy of the GNU General Public License is provided in the file gpl.txt. You
// may also obtain a copy of the GNU General Public License on the World Wide
// Web at http://www.gnu.org/licenses/gpl.html.
//
//******************************************************************************
/**
* Compute the sum of two double matrices.
* <P>
* Called with a 2-D grid of 2-D blocks. Each thread computes one element of the
* output matrix.
*
* @param a First input matrix.
* @param b Second input matrix.
* @param c Output matrix.
* @param rows Number of matrix rows.
* @param cols Number of matrix columns.
*
* @author Alan Kaminsky
* @version 05-Apr-2014
*/
extern "C" __global__ void addDoubleMatrices
(double **a,
double **b,
double **c,
int rows,
int cols)
{
int row = blockIdx.y*blockDim.y + threadIdx.y;
int col = blockIdx.x*blockDim.x + threadIdx.x;
if (row < rows && col < cols)
c[row][col] = a[row][col] + b[row][col];
}