<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<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-family: "Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, system-ui, Roboto, "Helvetica Neue", sans-serif; 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 style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" 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>
</body>
</html>