JSOC_Documentation: proj/globalhs/libs/pkbgn/predic.c Source File

00001 // $Header: /home/cvsuser/cvsroot/JSOC/proj/globalhs/libs/pkbgn/predic.c,v 1.1 2013/04/28 07:50:36 tplarson Exp $
00002 
00003 void predic(data,ndata,d,npoles,future,nfut)
00004 float data[],d[],future[];
00005 int ndata,npoles,nfut;
00006 {
00007     int k,j;
00008     float sum,discrp,*reg,*vector();
00009     void free_vector();
00010 
00011     reg=vector(1,npoles);
00012     for (j=1;j<=npoles;j++) reg[j]=data[ndata+1-j];
00013     for (j=1;j<=nfut;j++) {
00014         discrp=0.0;
00015         sum=discrp;
00016         for (k=1;k<=npoles;k++) sum += d[k]*reg[k];
00017         for (k=npoles;k>=2;k--) reg[k]=reg[k-1];
00018         future[j]=reg[1]=sum;
00019     }
00020     free_vector(reg,1,npoles);
00021 }