134135136137138139140141142143144
} } } // Compare Daubs to Daub4 for the case of 4 coefficients Daub4 db4 = new Daub4(); sbeps=1.e-15; for (i=0;i<N;i++) { for (j=0;j<N;j++) b[j]=a[j]=0.0; b[i]=a[i]=1.0/sqrt(N/2.0); wt1(a,-1,ds4);
3132333435363738394041
// Test wt1 System.out.println("Testing wt1"); Daub4 d4=new Daub4(); for (i=0;i<N;i++) { for (j=0;j<N;j++) a[j]=0.0; a[i]=1.0/sqrt(N/2.0); wt1(a,-1,d4); // for (j=0;j<N;j++) System.out.printf(setw(8) << a[j] << " ";
3233343536373839404142
// Test wtn System.out.println("Testing wtn"); Daub4 d4=new Daub4(); Ran myran = new Ran(17); nn[0]=N; nn[1]=M; for (i=0;i<NTRY;i++) { for (j=0;j<M*N;j++) a[j]=0.0;