public static ProcResult GetEnergy() { var strRet = String.Empty; int[] res = new int[_area / 1024]; strRet = strRet + _cudaArray.RunBlockAddInts_32_Kernel( destPtr: d_energyBlocks, srcPtr: d_energy, span: _span ); strRet = strRet + _cudaArray.CopyIntsFromDevice(res, d_energyBlocks, _area / 1024); float tot = res.Sum(); tot /= _area; var dRet = new Dictionary <string, object>(); dRet["Grid"] = res; dRet["Energy"] = tot; return(new ProcResult(data: dRet, err: strRet, stepsCompleted: 0, time: _stopwatch.ElapsedMilliseconds)); }
public static ProcResult ProcMarkBlocks(int steps) { var strRet = String.Empty; int[] res = new int[_area]; _stopwatch.Reset(); _stopwatch.Start(); for (var i = 0; i < steps; i++) { //uint[] rrands = IntArrayGen.RandUInts((int)DateTime.Now.Ticks, _blockCount); //strRet = strRet + _cudaArray.CopyUIntsToDevice(rrands, d_rands, _blockCount); strRet = strRet + _randoProcs.MakeRandomInts(d_indexRands, _blockCount); strRet = strRet + _gridProcs.Run_k_RandBlockPick( destPtr: d_grid, randPtr: d_indexRands, block_size: _block_size, blocks_per_span: _blocks_per_span ); } strRet = strRet + _cudaArray.CopyIntsFromDevice(res, d_grid, _area); _stopwatch.Stop(); //uint[] res2 = new uint[_blockCount]; //strRet = strRet + _cudaArray.CopyUIntsFromDevice(res2, d_rands, _blockCount); var dRet = new Dictionary <string, object>(); dRet["Grid"] = new SimGrid <int>(name: "Update", width: _span, height: _span, data: res); return(new ProcResult(data: dRet, err: strRet, stepsCompleted: steps, time: _stopwatch.ElapsedMilliseconds)); }
public static string Update(int[] results, int steps) { var strRet = String.Empty; for (int i = 0; i < steps; i++) { if (_backwards) { strRet = _gridProcs.Runk_Gol(d_In, d_Out, _span); if (!String.IsNullOrEmpty(strRet)) { return(strRet); } } else { strRet = _gridProcs.Runk_Gol(d_Out, d_In, _span); if (!String.IsNullOrEmpty(strRet)) { return(strRet); } } _backwards = !_backwards; } if (_backwards) { strRet = _cudaArray.CopyIntsFromDevice(results, d_In, _area); } else { strRet = _cudaArray.CopyIntsFromDevice(results, d_Out, _area); } return(strRet); }
static string TestCopyIntsDeviceToDevice() { string testName = "TestCopyIntsDeviceToDevice"; uint arrayLen = 1000; var alist = Enumerable.Range(4, (int)arrayLen).ToArray(); var aa = new CudaArray(); IntPtr devDataA = new System.IntPtr(); IntPtr devDataB = new System.IntPtr(); var retlist = new int[(int)arrayLen]; try { var res = aa.ResetDevice(); res = res + aa.MallocIntsOnDevice(ref devDataA, arrayLen); res = res + aa.MallocIntsOnDevice(ref devDataB, arrayLen); res = res + aa.CopyIntsToDevice(alist, devDataA, arrayLen); res = res + aa.CopyIntsDeviceToDevice(devDataB, devDataA, arrayLen); res = res + aa.CopyIntsFromDevice(retlist, devDataB, arrayLen); res = res + aa.ReleaseDevicePtr(devDataA); res = res + aa.ReleaseDevicePtr(devDataB); if (!alist.SequenceEqual(retlist)) { return(testName + " fail: sequences do not match"); } if (res != String.Empty) { return(testName + " fail: " + res); } return(testName + " pass"); } catch (Exception ex) { return(testName + " exception " + ex.Message); } finally { aa.ReleaseDevicePtr(devDataA); aa.ReleaseDevicePtr(devDataB); aa.ResetDevice(); } }
public static string Init(float[] temp_inputs, int[] flip_inputs, uint span, uint blockSize, int seed) { _span = span; _block_size = blockSize; _area = _span * _span; _blocks_per_span = span / blockSize; _blockCount = _blocks_per_span * _blocks_per_span; d_flipData = new IntPtr(); d_tempData = new IntPtr(); d_flipRands = new IntPtr(); d_indexRands = new IntPtr(); d_threshes = new IntPtr(); _cudaArray = new CudaArray(); _gridProcs = new GridProcs(); _randoProcs = new RandoProcs(); var strRet = _cudaArray.ResetDevice(); strRet = strRet + _randoProcs.MakeGenerator32(seed); strRet = strRet + _cudaArray.MallocFloatsOnDevice(ref d_tempData, _area); strRet = strRet + _cudaArray.MallocFloatsOnDevice(ref d_heatBlocks, _area / 1024); strRet = strRet + _cudaArray.MallocIntsOnDevice(ref d_flipData, _area); strRet = strRet + _cudaArray.MallocIntsOnDevice(ref d_indexRands, _blockCount); strRet = strRet + _cudaArray.MallocIntsOnDevice(ref d_flipRands, _blockCount); strRet = strRet + _cudaArray.MallocFloatsOnDevice(ref d_threshes, _allTempSteps); strRet = strRet + _cudaArray.CopyIntsToDevice(flip_inputs, d_flipData, _area); strRet = strRet + _cudaArray.CopyFloatsToDevice(temp_inputs, d_tempData, _area); var res9 = new int[_area]; strRet = strRet + _cudaArray.CopyIntsFromDevice(res9, d_flipData, _area); return(strRet); }
static string TestMakeRandomInts() { string testName = "TestMakeRandomInts"; uint arrayLen = 1000; int seed = 1234; var aa = new CudaArray(); var rdo = new RandoClr.RandoProcs(); IntPtr devRando = new IntPtr(); IntPtr devData = new IntPtr(); var retlist = new int[(int)arrayLen]; try { var res = aa.ResetDevice(); res = res + rdo.MakeGenerator32(ref devRando, seed); res = res + aa.MallocIntsOnDevice(ref devData, arrayLen); res = res + rdo.MakeRandomInts(devData, devRando, arrayLen); res = res + aa.CopyIntsFromDevice(retlist, devData, arrayLen); res = res + aa.ReleaseDevicePtr(devData); res = res + rdo.DestroyGenerator(devRando); if (res != String.Empty) { return(testName + " fail: " + res); } return(testName + " pass"); } catch { return(testName + " fail"); } finally { //rdo.DestroyGenerator(devRando); aa.ReleaseDevicePtr(devData); aa.ResetDevice(); } }
public static ProcResult UpdateH(int steps, float qRate, float filpEnergy, float beta) { var strRet = String.Empty; _stopwatch.Reset(); _stopwatch.Start(); for (var s = 0; s < steps; s++) { var res9 = new int[_area]; strRet = strRet + _cudaArray.CopyIntsFromDevice(res9, d_flipData, _area); var bbs = FloatFuncs.Betas(_tempSteps, beta); strRet = strRet + _cudaArray.CopyFloatsToDevice(bbs, d_threshes, _allTempSteps); strRet = strRet + _randoProcs.MakeRandomInts(d_indexRands, _blockCount); strRet = strRet + _randoProcs.MakeUniformRands(d_flipRands, _blockCount); strRet = strRet + _gridProcs.Run_k_ThermoIsing_bp( temp_data: d_tempData, flip_data: d_flipData, index_rands: d_indexRands, flip_rands: d_flipRands, threshes: d_threshes, flip_energy: filpEnergy, block_size: _block_size, blocks_per_span: _blocks_per_span, q_rate: qRate); } _stopwatch.Stop(); var dRet = new Dictionary <string, object>(); float[] bhts = new float[_area / 1024]; strRet = strRet + _cudaArray.RunBlockAddFloats_32_Kernel( destPtr: d_heatBlocks, srcPtr: d_tempData, span: _span ); strRet = strRet + _cudaArray.CopyFloatsFromDevice(bhts, d_heatBlocks, _area / 1024); float tot = bhts.Sum(); tot /= _area; dRet["TotalHeat"] = tot; var tres = new float[_area]; strRet = strRet + _cudaArray.CopyFloatsFromDevice(tres, d_tempData, _area); dRet["ThermGrid"] = new SimGrid <float>(name: "Therms", width: _span, height: _span, data: tres); var fres = new int[_area]; strRet = strRet + _cudaArray.CopyIntsFromDevice(fres, d_flipData, _area); dRet["FlipGrid"] = new SimGrid <int>(name: "Flips", width: _span, height: _span, data: fres); return(new ProcResult(data: dRet, err: strRet, stepsCompleted: steps, time: _stopwatch.ElapsedMilliseconds)); }