Почему этот CPP-файл выдает ошибку во время mex -largeArrayDims?

ProgramBox

Почему этот CPP-файл выдает ошибку во время mex -largeArrayDims?

Post author:admin
Запись опубликована:6 мая, 2023
Post category:Вопросы по программированию

#c #matlab #mex

Вопрос:

Я загрузил код:

 /*================================================================
a_times_b_cmplx.c = used by a couple of mex functions
provide Matrix vector multiplications,
and solve triangular systems
(sparse matrix and full vector)

CSC_CmplxVecMult_CAB_double, CSR_CmplxVecMult_CAB_double,
CSCsymm_CmplxVecMult_CAB_double added by Mirko Visontai (10/24/2003)

*=================================================================*/
# include "math.h"

///*c<-a'*b */
//void scalar_product(
//  const int m, const int k, /*nb_rows, nb_columns*/
//  const double *a,
//  const double *b,
//  double *c)
//{
//  int i;
//  double d;
//  d = 0;
//  for (i=0;i!=m;i  )
//    d =a[i]*b[i];
//  c[0] = d;
//}

/*C<-a*A*B C*/
void CSC_VecMult_CaABC_double(const int m, const int k, const double alpha, const double *val, const int *indx, const int *pntrb,const double *b, double *c)
{
  int i,j,jb,je;

  for (i=0;i!=k;i  ){
    jb = pntrb[i];
    je = pntrb[i 1];
    for (j=jb;j!=je;j  )
      c[indx[j]]  = alpha * b[i] * val[j];
  }
}

/*C<-a*A'*B C*/
void CSR_VecMult_CaABC_double(const int k,  const int m, const double alpha,const double *val, const mwIndex *indx, const mwIndex *pntrb,const double *b,  double *c)
{
  double t;
  const double *pval;
  int i,j,jb,je;

  pval = val;
  for (i=0;i!=m;i  ) {
    t = 0;
    jb = pntrb[i];
    je = pntrb[i 1];
    for (j=jb;j!=je;j  )
      t  = alpha * b[indx[j]] * (*pval  );
    c[i]  = t;
  }
}

/*C<-A*b */
void CSC_VecMult_CAB_double(const int m, const int k, /*nb_rows, nb_columns*/ const double *val, const int *indx, const int *pntrb, const double *b, double *c)
{
  int i,j,jb,je;
  double *pc=c;
  for (i=0;i!=m;i  )
    *pc   = 0;

  for (i=0;i!=k;i  ){
    jb = pntrb[i];
    je = pntrb[i 1];
    for (j=jb;j!=je;j  )
      c[indx[j]]  =  b[i] * val[j];
  }
}

/*C<-A*b (complex)*/
void CSC_CmplxVecMult_CAB_double(const int m,  const int k, const double *valr, const double *vali,const int *indx, const int *pntrb,const double *br, const double *bi,double *cr, double *ci )
{
  int i,j,jb,je;
  double *pcr=cr;
  double *pci=ci;
  for (i=0;i!=m;i  ){
    *pcr   = 0.0;
    *pci   = 0.0;
  }

  for (i=0;i!=k;i  ){
    jb = pntrb[i];
    je = pntrb[i 1];
    for (j=jb;j!=je;j  ){
      cr[indx[j]]  =  (br[i] * valr[j]) - (bi[i] * vali[j]);
      ci[indx[j]]  =  (br[i] * vali[j])   (bi[i] * valr[j]);
    }
  }
}

/*C<-A'*b 
 plus rapide que CSC_VecMult_CAB_double */
void CSR_VecMult_CAB_double(const int k, const int m, const double *val, const int *indx, const int *pntrb, const double *b, double *c)
{
  double t;
  const double *pval;
  double *pc=c;
  int i,j,jb,je;

  for (i=0;i!=m;i  )
    *pc   = 0;

  pval = val;
  for (i=0;i!=m;i  ) {
    t = 0;
    jb = pntrb[i];
    je = pntrb[i 1];
    for (j=jb;j!=je;j  )
      t  =  b[indx[j]] * (*pval  );
    c[i]  = t;
  }
}

/*C<-A'*b (complex) 
 plus rapide que CSC_VecMult_CAB_double */
void CSR_CmplxVecMult_CAB_double(const int k, const int m, const double *valr, const double *vali, const int *indx, const int *pntrb, const double *br, const double *bi, double *cr, double *ci)
{
  double tr, ti;
  const double *pvalr;
  const double *pvali;
  double *pcr=cr;
  double *pci=ci;
  int i,j,jb,je;

  for (i=0;i!=m;i  ){
    *pcr   = 0.0;
    *pci   = 0.0;
  }

  pvalr = valr;
  pvali = vali;
  for (i=0;i!=m;i  ) {
    tr = 0.0;
    ti = 0.0;
    jb = pntrb[i];
    je = pntrb[i 1];
    for (j=jb;j!=je;j  ){
      tr  =  (br[indx[j]] * (*pvalr)) - (bi[indx[j]] * (*pvali));
      ti  =  (br[indx[j]] * (*pvali  ))   (bi[indx[j]] * (*pvalr  ));
    }
    cr[i]  = tr;
    ci[i]  = ti;
  }
}

