public static cDot ( |
||
aStore | The first vector a. | |
aOffset | int | Offset in the vector a. |
aStride | int | The offset between two elements of the vector a. |
bStore | The second vector b. | |
bOffset | int | Offset in the vector a. |
bStride | int | The offset between two elements of the vector a. |
count | int | The number of elements to consider. |
return |
/// <summary> /// y = A * x + y, where x, y are vectors and A is a matrix. /// </summary> /// <param name="aStore">1 dim double array for A.</param> /// <param name="aOffset">Offset in the array for A.</param> /// <param name="aRowStride">Rows in A.</param> /// <param name="aColStride">Columns in A.</param> /// <param name="xStore">1 dim double array for x.</param> /// <param name="xOffset">Offset in the array for x.</param> /// <param name="xStride">Number of entries in x.</param> /// <param name="yStore">1 dim double array for y.</param> /// <param name="yOffset">Offset in the array for y.</param> /// <param name="yStride">Number of entries in y.</param> /// <param name="rows">Geometry information for the rows.</param> /// <param name="cols">Geometry information for the columns.</param> public static void cGemv( ScalarValue[] aStore, int aOffset, int aRowStride, int aColStride, ScalarValue[] xStore, int xOffset, int xStride, ScalarValue[] yStore, int yOffset, int yStride, int rows, int cols ) { int aIndex = aOffset; int yIndex = yOffset; for (int n = 0; n < rows; n++) { yStore[yIndex] += BlasL1.cDot(aStore, aIndex, aColStride, xStore, xOffset, xStride, cols); aIndex += aRowStride; yIndex += yStride; } }
/// <summary> /// Performs a matrix times matrix operation C = A * B with complex matrices. /// </summary> /// <param name="aStore">The 1-dimensional complex array for the matrix A.</param> /// <param name="aOffset">The offset in the array for A.</param> /// <param name="aRowStride">The difference for skipping one row.</param> /// <param name="aColStride">The difference for skipping one column.</param> /// <param name="bStore">The 1-dimensional complex array for the matrix B.</param> /// <param name="bOffset">The offset in the array for B.</param> /// <param name="bRowStride">The difference for skipping one row.</param> /// <param name="bColStride">The difference for skipping one column.</param> /// <param name="cStore">The 1-dimensional complex array for the matrix C.</param> /// <param name="cOffset">The offset in the array for C.</param> /// <param name="cRowStride">The difference for skipping one row.</param> /// <param name="cColStride">The difference for skipping one column.</param> /// <param name="rowsA">The rows to handle in the matrix A.</param> /// <param name="colsB">The coluumns to handle in the matrix B.</param> /// <param name="length">The columns of A / rows of B - or length for the multiplication.</param> public static void cGemm(ScalarValue[] aStore, int aOffset, int aRowStride, int aColStride, ScalarValue[] bStore, int bOffset, int bRowStride, int bColStride, ScalarValue[] cStore, int cOffset, int cRowStride, int cColStride, int rowsA, int colsB, int length) { int aIndex = aOffset; int cStart = cOffset; for (int n = 0; n < rowsA; n++) { int bIndex = bOffset; int cIndex = cStart; for (int m = 0; m < colsB; m++) { cStore[cIndex] = BlasL1.cDot(aStore, aIndex, aColStride, bStore, bIndex, bRowStride, length); bIndex += bColStride; cIndex += cColStride; } aIndex += aRowStride; cStart += cRowStride; } }