15、数据结构扩展以实现数据实例特定优化及GPU上信念传播优化

数据结构扩展以实现数据实例特定优化及GPU上信念传播优化

1. 直接求解器相关内容

在矩阵运算中，对于直接求解器，其代码有一些特定的特征和转换步骤。代码中数据被分割到两个不同的数组 M 和 M2 中，分别代表矩阵的下三角和上三角部分。 k 和 k2 定义了一种单射访问模式，并被置于内循环中，同时有相应的循环上下界。

以下是直接求解器的部分代码示例：

I = 1; k = 0;
for (g=0; g<glim; g++) {
    Temp = Intermediate[I];
    Temp = Temp / pivot[g];
    Intermediate[I] = Temp;
    for (m=0; m<mlim[i]; m++) {
        Intermediate[rowIdx[k]] -= Temp * M[k];
        k++;
    }
    I++;
}
I2 = Size; k2 = 0;
for (h=0; h<hlim; h++) {
    Temp = Intermediate[I2];
    for (n=0; n<nlim[h]; n++) {
        Temp -= M2[k2]*Intermediate[colIdx[k2]];
        k2++;
    }
    Intermediate[I2] = Temp;
    I2--;
}

在第一个循环中，变量