/*  ftuns.c    CCMATH mathematics library source code.
 *
 *  Copyright (C)  2000   Daniel A. Atkinson    All rights reserved.
 *  This code may be redistributed under the terms of the GNU library
 *  public license (LGPL). ( See the lgpl.license file for details.)
 * ------------------------------------------------------------------------
 */
struct complex {double re,im;};
ftuns(struct complex **pt,int n)
{ struct complex **p,**q; int j; double x,y,u,v,h=2.;
  p=pt+1; q=pt+n-1;
  if(n%2==0) n/=2; else n=n/2+1;
  for(j=1; j<n ;++j){
    x=((*p)->re+(*q)->re)/h; y=((*p)->im-(*q)->im)/h;
    u=((*p)->im+(*q)->im)/h; v=((*p)->re-(*q)->re)/h;
    (*p)->re=x; (*p++)->im=y;
    (*q)->re=u; (*q--)->im=v;
   }
}


syntax highlighted by Code2HTML, v. 0.9.1