当前时区为 UTC + 8 小时



发表新帖 回复这个主题  [ 13 篇帖子 ] 
作者 内容
1 楼 
 文章标题 : 大家来看看这个c程序的错误
帖子发表于 : 2010-05-19 22:11 

注册: 2010-05-19 22:01
帖子: 3
送出感谢: 0 次
接收感谢: 0 次
大家好!
这段时间想把跑一下splash2测试程序,但是我在编译其中的fft程序时
出现了如下错误,大家可以先忽略其中的MAIN_ENV宏定义。
其中的, error: expected ‘=’, ‘,’, ‘;’, ‘asm’ or ‘__attribute__’ before ‘struct’ 和 ‘struct GlobalMemory’ has no member named ‘totaltimes’错误是怎么回事? 应该不是程序的错误吧? 是不是由于ubuntu缺少什么包?
ubuntu : 9.10
gcc : 4.4.1
已安装build-essential

fft.c:
代码:
/*************************************************************************/
/*                                                                       */
/*  Copyright (c) 1994 Stanford University                               */
/*                                                                       */
/*  All rights reserved.                                                 */
/*                                                                       */
/*  Permission is given to use, copy, and modify this software for any   */
/*  non-commercial purpose as long as this copyright notice is not       */
/*  removed.  All other uses, including redistribution in whole or in    */
/*  part, are forbidden without prior written permission.                */
/*                                                                       */
/*  This software is provided with absolutely no warranty and no         */
/*  support.                                                             */
/*                                                                       */
/*************************************************************************/

/*************************************************************************/
/*                                                                       */
/*  Perform 1D fast Fourier transform using six-step FFT method          */
/*                                                                       */
/*  1) Performs staggered, blocked transposes for cache-line reuse       */
/*  2) Roots of unity rearranged and distributed for only local          */
/*     accesses during application of roots of unity                     */
/*  3) Small set of roots of unity elements replicated locally for       */
/*     1D FFTs (less than root N elements replicated at each node)       */
/*  4) Matrix data structures are padded to reduce cache mapping         */
/*     conflicts                                                         */
/*                                                                       */
/*  Command line options:                                                */
/*                                                                       */
/*  -mM : M = even integer; 2**M total complex data points transformed.  */
/*  -pP : P = number of processors; Must be a power of 2.                */
/*  -nN : N = number of cache lines.                                     */
/*  -lL : L = Log base 2 of cache line length in bytes.                  */
/*  -s  : Print individual processor timing statistics.                  */
/*  -t  : Perform FFT and inverse FFT.  Test output by comparing the     */
/*        integral of the original data to the integral of the data      */
/*        that results from performing the FFT and inverse FFT.          */
/*  -o  : Print out complex data points.                                 */
/*  -h  : Print out command line options.                                */
/*                                                                       */
/*  Note: This version works under both the FORK and SPROC models        */
/*                                                                       */
/*************************************************************************/

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define PAGE_SIZE               4096
#define NUM_CACHE_LINES        65536
#define LOG2_LINE_SIZE             4
#define PI                         3.1416
#define DEFAULT_M                 10
#define DEFAULT_P                  1

MAIN_ENV

#define SWAP_VALS(a,b) {double tmp; tmp=a; a=b; b=tmp;}

struct GlobalMemory {
  long id;
  LOCKDEC(idlock)
  BARDEC(start)
  long *transtimes;
  long *totaltimes;
  unsigned long starttime;
  unsigned long finishtime;
  unsigned long initdonetime;
} *Global;


long P = DEFAULT_P;
long M = DEFAULT_M;
long N;                  /* N = 2^M                                */
long rootN;              /* rootN = N^1/2                          */
double *x;              /* x is the original time-domain data     */
double *trans;          /* trans is used as scratch space         */
double *umain;          /* umain is roots of unity for 1D FFTs    */
double *umain2;         /* umain2 is entire roots of unity matrix */
long test_result = 0;
long doprint = 0;
long dostats = 0;
long transtime = 0;
long transtime2 = 0;
long avgtranstime = 0;
long avgcomptime = 0;
unsigned long transstart = 0;
unsigned long transend = 0;
long maxtotal=0;
long mintotal=0;
double maxfrac=0;
double minfrac=0;
double avgfractime=0;
long orig_num_lines = NUM_CACHE_LINES;     /* number of cache lines */
long num_cache_lines = NUM_CACHE_LINES;    /* number of cache lines */
long log2_line_size = LOG2_LINE_SIZE;
long line_size;
long rowsperproc;
double ck1;
double ck3;                        /* checksums for testing answer */
long pad_length;

void SlaveStart(void);
double TouchArray(double *x, double *scratch, double *u, double *upriv, long MyFirst, long MyLast);
double CheckSum(double *x);
void InitX(double *x);
void InitU(long N, double *u);
void InitU2(long N, double *u, long n1);
long BitReverse(long M, long k);
void FFT1D(long direction, long M, long N, double *x, double *scratch, double *upriv, double *umain2,
      long MyNum, long *l_transtime, long MyFirst, long MyLast, long pad_length, long test_result, long dostats);
