Thread: Error from Matlab to C (just one line code)

Hi people. I have en error translating a function from Matlab to C/C++ :

Code:

radius(floor(rows/2+1),floor(cols/2+1)) = 1;
 % Get rid of the 0 radius value at the 0
 % frequency point (now at top-left corner)
 % so that taking the log of the radius will 
 % not cause trouble.

What I did is:

Code:

radius[rows/2][cols/2]=1;

I m sure it s there the problem, because before this line I can print what I want and I can see it at run-time, but after this line I can not print nothing and the work is like blocked. Do you know the error?

Please show where you define the radius[][] array. Btw, this is the C++ forum. There's a separate forum for C questions.

Code:

double **createmat(int rows,int cols){

Code:

    double **m;

    m=new double*[rows];
     for(int i=0;i<rows;i++)
         m[i]=new double[cols];

    return m;
}

In main:

Code:

typedef double**matriz;

matriz radius,theta;

Both are built in a function

Hi, the problem is that before using it I used vectors of vectors and I had a problem about glibc. My teacher of University told me I have it because of vector of vectors utilization. So I used this form to do it. Also with that error we did not understand where were the error neither. Can I have an example of what you mean saying a smallest and simplest compilable program. Because to arrive where is that invocation, before there is a lot of code. Maybe I can send you the function where I have the problem.

All my project compile. But running it, it stopped in this line his work.This is the function where i have the problem:

Code:

void pre_filter_Computations(matriz radius,matriz theta,int cols,int rows){

  double X[rows][cols],Y[rows][cols];
  double x[cols],y[rows];
  double epsilon = 0.0001;

  for(int i=0;i<cols;i++){
      x[i]=((double)(i-cols)/2)/((double)cols/2);
  }


  for(int z=0;z<rows;z++){
       y[z]=-(((double)(z-rows)/2)/((double)rows/2));
   }


 for(int m=0;m<cols;m++){
  for(int n=0;n<rows;n++){
       X[m][n]=x[m];
       Y[m][n]=y[n];
  }
 }

      for(int a=0;a<rows;a++){
            for(int b=0;b<cols;b++){

                 X[a][b] = pow(X[a][b],2);
                 Y[a][b] = pow(Y[a][b],2);
                 X[a][b] = X[a][b] + Y[a][b];
                  radius[a][b] = sqrt(X[a][b]);

          }
      }
    printf("Hello\n"); // it print
    radius[rows/2][cols/2]=1;
   printf("Hello\n"); //doesn' print
      for(int a=0;a<rows;a++){
           for(int b=0;b<cols;b++){
               radius [a][b]= radius[a][b] + epsilon;
               theta[a][b] = atan2(Y[a][b],X[a][b])*180/PI;
           }
      }
}