static string TestcublasSgemm2()
{
string testName = "TestcublasSgemm2";
uint aw = 2;
uint bh = aw;
uint ah = 3;
uint bw = 3;
uint ch = ah;
uint cw = bw;
GpuMatrix gpuA;
GpuMatrix gpuB;
GpuMatrix gpuC;
var dataA = MatrixUtils.AA();
var dataB = MatrixUtils.BB();
var cRes = new float[ch * cw];
var cuby = new CublasClr.Cublas();
var aa = new CudaArray();
var res = aa.ResetDevice();
res = res + GpuMatrixOps.SetupGpuMatrix(
out gpuA,
new Matrix <float>(_rows: ch, _cols: cw,
host_data: ImmutableArray.Create(dataA),
matrixFormat: MatrixFormat.Column_Major));
res = res + GpuMatrixOps.SetupGpuMatrix(
out gpuB,
new Matrix <float>(_rows: bh, _cols: bw,
host_data: ImmutableArray.Create(dataB),
matrixFormat: MatrixFormat.Column_Major));
res = res + GpuMatrixOps.SetupGpuMatrix(
out gpuC,
new Matrix <float>(_rows: ch, _cols: cw,
host_data: ImmutableArray.Create(cRes),
matrixFormat: MatrixFormat.Column_Major));
IntPtr cublasHandle = new IntPtr();
res = res + cuby.MakeCublasHandle(ref cublasHandle);
GpuMatrix gpuProd;
res = res + GpuMatrixOps.Multiply(
gmOut: out gpuProd,
cublasHandle: cublasHandle,
gmA: gpuA, gmB: gpuB, gmC: gpuC);
GpuMatrix gpuSynched;
res = res + GpuMatrixOps.CopyToHost(out gpuSynched, gpuProd);
// GpuMatrixUtils.MatrixMult(C: cRes, A: dataA, B: dataB, wA: aw, hA: ah, wB: bw);
return(string.Empty);
}
static string TestGpuMatrix()
{
var data = MatrixUtils.MakeIdentiPoke(rows: 5, cols: 3);
GpuMatrix gpuSetup;
var res = GpuMatrixOps.SetupGpuMatrix(out gpuSetup,
new Matrix <float>(_rows: 5, _cols: 3,
host_data: ImmutableArray.Create(data),
matrixFormat: MatrixFormat.Column_Major));
GpuMatrix gpuReturned;
res = res + GpuMatrixOps.CopyToHost(out gpuReturned, gpuSetup);
GpuMatrix gpuCleared;
res = res + GpuMatrixOps.ClearOnDevice(out gpuCleared, gpuReturned);
var dout = gpuReturned.Matrix.Data;
return(String.Empty);
}
static string TestcublasSgemm1()
{
string testName = "TestcublasSgemm";
uint aw = 5;
uint bh = aw;
uint ah = 5;
uint bw = 5;
uint ch = ah;
uint cw = bw;
var cuby = new CublasClr.Cublas();
var aa = new CudaArray();
var res = aa.ResetDevice();
var dataA = MatrixUtils.MakeIdentity(rows: ah, cols: aw);
GpuMatrix gpuA;
res = res + GpuMatrixOps.SetupGpuMatrix(
out gpuA,
new Matrix <float>(_rows: ah, _cols: aw,
host_data: ImmutableArray.Create(dataA),
matrixFormat: MatrixFormat.Column_Major));
var dataB = MatrixUtils.MakeIdentiPoke(rows: bh, cols: bw);
GpuMatrix gpuB;
res = res + GpuMatrixOps.SetupGpuMatrix(
out gpuB,
new Matrix <float>(_rows: bh, _cols: bw,
host_data: ImmutableArray.Create(dataB),
matrixFormat: MatrixFormat.Column_Major));
var dataC = MatrixUtils.MakeZeroes(rows: bh, cols: bw);
GpuMatrix gpuC;
res = res + GpuMatrixOps.SetupGpuMatrix(
out gpuC,
new Matrix <float>(_rows: ch, _cols: cw,
host_data: ImmutableArray.Create(dataC),
matrixFormat: MatrixFormat.Column_Major));
IntPtr cublasHandle = new IntPtr();
res = res + cuby.MakeCublasHandle(ref cublasHandle);
GpuMatrix gpuProd;
res = res + GpuMatrixOps.Multiply(
gmOut: out gpuProd,
cublasHandle: cublasHandle,
gmA: gpuA, gmB: gpuB, gmC: gpuC);
GpuMatrix gpuSynched;
res = res + GpuMatrixOps.CopyToHost(out gpuSynched, gpuProd);
var cpuRes = new float[ah * bw];
MatrixUtils.RowMajorMatrixMult(C: cpuRes, A: dataA, B: dataB, wA: aw, hA: ah, wB: bw);
var cpuRes2 = new float[bh * aw];
MatrixUtils.RowMajorMatrixMult(C: cpuRes2, A: dataB, B: dataA, wA: bw, hA: bh, wB: aw);
return(res);
}