/* C<-A*b (A is symmetric) */
void CSRsymm_VecMult_CAB_double(const int k, const int m, const double *val, const mwIndex *indx, const mwIndex *pntrb, const double *b, double *c)
{
  const double *pval;
  double *pc=c;
  int i,j;
  int jj;
  int rpntrb, rpntre;
  int index, nvals;

  for (i=0;i!=m;i  )
    *pc   = 0;
  pval = val;
  for (j=0;j!=k;j  ){
    rpntrb = pntrb[j];
    rpntre = pntrb[j 1];
    for (jj=rpntrb;jj!=rpntre;jj  ) {
      index = indx[jj];
      if ( index == j ) {
        c[j]  =  b[j] * (*pval  );
        continue;
      }
      if ( index > j ) {
        c[index]  =  b[j] * (*pval);
        c[j]  =  b[index] * (*pval  );
      }
      else {
        pval  ;
      }
    }
  }
}

/* C<-A*b (A is symmetric and complex) */
void CSRsymm_CmplxVecMult_CAB_double(const int k, const int m, const double *valr, const double *vali, const int *indx, const int *pntrb, const double *br, const double *bi, double *cr, double *ci)
{
  const double *pvalr, *pvali;
  double *pcr=cr;
  double *pci=ci;
  int i,j;
  int jj;
  int rpntrb, rpntre;
  int index, nvals;

  for (i=0;i!=m;i  ){
    *pcr   = 0.0;
    *pci   = 0.0;
  }

  pvalr = valr;
  pvali = vali;
  for (j=0;j!=k;j  ){
    rpntrb = pntrb[j];
    rpntre = pntrb[j 1];
    for (jj=rpntrb;jj!=rpntre;jj  ) {
      index = indx[jj];
      if ( index == j ) {
        cr[j]  =  (br[j] * (*pvalr)) - (bi[j] * (*pvali));
        ci[j]  =  (br[j] * (*pvali  ))   (bi[j] * (*pvalr  ));
        continue;
      }
      if ( index > j ) {
        cr[index]  =  (br[j] * (*pvalr)) - (bi[j] * (*pvali));
        ci[index]  =  (br[j] * (*pvali))   (bi[j] * (*pvalr));

        cr[j]  =  (br[index] * (*pvalr)) - (bi[index] * (*pvali));
        ci[j]  =  (br[index] * (*pvali  ))   (bi[index] * (*pvalr  ));
      }
      else {
        pvalr  ;
        pvali  ;
      }
    }
  }
}

/*C<-AB; with Lower triangular A*/
void CSC_VecTriangSlvLD_CAB_double(const int m, const double *val, const int *indx, const int *pntrb, const double *b, double *c)
{
  int i, j, jb, je;
  double *pc=c;
  double z; 

  for (i=0;i!=m;i  ){
    *pc = b[i];
    pc  ;
  }

  pc=c;
  for (i=0;i!=m;i  ) {
    jb = pntrb[i];
    je = pntrb[i 1];
    z =  pc[i] / val[jb];
    pc[i] = z;
    for (j=jb 1;j<je;j  ) {
      c[indx[j]] -= z*val[j];
    }
  }
}

/*C<-AB; with Upper triangular A*/
void CSC_VecTriangSlvUD_CAB_double(const int m, const double *val, const int *indx, const int *pntrb, const double *b, double *c)
{
  int i, j, jb, je, index;
  double *pc=c;
  double z; 

  for (i=0;i!=m;i  ){
    *pc = b[i];
    pc  ;
  }

  pc=c;
  for (i=m-1;i!=-1;i--) {
    jb = pntrb[i];
    je = pntrb[i 1]-1;
    z = pc[i] /val[je];
    pc[i] = z;
    for (j=jb;j<je;j  ) {
      c[indx[j]] -= z * val[j];
    }
  }
}

/*C<-A'B; where A is upper (little slower than CSC)*/
void CSR_VecTriangSlvLD_CAB_double(const int m, const double *val, const int *indx, const int *pntrb, const double *b, double *c)
{
  int i, j, jb, je, index;
  double *pc=c;
  double z; 
  double valtmp;

  pc=c;
  for (i=0;i!=m;i  ) {
    z = 0;
    jb = pntrb[i];
    je = pntrb[i 1];
    for (j=jb;j<je;j  ) {
      index = indx[j];
      if ( index == i ) {
        valtmp = val[j];
      } else {
        z  = c[index] * val[j];
      }
    }
    pc[i] = (b[i] - z) / valtmp;
  }
}

/*C<-A'B; where A is lower  (little slower than CSC)*/
void CSR_VecTriangSlvUD_CAB_double(const int m, const double *val, const int *indx, const int *pntrb, const double *b, double *c)
{
  int i, j, jb, je, index;
  double *pc=c;
  double valtmp;
  double z; 

  pc=c;
  for (i=m-1;i!=-1; i--) {
    z = 0;
    jb = pntrb[i];
    je =  pntrb[i 1];
    for (j=jb 1; j<je; j  ) {
      z  = c[indx[j]] * val[j];
    }
    pc[i] = (b[i] - z) / val[jb];
  }
}

