<!DOCTYPE html>

<!-- BaNnErBlUrFlE-BoDy-start -->

<!-- Preheader Text : BEGIN -->

<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">

 You can use the PETSC_OPTIONS environment variable to specify options if you don't pass the command line arguments through. You can set -tao_grtol smaller to handle this difference in scales between the objective and the gradient, though applying

</div>

<!-- Preheader Text : END -->

<!-- Email Banner : BEGIN -->

<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerStart</div>

<!--[if ((ie)|(mso))]>

  <table border="0" cellspacing="0" cellpadding="0" width="100%" style="padding: 16px 0px 16px 0px; direction: ltr" ><tr><td>

    <table border="0" cellspacing="0" cellpadding="0" style="padding: 0px 10px 5px 6px; width: 100%; border-radius:4px; border-top:4px solid #90a4ae;background-color:#D0D8DC;"><tr><td valign="top">

      <table align="left" border="0" cellspacing="0" cellpadding="0" style="padding: 4px 8px 4px 8px">

        <tr><td style="color:#000000; font-family: 'Arial', sans-serif; font-weight:bold; font-size:14px; direction: ltr">

          This Message Is From an External Sender

        </td></tr>

        <tr><td style="color:#000000; font-weight:normal; font-family: 'Arial', sans-serif; font-size:12px; direction: ltr">

          This message came from outside your organization.

        </td></tr>

      </table>

    </td></tr></table>

  </td></tr></table>

<![endif]-->

<![if !((ie)|(mso))]>

  <div dir="ltr"  id="pfptBannerdyfqif0" style="all: revert !important; display:block !important; text-align: left !important; margin:16px 0px 16px 0px !important; padding:8px 16px 8px 16px !important; border-radius: 4px !important; min-width: 200px !important; background-color: #D0D8DC !important; background-color: #D0D8DC; border-top: 4px solid #90a4ae !important; border-top: 4px solid #90a4ae;">

    <div id="pfptBannerdyfqif0" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">

      <div id="pfptBannerdyfqif0" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-weight:bold !important; font-weight:bold; font-size:14px !important; line-height:18px !important; line-height:18px">

        This Message Is From an External Sender

      </div>

      <div id="pfptBannerdyfqif0" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-weight:normal; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-size:12px !important; line-height:18px !important; line-height:18px; margin-top:2px !important;">

This message came from outside your organization.

      </div>

    </div>

    <div style="clear: both !important; display: block !important; visibility: hidden !important; line-height: 0 !important; font-size: 0.01px !important; height: 0px"> </div>

  </div>

<![endif]>

<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerEnd</div>

<!-- Email Banner : END -->

<!-- BaNnErBlUrFlE-BoDy-end -->

<html>

<head><!-- BaNnErBlUrFlE-HeAdEr-start -->

<style>

  #pfptBannerdyfqif0 { all: revert !important; display: block !important; 

    visibility: visible !important; opacity: 1 !important; 

    background-color: #D0D8DC !important; 

    max-width: none !important; max-height: none !important }

  .pfptPrimaryButtondyfqif0:hover, .pfptPrimaryButtondyfqif0:focus {

    background-color: #b4c1c7 !important; }

  .pfptPrimaryButtondyfqif0:active {

    background-color: #90a4ae !important; }

</style>

<!-- BaNnErBlUrFlE-HeAdEr-end -->

You can set -tao_grtol smaller to handle this difference in scales between the objective and the gradient, though applying some nondimensionalization/choice of appropriate units is still recommended.

"Palmer, Bruce J via petsc-users" <petsc-users@mcs.anl.gov> writes:

> This is a fortran code that doesn’t make use of argc,argv (I tried running with the runtime options anyway, in case you implemented some magic I’m not familiar with, but didn’t see anything new in the output). I have a call to TaoView(tao, PETSC_VIEWER_STDOUT_SELF,ierr) in the code and it reports back

>

>

>

> Tao Object: 1 MPI process

>

>   type: cg

>

>     CG Type: prp

>

>     Gradient steps: 0

>

>     Reset steps: 0

>

>   TaoLineSearch Object: 1 MPI process

>

>     type: more-thuente

>

>     maximum function evaluations=30

>

>     tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9

>

>     total number of function evaluations=0

>

>     total number of gradient evaluations=0

>

>     total number of function/gradient evaluations=0

>

>     Termination reason: 0

>

>   convergence tolerances: gatol=1e-08,   steptol=0.,   gttol=0.

>

>   Residual in Function/Gradient:=7.54237e+75

>

>   Objective value=2.96082e+86

>

>   total number of iterations=0,                          (max: 100)

>

>   total number of function/gradient evaluations=1,      (max: 4000)

>

>   Solution converged:    ||g(X)||/|f(X)| <= grtol

>

>

>

> Bruce

>

> From: Barry Smith <bsmith@petsc.dev>

> Date: Wednesday, June 26, 2024 at 2:02 PM

> To: Palmer, Bruce J <Bruce.Palmer@pnnl.gov>

> Cc: petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov>

> Subject: Re: [petsc-users] Unconstrained optimization question

> Check twice before you click! This email originated from outside PNNL.

>

>

>   Please run with -tao_monitor -tao_converged_reason and see why it has stopped.

>

>   Barry

>

>

>

> On Jun 26, 2024, at 4:34 PM, Palmer, Bruce J via petsc-users <petsc-users@mcs.anl.gov> wrote:

>

> This Message Is From an External Sender

> This message came from outside your organization.

> Hi,

>

> I’m trying to do an unconstrained optimization on a molecular scale problem. Previously, I was looking at an artificial molecular problem where all parameters were of order 1 and so the objective function and variables were also in the range of 1 or at least within a few orders of magnitude of 1.

>

> More recently, I’ve been trying to apply this optimization to a real molecular system. Between Avogadro’s number (6.022e23) and Boltzmann’s constant (1.38e-16) combined with very small distances (1.0e-8 cm), etc. the objective function values and the values of the optimization variables have very large values (~1e86 and ~1e9, respectively). I’ve verified that the analytic gradients of the objective function that I’m calculating are correct by comparing them with numerical derivatives.

>

> I’ve tried using the LMVM and Conjugate Gradient optimizations, both of which worked previously, but I find that the optimization completes one objective function evaluation and then declares that the problem is converged and stops. I could find a set of units where everything is approximately 1 but I was hoping that there are some parameters I can set in the optimization that will get it moving again. Any suggestions?

>

> Bruce Palmer

</pre></body></html>