Реализация алгоритма умножения двух сжатых матриц по алгоритму Рейнбольдта - Местеньи на языке Си

Нужно реализовать алгортим множения двух сжатых матриц. Матрицы сжаты по кольцевой схеме Рейнбольдта - Местеньи. То есть матрица имеет 5 массивов: 1 - массив с ненулевыми элементами матрицы, 2 - массив с индексами следующих элементов в строке для каждого из элемента, 3 - массив с индексами следующих элементов в столбце для каждого из элемента, 4 - массив с индексами первых элементов в каждой строке (-1 если такого элемента нет), 5 - массив с индексами первых элементов в каждом столбце (-1 если такого элемента нет),