Research Group of Prof. Dr. M. Griebel
Institute for Numerical Simulation
maximize
next up previous
Next: Solution of the Navier-Stokes Up: Examples Previous: Calculation of Partial Derivatives


Adaptive Solution of a Poisson Equation

In this example a Poisson problem with Dirichlet and/or Neumann boundary conditions is solved. The solution is the spiky function from our first example. The non-homogeneous Dirichlet or Neumann-values are put on the right hand side to obtain a problem with homogeneous boundary conditions. The common solve-refine-cycle is employed to successively increase accuracy. Refinement is controlled by a threshold criterion. In former experiments we observed oscillating index sets, that means indices which where included in the adaptive basis, because the corresponding wavelet coefficient was large, were thrown out in the next cycle, and after that again included. To avoid this behaviour we mark all indices whose coefficients fulfilled the threshold criterion at least for one time. These are kept in the index set for all times.

The source to this example is Sources/Examples/PO.cc. You can compile it by make PO in the Examples directory.

see the code PO.cc

Visualize the results in MATLAB by (AWFD/Sources/Examples/po.m)

U=ReadUDF('../../Data/Test/U') ;
figure(1) ;
[x,y]=meshgrid(U.x{1},U.x{2}) ;
q=[1:8:1025] ; surf(x(q,q),y(q,q),U.a(q,q)) ; shading flat;
figure(2) ; G=ReadSparse('../../Data/Test/U.adp', 0 ) ;
m=2^G.L ;
g=sparse(G.index(:,1)+1 , m-(G.index(:,2)-1) +1 , G.val) ;
spy(g',0.1) ;
set(gca,'XTick',[],'YTick',[], 'XTickLabel',[], 'YTickLabel',[]) ;
xlabel('') ;

You should get the following two figures. Figure 2 shows the distribution of points for the adaptive grid of the final solution.

1806 1807


next up previous
Next: Solution of the Navier-Stokes Up: Examples Previous: Calculation of Partial Derivatives
koster 2003-07-29