JCuda:对数组进行并行算术运算?

时间:2011-01-27 22:30:56

标签: java cuda parallel-processing

我找不到最简单任务之一的JCuda示例:将两个值数组添加到第三个数组。

我总是只找到C / C ++示例,这些示例依赖于某些CUDA预编译器,因此对我来说无法使用。可怕的事情发生在那里,例如:

VecAdd<<<blocksPerGrid, threadsPerBlock>>>(d_A, d_B, d_C, N)

好吧,在JCuda图书馆中,与VecAdd没什么相似。

JCuda中有很多子库,例如JCudppJCublas,......我真的不知道哪里来查找和什么

请帮忙!

3 个答案:

答案 0 :(得分:2)

我设法组建了一个工作示例:

来源

import java.util.*;

import jcuda.*;
import static jcuda.jcublas.JCublas.*;

public class Main {

  private static final int VECTOR_SIZE = 10;

  static {
    Runtime.getRuntime().addShutdownHook(new Thread(new Runnable() {

      @Override
      public void run() {
        cublasShutdown();
      }

    }));

    cublasInit();
  }

  static void printArrays(float[] a1, float[] a2) {
    System.out.println("\t" + Arrays.toString(a1));
    System.out.println("\t" + Arrays.toString(a2));
    System.out.println();
  }

  public static void main(String[] args) {
    float[] a1 = new float[VECTOR_SIZE];
    float[] a2 = new float[VECTOR_SIZE];

    for (int i = 0; i < VECTOR_SIZE; ++i) {
      a1[i] = i + 1;
    }

    Pointer p1 = new Pointer();
    Pointer p2 = new Pointer();

    cublasAlloc(VECTOR_SIZE, Sizeof.FLOAT, p1);
    cublasAlloc(VECTOR_SIZE, Sizeof.FLOAT, p2);

    System.out.println("BEFORE:");
    printArrays(a1, a2);

    cublasSetVector(VECTOR_SIZE, Sizeof.FLOAT, Pointer.to(a1), 1, p1, 1);
    cublasSetVector(VECTOR_SIZE, Sizeof.FLOAT, Pointer.to(a2), 1, p2, 1);

    // THE ACTUAL OPERATION: MULTIPLY AND ADD
    cublasSaxpy(VECTOR_SIZE, 20f, p1, 1, p2, 1);

    cublasGetVector(VECTOR_SIZE, Sizeof.FLOAT, p2, 1, Pointer.to(a2), 1);

    System.out.println("AFTER:");
    printArrays(a1, a2);

    cublasFree(p1);
    cublasFree(p2);
  }

}

输出

BEFORE:
    [1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0]
    [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]

AFTER:
    [1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0]
    [20.0, 40.0, 60.0, 80.0, 100.0, 120.0, 140.0, 160.0, 180.0, 200.0]

注意

现在我将使用AMD Aparapi

答案 1 :(得分:1)

我会玩 JCublas 。多亏了你,你将能够在Java应用程序中使用CUBLAS,即Basic Linear Algebra Subprograms的NVIDIA CUDA实现。

以下页面提供的示例:http://www.jcuda.de/jcuda/jcudpp/JCudpp.html非常完整。

答案 2 :(得分:0)

您可以使用用C / C ++编写的VecAdd函数。似乎JCuda只能加载cubin文件,这是有道理的,因为cuda代码甚至不是由C编译器创建的。因此,您可以使用JCublas中的内置函数或在C中创建函数并创建cubin文件。此功能可以在JCuda中启动。