void TwiddleOneCol(long direction, long n1, long j, double *u, double *x, long pad_length);
void Scale(long n1, long N, double *x);
void Transpose(long n1, double *src, double *dest, long MyNum, long MyFirst, long MyLast, long pad_length);
void CopyColumn(long n1, double *src, double *dest);
void Reverse(long N, long M, double *x);
void FFT1DOnce(long direction, long M, long N, double *u, double *x);
void PrintArray(long N, double *x);
void printerr(char *s);
long log_2(long number);

void srand48(long int seedval);
double drand48(void);

int main(int argc, char *argv[])
{
  long i;
  long c;
  extern char *optarg;
  long m1;
  long factor;
  long pages;
  unsigned long start;

  CLOCK(start);

  while ((c = getopt(argc, argv, "p:m:n:l:stoh")) != -1) {
    switch(c) {
      case 'p': P = atoi(optarg);
                if (P < 1) {
                  printerr("P must be >= 1\n");
                  exit(-1);
                }
                if (log_2(P) == -1) {
                  printerr("P must be a power of 2\n");
                  exit(-1);
                }
           break; 
      case 'm': M = atoi(optarg);
                m1 = M/2;
                if (2*m1 != M) {
                  printerr("M must be even\n");
                  exit(-1);
                }
           break; 
      case 'n': num_cache_lines = atoi(optarg);
                orig_num_lines = num_cache_lines;
                if (num_cache_lines < 1) {
                  printerr("Number of cache lines must be >= 1\n");
                  exit(-1);
                }
           break; 
      case 'l': log2_line_size = atoi(optarg);
                if (log2_line_size < 0) {
                  printerr("Log base 2 of cache line length in bytes must be >= 0\n");
                  exit(-1);
                }
           break; 
      case 's': dostats = !dostats;
           break;
      case 't': test_result = !test_result;
           break;
      case 'o': doprint = !doprint;
           break;
      case 'h': printf("Usage: FFT <options>\n\n");
                printf("options:\n");
                printf("  -mM : M = even integer; 2**M total complex data points transformed.\n");
                printf("  -pP : P = number of processors; Must be a power of 2.\n");
                printf("  -nN : N = number of cache lines.\n");
                printf("  -lL : L = Log base 2 of cache line length in bytes.\n");
                printf("  -s  : Print individual processor timing statistics.\n");
                printf("  -t  : Perform FFT and inverse FFT.  Test output by comparing the\n");
                printf("        integral of the original data to the integral of the data that\n");
                printf("        results from performing the FFT and inverse FFT.\n");
                printf("  -o  : Print out complex data points.\n");
                printf("  -h  : Print out command line options.\n\n");
                printf("Default: FFT -m%1d -p%1d -n%1d -l%1d\n",
                       DEFAULT_M,DEFAULT_P,NUM_CACHE_LINES,LOG2_LINE_SIZE);
      exit(0);
           break;
    }
  }

//  MAIN_INITENV(,80000000);

  N = 1<<M;
  rootN = 1<<(M/2);
  rowsperproc = rootN/P;
  if (rowsperproc == 0) {
    printerr("Matrix not large enough. 2**(M/2) must be >= P\n");
    exit(-1);
  }

  line_size = 1 << log2_line_size;
  if (line_size < 2*sizeof(double)) {
    printf("WARNING: Each element is a complex double (%ld bytes)\n",2*sizeof(double));
    printf("  => Less than one element per cache line\n");
    printf("     Computing transpose blocking factor\n");
    factor = (2*sizeof(double)) / line_size;
    num_cache_lines = orig_num_lines / factor;
  } 
  if (line_size <= 2*sizeof(double)) {
    pad_length = 1;
  } else {
    pad_length = line_size / (2*sizeof(double));
  }

  if (rowsperproc * rootN * 2 * sizeof(double) >= PAGE_SIZE) {
    pages = (2 * pad_length * sizeof(double) * rowsperproc) / PAGE_SIZE;
    if (pages * PAGE_SIZE != 2 * pad_length * sizeof(double) * rowsperproc) {
      pages ++;
    }
    pad_length = (pages * PAGE_SIZE) / (2 * sizeof(double) * rowsperproc);
  } else {
    pad_length = (PAGE_SIZE - (rowsperproc * rootN * 2 * sizeof(double))) /

                 (2 * sizeof(double) * rowsperproc);
    if (pad_length * (2 * sizeof(double) * rowsperproc) !=
        (PAGE_SIZE - (rowsperproc * rootN * 2 * sizeof(double)))) {
      printerr("Padding algorithm unsuccessful\n");
      exit(-1);
    }
  }

  Global = (struct GlobalMemory *) G_MALLOC(sizeof(struct GlobalMemory));
  x = (double *) G_MALLOC(2*(N+rootN*pad_length)*sizeof(double)+PAGE_SIZE);
  trans = (double *) G_MALLOC(2*(N+rootN*pad_length)*sizeof(double)+PAGE_SIZE);
  umain = (double *) G_MALLOC(2*rootN*sizeof(double)); 
  umain2 = (double *) G_MALLOC(2*(N+rootN*pad_length)*sizeof(double)+PAGE_SIZE);

  Global->transtimes = (long *) G_MALLOC(P*sizeof(long)); 
  Global->totaltimes = (long *) G_MALLOC(P*sizeof(long)); 
  if (Global == NULL) {
    printerr("Could not malloc memory for Global\n");
    exit(-1);
  } else if (x == NULL) {
    printerr("Could not malloc memory for x\n");
    exit(-1);
  } else if (trans == NULL) {
    printerr("Could not malloc memory for trans\n");
    exit(-1);
  } else if (umain == NULL) {
    printerr("Could not malloc memory for umain\n");
    exit(-1);
  } else if (umain2 == NULL) {
    printerr("Could not malloc memory for umain2\n");
    exit(-1);
  }

  x = (double *) (((unsigned long) x) + PAGE_SIZE - ((unsigned long) x) % PAGE_SIZE);
  trans = (double *) (((unsigned long) trans) + PAGE_SIZE - ((unsigned long) trans) % PAGE_SIZE);
  umain2 = (double *) (((unsigned long) umain2) + PAGE_SIZE - ((unsigned long) umain2) % PAGE_SIZE);

/* In order to optimize data distribution, the data structures x, trans,
   and umain2 have been aligned so that each begins on a page boundary.
   This ensures that the amount of padding calculated by the program is
   such that each processor's partition ends on a page boundary, thus
   ensuring that all data from these structures that are needed by a
   processor can be allocated to its local memory */

/* POSSIBLE ENHANCEMENT:  Here is where one might distribute the x,
   trans, and umain2 data structures across physically distributed
   memories as desired.
   
   One way to place data is as follows:

   double *base;
   long i;

   i = ((N/P)+(rootN/P)*pad_length)*2;
   base = &(x[0]);
   for (j=0;j<P;j++) {
    Place all addresses x such that (base <= x < base+i) on node j
    base += i;
   }

   The trans and umain2 data structures can be placed in a similar manner.

   */

  printf("\n");
  printf("FFT with Blocking Transpose\n");
  printf("   %ld Complex Doubles\n",N);
  printf("   %ld Processors\n",P);
  if (num_cache_lines != orig_num_lines) {
    printf("   %ld Cache lines\n",orig_num_lines);
    printf("   %ld Cache lines for blocking transpose\n",num_cache_lines);
  } else {
    printf("   %ld Cache lines\n",num_cache_lines);
  }
  printf("   %d Byte line size\n",(1 << log2_line_size));
  printf("   %d Bytes per page\n",PAGE_SIZE);
  printf("\n");

  BARINIT(Global->start, P);
  LOCKINIT(Global->idlock);
  Global->id = 0;
  InitX(x);                  /* place random values in x */

  if (test_result) {
    ck1 = CheckSum(x);
  }
  if (doprint) {
    printf("Original data values:\n");
    PrintArray(N, x);
  }

  InitU(N,umain);               /* initialize u arrays*/
  InitU2(N,umain2,rootN);

  /* fire off P processes */

  CREATE(SlaveStart, P);
  WAIT_FOR_END(P);

  if (doprint) {
    if (test_result) {
      printf("Data values after inverse FFT:\n");
    } else {
      printf("Data values after FFT:\n");
    }
    PrintArray(N, x);
  }

  transtime = Global->transtimes[0];
  printf("\n");
  printf("                 PROCESS STATISTICS\n");
  printf("            Computation      Transpose     Transpose\n");
  printf(" Proc          Time            Time        Fraction\n");
  printf("    0        %10ld     %10ld      %8.5f\n",
         Global->totaltimes[0],Global->transtimes[0],
         ((double)Global->transtimes[0])/Global->totaltimes[0]);
  if (dostats) {
    transtime2 = Global->transtimes[0];
    avgtranstime = Global->transtimes[0];
    avgcomptime = Global->totaltimes[0];
    maxtotal = Global->totaltimes[0];
    mintotal = Global->totaltimes[0];
    maxfrac = ((double)Global->transtimes[0])/Global->totaltimes[0];
    minfrac = ((double)Global->transtimes[0])/Global->totaltimes[0];
    avgfractime = ((double)Global->transtimes[0])/Global->totaltimes[0];
    for (i=1;i<P;i++) {
      if (Global->transtimes[i] > transtime) {
        transtime = Global->transtimes[i];
      }
      if (Global->transtimes[i] < transtime2) {
        transtime2 = Global->transtimes[i];
      }
      if (Global->totaltimes[i] > maxtotal) {
        maxtotal = Global->totaltimes[i];
      }
      if (Global->totaltimes[i] < mintotal) {
        mintotal = Global->totaltimes[i];
      }
      if (((double)Global->transtimes[i])/Global->totaltimes[i] > maxfrac) {
        maxfrac = ((double)Global->transtimes[i])/Global->totaltimes[i];
      }
      if (((double)Global->transtimes[i])/Global->totaltimes[i] < minfrac) {
        minfrac = ((double)Global->transtimes[i])/Global->totaltimes[i];
      }
      printf("  %3ld        %10ld     %10ld      %8.5f\n",
             i,Global->totaltimes[i],Global->transtimes[i],
             ((double)Global->transtimes[i])/Global->totaltimes[i]);
      avgtranstime += Global->transtimes[i];
      avgcomptime += Global->totaltimes[i];
      avgfractime += ((double)Global->transtimes[i])/Global->totaltimes[i];
    }
    printf("  Avg        %10.0f     %10.0f      %8.5f\n",
           ((double) avgcomptime)/P,((double) avgtranstime)/P,avgfractime/P);
    printf("  Max        %10ld     %10ld      %8.5f\n",
      maxtotal,transtime,maxfrac);
    printf("  Min        %10ld     %10ld      %8.5f\n",
      mintotal,transtime2,minfrac);
  }
  Global->starttime = start;
  printf("\n");
  printf("                 TIMING INFORMATION\n");
  printf("Start time                        : %16lu\n",
     Global->starttime);
  printf("Initialization finish time        : %16lu\n",
     Global->initdonetime);
  printf("Overall finish time               : %16lu\n",
     Global->finishtime);
  printf("Total time with initialization    : %16lu\n",
     Global->finishtime-Global->starttime);
  printf("Total time without initialization : %16lu\n",
     Global->finishtime-Global->initdonetime);
  printf("Overall transpose time            : %16ld\n",
         transtime);
  printf("Overall transpose fraction        : %16.5f\n",
         ((double) transtime)/(Global->finishtime-Global->initdonetime));
  printf("\n");

  if (test_result) {
    ck3 = CheckSum(x);
    printf("              INVERSE FFT TEST RESULTS\n");
    printf("Checksum difference is %.3f (%.3f, %.3f)\n",
      ck1-ck3, ck1, ck3);
    if (fabs(ck1-ck3) < 0.001) {
      printf("TEST PASSED\n");
    } else {
      printf("TEST FAILED\n");
    }
  }

  MAIN_END;
}


