Changeset 550


Ignore:
Timestamp:
Oct 24, 2005 5:54:44 PM (14 years ago)
Author:
andreasw
Message:

astyle

Location:
branches/dev-penalty/Algorithm
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/dev-penalty/Algorithm/IpCGPenaltyLSAcceptor.cpp

    r549 r550  
    8989      // ToDo we need to make this cleaner!
    9090      if (IpData().HaveCgPenDeltas()) {
    91         // use the fast step
    92         reference_direct_deriv_penalty_function_ =
    93           IpCq().curr_fast_direct_deriv_penalty_function();
    94         // Overwrite the original step (hmm.....)
    95         SmartPtr<const IteratesVector> delta_cgpen = IpData().delta_cgpen();
    96         IpData().set_delta(delta_cgpen);
     91        // use the fast step
     92        reference_direct_deriv_penalty_function_ =
     93          IpCq().curr_fast_direct_deriv_penalty_function();
     94        // Overwrite the original step (hmm.....)
     95        SmartPtr<const IteratesVector> delta_cgpen = IpData().delta_cgpen();
     96        IpData().set_delta(delta_cgpen);
    9797      }
    9898      else {
    99         reference_direct_deriv_penalty_function_ =
    100           IpCq().curr_direct_deriv_penalty_function();
     99        reference_direct_deriv_penalty_function_ =
     100          IpCq().curr_direct_deriv_penalty_function();
    101101      }
    102102    }
     
    125125    if (Jnlst().ProduceOutput(J_DETAILED, J_LINE_SEARCH)) {
    126126      Jnlst().Printf(J_DETAILED, J_LINE_SEARCH,
    127                      "curr_barr  = %23.16e curr_inf  = %23.16e\n",
    128                      IpCq().curr_barrier_obj(),
    129                      IpCq().curr_primal_infeasibility(NORM_2));
     127                     "curr_barr  = %23.16e curr_inf  = %23.16e\n",
     128                     IpCq().curr_barrier_obj(),
     129                     IpCq().curr_primal_infeasibility(NORM_2));
    130130      Jnlst().Printf(J_DETAILED, J_LINE_SEARCH,
    131                      "trial_barr = %23.16e trial_inf = %23.16e\n",
    132                      IpCq().trial_barrier_obj(),
    133                      IpCq().trial_primal_infeasibility(NORM_2));
     131                     "trial_barr = %23.16e trial_inf = %23.16e\n",
     132                     IpCq().trial_barrier_obj(),
     133                     IpCq().trial_primal_infeasibility(NORM_2));
    134134    }
    135135
     
    224224
    225225    Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    226                    "Starting tests for penalty parameter update:\n");
     226                   "Starting tests for penalty parameter update:\n");
    227227    // Lifeng: Should we use the new infeasibility here?
    228228    Number curr_inf = IpCq().curr_primal_infeasibility(NORM_2);
    229229    Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    230                    "current infeasibility = %8.2\n", curr_inf);
     230                   "current infeasibility = %8.2\n", curr_inf);
    231231    bool increase = (curr_inf >= penalty_update_infeasibility_tol_);
    232232    if (!increase) {
     
    236236    if (increase) {
    237237      Number max_step = Max(IpData().delta()->x()->Amax(),
    238                             IpData().delta()->s()->Amax());
    239       Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    240                      "Max norm of step = %8.2\n", max_step);
     238                            IpData().delta()->s()->Amax());
     239      Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
     240                     "Max norm of step = %8.2\n", max_step);
    241241      increase = (max_step <= curr_eta_);
    242242      if (!increase) {
     
    253253      if (IpNLP().x_L()->Dim()>0) {
    254254        SmartPtr<const Vector> compl_x_L = IpCq().curr_compl_x_L();
    255         min_compl = Min(min_compl, compl_x_L->Min());
    256         max_compl = Max(max_compl, compl_x_L->Max());
     255        min_compl = Min(min_compl, compl_x_L->Min());
     256        max_compl = Max(max_compl, compl_x_L->Max());
    257257      }
    258258      if (IpNLP().x_U()->Dim()>0) {
    259259        SmartPtr<const Vector> compl_x_U = IpCq().curr_compl_x_U();
    260         min_compl = Min(min_compl, compl_x_U->Min());
    261         max_compl = Max(max_compl, compl_x_U->Max());
     260        min_compl = Min(min_compl, compl_x_U->Min());
     261        max_compl = Max(max_compl, compl_x_U->Max());
    262262      }
    263263      if (IpNLP().d_L()->Dim()>0) {
    264264        SmartPtr<const Vector> compl_s_L = IpCq().curr_compl_s_L();
    265         min_compl = Min(min_compl, compl_s_L->Min());
    266         max_compl = Max(max_compl, compl_s_L->Max());
     265        min_compl = Min(min_compl, compl_s_L->Min());
     266        max_compl = Max(max_compl, compl_s_L->Max());
    267267      }
    268268      if (IpNLP().d_U()->Dim()>0) {
    269269        SmartPtr<const Vector> compl_s_U = IpCq().curr_compl_s_U();
    270         min_compl = Min(min_compl, compl_s_U->Min());
    271         max_compl = Max(max_compl, compl_s_U->Max());
    272       }
    273       Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    274                       "Minimal compl = %8.2\n", min_compl);
    275       Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    276                       "Maximal compl = %8.2\n", max_compl);
     270        min_compl = Min(min_compl, compl_s_U->Min());
     271        max_compl = Max(max_compl, compl_s_U->Max());
     272      }
     273      Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
     274                     "Minimal compl = %8.2\n", min_compl);
     275      Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
     276                     "Maximal compl = %8.2\n", max_compl);
    277277      increase = (min_compl >= mu*penalty_update_compl_tol_ &&
    278                   max_compl <= mu/penalty_update_compl_tol_);
     278                  max_compl <= mu/penalty_update_compl_tol_);
    279279      if (!increase) {
    280280        info_alpha_primal_char='c';
     
    291291      Number omega_test = vec->Amax();
    292292      Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    293                       "omega_test for c = %8.2\n", omega_test);
     293                     "omega_test for c = %8.2\n", omega_test);
    294294      increase = (omega_test < curr_eta_);
    295295      if (increase) {
     
    298298                           -1./IpCq().curr_cg_pert_fact(), *IpCq().curr_d_minus_s(),
    299299                           1.);
    300         omega_test = vec->Amax();
    301           Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    302                         "omega_test for d = %8.2\n", omega_test);
    303           increase = (omega_test < curr_eta_);
     300        omega_test = vec->Amax();
     301        Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
     302                      "omega_test for d = %8.2\n", omega_test);
     303        increase = (omega_test < curr_eta_);
    304304      }
    305305      if (!increase) {
     
    315315      curr_eta_ = Max(eta_min_, curr_eta_/2.);
    316316      Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    317                      "Updating eta to = %8.2\n", curr_eta_);
     317                     "Updating eta to = %8.2\n", curr_eta_);
    318318      Number penalty = IpData().curr_penalty();
    319319      Number y_full_step_max;
  • branches/dev-penalty/Algorithm/IpCGSearchDirCalc.cpp

    r549 r550  
    138138
    139139    Number hat_y_nrm = sqrt(pow(y_c->Nrm2(), 2.) +
    140                             pow(y_d->Nrm2(), 2.) +
    141                             2.*y_c->Dot(*delta_y_c) +
    142                             2.*y_d->Dot(*delta_y_d) +
    143                             pow(delta_y_c->Nrm2(), 2.) +
    144                             pow(delta_y_d->Nrm2(), 2.));
     140                            pow(y_d->Nrm2(), 2.) +
     141                            2.*y_c->Dot(*delta_y_c) +
     142                            2.*y_d->Dot(*delta_y_d) +
     143                            pow(delta_y_c->Nrm2(), 2.) +
     144                            pow(delta_y_d->Nrm2(), 2.));
    145145    Number diff_y = sqrt(pow(delta_fast_y_c->Nrm2(), 2.) +
    146                         pow(delta_fast_y_d->Nrm2(), 2.) -
    147                         2.*delta_y_c->Dot(*delta_fast_y_c) -
    148                         2.*delta_y_d->Dot(*delta_fast_y_d) +
    149                         pow(delta_y_c->Nrm2(), 2.) +
    150                         pow(delta_y_d->Nrm2(), 2.));
     146                        pow(delta_fast_y_d->Nrm2(), 2.) -
     147                        2.*delta_y_c->Dot(*delta_fast_y_c) -
     148                        2.*delta_y_d->Dot(*delta_fast_y_d) +
     149                        pow(delta_y_c->Nrm2(), 2.) +
     150                        pow(delta_y_d->Nrm2(), 2.));
    151151
    152152    Number kappa_dis_ = 1e1;
     
    155155      IpData().Append_info_string("G");
    156156    }
    157  
     157
    158158    if (keep_fast_delta) {
    159159      // For now, I just check if the directional derivative for the
     
    162162      Number fast_direct_deriv = IpCq().curr_fast_direct_deriv_penalty_function();
    163163      Jnlst().Printf(J_MOREDETAILED, J_LINE_SEARCH,
    164                      "direct_deriv = %23.16e  fast_direct_deriv = %23.15e\n",
    165                      direct_deriv, fast_direct_deriv);
     164                     "direct_deriv = %23.16e  fast_direct_deriv = %23.15e\n",
     165                     direct_deriv, fast_direct_deriv);
    166166      Number need_name_ = 1e-1;
    167167      if (fast_direct_deriv > need_name_*direct_deriv) {
    168         // Discard the fast direction
    169         //delta_fast = NULL;
    170         //IpData().set_delta_cgpen(delta_fast);
    171         keep_fast_delta = false;
    172         IpData().Append_info_string("g");
     168        // Discard the fast direction
     169        //delta_fast = NULL;
     170        //IpData().set_delta_cgpen(delta_fast);
     171        keep_fast_delta = false;
     172        IpData().Append_info_string("g");
    173173      }
    174174    }
  • branches/dev-penalty/Algorithm/IpIpoptData.hpp

    r549 r550  
    230230      have_cgpen_deltas_ = have_cgpen_deltas;
    231231    }
    232    
     232
    233233
    234234    /** @name Public Methods for updating iterates */
Note: See TracChangeset for help on using the changeset viewer.