读简单的源码1-6.c
最近还是经常读gawk1.01,grep的代码，感觉grep的代码还是比较好读，毕竟比较简单，如果结合《精通正则表达式》来读，就基本上能理解了。
但gawk就复杂了，要懂bison，如果不懂这个，你C语言来通，还是不能串起来。所以要补充编译原理的知识。
最近还想读linux0.11的源码。感觉还是要把汇编搞通，才是。
于是，手头就找些简单的代码来读。现在感觉，如果不看到代码，心里就有些慌。使用软件时，也希望能是开源的。这样，就更有控制的感觉。
今天读了一段代码。感觉很是难懂。关键是其中：
typedef struct
{
  ElemType *base; /* 数组元素基址，由InitArray分配 */
  int dim; /* 数组维数 */
  int *bounds; /* 数组维界基址，由InitArray分配 */
  int *constants; /* 数组映象函数常量基址，由InitArray分配 */
}Array;

对constants总是不能理解。作者给它赋值的过程也很怪。
  (*A).constants[dim-1]=1;
  for(i=dim-2;i>=0;--i)
    (*A).constants[i]=(*A).bounds[i+1]*(*A).constants[i+1];
  return OK;
作者实现了数组的操作，感觉写得还是很流畅的。作者的演示程序是三维数组，我改成了二维数组。当然先要编译后执行一次。结果如下：
 ./a.out
A.bounds=3 4
A.constants=4 1
3ҳ40оԪ:
A[0][0]= 0
A[0][1]= 1
A[0][2]= 2
A[0][3]= 3

A[1][0]=10
A[1][1]=11
A[1][2]=12
A[1][3]=13

A[2][0]=20
A[2][1]=21
A[2][2]=22
A[2][3]=23

A.base=
   0   1   2   3
  10  11  12  13
  20  21  22  23

