<div dir="auto">Thanks!</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 2, 2020, 20:51 Zhang, Hong via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div style="word-wrap:break-word;line-break:after-white-space">
I checked the results again and can confirm that it is a bug on our side. A merge request has been created to fix it: 
<div><a href="https://gitlab.com/petsc/petsc/-/merge_requests/2830" target="_blank" rel="noreferrer">https://gitlab.com/petsc/petsc/-/merge_requests/2830</a></div>
<div><br>
</div>
<div>Thanks,</div>
<div>Hong (Mr.)<br>
<div><br>
<blockquote type="cite">
<div>On Jun 2, 2020, at 6:26 PM, Salazar De Troya, Miguel <<a href="mailto:salazardetro1@llnl.gov" target="_blank" rel="noreferrer">salazardetro1@llnl.gov</a>> wrote:</div>
<br>
<div>
<div style="font-family:Verdana;font-size:14px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<div style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
Hong,<u></u><u></u></div>
<div style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
<div style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
That is not the correct result, however, we can obtain the correct result if we comment ts.setThetaEndpoint(True). Why is it so?<u></u><u></u></div>
<div style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
<div style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
Thanks<u></u><u></u></div>
<div style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
Miguel<u></u><u></u></div>
<div style="margin:0in 0in 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
<div style="border-style:solid none none;border-top-width:1pt;border-top-color:rgb(181,196,223);padding:3pt 0in 0in">
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<b><span style="font-size:12pt">From:<span> </span></span></b><span style="font-size:12pt">"Zhang, Hong" <<a href="mailto:hongzhang@anl.gov" target="_blank" rel="noreferrer">hongzhang@anl.gov</a>><br>
<b>Date:<span> </span></b>Tuesday, June 2, 2020 at 3:31 PM<br>
<b>To:<span> </span></b>"Salazar De Troya, Miguel" <<a href="mailto:salazardetro1@llnl.gov" target="_blank" rel="noreferrer">salazardetro1@llnl.gov</a>><br>
<b>Cc:<span> </span></b>"<a href="mailto:petsc-users@mcs.anl.gov" target="_blank" rel="noreferrer">petsc-users@mcs.anl.gov</a>" <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank" rel="noreferrer">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:<span> </span></b>Re: [petsc-users] TSAdjoint not working correctly when using TSThetaSetEndpoint(true)<u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Miguel,<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
</div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
After I uncommented ts.setThetaEndpoint(True) and removed ts.setProblemType(False.ProblemType.LINEAR) in your code, I got the following result:<u></u><u></u></div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
hongzhang@Hongs-MacBook-Pro$ python3 petsc_question_2.py<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Cost function<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Vec Object: 1 MPI processes<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
  type: seq<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
127.781<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Exact value: 127.78112<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Numerical sensitivity w.r.t. a<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
12.778122049945212<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Real sensitivity w.r.t a<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
12.7781121978613<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Numerical sensitivity w.r.t. b<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
83.8907580310942<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Real sensitivity w.r.t b<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
211.67168296791954<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
-ts_type cn gives the same result. What was the problem you encountered?<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Thanks,<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Hong (Mr.)<u></u><u></u></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<br>
<br>
<u></u><u></u></div>
<blockquote style="margin-top:5pt;margin-bottom:5pt" type="cite">
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
On Jun 2, 2020, at 2:36 PM, Salazar De Troya, Miguel via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" style="color:blue;text-decoration:underline" target="_blank" rel="noreferrer">petsc-users@mcs.anl.gov</a>> wrote:<u></u><u></u></div>
</div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<u></u> <u></u></div>
<div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Hello,<span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
 <span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
I am attaching a simple example that uses TSAdjoint to calculate the sensitivity of a 1 degree of freedom ODE. When using theta methods, it returns the right sensitivity as given by the analytical solution. When I set `ts.setThetaEndpoint(True)`, it does not
 work. Is there a theoretical reason why? I was using this option because using the options<span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
 <span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
ts.setType(ts.Type.THETA)<span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
ts.setTheta(0.5) <span> </span><span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
ts.setThetaEndpoint(True)<span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
 <span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
is equivalent to using <a href="http://ts.Type.CN" target="_blank" rel="noreferrer">ts.Type.CN</a><span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
 <span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Thanks<span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
Miguel<span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
 <span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<span lang="ES" style="font-size:9pt;font-family:Consolas">Miguel A. Salazar de Troya</span><span style="font-size:12pt"><u></u><u></u></span></div>
</div>
<div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<span style="font-size:9pt;font-family:Consolas">Postdoctoral Researcher, Lawrence Livermore National Laboratory</span><span style="font-size:12pt"><u></u><u></u></span></div>
</div>
</div>
<div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<span style="font-size:9pt;font-family:Consolas">B141</span><span style="font-size:12pt"><u></u><u></u></span></div>
</div>
</div>
<div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<span style="font-size:9pt;font-family:Consolas">Rm: 1085-5</span><span style="font-size:12pt"><u></u><u></u></span></div>
</div>
</div>
<div>
<div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<span style="font-size:9pt;font-family:Consolas">Ph: 1(925) 422-6411</span><span style="font-size:12pt"><u></u><u></u></span></div>
</div>
</div>
<div style="margin:0in 0in 0.0001pt 0.5in;font-size:11pt;font-family:Calibri,sans-serif">
<petsc_question_2.py></div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>

</blockquote></div>