void SlaveStart()
{
  long i;
  long MyNum;
  double *upriv;
  long initdone;
  long finish;
  long l_transtime=0;
  long MyFirst;
  long MyLast;

  LOCK(Global->idlock);
    MyNum = Global->id;
    Global->id++;
  UNLOCK(Global->idlock);

  BARINCLUDE(Global->start);

/* POSSIBLE ENHANCEMENT:  Here is where one might pin processes to
   processors to avoid migration */

  BARRIER(Global->start, P);

  upriv = (double *) malloc(2*(rootN-1)*sizeof(double)); 
  if (upriv == NULL) {
    fprintf(stderr,"Proc %ld could not malloc memory for upriv\n",MyNum);
    exit(-1);
  }
  for (i=0;i<2*(rootN-1);i++) {
    upriv[i] = umain[i];
  }   

  MyFirst = rootN*MyNum/P;
  MyLast = rootN*(MyNum+1)/P;

  TouchArray(x, trans, umain2, upriv, MyFirst, MyLast);

  BARRIER(Global->start, P);

/* POSSIBLE ENHANCEMENT:  Here is where one might reset the
   statistics that one is measuring about the parallel execution */

  if ((MyNum == 0) || (dostats)) {
    CLOCK(initdone);
  }

  /* perform forward FFT */
  FFT1D(1, M, N, x, trans, upriv, umain2, MyNum, &l_transtime, MyFirst,
   MyLast, pad_length, test_result, dostats);

  /* perform backward FFT */
  if (test_result) {
    FFT1D(-1, M, N, x, trans, upriv, umain2, MyNum, &l_transtime, MyFirst,
     MyLast, pad_length, test_result, dostats);
  } 

  if ((MyNum == 0) || (dostats)) {
    CLOCK(finish);
    Global->transtimes[MyNum] = l_transtime;
    Global->totaltimes[MyNum] = finish-initdone;
  }
  if (MyNum == 0) {
    Global->finishtime = finish;
    Global->initdonetime = initdone;
  }
}


