<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, May 11, 2021 at 4:03 PM Frederico Teixeira <teixeira@zmt.swiss> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Hi Mark,<br><br>I think so. The complex-valued linear system Ax = b, with A=T+Wi, u=x+yi, b=p+qi (i is the unit complex) can be transformed into a 2x2 real-valued system:<br>[ T  -W ] (x) = (p)<br>[ W   T ] (y)    (q)<br></div></div></blockquote><div><br></div><div>Yes, this is ERF. </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>Perhaps a bit more background can clarify things. System A stems from a Poisson equation with complex coefficients. There are ~130 piecewise constant coefficients ranging across 5 orders of magnitude. </div></div></blockquote><div><br></div><div>I assume these are big jumps, are they?</div><div><br></div><div>You might be able to use GAMG or Hypre and tune the threshold to "see" these jumps if they are big enough.  In general large random coefficient jumps are hard with generic AMG. You have 130 patches so a coarse grid could in theory resolve this structure with a not too large "grid"</div><div><br></div><div>This may not be clear, but I would start with a real valued version of this and see if you can tune Hypre or GAMG to work OK. This could be done with, as I said, tuning the threshold parameter so that coarsening "sees" your jumps and respects them.</div><div><br></div><div>Next, use a coarse grid solve that is big enough to capture the geometry of your jumps, approximately, where an exact solver is used and it does not care about jumps.</div><div><br></div><div>As far as ERF, if you can translate it into a (complex) PetscScalar matrix somehow then you can use GAMG. I don't know about Hypre.</div><div><br></div><div>Or, I have used ERF with some success. A long time ago. You want to interleave the real and complex parts. How big is |W| / |T| ? If this is small it should work OK, but it is hard to be sure.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Most of the KSPs and PCs I managed to test take between 6000-9000 iterations to reach 10^-9. I hope to get some speed up with preconditioners that exploit the structure above.<br><br>Regards,<br> Frederico.<br><br><div id="gmail-m_-6225517201256615678e1bde0fd-f1a5-4833-935b-0cc70f26a0a4"><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div></div><div style="text-align:left"><span style="font-size:x-small"></span><span style="font-size:x-small"><span style="font-size:x-small"></span></span><div style="text-align:left"><span style="font-size:11pt;font-family:arial,helvetica,sans-serif;color:rgb(0,0,0);background-color:rgb(255,255,255)">Dr. Frederico Teixeira</span></div><div style="text-align:left"><span style="font-size:11pt;font-family:arial,helvetica,sans-serif;background-color:rgb(255,255,255);color:rgb(0,0,0)">Computational Modeler and Software Developer, <a href="http://www.itis.swiss" style="background-color:rgb(255,255,255);color:rgb(0,0,0)" rel="nofollow%20noopener%20noreferrer%20nofollow%20noopener%20noreferrer%20nofollow%20noopener%20noreferrer%20nofollow%20noopener%20noreferrer nofollow noopener noreferrer nofollow noopener noreferrer" target="_blank">ZMT</a> (member of <a href="https://www.z43.swiss" style="background-color:rgb(255,255,255);color:rgb(0,0,0)" rel="nofollow%20noopener%20noreferrer%20nofollow%20noopener%20noreferrer%20nofollow%20noopener%20noreferrer%20nofollow%20noopener%20noreferrer nofollow noopener noreferrer nofollow noopener noreferrer" target="_blank">Zurich43</a>)</span></div><div style="text-align:left"><span style="font-size:11pt;font-family:arial,helvetica,sans-serif;background-color:rgb(255,255,255);color:rgb(0,0,0)"><br></span></div><div style="text-align:left">







<p><span style="font-size:11pt;font-family:arial,helvetica,sans-serif;background-color:rgb(255,255,255);color:rgb(0,0,0)">P +41 44 245 9698<br>
Zeughausstrasse 43, 8004 Zurich, Switzerland</span></p></div></div></div><br><hr id="gmail-m_-6225517201256615678zwchr"><div><b>From: </b>"Mark Adams" <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>><br><b>To: </b>"Frederico Teixeira" <teixeira@zmt.swiss><br><b>Cc: </b>"petsc-users" <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br><b>Sent: </b>Tuesday, May 11, 2021 5:16:35 PM<br><b>Subject: </b>Re: [petsc-users] Binary format in real vs. complex scalar type configurations<br></div><br><div><div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, May 11, 2021 at 9:30 AM Frederico Teixeira <teixeira@zmt.swiss> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><div><div><div><div>Dear fellows,<br></div></div></div></div><div><div style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255)"><br><div>I hope this message finds you safe and well.</div><br><div>I have<span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);float:none;display:inline"> </span><span style="font-size:12pt">a complex-valued matrix and its real/imaginary components <span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);float:none;display:inline">in </span><span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);float:none;display:inline">binary format</span></span><span style="font-size:12pt">. They were extracted from a solver that only works with "scalar-type=complex" configuration.</span></div><div>I am getting weird results when I load them into a small test program that's configured with "scalar-type=real", but I believe this is expected.</div><div>At the end of the day, I would like to have both real and imaginary components as real-valued matrices.</div><div>Is it possible to do it? I want to<span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);float:none;display:inline"> test preconditioners that are tailored for this sort of problem.</span></div></div></div></div></div></blockquote><br><div>Do you mean you want what is called "equivalent real form" where the real and complex parts are stored as type 'double' for example and operations like multiply take two pairs of doubles, do a complex multiply manually, and return a real/complex pair of doubles?</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><div><div style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-size:16px;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255)"><br><div>Regards,</div><div> Frederico.</div></div><span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif;font-style:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:11pt;background-color:rgb(255,255,255)"></span></div></div></div></blockquote></div></div><br></div></div></div></blockquote></div></div>