WTDaub.cc

00001 // Example for non-adaptive wavelet transforms
00002 
00003 # include "UniformData.hpp"
00004 # include "Function.hpp"
00005 
00006 int debugRefine ;
00007 
00008 int main() {
00009  Wavelets WC("Daubechies",3) ;
00010  
00011  int L=7 ;
00012  int i,a[2]={0},e[2]={1<<L,1<<L} ;
00013  UniformData<2> A(a,e,&WC), B(a,e,&WC) ;
00014 
00015  // read nodal values
00016  A.ReadUDF("../../Data/Test/A") ;
00017 
00018  // compute the coefficients for the isotropic Meyer-wavelets
00019  B.MRATransform(&A) ;
00020  B.WriteUDF("../../Data/Test/Bd1") ;
00021 
00022  // tensor product qudrature for the scaling function coefficients
00023  for (i=0; i<2; i++) A.ApplyOp(A.Ext.BC[i], &A , i, SHORTQUAD) ;
00024  B.MRATransform(&A) ;
00025  B.WriteUDF("../../Data/Test/Bd2") ; 
00026 
00027 }

Generated at Mon Aug 19 10:02:32 2002 for AWFD by