Оптимизация произведения матриц через указатели. Как избавится от счетчиков?

@DioxiT · Регистрация: 11.12.2018

Author24 — интернет-сервис помощи студентам

Разрабатываю оператор для класса матриц, позволяющий умножать матрицы друг на друга. Есть рабочий вариант, однако преподаватель попросил убрать в нем счетчики в циклах и использовать для этого сами указатели. Однако как это реализовать идей уже нет.
Информация об операторе:

Кликните здесь для просмотра всего текста

В классе имеется информация о количестве строк и столбцов матрицы:

C++

1	unsigned int rows, cols;

И указатель на первый элемент матрицы:

C++

1	double* elems;

Метод возвращающий указатель на первый элемент указанной в аргументе строки:

C++

1	double* RowPtr(unsigned row = 1);

Сам оператор умножения:

C++

Matrix operator*(const Matrix &B) const //произведение двух матриц (ЛОГИНОВ)
    {
        if (elems == NULL && B.elems == NULL) //предпроверки на пустую матрицу, сравнение строк и столбцов
        {
            EC = 1;
            return Matrix();
        }
        if (cols != B.rows)
        {
            EC = 7;
            return Matrix();
        }
        Matrix C = Matrix(rows, B.cols); //создаем матрицу для записи результатов
        double* colsA = elems; double* rowsA = elems; double* ptrC = C.elems;
        double* colsB = B.elems; double* rowsB = B.elems;
    
        for (unsigned j = 1; j <= rows; j++)
        {
            rowsA = RowPtr(j); colsB = B.elems;
                for (unsigned i = 1; i <= B.cols; i++, ptrC++, colsB++)
                {
                    colsA = rowsA; rowsB = colsB;
                    for (unsigned m = 1; m <= cols; m++, colsA++, rowsB = B.RowPtr(m) + (i - 1))
                        *ptrC += *colsA*(*rowsB);
                }
        }
        return C;
    }

Собственно, как избавиться от j,i,m не потеряв работоспособность не понимаю. Т.к. счетчики увеличиваются на единицу так же как и указатели, пробовал условия цикла ставить через указатели, например внутренний цикл с m:

C++

for (unsigned i = 1; i <= B.cols; i++, ptrC++, colsB++)
                {
                    colsA = rowsA; rowsB = colsB;
                    for (; colsA <=elems + cols;colsA++, rowsB = B.RowPtr(m) + (i - 1))
                        *ptrC += *colsA*(*rowsB);
                }

Первую строку умножает, естественно отлично, но на второй сразу выходит из цикла, т.к. номер первого элемента второй строки уже сразу за границами elems+cols (т.е. последнего элемента первой строки). Получается нужно условие делать каким то динамическим для каждой строки, только как?!!
P.S. Надеюсь раздел выбрал нужный, заранее извиняюсь за возможную ошибку.

	23.12.2019, 21:00