<div dir="ltr">Jed creates the LocalToGlobal that does this elimination in plexsfc.c<div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jul 18, 2023 at 12:07 PM Barry Smith <<a href="mailto:bsmith@petsc.dev">bsmith@petsc.dev</a>> 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><br></div>   They are never really "eliminated" because extra copies in the global vector never exist. <br><div><br><blockquote type="cite"><div>On Jul 18, 2023, at 12:03 PM, Karthikeyan Chockalingam - STFC UKRI <<a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" target="_blank">karthikeyan.chockalingam@stfc.ac.uk</a>> wrote:</div><br><div><div style="font-family:Helvetica;font-size:18px;font-style:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Thank you, Barry. I am using the MPIAIJ format for a Finite Element application. So, I am trying to understand what is implemented in DMDA to eliminate those extra nodes.<u></u><u></u></span></div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"><u></u> <u></u></span></div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Best,<u></u><u></u></span></div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Karthik.<u></u><u></u></span></div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"><u></u> <u></u></span></div><div id="m_5778741566562819247mail-editor-reference-message-container"><div><div style="border-width:1pt medium medium;border-style:solid none none;border-color:rgb(181,196,223) currentcolor currentcolor;padding:3pt 0cm 0cm"><p class="MsoNormal" style="margin:0cm 0cm 12pt;font-size:10pt;font-family:Calibri,sans-serif"><b><span style="font-size:12pt">From:<span> </span></span></b><span style="font-size:12pt">Barry Smith <<a href="mailto:bsmith@petsc.dev" style="color:blue;text-decoration:underline" target="_blank">bsmith@petsc.dev</a>><br><b>Date:<span> </span></b>Tuesday, 18 July 2023 at 16:58<br><b>To:<span> </span></b>Chockalingam, Karthikeyan (STFC,DL,HC) <<a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" style="color:blue;text-decoration:underline" target="_blank">karthikeyan.chockalingam@stfc.ac.uk</a>><br><b>Cc:<span> </span></b>Matthew Knepley <<a href="mailto:knepley@gmail.com" style="color:blue;text-decoration:underline" target="_blank">knepley@gmail.com</a>>,<span> </span><a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank">petsc-users@mcs.anl.gov</a>><br><b>Subject:<span> </span></b>Re: [petsc-users] periodic boundary conditions<u></u><u></u></span></p></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"><u></u> <u></u></span></div></div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">   If you are using DMDA with periodic boundary conditions for example only one "copy" of such nodes exists in the global vector (the vector the solvers see) so one does not need to eliminate extra ones <u></u><u></u></span></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"><br><br><u></u><u></u></span></div><blockquote style="margin-top:5pt;margin-bottom:5pt"><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">On Jul 18, 2023, at 11:51 AM, Karthikeyan Chockalingam - STFC UKRI via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<u></u><u></u></span></div></div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"><u></u> <u></u></span></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Yes, I clearly understand I need to eliminate one set of periodic nodes. I was hoping to use<span> </span>x = P x’ to eliminate one set. It is a kind of mapping.</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Sorry, I am not sure if it is the<span> </span>LocalToGlobal mapping you are referring to. Is there an example or reference to show how<span> </span>the<span> </span>LocalToGlobal mapping is being used to impose PBC?</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Best,</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Karthik.</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div id="m_5778741566562819247mail-editor-reference-message-container"><div><div style="border-width:1pt medium medium;border-style:solid none none;padding:3pt 0cm 0cm;border-color:currentcolor"><p class="MsoNormal" style="margin:0cm 0cm 12pt;font-size:10pt;font-family:Calibri,sans-serif"><b><span style="font-size:12pt">From:<span> </span></span></b><span style="font-size:12pt">Matthew Knepley <<a href="mailto:knepley@gmail.com" style="color:blue;text-decoration:underline" target="_blank">knepley@gmail.com</a>><br><b>Date:<span> </span></b>Tuesday, 18 July 2023 at 16:38<br><b>To:<span> </span></b>Chockalingam, Karthikeyan (STFC,DL,HC) <<a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" style="color:blue;text-decoration:underline" target="_blank">karthikeyan.chockalingam@stfc.ac.uk</a>><br><b>Cc:<span> </span></b><a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank">petsc-users@mcs.anl.gov</a>><br><b>Subject:<span> </span></b>Re: [petsc-users] periodic boundary conditions</span><u></u><u></u></p></div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">On Tue, Jul 18, 2023 at 11:18 AM Karthikeyan Chockalingam - STFC UKRI <</span><a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:11pt">karthikeyan.chockalingam@stfc.ac.uk</span></a><span style="font-size:11pt">> wrote:</span><u></u><u></u></div></div></div><div><blockquote style="border-width:medium medium medium 1pt;border-style:none none none solid;padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt;border-color:currentcolor currentcolor currentcolor rgb(204,204,204)"><div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Thanks Matt.</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">The mesh is structured (rectilinear), so it is periodic in that sense.</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Can you please explain how I can impose it strongly?</span><u></u><u></u></div></div></div></div></div></blockquote><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Strongly means make those variables equal in a pointwise sense. We do this in the</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">LocalToGlobal mapping, so one set is eliminated in the global problem.</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">  Thanks,</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">    Matt</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><blockquote style="border-width:medium medium medium 1pt;border-style:none none none solid;padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt;border-color:currentcolor currentcolor currentcolor rgb(204,204,204)"><div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">My initial thought was to come up with a relation between the periodic nodes:</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">x = P x’</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Say for 1-D problem with two elements</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">(1)-------------(2)------------(3)</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">P = [1 0, 0 1, 1 0]</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">x = [x1 x2 x3]</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">x’ = [x1 x2]</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">and solve</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">[P^T A P] x’ = P^T b</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">I don’t think [P^T A P] is deterministic.</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Kind regards,</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Karthik.</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div id="m_5778741566562819247m_-6387936294348734192mail-editor-reference-message-container"><div><div style="border-width:1pt medium medium;border-style:solid none none;padding:3pt 0cm 0cm;border-color:currentcolor"><p class="MsoNormal" style="margin:0cm 0cm 12pt;font-size:10pt;font-family:Calibri,sans-serif"><b><span style="font-size:12pt">From:<span> </span></span></b><span style="font-size:12pt">Matthew Knepley <</span><a href="mailto:knepley@gmail.com" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:12pt">knepley@gmail.com</span></a><span style="font-size:12pt">><br><b>Date:<span> </span></b>Tuesday, 18 July 2023 at 14:31<br><b>To:<span> </span></b>Chockalingam, Karthikeyan (STFC,DL,HC) <</span><a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:12pt">karthikeyan.chockalingam@stfc.ac.uk</span></a><span style="font-size:12pt">><br><b>Cc:<span> </span></b></span><a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:12pt">petsc-users@mcs.anl.gov</span></a><span><span style="font-size:12pt"> </span></span><span style="font-size:12pt"><</span><a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:12pt">petsc-users@mcs.anl.gov</span></a><span style="font-size:12pt">><br><b>Subject:<span> </span></b>Re: [petsc-users] periodic boundary conditions</span><u></u><u></u></p></div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">On Tue, Jul 18, 2023 at 9:02 AM Karthikeyan Chockalingam - STFC UKRI via petsc-users <</span><a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:11pt">petsc-users@mcs.anl.gov</span></a><span style="font-size:11pt">> wrote:</span><u></u><u></u></div></div></div><div><blockquote style="border-width:medium medium medium 1pt;border-style:none none none solid;padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt;border-color:currentcolor currentcolor currentcolor rgb(204,204,204)"><div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Hello,</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">This is exactly not a PETSc question. I am solving a Poisson equation using finite elements. I would like to impose PBC. I am thinking of using the Lagrange multiplier method to impose them as constraints. Or do you think I could take an alternative approach?</span><u></u><u></u></div></div></div></div></div></blockquote><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">There are several options:</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">1) Just make a periodic mesh. This is what Plex does by default.</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">2) Impose the conditions strongly. This is what is done if you create the ZBox shape in Plex.</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">3) Impose the conditions weakly. This is what you are doing with Lagrange multipliers. You could</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">also do a Nitsche boundary condition for this.</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Since the constraint is so simple, I do not see an advantage to imposing it weakly.</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">  Thanks,</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">     Matt</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><blockquote style="border-width:medium medium medium 1pt;border-style:none none none solid;padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt;border-color:currentcolor currentcolor currentcolor rgb(204,204,204)"><div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Thank you for your help.</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Kind regards,</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">Karthik.</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">-- </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><b><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:rgb(45,45,98)">Dr. Karthik Chockalingam</span></b><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:rgb(65,65,65)">High Performance Software Engineering Group</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:rgb(65,65,65)">Hartree Centre | Science and Technology Facilities Council</span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:rgb(5,99,193)">karthikeyan.chockalingam@stfc.ac.uk</span></a><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:12pt;font-family:"Times New Roman",serif"> </span><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:rgb(65,65,65)"> </span><u></u><u></u></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><b><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:rgb(65,65,65)"> </span></b><span style="font-size:12pt;font-family:"Times New Roman",serif;color:rgb(31,73,125)"><image001.png></span><u></u><u></u></div></div></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div></div></div></blockquote></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"><br clear="all"></span><u></u><u></u></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span><span style="font-size:11pt">--</span></span><u></u><u></u></div></div><div><div><div><div><div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><a href="http://www.cse.buffalo.edu/~knepley/" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:11pt">https://www.cse.buffalo.edu/~knepley/</span></a><u></u><u></u></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></blockquote></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"><br clear="all"></span><u></u><u></u></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span><span style="font-size:11pt">--</span></span><u></u><u></u></div></div><div><div><div><div><div><div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><span style="font-size:11pt"> </span><u></u><u></u></div></div></div><div><div><div style="margin:0cm;font-size:10pt;font-family:Calibri,sans-serif"><a href="http://www.cse.buffalo.edu/~knepley/" style="color:blue;text-decoration:underline" target="_blank"><span style="font-size:11pt">https://www.cse.buffalo.edu/~knepley/</span></a></div></div></div></div></div></div></div></div></div></div></div></div></div></blockquote></div></div></div></div></div></blockquote></div><br></div></blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div>