// fourier.c
# include "alfa.h"

// Per numeri complessi.
nc Valic (nc *A, int N, nc z)
{if (N==1) return A[0];
if (N==2) return Ncadd(A[0],Ncmolt(z,A[1]));
int M=N/2; int k; nc z2;
nc a1[M],a2[M]; nc *X,*X1,*X2;
for (k=0,X=A,X1=a1,X2=a2;k<M;k++)
{*(X1++)=*(X++); *(X2++)=*(X++);}
z2=Ncmolt(z,z);
return Ncadd(Valic(a1,M,z2),
  Ncmolt(z,Valic(a2,M,z2)));}

// N deve essere una poteza di 2. I coefficienti
// sono elencati a partire dal termine costante.
double Valicreale  (Double A, int N, double x)
{if (N==1) return A[0];
if (N==2) return A[0]+x*A[1];
int M=N/2; int k; double x2;
double a1[M],a2[M]; Double X,X1,X2;
for (k=0,X=A,X1=a1,X2=a2;k<M;k++)
{*(X1++)=*(X++); *(X2++)=*(X++);}
x2=x*x;
return Valicreale (a1,M,x2)
  +x*Valicreale(a2,M,x2);}