pniedzwiedzinski / cpp-cuda-linear-regression Goto Github PK

__global__ void transpose(matrix* A, matrix* Return) {
  int y = blockIdx.y * blockDim.y + threadIdx.y;
  int x = blockIdx.x * blockDim.x + threadIdx.x;

  Return[y][x] = A[x][y];
}

__global__ void dot(matrix* A, matrix* B, matrix* Return) {
  int row = blockIdx;
  int col = threadIdx.x;

  for i in range(A.x.length())
    sum += A[row][i] * B[i][col];

  Return[row][col] = sum;
}

[[a, b],
 [c, d]]

[[a, b, c],
 [d, e, f],
 [g, h, i]]

float determinant(matrix* A) {
  if matrix.dim = 2 {
    return A[0][0] * A[1][1] - A[1][0] * A[0][1];
  } else {
    // split???
}