Zaimplementowałem jedno mnożenie macierzy za pomocąboost::numeric::ublas::matrix (widziećmój pełny, działający kod doładowania)

Result result = read ();

boost::numeric::ublas::matrix<int> C;
C = boost::numeric::ublas::prod(result.A, result.B);

i drugi ze standardowym algorytmem (patrzpełny standardowy kod):

vector< vector<int> > ijkalgorithm(vector< vector<int> > A, 
                                    vector< vector<int> > B) {
    int n = A.size();

    // initialise C with 0s
    vector<int> tmp(n, 0);
    vector< vector<int> > C(n, tmp);

    for (int i = 0; i < n; i++) {
        for (int k = 0; k < n; k++) {
            for (int j = 0; j < n; j++) {
                C[i][j] += A[i][k] * B[k][j];
            }
        }
    }
    return C;
}

W ten sposób testuję prędkość:

time boostImplementation.out > boostResult.txt
diff boostResult.txt correctResult.txt

time simpleImplementation.out > simpleResult.txt
diff simpleResult.txt correctResult.txt

Oba programy odczytują na stałe zapisany plik tekstowy zawierający dwie macierze 2000 x 2000. Oba programy zostały skompilowane z tymi flagami:

g++ -std=c++98 -Wall -O3 -g $(PROBLEM).cpp -o $(PROBLEM).out -pedantic

mam15 sekund dla mojej realizacji i więcej4 minuty na wdrożenie doładowania!

edytuj: po skompilowaniu z

g++ -std=c++98 -Wall -pedantic -O3 -D NDEBUG -DBOOST_UBLAS_NDEBUG library-boost.cpp -o library-boost.out

mam28,19 sekundy dla algorytmu ikj i60,99 sekund dla wzmocnienia. Więc Boost jest nadal znacznie wolniejszy.

Dlaczego wzmocnienie jest znacznie wolniejsze niż moje wdrożenie?