----
4,
1,
读代码一定要编译，调试，因此，除非你能在头脑中很好的模拟代码的执行，否则，读代码真是要在电脑上读。读别人的代码，感觉自己的水平真是差呀。
现在把代码放后面。
#include<stdio.h> /* EOF(=^Z或F6),NULL */
#include<math.h> /* floor(),ceil(),abs() */
#include<stdarg.h>
#define OK 1
#define ERROR 0
#define MAX_ARRAY_DIM 8
#include<stdlib.h>
#define OVERFLOW 88888
#define UNDERFLOW 999999
typedef int ElemType;
typedef int Status; /* Status是函数的类型,其值是函数结果状态代码，如OK等 */
typedef struct
{
  ElemType *base; /* 数组元素基址，由InitArray分配 */
  int dim; /* 数组维数 */
  int *bounds; /* 数组维界基址，由InitArray分配 */
  int *constants; /* 数组映象函数常量基址，由InitArray分配 */
}Array;
/* 顺序存储数组的基本操作*/
Status InitArray(Array *A,int dim,...)
{ /* 若维数dim和各维长度合法，则构造相应的数组A，并返回OK */
  int elemtotal=1,i; /* elemtotal是元素总值 */
  va_list ap;
  if(dim<1||dim>MAX_ARRAY_DIM)
    return ERROR;
  (*A).dim=dim;
  (*A).bounds=(int *)malloc(dim*sizeof(int));
  if(!(*A).bounds)
    exit(OVERFLOW);
  va_start(ap,dim);
  for(i=0;i<dim;++i)
  {
    (*A).bounds[i]=va_arg(ap,int);
    if((*A).bounds[i]<0)
      return UNDERFLOW; /* 在math.h中定义为4 */
    elemtotal*=(*A).bounds[i];
  }
  va_end(ap);
  (*A).base=(ElemType *)malloc(elemtotal*sizeof(ElemType));
  if(!(*A).base)
    exit(OVERFLOW);
  (*A).constants=(int *)malloc(dim*sizeof(int));
  if(!(*A).constants)
    exit(OVERFLOW);
  (*A).constants[dim-1]=1;
  for(i=dim-2;i>=0;--i)
    (*A).constants[i]=(*A).bounds[i+1]*(*A).constants[i+1];
  return OK;
}
Status DestroyArray(Array *A)
{ /* 销毁数组A */
  if((*A).base)
  {
    free((*A).base);
    (*A).base=NULL;
  }
  else
    return ERROR;
  if((*A).bounds)
  {
    free((*A).bounds);
    (*A).bounds=NULL;
  }
  else
    return ERROR;
  if((*A).constants)
  {
    free((*A).constants);
    (*A).constants=NULL;
  }
  else
    return ERROR;
  return OK;
}
Status Locate(Array A,va_list ap,int *off) /* Value()、Assign()调用此函数 */
{ /* 若ap指示的各下标值合法，则求出该元素在A中的相对地址off */
  int i,ind;
  *off=0;
  for(i=0;i<A.dim;i++)
  {
    ind=va_arg(ap,int);
    if(ind<0||ind>=A.bounds[i])
      return OVERFLOW;
    *off+=A.constants[i]*ind;
  }
  return OK;
}
Status Value(ElemType *e,Array A,...)
{ /* ...依次为各维的下标值，若各下标合法，则e被赋值为A的相应的元素值 */
  va_list ap;
  Status result;
  int off;
  va_start(ap,A);
  if((result=Locate(A,ap,&off))==OVERFLOW) /* 调用Locate() */
    return result;
  *e=*(A.base+off);
  return OK;
}
Status Assign(Array *A,ElemType e,...)
{ /* ...依次为各维的下标值，若各下标合法，则将e的值赋给A的指定的元素 */
  va_list ap;
  Status result;
  int off;
  va_start(ap,e);
  if((result=Locate(*A,ap,&off))==OVERFLOW) /* 调用Locate() */
    return result;
  *((*A).base+off)=e;
  return OK;
}
void main()
{
  Array A;
  //int i,j,k,*p,dim=3,bound1=3,bound2=4,bound3=2; /* a[3][4][2]数组 */
  int i,j,k,*p,dim=2,bound1=3,bound2=4; /* a[3][4][2]数组 */
  ElemType e,*p1;
//  InitArray(&A,dim,bound1,bound2,bound3); /* 构造3＊4＊2的3维数组A */
  InitArray(&A,dim,bound1,bound2); /* 构造3＊4＊2的3维数组A */
  p=A.bounds;
  printf("A.bounds=");
  for(i=0;i<dim;i++) /* 顺序输出A.bounds */
    printf("%d ",*(p+i));
  p=A.constants;
  printf("\nA.constants=");
  for(i=0;i<dim;i++) /* 顺序输出A.constants */
    printf("%d ",*(p+i));
  //printf("\n%d页%d行%d列矩阵元素如下:\n",bound1,bound2,bound3);
  printf("\n%d页%d行%d列矩阵元素如下:\n",bound1,bound2);
  for(i=0;i<bound1;i++)
  {
    for(j=0;j<bound2;j++)
    {
      //for(k=0;k<bound3;k++)
      {
        //Assign(&A,i*100+j*10+k,i,j,k); /* 将i*100+j*10+k赋值给A[i][j][k] */
        Assign(&A,i*10+j*1,i,j); /* 将i*100+j*10+k赋值给A[i][j][k] */
        Value(&e,A,i,j); /* 将A[i][j][k]的值赋给e */
        printf("A[%d][%d]=%2d ",i,j,e); /* 输出A[i][j][k] */
      }
      printf("\n");
    }
    printf("\n");
  }
  p1=A.base;
  printf("A.base=\n");
  for(i=0;i<bound1*bound2;i++) /* 顺序输出A.base */
  {
    printf("%4d",*(p1+i));
    if(i%(bound2)==bound2-1)
      printf("\n");
  }
  printf("\n----\n");
  for(i=0;i<A.dim;i++){
      printf("%d,\n",A.constants[i]);
  }
  DestroyArray(&A);
}