There is a C # code, and I need to redo it so that everything is counted on a video card:

public static double[] SATL(double[] data) { double[] matrix = { 0.0982862174, ..., -0.0229204861, ..., 0.0161380976 }; int dataLength = data.Length; int matrixLength = matrix.Length; int start = matrixLength + 1; var newData = new double[dataLength]; if (dataLength <= matrixLength) { return null; } for (int i = matrixLength; i < dataLength; i++) { int counter = i - matrixLength; for (int j = 0; j < matrixLength; j++) { newData[i] += matrix[j] * data[counter++]; } } return newData; }

The problem is that I still can not figure out how to run through the second array. I don’t understand how indexing happens in Cudafy . This is what I did, but it does not work correctly:

 public static double[] FATLCuda(double[] data, double[] matrix) { CudafyModule km = CudafyTranslator.Cudafy(); GPGPU gpu = CudafyHost.GetDevice(CudafyModes.Target, CudafyModes.DeviceId); gpu.LoadModule(km); ... gpu.Launch().add(dev_data, dev_matrix, dev_newdata); } [Cudafy] public static void add(GThread thread, double[] data, double[] matrix, double[] newdata) { int tid = thread.blockIdx.x; if (tid < data.Length) { int jid = 0; int counter = tid - matrix.Length; if (jid < matrix.Length) { newdata[tid] += matrix[jid] * data[counter++]; jid++; } tid++; } }
  • one
    To get an answer faster, try to formalize the problem. Of course, the reward stimulates the responders, but: not everyone will understand the C # code; another person who has encountered a problem cannot find it by searching. - Cyrus
  • Nothing is clear: | What is the "second array"? matrix? What does it mean: "can not run"? Go through the loop? - Mstislav Pavlov

1 answer 1

Not so long ago in .net added support for SIMD . This namespace allows you to use hardware acceleration. To use, you must have RyuJIT compiler, .NET 4.6 and System.Numerics.Vectors, which is put through Nuget.

An example of the simplest program using this namespace:

 using System; using System.Numerics; class Program { static void Main(string[] args) { const Int32 N = 8; Single[] a = { 41982.0F, 81.5091F, 3.14F, 42.666F, 54776.45F, 342.4556F, 6756.2344F, 4563.789F }; Single[] b = { 85989.111F, 156.5091F, 3.14F, 42.666F, 1006.45F, 9999.4546F, 0.2344F, 7893.789F }; Single[] c = new Single[N]; for (int i = 0; i < N; i += Vector<Single>.Count) // Count возвращает 16 для char, 4 для float, 2 для double и т.п. { var aSimd = new Vector<Single>(a, i); // создать экземпляр со смещением i var bSimd = new Vector<Single>(b, i); Vector<Single> cSimd = aSimd + bSimd; // или так Vector<Single> c_simd = Vector.Add(b_simd, a_simd); cSimd.CopyTo(c, i); //копировать в массив со смещением } for (int i = 0; i < a.Length; i++) { Console.WriteLine(c[i]); } Console.ReadKey(); } }

In other words, using this namespace, you can make everything much simpler and more beautiful.