<div dir="ltr"><div><span style="font-family:arial,sans-serif;font-size:13px">Thanks Jed for the prompt reply.</span></div><div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><span style="font-family:arial,sans-serif;font-size:13px">> How do you distinguish "Lagrange" versus "Spectral" elements?  Are you </span><span style="font-family:arial,sans-serif;font-size:13px">referring to nodal versus modal bases, choice of collocation quadrature, </span><span style="font-family:arial,sans-serif;font-size:13px">or something else?</span><br>
<div><span style="font-family:arial,sans-serif">By Spectral Elements I mean spectral quadrilateral/hexahedral elements based on tensor product lagrangian polynomials on Gauss Lobatto Legendre points. </span><br></div><div>
<font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif"> ></font>You could reorder your equations, but multicolor GS is not a very good or representative algorithm on cache-based architectures, due to its poor cache reuse.  I suggest just using standard GS smoothers (-pc_type sor with default relaxation parameter of 1.0).</div>
<div>I plan to implement multicolor GS precisely to demonstrate its poor performance as compared to other iterative and MG schemes, because in the Phase Field community multicolor GS is still quite popular and lingers around as a solver. The main point of this work is to clearly demonstrate the ill-suitedness of GS for  these coupled transport problems. </div>
<div>So just wondering if there are any related examples showing multicolor GS as a solver. Also, since you mentioned, are there any references which demonstrate the poor cache reuse of multicolor GS or is it too obvious?... just curious.</div>
<div><br></div><div><div>>Do you already have spectral element code or are you planning to write one from scratch?  There are several unstructured FEM examples using DMPlex, though they do not explicitly expose the tensor product that is</div>
<div>important for efficient high-order methods.  How would you plan to precondition these things?</div></div><div><br></div><div>Sorry I forgot to mention..... I am only interested in structured quad/hex elements. I have my old implementations of higher order Lagrange elements and also used deal.ii's Spectral elements.... but for this work I will more or less write one from scratch. So any pointers to efficient tensor grid FEM implementation will really help me. </div>
<div><br></div><div>And thanks for your paper reference.</div><div><br></div><div>> That is just -ksp_type richardson.  Or are you talking about nonlinear multigrid (-snes_type fas)?</div><div>Cool..... -snes_type fas is what I was looking for. </div>
<div><br></div><div>Thanks,</div><div>Shiva</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Oct 17, 2013 at 10:50 AM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">Shiva Rudraraju <<a href="mailto:rudraa@umich.edu">rudraa@umich.edu</a>> writes:<br>
<br>
> I am planning to study the effects of geometric multigrid method both as a<br>
> Solver and Preconditioner on a system of transport problems involving<br>
> coupled diffusion + finite strain mechanics. Part of that study will be to<br>
> evaluate the performance of FD vs higher order FEM (Both Lagrange and<br>
> Spectral elements).<br>
<br>
</div>How do you distinguish "Lagrange" versus "Spectral" elements?  Are you<br>
referring to nodal versus modal bases, choice of collocation quadrature,<br>
or something else?<br>
<div class="im"><br>
> So in this context I have a few questions:<br>
> 1)  To motivate the use of MG schemes in this coupled system of equations,<br>
> I plan to initially implement a Red-Black Gauss Siedel solver for FD to<br>
> serve as a benchmark. Is there a good example related to the Red-Black<br>
> Gauss Siedel implementation?.<br>
<br>
</div>You could reorder your equations, but multicolor GS is not a very good<br>
or representative algorithm on cache-based architectures, due to its<br>
poor cache reuse.  I suggest just using standard GS smoothers (-pc_type<br>
sor with default relaxation parameter of 1.0).<br>
<div class="im"><br>
> 2) I plan to implement a higher order Lagrange and Spectral element FEM<br>
> code. I found one related example in<br>
> petsc/src/ksp/ksp/examples/tutorials/ex59.c. Are there any better examples<br>
> of general FEM code implementations.<br>
<br>
</div>Do you already have spectral element code or are you planning to write<br>
one from scratch?  There are several unstructured FEM examples using<br>
DMPlex, though they do not explicitly expose the tensor product that is<br>
important for efficient high-order methods.  How would you plan to<br>
precondition these things?<br>
<br>
I did some work on this using embedded low-order discretizations (see<br>
below) and the library used in that paper is available.  The paper used<br>
algebraic multigrid, but geometric and semi-geometric methods would be<br>
better in many settings.  (I'm hoping to merge the good attributes of<br>
that library into PETSc using the new DMPlex.)<br>
<div class="im"><br>
> 3) MG is often used as a preconditioner, but are there any examples showing<br>
> its use as a solver?.<br>
<br>
</div>That is just -ksp_type richardson.  Or are you talking about nonlinear<br>
multigrid (-snes_type fas)?<br>
<div class="im"><br>
> I found a few examples here<br>
> <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/concepts/multigrid.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/concepts/multigrid.html</a>,<br>
> but nothing in the current version of those codes refer to any MG<br>
> datastructures like DMMG which were part of the earlier versions.<br>
<br>
</div>The DMMG object was a broken abstraction; its functionality was rolled into DM and SNES.<br>
<br>
<br>
<br>
@article {brown2010ens,<br>
  author = {Brown, Jed},<br>
  affiliation = {Versuchsanstalt für Wasserbau, Hydrologie und Glaziologie (VAW), ETH Zürich, Zürich, Switzerland},<br>
  title = {Efficient Nonlinear Solvers for Nodal High-Order Finite Elements in {3D}},<br>
  journal = {Journal of Scientific Computing},<br>
  publisher = {Springer Netherlands},<br>
  issn = {0885-7474},<br>
  keyword = {Mathematics and Statistics},<br>
  pages = {48-63},<br>
  volume = {45},<br>
  issue = {1},<br>
  url = {<a href="http://dx.doi.org/10.1007/s10915-010-9396-8" target="_blank">http://dx.doi.org/10.1007/s10915-010-9396-8</a>},<br>
  doi = {10.1007/s10915-010-9396-8},<br>
  year = {2010}<br>
}<br>
</blockquote></div><br></div>