Как получить итераторы, если тип шаблона сам шаблон? - PullRequest
Новая
       21

Как получить итераторы, если тип шаблона сам шаблон?

2 голосов
/ 13 декабря 2010

Мой первый шаблонный класс попробовал упростить преобразование из некоторого надбавленного сжатого_матрицы в формат PETSc Sparse AIJ, но у меня возникли проблемы с итератором typedef в самом начале. Таким образом, T - сжатая_матрица, наивное использование приведенного ниже кода действительно приводит к ошибкам времени компиляции, как следует вместо этого определять эти итераторы. Я должен прочитать немного больше об этих шаблонах ... 

template <class T, class T1>
int converter<T, T1>::convertMBo2Pe( const T& A,
                           T1 A_ ){
    PetscErrorCode ierr;
    int cntNnz = 0;
    typedef T::iterator1 i1_t;
    typedef T::iterator2 i2_t;
    //int nnz[ooFelieMatrix.size1()];
    int nnz[A.size1()];
    unsigned ind=0; 
    //get information about the matrix

    double* vals = NULL;
    for (i1_t i1 = A.begin1(); i1 != A.end1(); ++i1)
    {
      nnz[ind] = distance(i1.begin(), i1.end()); 
      ind++; 
    }
    // create the matrix depending 
    // on the values of the nonzeros
    // on each row
    ierr = MatCreateSeqAIJ( PETSC_COMM_SELF, A.size1(), 
                                                  A.size2(), cntNnz, nnz, 
                                                  A_ );
    PetscInt rInd = 0, cInd=0; 
    PetscInt* rCount, dummy; 
    rCount = &dummy;
    // pointer to values in a row
    PetscScalar*   valsOfRowI = NULL;
    PetscInt*  colIndexOfRowI = NULL;
    PetscInt rC = 1;
    for(i1_t i1 = A.begin1(); i1 != A.end1(); ++i1)
    {
      // allocate space for the values of row I
      valsOfRowI     = new PetscScalar[nnz[rInd]];
      colIndexOfRowI = new PetscInt[nnz[rInd]];
      for(i2_t i2 = i1.begin(); i2 != i1.end(); ++i2)
      {
         colIndexOfRowI[cInd] = i2.index2();
         valsOfRowI[cInd]     = *i2;
         cInd++;
      }
      // setting one row each time
      *rCount = rInd;
      MatSetValues( A_, rC, rCount, nnz[rInd], 
                                  colIndexOfRowI, valsOfRowI,
                  INSERT_VALUES );
      // delete
      delete [] valsOfRowI;
      delete [] colIndexOfRowI;
      rInd++; cInd = 0;
    }
    //
    MatAssemblyBegin( A_, MAT_FINAL_ASSEMBLY );
    MatAssemblyEnd( A_, MAT_FINAL_ASSEMBLY );
    // return 
    return 0;   
}

1 Ответ

1 голос
/ 13 декабря 2010

Вы должны явно указать компилятору, что T::iterator1 - это имя типа (в отличие от статической переменной или чего-то подобного).

typedef typename T::iterator1 i1_t;

...