/*C<-A*B, where A is (m,k) and B is (k,n)*/
void CSC_MatMult_CAB_double(const int m, const int n, const int k, const double *val, const int *indx, const int *pntrb, const double *b, const int ldb, double *c, const int ldc)
{
  int i,j,jb,je;
  double *pc=c;
  int l;

  for (l=0;l!=n;l  )
    for (i=0;i!=m;i  )
      *pc   = 0;

  for (l=0;l!=n;l  ) {
    for (i=0;i!=k;i  ){
      jb = pntrb[i];
      je = pntrb[i 1];
      for (j=jb;j!=je;j  )
        c[indx[j]]  =  b[i] * val[j];
    }
    /*c  = ldc; b  = ldb;  */
    c  = m; b  = m;
  }
}

Я использовал mex -largeArrayDims thisfile.cpp с выбором компилятора Visual C 2010. Однако он выдает мне эти ошибки:

mex -largeArrayDims a_times_b_cmplx.cpp

a_times_b_cmplx.cpp(45) : ошибка C4430: отсутствует спецификатор типа — предполагается значение int. Примечание: C не поддерживает default-int
a_times_b_cmplx.cpp(45) : ошибка C2143: синтаксическая ошибка: отсутствует ',' перед '*'

a_times_b_cmplx.cpp(54) : ошибка C2065: 'pntrb' : необъявленный идентификатор
a_times_b_cmplx.cpp(55) : ошибка C2065: 'pntrb' : необъявленный идентификатор
a_times_b_cmplx.cpp(57) : ошибка C2065: 'b' : необъявленный идентификатор
a_times_b_cmplx.cpp(57) : ошибка C2065: 'indx' : необъявленный идентификатор
a_times_b_cmplx.cpp(58) : ошибка C2065: 'c' : необъявленный идентификатор
a_times_b_cmplx.cpp(156) : ошибка C4430: отсутствует спецификатор типа — предполагается значение int. Примечание: C не поддерживает default-int
a_times_b_cmplx.cpp(156) : ошибка C2143: синтаксическая ошибка: отсутствует ',' перед '*'

a_times_b_cmplx.cpp(159) : ошибка C2065: 'c' : необъявленный идентификатор
a_times_b_cmplx.cpp(169) : ошибка C2065: 'pntrb' : необъявленный идентификатор
a_times_b_cmplx.cpp(170) : ошибка C2065: 'pntrb' : необъявленный идентификатор
a_times_b_cmplx.cpp(172) : ошибка C2065: 'indx' : необъявленный идентификатор
a_times_b_cmplx.cpp(174) : ошибка C2065: 'c' : необъявленный идентификатор
a_times_b_cmplx.cpp(174) : ошибка C2065: 'b' : необъявленный идентификатор
a_times_b_cmplx.cpp(178) : ошибка C2065: 'c' : необъявленный идентификатор
a_times_b_cmplx.cpp(178) : ошибка C2065: 'b' : необъявленный идентификатор
a_times_b_cmplx.cpp(180) : ошибка C2065: 'c' : необъявленный идентификатор
a_times_b_cmplx.cpp(180) : ошибка C2065: 'b' : необъявленный идентификатор

C:PROGRA~1MATLABR2010ABINMEX.PL : Ошибка: не удалось выполнить компиляцию 'a_times_b_cmplx.cpp' .

Я не знаю, почему я получаю эти ошибки. Это готовый лицензированный код, и я должен иметь возможность его редактировать. Можете ли вы назвать мне причину и показать способ исправить это. Я действительно не знаю, проблема в компиляторе или во мне.

Я думаю, это может быть из-за mwIndex . Но как я могу это решить?

1. Не хочу показаться грубым, но когда вы просматривали внешний вид этого вопроса, действительно ли выглядело так, будто весь этот код был правильно отформатирован до того, как вы нажали кнопку отправки? Вот это да.

2. Извините, я разбираюсь с этим часами, и это просто сводит меня с ума, я только хотел немного поторопиться.

3. Если сделать ваш вопрос совершенно болезненным для чтения (и длинным ), вы не получите ответы быстрее … ;-/

4. @ildjarn: Думаю, на этот раз мне повезло, и я решил проблему самостоятельно…

Ответ №1:

Я думаю, что ваш компилятор не может найти mwIndex тип (как вы подозреваете).

Быстрый поиск в Google по некоторой документации matlab показывает, что это mwIndex объявлено в matrix.h заголовочном файле.

Итак, попробуйте добавить

 #include <matrix.h>

…к вашему коду. Я подозреваю, что большинство этих сообщений об ошибках будут отправлены, хотя кто знает, скрываются ли за этим какие-либо другие проблемы.

Вопрос:

Комментарии:

Ответ №1:

Вам также может понравиться

получение максимального значения столбца метки времени в качестве подзапроса в улье

Проверьте, находится ли каждый элемент 1-го массива numpy в соответствующей строке 2-го массива numpy

Не удается открыть выпадающий список