0

我正在尝试制作子矩阵以将两个矩阵相乘。我已经获得了伪代码,因此所有试图找到在 java 中实现它的正确方法的尝试,我仍然停留在基础知识上。这个伪代码的 java 代码的整体结构是什么样的?

SQUARE-MATRIX-MULTITPLY-RECURSIVE(A,B)

n = A.rows
let C be a n x n matrix 
if n == 1  
   c11 = a11 * b11
else partition A, B and C
     C11 = SQURAE-MATRIX-MULTIPLY.RECURSIVE(A11,B11)
          + SQUARE-MARTIX-MULTIPLY-RECURSIVE(A12,B21)
     C12 = SQUARE-MATRIX-MULTIPLY-RECURSUÌVE(A11,B12)
          + SQUARE-MATRIX-MULTIPLY-RECURSIVE(A11,B12)
     C21 = SQUARE-MATRIX-MULTIPLY-RECURSUÌVE(A12,B22)
          + SQUARE-MATRIX-MULTIPLY-RECURSUÌVE(A22,B21)
     C22 = SQUARE-MATRIX-MULTIPLY-RECURSUÌVE()A21,B12)
          + SQUARE-MATRIX-MULTIPLY-RECURSUÌVE(A22,B22)
     return C;

由于矩阵 A 和 B 是 nxn 矩阵,它们可以用大小为 n^2 的规则数组来表示。使用 a[in + j] 在数组 a 中访问的矩阵(0 索引)时代的系数 Cij。

我的任务是实现两个版本:一个将子矩阵复制到新数组中,另一个使用两个索引变量来设置正在使用的矩阵的哪一部分

4

1 回答 1

0

这是SparceMatrix乘法的实现

稀疏表示的向量类:

public class SparseVector {
    private HashST<Integer, Double> st;

    public SparseVector() {
        st = new HashST<Integer, Double>();
    }

    public int size() {
        return st.size();
    }

    public void put(int i, double x) {
        st.put(i, x);
    }

    public double get(int i) {
        if (!st.contains(i))
            return 0.0;
        else
            return st.get(i);
    }

    public double dot(double[] that) {
        double sum = 0.0;
        for (int i : st.keys())
            sum += that[i] * this.get(i);
        return sum;
    }
}

主要方法/调用代码:

SparseVector[] a;
a = new SparseVector[N];
double[] x = new double[N];
double[] b = new double[N]; // Initialize a[] and x[]. ...
for (int i = 0; i < N; i++)
    b[i] = a[i].dot(x);
于 2014-03-11T12:00:03.190 回答