Definition at line 27 of file MakeGaussian1d.h. Referenced by makeFilter(). 00028 { 00029 if (!order) 00030 return H0; 00031 00032 double a = 1./(sigma*sigma); 00033 double ax = a*x; 00034 double Hnmin1 = H0; 00035 double Hn = -ax*H0; 00036 00037 int n; 00038 for( n = 2; n <= order; n++ ) { 00039 double Hnmin2 = Hnmin1; 00040 Hnmin1 = Hn; 00041 Hn = a*(1-n)*Hnmin2 - ax*Hnmin1; 00042 } 00043 00044 return (order%2) ? -Hn : Hn; 00045 }
|