<div class="gmail_quote">On Fri, Jan 20, 2012 at 08:19, Alexander Grayver <span dir="ltr"><<a href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#ffffff" text="#000000">My initial motivation to try gamg was that paper:<br>
    <span>W.A. <em>Mulder</em>, A <em>multigrid solver for
        3D</em> electromagnetic diffusion<br></span></div></blockquote><div><br></div><div>This paper does entirely geometric coarse grids. You could use PCMG, but you would have to write your own interpolation and restriction operators (because even in the structured case, DMDA does not have specific support for edge discretizations; you can use DMDA, but it won't automatically build correct interpolants).</div>
<div><br></div><div>The smoother in this paper is node-centered overlapping relaxation. It solves overlapping 6x6 problems (all edges connected to that node). I don't know of a reliable way to extract these node connections from the matrix alone, so you would have to implement this.</div>
<div><br></div><div>If you implement these two ingredients (interpolation and node-based relaxation), then you can run the geometric multigrid method in this paper.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#ffffff" text="#000000"><span>
      <br>
      Where it seems to work fine (on uniform grids only at least).<br>
    </span>I use MUMPS currently and it is very robust. It properly
    solves systems with extremely low \omega and \sigma for any
    stretched (which would 'kill' any multigrid I guess) grids.<br>
    I don't expect iterative solvers to be that robust, but I would like
    to solve even simple models with ~uniform grid of the order 10^7<br>
    <br>
    So, I will try to look at HYPRE's approach for Maxwell equations. </div></blockquote></div><br>