Method transposeMultiply

JSAT/src/jsat/linear/SparseMatrix.java:308–334

(Matrix B, Matrix C)

    }

    @Override
    public void transposeMultiply(Matrix B, Matrix C)
    {
        if(this.rows() != B.rows())//Normaly it is A_cols == B_rows, but we are doint A'*B, not A*B
            throw new ArithmeticException("Matrix dimensions do not agree");
        else if(this.cols() != C.rows() || B.cols() != C.cols())
            throw new ArithmeticException("Destination matrix does not have matching dimensions");
        final SparseMatrix A = this;
        ///Should choose step size such that 2*NB2^2 * dataTypeSize <= CacheSize
        
        final int kLimit = this.rows();

        for (int k = 0; k < kLimit; k++)
        {
            Vec bRow_k = B.getRowView(k);
            Vec aRow_k = A.getRowView(k);

            for (IndexValue iv : aRow_k)//iterating over "i"
            {

                Vec cRow_i = C.getRowView(iv.getIndex());
                double a = iv.getValue();//A.get(k, i);

                cRow_i.mutableAdd(a, bRow_k);
            }
        }
    }

    @Override
    public void transposeMultiply(final Matrix B, final Matrix C, ExecutorService threadPool)