<html>
<head>
<!-- Template generated by Exclaimer Mail Disclaimers on 05:00:52 Friday, 11 September 2020 -->
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} p
{margin-top:0;
margin-bottom:0}--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p></p>
<p>Sure, that was the advise 9 years ago in the ancient thread. It's not a big problem.</p>
<p><br>
</p>
<p>Chris<br>
</p>
<div dir="ltr" style="color: rgb(33, 33, 33);"></div>
<br>
<font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;">dr. ir.</font><font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;">
</font><font style="font-family: Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;">Christiaan</font><font style="font-family: Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;">
</font><font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;">Klaij</font><font size="3" face="Segoe UI"></font><font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;"> |
</font><font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;">Senior Researcher</font><font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;"> |
</font><font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;font-weight:bold;">Research & Development</font><br>
<font style="font-family:Verdana;font-size:9pt;color:#7D7D7D;font-weight:bold;">MARIN</font><font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;"> | T
</font><font style="font-family:Verdana;font-size:8pt;color:#7D7D7D;">+31 317 49 33 44</font><font color="#7d7d7d" size="2" face="Segoe UI"></font><font style="font-family:Verdana;font-size:8pt;"> |
</font><span style="font-family:Verdana;font-size:8pt;"><a href="mailto:C.Klaij@marin.nl" title="" target="" style="font-family:Verdana;font-size:8pt;"><span style="font-family:Verdana; font-size:8pt;">C.Klaij@marin.nl</span></a></span><font color="#7d7d7d" size="2" face="Segoe UI"></font><font style="font-family:Verdana;font-size:8pt;"> |
</font><span style="font-family:Verdana;font-size:8pt;"><a href="http://www.marin.nl" title="" target="" style="font-family:Verdana;font-size:8pt;"><span style="font-family:Verdana; font-size:8pt;">www.marin.nl</span></a></span><br>
<br>
<a href="https://www.linkedin.com/company/marin" target=""><img width="16" height="16" style="border: 0px Solid ; " src="cid:imagefd46c4.PNG@806d4b51.459cf191" alt="LinkedIn"></a> <a href="http://www.youtube.com/marinmultimedia" target=""><img width="16" height="16" style="border: 0px Solid ; " src="cid:image97655b.PNG@14a86986.4b9c8c49" alt="YouTube"></a> <a href="https://twitter.com/MARIN_nieuws" target=""><img width="16" height="16" style="border: 0px Solid ; " src="cid:imagedce627.PNG@133da3ea.4c916ca8" alt="Twitter"></a> <a href="https://www.facebook.com/marin.wageningen" target=""><img width="16" height="16" style="border: 0px Solid ; " src="cid:imagef2324a.PNG@58eef720.49bd5877" alt="Facebook"></a><br>
<font style="FONT-FAMILY: Verdana; FONT-SIZE: 8pt" color="#7d7d7d" size="2" face="Verdana">MARIN news:
<a href="https://www.marin.nl/news/new-publication-on-noise-measurements-on-a-cavitating-propeller" title="" target="" style="">
New publication on noise measurements of a cavitating propeller</a></font>
<p></p>
<p id="undefined"></p>
<p></p>
<div dir="ltr" style="color: rgb(33, 33, 33);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Zhang, Hong <hzhang@mcs.anl.gov><br>
<b>Sent:</b> Friday, September 11, 2020 4:57 PM<br>
<b>To:</b> Klaij, Christiaan; petsc-users@mcs.anl.gov<br>
<b>Subject:</b> Re: [petsc-users] How to activate the modified Gram-Schmidt orthogonalization process in Fortran?</font>
<div> </div>
</div>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Sorry, we have not done it. Can you use
<pre>PetscOptionsSetValue("-ksp_gmres_modifiedgramschmidt", "1")</pre>
<pre>for now? We'll try to add <span style="color:rgb(32,31,30); font-size:14.6667px; background-color:rgb(255,255,255); display:inline!important">the fortran binding later.</span></pre>
<pre>Hong</pre>
</div>
<div id="appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> petsc-users <petsc-users-bounces@mcs.anl.gov> on behalf of Klaij, Christiaan <C.Klaij@marin.nl><br>
<b>Sent:</b> Friday, September 11, 2020 7:50 AM<br>
<b>To:</b> petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Subject:</b> Re: [petsc-users] How to activate the modified Gram-Schmidt orthogonalization process in Fortran?</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="PlainText"><br>
Make me feel ancient. Would be nice to have the fortran binding though...<br>
<br>
Chris<br>
<br>
> ----------------------------------------------------------------------<br>
><br>
> Message: 1<br>
> Date: Thu, 10 Sep 2020 19:41:30 -0600<br>
> From: Zhuo Chen <chenzhuotj@gmail.com><br>
> To: "Zhang, Hong" <hzhang@mcs.anl.gov><br>
> Cc: "petsc-users@mcs.anl.gov" <petsc-users@mcs.anl.gov><br>
> Subject: Re: [petsc-users] How to activate the modified Gram-Schmidt<br>
> orthogonalization process in Fortran?<br>
> Message-ID:<br>
> <CAHjeDn8HmgxttujK6YP0qfgPcmQMbVNa0UTFN4Wd_Au-0UjDAg@mail.gmail.com><br>
> Content-Type: text/plain; charset="utf-8"<br>
><br>
> Hi Hong,<br>
><br>
> According to that very old thread, KSPGMRESSetOrthogonalization was not<br>
> implemented in Fortran. I did as you suggested and the compiler will tell<br>
> me<br>
><br>
> undefined reference to `kspgmressetorthogonalization_'<br>
><br>
> I think I will use the -ksp_gmres_modifiedgramschmidt method. Thank you so<br>
> much!<br>
><br>
<br>
dr. ir. Christiaan Klaij | Senior Researcher | Research & Development<br>
MARIN | T +31 317 49 33 44 | <a href="mailto:C.Klaij@marin.nl">mailto:C.Klaij@marin.nl</a> |
<a href="http://www.marin.nl">http://www.marin.nl</a><br>
<br>
<br>
> On Thu, Sep 10, 2020 at 7:32 PM Zhang, Hong <hzhang@mcs.anl.gov> wrote:<br>
><br>
> > Zhuo,<br>
> > Call<br>
> > KSPSetType(ksp,KSPGMRES);<br>
> ><br>
> > KSPGMRESSetOrthogonalization(ksp,KSPGMRESModifiedGramSchmidtOrthogonalization);<br>
> > Hong<br>
> ><br>
> > ------------------------------<br>
> > *From:* Zhuo Chen <chenzhuotj@gmail.com><br>
> > *Sent:* Thursday, September 10, 2020 8:17 PM<br>
> > *To:* Zhang, Hong <hzhang@mcs.anl.gov><br>
> > *Cc:* petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
> > *Subject:* Re: [petsc-users] How to activate the modified Gram-Schmidt<br>
> > orthogonalization process in Fortran?<br>
> ><br>
> > Hi Hong,<br>
> ><br>
> > Thank you very much for your help.<br>
> ><br>
> > It seems that if I simply append -ksp_gmres_modifiedgramschmidt the<br>
> > warning goes away. However<br>
> > KSPGMRESSetOrthogonalization(ksp,KSPGMRESModifiedGramSchmidtOrthogonalization,ierr)<br>
> > has another issue.<br>
> ><br>
> > Error: Symbol ?kspgmresmodifiedgramschmidtorthogonalization? at (1) has no<br>
> > IMPLICIT type<br>
> ><br>
> > Is it because the argument is too long? I am using gcc 8.4.0 instead of<br>
> > ifort<br>
> ><br>
> > On Thu, Sep 10, 2020 at 7:08 PM Zhang, Hong <hzhang@mcs.anl.gov> wrote:<br>
> ><br>
> > Zhuo,<br>
> > Run your code with option '-ksp_gmres_modifiedgramschmidt'. For example,<br>
> > petsc/src/ksp/ksp/tutorials<br>
> > mpiexec -n 2 ./ex2 -ksp_view -ksp_gmres_modifiedgramschmidt<br>
> > KSP Object: 2 MPI processes<br>
> > type: gmres<br>
> > restart=30, using Modified Gram-Schmidt Orthogonalization<br>
> > happy breakdown tolerance 1e-30<br>
> > maximum iterations=10000, initial guess is zero<br>
> > tolerances: relative=0.000138889, absolute=1e-50, divergence=10000.<br>
> > left preconditioning<br>
> > using PRECONDITIONED norm type for convergence test<br>
> > PC Object: 2 MPI processes<br>
> > type: bjacobi<br>
> > ...<br>
> ><br>
> > You can<br>
> > call KSPGMRESSetOrthogonalization(ksp,KSPGMRESModifiedGramSchmidtOrthogonalization)<br>
> > in your program.<br>
> ><br>
> > Hong<br>
> ><br>
> > ------------------------------<br>
> > *From:* petsc-users <petsc-users-bounces@mcs.anl.gov> on behalf of Zhuo<br>
> > Chen <chenzhuotj@gmail.com><br>
> > *Sent:* Thursday, September 10, 2020 7:52 PM<br>
> > *To:* petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
> > *Subject:* [petsc-users] How to activate the modified Gram-Schmidt<br>
> > orthogonalization process in Fortran?<br>
> ><br>
> > Dear Petsc users,<br>
> ><br>
> > I found an ancient thread discussing this problem.<br>
> ><br>
> > <a href="https://lists.mcs.anl.gov/pipermail/petsc-users/2011-October/010607.html">
https://lists.mcs.anl.gov/pipermail/petsc-users/2011-October/010607.html</a><br>
> ><br>
> > However, when I add<br>
> ><br>
> > call KSPSetType(ksp,KSPGMRES,ierr);CHKERRQ(ierr)<br>
> > call<br>
> > PetscOptionsSetValue(PETSC_NULL_OPTIONS,'-ksp_gmres_modifiedgramschmidt','1',ierr);CHKERRQ(ierr)<br>
> ><br>
> > the program will tell me<br>
> ><br>
> > WARNING! There are options you set that were not used!<br>
> > WARNING! could be spelling mistake, etc!<br>
> > There is one unused database option. It is:<br>
> > Option left: name:-ksp_gmres_modifiedgramschmidt value: 1<br>
> ><br>
> > I would like to know the most correct way to activate the modified<br>
> > Gram-Schmidt orthogonalization process in Fortran. Thank you very much!<br>
> ><br>
> > Best regards.<br>
> ><br>
> ><br>
</div>
</span></font></div>
<br>
<br>
<font style="background-color:#ffffff">
<p align="center"><font style="background-color:#ffffff"><b>Help us improve the spam filter. If this message contains SPAM, click
<a href="https://www.mailcontrol.com/sr/k4QvsNfY0JbGX2PQPOmvUmuFhuAgpUQ2aIYGFgBh1XE-_-bwIG1ztLM6HR4_1LKTfIrk20xxkqZfcZxYauAECA==">
here</a> to report. Thank you, MARIN Support Group</b></font></p>
</font></div>
</div>
</EXCLAIMER:ORIGINAL_BODY_AFTER>
<p></p>
</EXCLAIMER:TEXT_ONLY></EXCLAIMER:AD></EXCLAIMER:AD>
</body>
</html>