|
00001 #include "AdaptiveData.hpp" 00002 int debugRefine ; 00003 int main() { 00004 Wavelets WC("Interpolet",4) ; 00005 // boundary conditions for A and B: 00006 int i, bc1[2][2], bc2[2][2] ; 00007 for (i=0; i<2; i++) { 00008 bc1[i][0]=bc1[i][1]= -1; // no BC for A 00009 bc2[i][0]=bc2[i][1]= 0 ; // homogeneous Dirichlet for B 00010 } 00011 00012 AdaptiveGrid<2> G(&WC) ; 00013 AdaptiveData<2> A(&G), B(&G), C(&G) ; 00014 G.SetPeriodicConditions(bc1) ; 00015 00016 ConstFunction F ; 00017 double p=1.0 ; 00018 F.SetParameters(&p) ; 00019 // get multiscale represenation of f(x,y)=1 00020 A.SetFunction(&F) ; 00021 00022 // set multiscale representation of g(x,y)=0 00023 B.Set(0.0) ; 00024 B.SetBoundaryConditions(bc2) ; 00025 00026 // try to add these 00027 C.Add(&A,&B) ; 00028 } 00029