[file head]:
#include"tch.h"
#define send(x, y, n); {{int k; for(k=0; k<n; k++) (y)[k]=(x)[k];}}

int **int_matrix_alloc(int row, int col)
{
int i;
int **m;

m=(int **)malloc(row*sizeof(int*));

for(i=0; i<row; i++)
m[i]=(int *)malloc(col*sizeof(int));

return m;
}
/*
void free_int_matrix(int **m, int row)
{
int i;

for(i=0; i<row; i++) free((int*) m[i]);
free((int*) m);
}
*/
void fixed_spline_smoothing(double data[], int n, int iter)
{
int i, j;
double *data1;

data1=(double *)malloc(n*sizeof(double));
for(i=0; i<iter; i++)
{
send(data, data1, n);
for(j=1; j<n-1; j++)
data[j]=(data1[j-1]+2.0*data1[j]+data1[j+1])/4.0;
}
free(data1);
}

double **db_matrix(int row, int col)
{
int i;
double **m;

m=(double **)malloc(row*sizeof(double *));
for(i=0; i<row; i++)
m[i]=(double *)malloc(col*sizeof(double));
return m;
}

void free_db_matrix(double **m, int row)
{
int i;

for(i=0; i<row; i++)
free((double *)m[i]);
free((double *)m);
}

... ...