double TouchArray(double *x, double *scratch, double *u, double *upriv, long MyFirst, long MyLast)
{
  long i,j,k;
  double tot = 0.0;

  /* touch my data */
  for (j=0;j<2*(rootN-1);j++) {
    tot += upriv[j];
  }   
  for (j=MyFirst; j<MyLast; j++) {
    k = j * (rootN + pad_length);
    for (i=0;i<rootN;i++) {
      tot += x[2*(k+i)] + x[2*(k+i)+1] +
             scratch[2*(k+i)] + scratch[2*(k+i)+1] +
        u[2*(k+i)] + u[2*(k+i)+1];
    }
  } 
  return tot;
}


double CheckSum(double *x)
{
  long i,j,k;
  double cks;

  cks = 0.0;
  for (j=0; j<rootN; j++) {
    k = j * (rootN + pad_length);
    for (i=0;i<rootN;i++) {
      cks += x[2*(k+i)] + x[2*(k+i)+1];
    }
  }

  return(cks);
}


void InitX(double *x)
{
  long i,j,k;

  srand48(0);
  for (j=0; j<rootN; j++) {
    k = j * (rootN + pad_length);
    for (i=0;i<rootN;i++) {
      x[2*(k+i)] = drand48();
      x[2*(k+i)+1] = drand48();
    }
  }
}


