c++ умножение разреженных матриц

[world12_tk]

Здравствуйте уважаемые форумчанины....

Возникла следующая проблема....

есть односвязный список, который хранит номер ячейки, в котором находится не нулевой элемент и само значение:

template <class T>
struct SparseMatrixStorage
{
        long indexI;
        long indexJ;
        T info;
        SparseMatrixStorage *next;
};

Есть класс SparseMatrix, который работает с данной структурой.
В классе есть метод умножение двух матриц.
Она у меня не работает.. Вылетают разные ошибки....
Не могли бы вы мне помочь написать нормальный метод??????

template <class T>
class SparseMatrix : BaseMatrix<T>
{
        public:
                SparseMatrix(long lengthI, long lengthJ): BaseMatrix<T>(lengthI, lengthJ, false)
                {
                        head = NULL;
                        count = 0;
                };
                SparseMatrix():BaseMatrix<T>()
                {
                        head = NULL;
                        count = 0;
                };
                void operator = (T ** matrix)
                {
                        SparseMatrixStorage<T> *temp;
                        long dimension = 0;
                        for (long i = 0; i < GetCountI(); i++)
                                for (long j = 0; j < GetCountJ(); j++)
                                {
                                        if (matrix[i][j] != 0)
                                        {
                                                SparseMatrixStorage<T> *elem = new SparseMatrixStorage<T>;
                                                if(GetHead() == NULL)
                                                        SetHead(elem);
                                                else
                                                        SetPrevAddress(elem);
                                                try
                                                {
                                                        elem->info = matrix[i][j];
                                                        elem->indexI = i;
                                                        elem->indexJ = j;
                                                        elem->next = NULL;
                                                        SetPrev(elem);
                                                        dimension ++ ;
                                                        //......
                                                }
                                                catch(...)
                                                {
                                                        ShowMessage("Неправильно заполнена матрица");
                                                        break;
                                                }
                                        }
                                }
                        SetCount(dimension);
                };

[world12_tk]

      void operator * (T ** matrix)
                {
                        SparseMatrix result;
                        for (long i = 0; i < GetCountI(); i++)
                                for (long j = 0; j < GetCountJ(); j++)
                                {
                                        if (matrix[j][i] != 0)
                                        {
                                                SparseMatrixStorage<T> *elem = new SparseMatrixStorage<T>;
                                                if(head == NULL)
                                                        head = elem;
                                                else
                                                        prev->next = elem;
                                                try
                                                {
                                                        elem->info = matrix[j][i];
                                                        elem->indexI = j;
                                                        elem->indexJ = i;
                                                        elem->next = NULL;
                                                        prev = elem;
                                                        count ++ ;
                                                }
                                                catch(...)
                                                {
                                                        ShowMessage("Неправильно заполнена матрица");
                                                        break;
                                                }
                                        }
                                }
                        SparseMatrixStorage<T> *tempHead = GetHead();
                        SparseMatrixStorage<T> *elem = new SparseMatrixStorage<T>;
 
                        for (long i = 0; i < GetCountJ(); i++)
                        {
                                if (result.head == NULL)
                                        result.head = elem;
                                else
                                        result.prev->next = elem;
                                SparseMatrixStorage<T> *temp = tempHead;
                                SparseMatrixStorage<T> *temp2Head = head;
 
                                if (temp != NULL)
                                {
                                        for (long j = 0; j < GetCountJ(); j++)
                                        {
                                                long number = 0;
                                                SparseMatrixStorage<T> *temp2 = temp2Head;
                                                if (temp2 != NULL)
                                                {
                                                        for (long k = 0; k < GetCountI(); k++)
                                                        {
                                                                if ((temp != NULL) && (temp->indexI == j) && (temp2->indexJ == k))
                                                                {
                                                                        number += temp->info * temp2->info;
                                                                        temp = temp->next;
                                                                        temp2 = temp2->next;
                                                                }
 
                                                        }
                                                        temp2Head = temp2Head->next;
                                                        if (number != 0)
                                                        {
                                                                elem->indexI = i;
                                                                elem->indexJ = j;
                                                                elem->info = number;
                                                                result.prev = elem;
                                                                result.count++;
                                                                SparseMatrixStorage<T> *elem = new SparseMatrixStorage<T>;
                                                        }
                                                }
                                        }
                                        tempHead = tempHead->next;
                                }
 
 
                        }

[world12_tk]

                      
                        SetHead(result.head);
                        SetPrev(result.prev);
                        DeleteStorage();
                        SetCountI(GetCountJ());
                };
                private:
                void DeleteStorage()
                {
                        while(head != NULL)
                        {
                                SparseMatrixStorage<T> *pTemp = head;
                                head = head->next;
                                delete pTemp;
                        }
                };
                SparseMatrixStorage<T> * head;
                SparseMatrixStorage<T> * prev;
                long count;
};

Вызываю я так:

        SparseMatrix<long double> sparseMatrix (column,row);
        sparseMatrix = realMatrix;
        sparseMatrix * realMatrix2;
        column = row;
        OutputPrint(column,row,sparseMatrix.Print());

Заранее спасибо.....