int main (int argc, char *argv[]) { int nmax, l, lc, n, m; Cheby12WorkSpace *wc=0; double *ptr,*phi; if (argc <3) { printf("chebyshev12 requires nmax as input argument followed by phi values\n"); return 0; } nmax=atoi(argv[1]); l=argc-2; if ((phi=(double*)calloc(l,sizeof(double)))==0) { printf("chebyshev12 could not allocate memory for %d phi-values\n",l); return 0; } if ((wc=chebyshev12_alloc(nmax,l))==0) { printf("chebyshev12 could not allocate memory for n=%d\n and l=%d\n",nmax, l); return 0; } for (n=0; n<l; n++) { phi[n]=atof(argv[n+2]); } chebyshev12(phi,wc); ptr=wc->t; for (lc=0; lc<l; lc++) { printf("pt %d:\n",lc); for (m=-nmax; m<=nmax; m++) { printf("T[%2d](%7.4f)=%7.4f\n",m,phi[lc],*ptr++); } } chebyshev12_free(wc); free(phi); return 1; }
static void allocMTXChdata (MTXCh *x) { x->phi=(double*)calloc(x->l,sizeof(double)); x->wc=chebyshev12_alloc(x->nmax,x->l); x->list_ch=(t_atom*)calloc(x->l*(2*x->nmax+1)+2,sizeof(t_atom)); }