void InitU(long N, double *u)
{
  long q;
  long j;
  long base;
  long n1;

  for (q=0; 1<<q<N; q++) { 
    n1 = 1<<q;
    base = n1-1;
    for (j=0; j<n1; j++) {
      if (base+j > rootN-1) {
   return;
      }
      u[2*(base+j)] = cos(2.0*PI*j/(2*n1));
      u[2*(base+j)+1] = -sin(2.0*PI*j/(2*n1));
    }
  }
}


void InitU2(long N, double *u, long n1)
{
  long i,j,k;

  for (j=0; j<n1; j++) { 
    k = j*(rootN+pad_length);
    for (i=0; i<n1; i++) { 
      u[2*(k+i)] = cos(2.0*PI*i*j/(N));
      u[2*(k+i)+1] = -sin(2.0*PI*i*j/(N));
    }
  }
}


long BitReverse(long M, long k)
{
  long i;
  long j;
  long tmp;

  j = 0;
  tmp = k;
  for (i=0; i<M; i++) {
    j = 2*j + (tmp&0x1);
    tmp = tmp>>1;
  }
  return(j);
}


void FFT1D(long direction, long M, long N, double *x, double *scratch, double *upriv, double *umain2,
           long MyNum, long *l_transtime, long MyFirst, long MyLast, long pad_length, long test_result, long dostats)
{
  long j;
  long m1;
  long n1;
  unsigned long clocktime1;
  unsigned long clocktime2;

  m1 = M/2;
  n1 = 1<<m1;

  BARRIER(Global->start, P);

  if ((MyNum == 0) || (dostats)) {
    CLOCK(clocktime1);
  }

  /* transpose from x into scratch */
  Transpose(n1, x, scratch, MyNum, MyFirst, MyLast, pad_length);
 
  if ((MyNum == 0) || (dostats)) {
    CLOCK(clocktime2);
    *l_transtime += (clocktime2-clocktime1);
  }

  /* do n1 1D FFTs on columns */
  for (j=MyFirst; j<MyLast; j++) {
    FFT1DOnce(direction, m1, n1, upriv, &scratch[2*j*(n1+pad_length)]);
    TwiddleOneCol(direction, n1, j, umain2, &scratch[2*j*(n1+pad_length)], pad_length);
  } 

  BARRIER(Global->start, P);

  if ((MyNum == 0) || (dostats)) {
    CLOCK(clocktime1);
  }
  /* transpose */
  Transpose(n1, scratch, x, MyNum, MyFirst, MyLast, pad_length);

  if ((MyNum == 0) || (dostats)) {
    CLOCK(clocktime2);
    *l_transtime += (clocktime2-clocktime1);
  }

  /* do n1 1D FFTs on columns again */
  for (j=MyFirst; j<MyLast; j++) {
    FFT1DOnce(direction, m1, n1, upriv, &x[2*j*(n1+pad_length)]);
    if (direction == -1)
      Scale(n1, N, &x[2*j*(n1+pad_length)]);
  }

  BARRIER(Global->start, P);

  if ((MyNum == 0) || (dostats)) {
    CLOCK(clocktime1);
  }

  /* transpose back */
  Transpose(n1, x, scratch, MyNum, MyFirst, MyLast, pad_length);

  if ((MyNum == 0) || (dostats)) {
    CLOCK(clocktime2);
    *l_transtime += (clocktime2-clocktime1);
  }

  BARRIER(Global->start, P);

  /* copy columns from scratch to x */
  if ((test_result) || (doprint)) { 
    for (j=MyFirst; j<MyLast; j++) {
      CopyColumn(n1, &scratch[2*j*(n1+pad_length)], &x[2*j*(n1+pad_length)]);
    } 
  } 

  BARRIER(Global->start, P);
}


