Question

我在java中使用Jama进行矩阵运算，但由于它们没有稀疏矩阵，所以我开始使用Parallel Cold Library (PColt)。它是Colt的多线程版本。我试图将两个平方矩阵乘以AxB（它是矩阵 - 矩阵乘法，非元素（或标量）乘法），大小（NxN）。我找不到一个提供的PColt中矩阵 - 矩阵乘法的方法（我不是我想要元素乘法），以便我将方法编码如下。当我乘以N = 1000两个矩阵时，需要将近5分钟才能完成。如果有人知道如何在pcolt中乘以两个矩阵，那将非常感激。或者，如果您发现代码中存在任何不必要的错误并且复杂性很高，请通知我。我的方法如下：

    /**
 * Linear algebraic matrix-matrix multiplication; (new)A = A x B. A[i,j] =
 * Sum(A[i,k] * B[k,j]), k=0..n-1. Matrix shapes: A(m x n), B(n x p), A(m x
 * p).
 * 
 * @param matrix1
 *            first matrix
 * @param matrix2
 *            second matrix
 * @return changes matrix1 with new values matrix-matrix multiplication
 */

public static FloatMatrix2D multiplyMatrix(FloatMatrix2D matrix1,
        FloatMatrix2D matrix2) {
    // create a matrix same size of input matrix
    FloatMatrix2D resultMatrix = matrix1.like();
    // matrix-matrix multiplication row of first matrix must be equal to
    // column of second matrix
    if (matrix1.rows() == matrix2.columns()) {
        for (int i = 0; i < matrix1.rows(); i++) {
            for (int j = 0; j < matrix2.columns(); j++) {
                FloatMatrix1D rowVec = getRow(matrix1, i).copy();
                FloatMatrix1D colVec = getCol(matrix2, j).copy();
                // first multiplies each row with each column and then sum
                // up the result and assign the result to value of matrix
                float multOfVects = rowVec.assign(colVec,
                        FloatFunctions.mult).aggregate(FloatFunctions.plus,
                        FloatFunctions.identity);
                // set sum of vectors to the new matrix
                resultMatrix.setQuick(i, j, multOfVects);

            }
            System.out.println(" i th row "+ i);
        }
    } else {
        System.err
                .println("Row size of first matrix must be equal to Column size of second matrix: "
                        + matrix1 + " != " + matrix2);
    }

    return resultMatrix;
}

////解决方案...... Okie Dokie，我得到了解决方案。实际上忘记了上面的代码。 PColt提供矩阵 - 矩阵乘法，但方法名称令人困惑。

无论如何，为了多个两个矩阵，请使用以下方法：

public DoubleMatrix2D zMult（DoubleMatrix2D B， DoubleMatrix2D C）线性代数矩阵 - 矩阵乘法; C = A x B;相当于A.zMult（B，C，1,0，false，false）

在此关注参数的顺序，因为结果被保存到最后的参数，可能 C ....这真的花了我很多时间：|

Answer 1

矩阵乘法方法在DoubleAlgebra类中。

DenseDoubleAlgebra algebra = new DenseDoubleAlgebra();
DoubleMatrix2D productMatrix = algebra.mult(aMatrix, anotherMatrix);

Java矩阵运算，并行colt矩阵 - 矩阵乘法

1 个答案: