Download User Manual Adaptive Multi-Resolution Wavelet

User Manual
initial guess (m = 0): um
G≥m
k and
h
i
while m = 0 or m > 1 and G≥m 6= G≥m−1 or kf J − LuJ≥ k∞ > δǫ
m=m+1
perform forward wavelet transform for each component of um
k
for all levels j = J : −1 : 1
create a mask M for |dµ,j
l | ≥ ǫ
end
extend the mask M with adjacent wavelets
perform the reconstruction check procedure
construct G≥m+1
if G≥m+1 6= G≥m
m+1
interpolate um
k to G≥
end if
solve Lu = f using Local Multilevel Elliptic Solver.
end
Algorithm 1. Global Elliptic Solver.
initial conditions (n = 0): unk and G≥n
while tn < tmax
tn+1 = tn + ∆t
integrate the system of equations using Krylov time integration to obtain un+1
k
n+1
perform forward wavelet transform for each component of uk
for all levels j = J : −1 : 1
create a mask M for |dµ,j
l | ≥ ǫ
end
extend the mask M with adjacent wavelets
perform the reconstruction check procedure
construct G≥n+1
if G≥n+1 6= G≥n
interpolate un+1
to G≥n+1
k
end if
n=n+1
end
Algorithm 2. Time Evolution Solver.
4