void TwiddleOneCol(long direction, long n1, long j, double *u, double *x, long pad_length)
{
  long i;
  double omega_r;
  double omega_c;
  double x_r;
  double x_c;

  for (i=0; i<n1; i++) {
    omega_r = u[2*(j*(n1+pad_length)+i)];
    omega_c = direction*u[2*(j*(n1+pad_length)+i)+1]; 
    x_r = x[2*i];
    x_c = x[2*i+1];
    x[2*i] = omega_r*x_r - omega_c*x_c;
    x[2*i+1] = omega_r*x_c + omega_c*x_r;
  }
}


void Scale(long n1, long N, double *x)
{
  long i;

  for (i=0; i<n1; i++) {
    x[2*i] /= N;
    x[2*i+1] /= N;
  }
}


void Transpose(long n1, double *src, double *dest, long MyNum, long MyFirst, long MyLast, long pad_length)
{
  long i;
  long j;
  long k;
  long l;
  long m;
  long blksize;
  long numblks;
  long firstfirst;
  long h_off;
  long v_off;
  long v;
  long h;
  long n1p;
  long row_count;

  blksize = MyLast-MyFirst;
  numblks = (2*blksize)/num_cache_lines;
  if (numblks * num_cache_lines != 2 * blksize) {
    numblks ++;
  }
  blksize = blksize / numblks;
  firstfirst = MyFirst;
  row_count = n1/P;
  n1p = n1+pad_length;
  for (l=MyNum+1;l<P;l++) {
    v_off = l*row_count;
    for (k=0; k<numblks; k++) {
      h_off = firstfirst;
      for (m=0; m<numblks; m++) {
        for (i=0; i<blksize; i++) {
     v = v_off + i;
          for (j=0; j<blksize; j++) {
       h = h_off + j;
            dest[2*(h*n1p+v)] = src[2*(v*n1p+h)];
            dest[2*(h*n1p+v)+1] = src[2*(v*n1p+h)+1];
          }
        }
   h_off += blksize;
      }
      v_off+=blksize;
    }
  }

  for (l=0;l<MyNum;l++) {
    v_off = l*row_count;
    for (k=0; k<numblks; k++) {
      h_off = firstfirst;
      for (m=0; m<numblks; m++) {
        for (i=0; i<blksize; i++) {
     v = v_off + i;
          for (j=0; j<blksize; j++) {
            h = h_off + j;
            dest[2*(h*n1p+v)] = src[2*(v*n1p+h)];
            dest[2*(h*n1p+v)+1] = src[2*(v*n1p+h)+1];
          }
        }
   h_off += blksize;
      }
      v_off+=blksize;
    }
  }

  v_off = MyNum*row_count;
  for (k=0; k<numblks; k++) {
    h_off = firstfirst;
    for (m=0; m<numblks; m++) {
      for (i=0; i<blksize; i++) {
        v = v_off + i;
        for (j=0; j<blksize; j++) {
          h = h_off + j;
          dest[2*(h*n1p+v)] = src[2*(v*n1p+h)];
          dest[2*(h*n1p+v)+1] = src[2*(v*n1p+h)+1];
   }
      }
      h_off += blksize;
    }
    v_off+=blksize;
  }
}


void CopyColumn(long n1, double *src, double *dest)
{
  long i;

  for (i=0; i<n1; i++) {
    dest[2*i] = src[2*i];
    dest[2*i+1] = src[2*i+1];
  }
}


void Reverse(long N, long M, double *x)
{
  long j, k;

  for (k=0; k<N; k++) {
    j = BitReverse(M, k);
    if (j > k) {
      SWAP_VALS(x[2*j], x[2*k]);
      SWAP_VALS(x[2*j+1], x[2*k+1]);
    }
  }
}


void FFT1DOnce(long direction, long M, long N, double *u, double *x)
{
  long j;
  long k;
  long q;
  long L;
  long r;
  long Lstar;
  double *u1;
  double *x1;
  double *x2;
  double omega_r;
  double omega_c;
  double tau_r;
  double tau_c;
  double x_r;
  double x_c;

  Reverse(N, M, x);

  for (q=1; q<=M; q++) {
    L = 1<<q; r = N/L; Lstar = L/2;
    u1 = &u[2*(Lstar-1)];
    for (k=0; k<r; k++) {
      x1 = &x[2*(k*L)];
      x2 = &x[2*(k*L+Lstar)];
      for (j=0; j<Lstar; j++) {
   omega_r = u1[2*j];
        omega_c = direction*u1[2*j+1];
   x_r = x2[2*j];
        x_c = x2[2*j+1];
   tau_r = omega_r*x_r - omega_c*x_c;
   tau_c = omega_r*x_c + omega_c*x_r;
   x_r = x1[2*j];
        x_c = x1[2*j+1];
   x2[2*j] = x_r - tau_r;
   x2[2*j+1] = x_c - tau_c;
   x1[2*j] = x_r + tau_r;
   x1[2*j+1] = x_c + tau_c;
      }
    }
  }
}


void PrintArray(long N, double *x)
{
  long i, j, k;

  for (i=0; i<rootN; i++) {
    k = i*(rootN+pad_length);
    for (j=0; j<rootN; j++) {
      printf(" %4.2f %4.2f", x[2*(k+j)], x[2*(k+j)+1]);
      if (i*rootN+j != N-1) {
        printf(",");
      }
      if ((i*rootN+j+1) % 8 == 0) {
        printf("\n");
      }
    }
  }
  printf("\n");
  printf("\n");
}


void printerr(char *s)
{
  fprintf(stderr,"ERROR: %s\n",s);
}


long log_2(long number)
{
  long cumulative = 1, out = 0, done = 0;

  while ((cumulative < number) && (!done) && (out < 50)) {
    if (cumulative == number) {
      done = 1;
    } else {
      cumulative = cumulative * 2;
      out ++;
    }
  }

  if (cumulative == number) {
    return(out);
  } else {
    return(-1);
  }
}



出现了如下错误,大家可以先忽略其中的MAIN_ENV宏定义。
其中的, error: expected ‘=’, ‘,’, ‘;’, ‘asm’ or ‘__attribute__’ before ‘struct’ 和 ‘struct GlobalMemory’ has no member named ‘totaltimes’错误是怎么回事? 应该不是程序的错误吧? 是不是由于ubuntu缺少什么包?
ubuntu : 9.10
gcc : 4.4.1
已安装build-essential

代码:
tingting@ubuntu:~/program/splash2/codes/kernels/fft$ gcc fft.c
fft.c:60: error: expected ‘=’, ‘,’, ‘;’, ‘asm’ or ‘__attribute__’ before ‘struct’
fft.c: In function ‘main’:
fft.c:206: warning: format ‘%ld’ expects type ‘long int’, but argument 2 has type ‘unsigned int’
fft.c:241: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:242: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:305: error: ‘struct GlobalMemory’ has no member named ‘start’
fft.c:306: error: ‘struct GlobalMemory’ has no member named ‘idlock’
fft.c:307: error: ‘struct GlobalMemory’ has no member named ‘id’
fft.c:335: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:341: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:341: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:342: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:342: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:344: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:345: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:346: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:347: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:348: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:349: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:349: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:350: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:350: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:351: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:351: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:353: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:354: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:356: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:357: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:359: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:360: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:362: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:363: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:365: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:365: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:366: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:366: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:368: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:368: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:369: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:369: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:372: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:372: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:373: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:373: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:374: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:375: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:376: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:376: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:385: error: ‘struct GlobalMemory’ has no member named ‘starttime’
fft.c:389: error: ‘struct GlobalMemory’ has no member named ‘starttime’
fft.c:391: error: ‘struct GlobalMemory’ has no member named ‘initdonetime’
fft.c:393: error: ‘struct GlobalMemory’ has no member named ‘finishtime’
fft.c:395: error: ‘struct GlobalMemory’ has no member named ‘finishtime’
fft.c:395: error: ‘struct GlobalMemory’ has no member named ‘starttime’
fft.c:397: error: ‘struct GlobalMemory’ has no member named ‘finishtime’
fft.c:397: error: ‘struct GlobalMemory’ has no member named ‘initdonetime’
fft.c:401: error: ‘struct GlobalMemory’ has no member named ‘finishtime’
fft.c:401: error: ‘struct GlobalMemory’ has no member named ‘initdonetime’
fft.c:416: error: ‘MAIN_END’ undeclared (first use in this function)
fft.c:416: error: (Each undeclared identifier is reported only once
fft.c:416: error: for each function it appears in.)
fft.c: In function ‘SlaveStart’:
fft.c:431: error: ‘struct GlobalMemory’ has no member named ‘idlock’
fft.c:432: error: ‘struct GlobalMemory’ has no member named ‘id’
fft.c:433: error: ‘struct GlobalMemory’ has no member named ‘id’
fft.c:434: error: ‘struct GlobalMemory’ has no member named ‘idlock’
fft.c:436: error: ‘struct GlobalMemory’ has no member named ‘start’
fft.c:441: error: ‘struct GlobalMemory’ has no member named ‘start’
fft.c:457: error: ‘struct GlobalMemory’ has no member named ‘start’
fft.c:478: error: ‘struct GlobalMemory’ has no member named ‘transtimes’
fft.c:479: error: ‘struct GlobalMemory’ has no member named ‘totaltimes’
fft.c:482: error: ‘struct GlobalMemory’ has no member named ‘finishtime’
fft.c:483: error: ‘struct GlobalMemory’ has no member named ‘initdonetime’
fft.c: In function ‘FFT1D’:
fft.c:604: error: ‘struct GlobalMemory’ has no member named ‘start’
fft.c:624: error: ‘struct GlobalMemory’ has no member named ‘start’
fft.c:644: error: ‘struct GlobalMemory’ has no member named ‘start’
fft.c:658: error: ‘struct GlobalMemory’ has no member named ‘start’
fft.c:667: error: ‘struct GlobalMemory’ has no member named ‘start’


各位大虾帮帮忙吧! 小女子不胜感激


最后由 parallel 编辑于 2010-05-19 22:18,总共编辑了 1 次

页首
 用户资料  
 
2 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-19 22:14 
头像

注册: 2007-09-09 15:41
帖子: 218
地址: Beijing
送出感谢: 0 次
接收感谢: 1
代码:
apt-get install build-essential

?


页首
 用户资料  
 
3 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-19 22:16 

注册: 2010-05-19 22:01
帖子: 3
送出感谢: 0 次
接收感谢: 0 次
代码:
apt-get install build-essential


这个我已经安装了


页首
 用户资料  
 
4 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-19 22:19 
头像

注册: 2006-07-02 11:16
帖子: 12522
地址: 廣州
送出感谢: 0 次
接收感谢: 8
gcc -E 展开宏看那个 struct 前是什么东西


_________________
^_^ ~~~
要理解递归,首先要理解递归。

地球人都知道,理论上,理论跟实际是没有差别的,但实际上,理论跟实际的差别是相当大滴。


页首
 用户资料  
 
5 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-19 22:21 
头像

注册: 2008-03-25 15:49
帖子: 25877
地址: 谁知道?
送出感谢: 8
接收感谢: 10
这是人看的么。。


页首
 用户资料  
 
6 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-19 22:23 

注册: 2010-05-19 22:01
帖子: 3
送出感谢: 0 次
接收感谢: 0 次
pocoyo 写道:
这是人看的么。。



不好意思 程序帖的长了点


页首
 用户资料  
 
7 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-19 22:29 
头像

注册: 2008-07-27 8:51
帖子: 717
地址: 广西玉林|广东深圳
系统: Ubuntu/Windows
送出感谢: 1
接收感谢: 3
parallel 写道:
各位大虾帮帮忙吧! 小女子不胜感激


请问小女子芳龄几何? :em09


_________________
走过去了也便有了路
http://www.ptpt52.com/


页首
 用户资料  
 
8 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-20 23:38 
头像

注册: 2007-09-09 15:41
帖子: 218
地址: Beijing
送出感谢: 0 次
接收感谢: 1
代码:
  LOCKDEC(idlock)
  BARDEC(start)

这两个一定是宏,在你某个头文件里把,没有包含进来,宏无法展开


页首
 用户资料  
 
9 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-20 23:57 
头像

注册: 2007-09-09 15:41
帖子: 218
地址: Beijing
送出感谢: 0 次
接收感谢: 1
和 M4 宏 有关? 不清楚了
http://www.gnu.org/software/m4/manual/m4.html


页首
 用户资料  
 
10 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2010-05-27 20:59 

注册: 2006-11-08 16:33
帖子: 286
送出感谢: 0 次
接收感谢: 0 次
我也是闲的,下载下来看了一眼
同ls,发现不能只用gcc,还要用m4
Lz看一眼makefile就明白了


页首
 用户资料  
 
11 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2013-05-13 12:15 

注册: 2013-05-13 12:10
帖子: 1
系统: UBUNTU
送出感谢: 0 次
接收感谢: 0 次
这是splash2的FFT么?


页首
 用户资料  
 
12 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2013-05-13 13:13 

注册: 2009-11-07 19:24
帖子: 451
送出感谢: 0 次
接收感谢: 1
:em20 从来不会用指针。


页首
 用户资料  
 
13 楼 
 文章标题 : Re: 大家来看看这个c程序的错误
帖子发表于 : 2013-06-11 16:45 

注册: 2013-06-02 15:44
帖子: 13
系统: ubuntu12.04 lst
送出感谢: 0 次
接收感谢: 1
你的GlobalMemory结构体里面右面所有的定义都没有,你检查一下你的这个结构体。把没有的成员加上就可以了。


页首
 用户资料  
 
显示帖子 :  排序  
发表新帖 回复这个主题  [ 13 篇帖子 ] 

当前时区为 UTC + 8 小时


在线用户

正在浏览此版面的用户:没有注册用户 和 2 位游客


不能 在这个版面发表主题
不能 在这个版面回复主题
不能 在这个版面编辑帖子
不能 在这个版面删除帖子
不能 在这个版面提交附件

前往 :  
本站点为公益性站点,用于推广开源自由软件,由 DiaHosting VPSBudgetVM VPS 提供服务。
我们认为:软件应可免费取得,软件工具在各种语言环境下皆可使用,且不会有任何功能上的差异;
人们应有定制和修改软件的自由,且方式不受限制,只要他们自认为合适。

Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group
简体中文语系由 王笑宇 翻译