Basic HTML version of Foils prepared March 15 00

Foil 22 Multigrid Algorithm: procedure MG(level, A, u, f)

From Collection of Extra Foils for CPS615 PDE Iterative Solution Discussion CPS615 Spring Semester 00 -- March 00. by Geoffrey C. Fox
if level = coarsest then
  • solve coarsest grid equation by another method "exactly"
else
  • smooth Alevel u = f (m1 times)
  • Compute residual r = f - Alevel u
  • Restrict F = R r ( R is Restriction Operator)
  • Call MG( level + 1, A(level+1), V, F) (mc times)
  • Interpolate v = P V (Interpolate new solution at this level)
  • correct unew = u + v
  • smooth Aunew = f (m2 times) and
  • set u = unew
endif
endprocedure
Alevel uexact = f
uexact = u + v
Alevel v = r = f - Alevel u

© Northeast Parallel Architectures Center, Syracuse University, npac@npac.syr.edu

If you have any comments about this server, send e-mail to webmaster@npac.syr.edu.

Page produced by wwwfoil on Mon Mar 20 2000