问并行化STXXL容器

EN

#include <iostream>
#include <stxxl/vector>

#define FEATURE_LENGTH 128
#define PARTITION_SIZE 16

typedef stxxl::VECTOR_GENERATOR<int>::result vec;

using namespace std;

struct MyMatrix
{
    float A[FEATURE_LENGTH*PARTITION_SIZE];
};

int main()
{
    const int N = (1 << 20);
    stxxl::vector<MyMatrix> C;
    for (int n = 0; n < N; ++n) {
        MyMatrix A;
        for (size_t j = 0; j < PARTITION_SIZE; ++j) {
            for (size_t i = 0; i < FEATURE_LENGTH; ++i) {
                A.A[j*FEATURE_LENGTH + i] = n;
            }
        }
        if (n == 0) {
            std::cout << "sizeof(A) = " << sizeof(A) << std::endl;
        }
        C.push_back(A);
    }

    #pragma omp parallel for shared(C)
    for (unsigned long long m = 0; m < C.size(); ++m) {
        for (int j = 0; j < PARTITION_SIZE; ++j) {
            MyMatrix& C_m = C[m];
            float sum = 0;
            for (int i = 0; i < FEATURE_LENGTH; ++i) {
                sum += C_m.A[j*FEATURE_LENGTH + i];
            }
            if (sum != 0) {
                for (int i = 0; i < FEATURE_LENGTH; ++i) {
                    C_m.A[j*FEATURE_LENGTH + i] /= sum;
                }
            }
        }
    }
    return 0;
}