Creating a Grid Integrated Plant Design Problem using IDAES Surrogates
Import necessary packages
To begin, we need to import necessary packages. The primary packages required are: * Pyomo : Modeling the optimization problem * IDAES: Surrogate modeling capabilities * Tensorflow: load trained Keras surrogates that represent market simulations
[1]:
import json
import pandas as pd
import tensorflow
#import pyomo and its objects
from pyomo.environ import (
ConcreteModel,
Var,
Constraint,
Expression,
NonNegativeReals,
Block,
value,
Objective,
)
import pyomo.environ as pyo
# Import IDAES components
from idaes.core.util import to_json, from_json
from idaes.core.solvers import get_solver
#Import IDAES SurrogateBlock
from idaes.core.surrogate.keras_surrogate import KerasSurrogate
from idaes.core.surrogate.sampling.scaling import OffsetScaler
from idaes.core.surrogate.surrogate_block import SurrogateBlock
# Import simple rankine cycle flowsheet
import dispatches.case_studies.simple_rankine_cycle.simple_rankine_cycle as src
Load Keras Surrogates
The first step towards assembling our optimization model is to load the pre-trained Keras neural networks (available in this DISPATCHES repository) as well as associated meta-data saved in corresponding .json files. This meta data contains input and output labels, scaling information, and the input bounds used during training.
The keras models themselves represent the results of market simulations using the Prescient Production Cost Modeling open-source package. We previously trained neural network models to predict the following outputs based on market inputs: * Annual Revenue [MM$]: (the total amount of revenue our generator earned per year) * Annual Number of Generator Startups * Annual Zone Output [hours]: (the amount of time our generator spent at different power output intervals). For this neural network, we trained 11 different zone outputs (off,0-10% of pmax,10-20% of pmax,…,90-100% of pmax)
For market inputs, our trianing data consisted of 8 primary attributes: * p_max (the nameplate plant capacity) * p_min_multipler (the minimum operating output as a fration of p_max) * ramp_rate_multiplier * minimum_up_time * minimum_dn_multiplier * marginal_cost * fixed_run_cost * startup_profile
[2]:
# load training meta-data for each surrogate. this includes things like scaling
# information, labels, and input bounds
with open("keras_models/models/training_parameters_revenue.json", 'rb') as f:
rev_data = json.load(f)
with open("keras_models/models/training_parameters_zones.json", 'rb') as f:
zone_data = json.load(f)
with open("keras_models/models/training_parameters_nstartups.json", 'rb') as f:
nstartups_data = json.load(f)
# unpack scaling information for revenue
xm_rev = rev_data['xm_inputs']
xstd_rev = rev_data['xstd_inputs']
zm_rev = rev_data['zm_revenue']
zstd_rev = rev_data['zstd_revenue']
# unpack scaling information for number of startups
xm_nstart = nstartups_data['xm_inputs']
xstd_nstart = nstartups_data['xstd_inputs']
zm_nstartups = nstartups_data['zm_nstartups']
zstd_nstartups = nstartups_data['zstd_nstartups']
# unpack scaling information for duration at different operating outputs
xm_zones = zone_data['xm_inputs']
xstd_zones = zone_data['xstd_inputs']
zm_zones = zone_data['zm_zones']
zstd_zones = zone_data['zstd_zones']
# load surrogates from keras directories
keras_revenue = tensorflow.keras.models.load_model('keras_models/models/keras_revenue', compile=False)
keras_nstartups = tensorflow.keras.models.load_model('keras_models/models/keras_nstartups', compile=False)
keras_zones = tensorflow.keras.models.load_model('keras_models/models/keras_zones', compile=False)
WARNING:tensorflow:SavedModel saved prior to TF 2.5 detected when loading Keras model. Please ensure that you are saving the model with model.save() or tf.keras.models.save_model(), *NOT* tf.saved_model.save(). To confirm, there should be a file named "keras_metadata.pb" in the SavedModel directory.
WARNING:tensorflow:SavedModel saved prior to TF 2.5 detected when loading Keras model. Please ensure that you are saving the model with model.save() or tf.keras.models.save_model(), *NOT* tf.saved_model.save(). To confirm, there should be a file named "keras_metadata.pb" in the SavedModel directory.
WARNING:tensorflow:SavedModel saved prior to TF 2.5 detected when loading Keras model. Please ensure that you are saving the model with model.save() or tf.keras.models.save_model(), *NOT* tf.saved_model.save(). To confirm, there should be a file named "keras_metadata.pb" in the SavedModel directory.
Create IDAES Keras Surrogate Objects
Now that we have our keras models loaded we can create KerasSurrogate objects using IDAES. These objects wrap our trained Keras Sequential models and store the additional training information we loaded such as scaling. The KerasSurrogate object provides a consistent interface for handling surrogate models in IDAES.
Once we have a KerasSurrogate, we can use it to build a mathematical model on a SurrogateBlock. This is shown later in this notebook.
[3]:
# revenue Keras surrogate
# load extra meta-data for the revenue surrogate
# input_labels = list('x'+str(i) for i in range(len(xm_rev)))
# output_labels = ["z_rev"]
input_labels = rev_data['input_labels']
output_labels = rev_data['output_labels']
xmin = rev_data['xmin']
xmax = rev_data['xmax']
# create scaling objects for the neural network input and outputs
inputs_scaler = OffsetScaler(
expected_columns=input_labels,
offset_series = pd.Series(dict(zip(input_labels, xm_rev))),
factor_series = pd.Series(dict(zip(input_labels, xstd_rev))),
)
outputs_scaler = OffsetScaler(
expected_columns=output_labels,
offset_series = pd.Series(dict(zip(output_labels, [zm_rev]))),
factor_series = pd.Series(dict(zip(output_labels, [zstd_rev]))),
)
# store input bounds in a dictionary where keys correspond to input labels and values are 2-length tuples
# of the lower and upper bound respectively
input_bounds={input_labels[i]: (xmin[i], xmax[i]) for i in range(len(input_labels))}
# create the KerasSurrogate object using the keras neural network and associated training data
keras_revenue_surrogate = KerasSurrogate(
keras_model=keras_revenue,
input_labels=input_labels,
output_labels=output_labels,
input_bounds=input_bounds,
input_scaler=inputs_scaler,
output_scaler=outputs_scaler,
)
# nstartups Keras surrogate
input_labels = nstartups_data['input_labels']
output_labels = nstartups_data['output_labels']
xmin = nstartups_data['xmin']
xmax = nstartups_data['xmax']
inputs_scaler = OffsetScaler(
expected_columns=input_labels,
offset_series = pd.Series(dict(zip(input_labels, xm_nstart))),
factor_series = pd.Series(dict(zip(input_labels, xstd_nstart))),
)
outputs_scaler = OffsetScaler(
expected_columns=output_labels,
offset_series = pd.Series(dict(zip(output_labels,[zm_nstartups]))),
factor_series = pd.Series(dict(zip(output_labels,[zstd_nstartups]))),
)
input_bounds={input_labels[i]: (xmin[i], xmax[i]) for i in range(len(input_labels))}
keras_nstartups_surrogate = KerasSurrogate(
keras_model=keras_nstartups,
input_labels=input_labels,
output_labels=output_labels,
input_bounds=input_bounds,
input_scaler=inputs_scaler,
output_scaler=outputs_scaler,
)
# zone hour power output Keras surrogate
input_labels = zone_data['input_labels']
output_labels = zone_data['output_labels']
xmin = zone_data['xmin']
xmax = zone_data['xmax']
inputs_scaler = OffsetScaler(
expected_columns=input_labels,
offset_series = pd.Series(dict(zip(input_labels, xm_zones))),
factor_series = pd.Series(dict(zip(input_labels, xstd_zones))),
)
outputs_scaler = OffsetScaler(
expected_columns=output_labels,
offset_series = pd.Series(dict(zip(output_labels, zm_zones))),
factor_series = pd.Series(dict(zip(output_labels, zstd_zones))),
)
input_bounds={input_labels[i]: (xmin[i], xmax[i]) for i in range(len(input_labels))}
keras_zones_surrogate = KerasSurrogate(
keras_model=keras_zones,
input_labels=input_labels,
output_labels=output_labels,
input_bounds=input_bounds,
input_scaler=inputs_scaler,
output_scaler=outputs_scaler,
)
Building the Rankine Cycle Model
We now build the rankine cycle model available from DISPATCHES. This cell creates a pyomo model and uses rankine cycle functions to construct a design flowsheet. More specifically, we construct a high level flowsheet to represent the plant design. We will later construct additional flowsheets that represent plant operation for different power outputs.
[4]:
# rankine cycle parameters
heat_recovery = True
calc_boiler_eff = True
capital_payment_years = 5
plant_lifetime = 20
coal_price = 30 #$/ton
m = ConcreteModel()
# Create capex plant
m.cap_fs = src.create_model(
heat_recovery=heat_recovery,
capital_fs=True,
calc_boiler_eff=False,
)
src.set_inputs(m.cap_fs)
src.initialize_model(m.cap_fs)
src.close_flowsheet_loop(m.cap_fs)
src.add_capital_cost(m.cap_fs)
# capital cost (M$/yr)
cap_expr = m.cap_fs.fs.capital_cost / capital_payment_years
2023-05-19 16:43:10 [INFO] idaes.init.cap_fs.fs.boiler.control_volume: Initialization Complete
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.boiler: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.turbine: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.pre_condenser.control_volume: Initialization Complete
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.pre_condenser: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.condenser.control_volume: Initialization Complete
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.condenser: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.bfw_pump.control_volume: Initialization Complete
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.bfw_pump: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.feed_water_heater.control_volume: Initialization Complete
2023-05-19 16:43:11 [INFO] idaes.init.cap_fs.fs.feed_water_heater: Initialization Complete: optimal - Optimal Solution Found
Ipopt 3.13.2: nlp_scaling_method=gradient-based
tol=1e-06
******************************************************************************
This program contains Ipopt, a library for large-scale nonlinear optimization.
Ipopt is released as open source code under the Eclipse Public License (EPL).
For more information visit http://projects.coin-or.org/Ipopt
This version of Ipopt was compiled from source code available at
https://github.com/IDAES/Ipopt as part of the Institute for the Design of
Advanced Energy Systems Process Systems Engineering Framework (IDAES PSE
Framework) Copyright (c) 2018-2019. See https://github.com/IDAES/idaes-pse.
This version of Ipopt was compiled using HSL, a collection of Fortran codes
for large-scale scientific computation. All technical papers, sales and
publicity material resulting from use of the HSL codes within IPOPT must
contain the following acknowledgement:
HSL, a collection of Fortran codes for large-scale scientific
computation. See http://www.hsl.rl.ac.uk.
******************************************************************************
This is Ipopt version 3.13.2, running with linear solver ma27.
Number of nonzeros in equality constraint Jacobian...: 107
Number of nonzeros in inequality constraint Jacobian.: 0
Number of nonzeros in Lagrangian Hessian.............: 27
Total number of variables............................: 44
variables with only lower bounds: 0
variables with lower and upper bounds: 21
variables with only upper bounds: 0
Total number of equality constraints.................: 44
Total number of inequality constraints...............: 0
inequality constraints with only lower bounds: 0
inequality constraints with lower and upper bounds: 0
inequality constraints with only upper bounds: 0
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 0.0000000e+00 1.74e+08 0.00e+00 -1.0 0.00e+00 - 0.00e+00 0.00e+00 0
1 0.0000000e+00 5.96e-08 1.83e+00 -1.0 1.74e+08 - 5.66e-01 1.00e+00h 1
Number of Iterations....: 1
(scaled) (unscaled)
Objective...............: 0.0000000000000000e+00 0.0000000000000000e+00
Dual infeasibility......: 0.0000000000000000e+00 0.0000000000000000e+00
Constraint violation....: 3.8518889586569750e-09 5.9604644775390625e-08
Complementarity.........: 0.0000000000000000e+00 0.0000000000000000e+00
Overall NLP error.......: 3.8518889586569750e-09 5.9604644775390625e-08
Number of objective function evaluations = 2
Number of objective gradient evaluations = 2
Number of equality constraint evaluations = 2
Number of inequality constraint evaluations = 0
Number of equality constraint Jacobian evaluations = 2
Number of inequality constraint Jacobian evaluations = 0
Number of Lagrangian Hessian evaluations = 1
Total CPU secs in IPOPT (w/o function evaluations) = 0.018
Total CPU secs in NLP function evaluations = 0.004
EXIT: Optimal Solution Found.
Build the grid market surrogates
We now create variables to represent the previously defined market inputs such as name-plate capacity, ramp rate, marginal cost, etc…
We use the IDAES SurrogateBlock to build the neural network constraints by using the previously creates KerasSurrogate objects. Notably, the build function automatically performs the input and output scaling because we provided the scaling information to each KerasSurrogate.
[5]:
########################################
# Surrogate market inputs
########################################
m.pmax = Var(within=NonNegativeReals, bounds=(175, 450), initialize=400) #MW
m.pmin_multi = Var(within=NonNegativeReals, bounds=(0.15, 0.45), initialize=0.3)
m.ramp_multi = Var(within=NonNegativeReals, bounds=(0.5, 1.0), initialize=0.75)
m.min_up_time = Var(within=NonNegativeReals, bounds=(1.0, 16.0), initialize=4.0)
m.min_dn_multi = Var(within=NonNegativeReals, bounds=(0.5, 2.0), initialize=1.0)
m.marg_cst = Var(within=NonNegativeReals, bounds=(5, 30), initialize=15)
m.no_load_cst = Var(within=NonNegativeReals, bounds=(0, 2.5), initialize=1)
m.startup_cst = Var(within=NonNegativeReals, bounds=(0, 136), initialize=75)
m.pmax_con = Constraint(expr=m.pmax == m.cap_fs.fs.net_cycle_power_output * 1e-6)
# Actual generator values for minimum operating output, minimum down time, and true ramp rate
m.pmin = Expression(expr=m.pmin_multi * m.pmax) # MW
m.min_dn_time = Expression(expr=m.min_dn_multi * m.min_up_time) # hours
m.ramp_rate = Expression(expr=m.ramp_multi * (m.pmax - m.pmin)) # MW/hr
[6]:
######################################
# Revenue surrogate
######################################
m.rev_surrogate = Var()
m.keras_revenue_surrogate = SurrogateBlock()
m.keras_revenue_surrogate.build_model(
keras_revenue_surrogate,
input_vars=[
m.pmax, m.pmin_multi, m.ramp_multi, m.min_up_time,
m.min_dn_multi, m.marg_cst, m.no_load_cst, m.startup_cst,
],
output_vars = [m.rev_surrogate],
formulation=KerasSurrogate.Formulation.REDUCED_SPACE,
)
# this is a smooth-max; it sets negative revenue to zero
m.revenue = Expression(expr=0.5*pyo.sqrt(m.rev_surrogate**2 + 0.001**2) + 0.5*m.rev_surrogate)
[7]:
#######################################
#nstartups surrogate
#######################################
m.nstartups_surrogate = Var()
m.keras_nstartups_surrogate = SurrogateBlock()
m.keras_nstartups_surrogate.build_model(
keras_nstartups_surrogate,
input_vars=[
m.pmax, m.pmin_multi, m.ramp_multi, m.min_up_time,
m.min_dn_multi, m.marg_cst, m.no_load_cst, m.startup_cst,
],
output_vars = [m.rev_surrogate],
formulation=KerasSurrogate.Formulation.REDUCED_SPACE,
)
m.nstartups = Expression(
expr=0.5*pyo.sqrt(m.nstartups_surrogate**2 + 0.001**2) + 0.5*m.nstartups_surrogate
)
[8]:
############################################
# zone surrogates
############################################
m.zone_hours_surrogate = Var(range(11))
m.keras_zones_surrogate = SurrogateBlock()
m.keras_zones_surrogate.build_model(
keras_zones_surrogate,
input_vars=[
m.pmax,m.pmin_multi,m.ramp_multi,m.min_up_time,
m.min_dn_multi,m.marg_cst,m.no_load_cst,m.startup_cst
],
output_vars = m.zone_hours_surrogate,
formulation=KerasSurrogate.Formulation.REDUCED_SPACE,
)
Build the Operation Flowsheets
We now create a flowsheet for each possible operating zone. Each flowsheet uses the market decisions defined above and we use the zone surrogates to predict the number of hours spent at each flowsheet’s powr output.
[9]:
off_fs = Block()
off_fs.fs = Block()
off_fs.fs.operating_cost = m.no_load_cst * m.pmax
off_fs.zone_hours = Expression(
expr=0.5*pyo.sqrt(m.zone_hours_surrogate[0]**2 + 0.001**2) + 0.5*m.zone_hours_surrogate[0]
)
setattr(m, 'zone_{}'.format('off'), off_fs)
# Denote the scaled power output for each of the 10 zones (0 corresponds to pmin, 1.0 corresponds to pmax)
zone_outputs = [0.0, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, 1.0]
# Create a surrogate flowsheet for each operating zone
op_zones = []
init_flag = 0
for (i, zone_output) in enumerate(zone_outputs):
print("Creating instance ", i)
op_fs = src.create_model(
heat_recovery=heat_recovery,
capital_fs=False,
calc_boiler_eff=calc_boiler_eff,
)
# Set model inputs for the capex and opex plant
src.set_inputs(op_fs)
# Fix the p_max of op_fs to p of cap_fs for initialization
op_fs.fs.net_power_max.fix(value(m.cap_fs.fs.net_cycle_power_output))
#initialize with json. this speeds up model instantiation. it writes a json file \
#for the first flowsheet which is used to initialize the next flowsheets
if init_flag == 0:
# Initialize the opex plant
src.initialize_model(op_fs)
# save model state after initializing the first instance
init_model = to_json(op_fs.fs, return_dict=True)
init_flag = 1
else:
# Initialize the capex and opex plant
from_json(op_fs.fs, sd=init_model)
# Closing the loop in the flowsheet
src.close_flowsheet_loop(op_fs)
src.add_operating_cost(op_fs, coal_price=coal_price)
# Unfix op_fs p_max and set constraint linking that to cap_fs p_max
op_fs.fs.net_power_max.unfix()
op_fs.fs.eq_p_max = Constraint(
expr=op_fs.fs.net_power_max == m.cap_fs.fs.net_cycle_power_output * 1e-6
)
# Fix zone power output
op_fs.fs.eq_fix_power = Constraint(
expr=op_fs.fs.net_cycle_power_output * 1e-6 == zone_output * (m.pmax-m.pmin) + m.pmin
)
# smooth max on zone hours (avoids negative hours)
op_fs.zone_hours = Expression(
expr=0.5*pyo.sqrt(m.zone_hours_surrogate[i+1]**2 + 0.001**2) + 0.5*m.zone_hours_surrogate[i+1]
)
# unfix the boiler flow rate
op_fs.fs.boiler.inlet.flow_mol[0].setlb(0.01)
op_fs.fs.boiler.inlet.flow_mol[0].unfix()
setattr(m, 'zone_{}'.format(i), op_fs)
op_zones.append(op_fs)
# scale zone hours such that they add up to 8736 (if the surrogate is good, the unscaled will be pretty close to this)
m.zone_total_hours = sum(op_zones[i].zone_hours for i in range(len(op_zones))) + off_fs.zone_hours
for op_fs in op_zones:
op_fs.scaled_zone_hours = Var(within=NonNegativeReals, bounds=(0, 8736), initialize=100)
# NOTE: scaled_hours_i = surrogate_i * 8736 / surrogate_total
op_fs.con_scale_zone_hours = Constraint(
expr=op_fs.scaled_zone_hours * m.zone_total_hours == op_fs.zone_hours * 8736
)
off_fs.scaled_zone_hours = Var(within=NonNegativeReals, bounds=(0, 8736), initialize=100)
off_fs.con_scale_zone_hours = Constraint(
expr=off_fs.scaled_zone_hours * m.zone_total_hours == off_fs.zone_hours * 8736
)
#operating cost in $MM (million dollars)
m.op_expr = sum(
op_zones[i].scaled_zone_hours * op_zones[i].fs.operating_cost
for i in range(len(op_zones))
) * 1e-6 + off_fs.scaled_zone_hours * off_fs.fs.operating_cost * 1e-6
#startup cost in MM$
m.startup_expr = m.startup_cst * m.nstartups * m.pmax * 1e-6 #MM$
# set zone flowsheets to pyomo model
m.op_zones = op_zones
# Piecewise cost limits, connect marginal cost to operating cost. We say marginal cost is the average operating cost
m.connect_mrg_cost = Constraint(
expr=m.marg_cst == 0.5*(op_zones[0].fs.operating_cost/m.pmin + op_zones[-1].fs.operating_cost/m.pmax)
)
# Expression for total cap and op cost - $
m.total_cost = Expression(
expr=plant_lifetime * (m.op_expr + m.startup_expr) + capital_payment_years*cap_expr
)
# Expression for total revenue
m.total_revenue = Expression(expr=plant_lifetime*m.revenue)
# Objective $
m.obj = Objective(expr=-(m.total_revenue - m.total_cost))
# Unfixing the boiler inlet flowrate for capex plant
m.cap_fs.fs.boiler.inlet.flow_mol[0].unfix()
# Setting bounds for the capex plant flowrate
m.cap_fs.fs.boiler.inlet.flow_mol[0].setlb(0.01)
# Setting bounds for net cycle power output for the capex plant
p_lower_bound=10
p_upper_bound=500
m.cap_fs.fs.eq_min_power = Constraint(
expr=m.cap_fs.fs.net_cycle_power_output >= p_lower_bound*1e6)
m.cap_fs.fs.eq_max_power = Constraint(
expr=m.cap_fs.fs.net_cycle_power_output <= p_upper_bound*1e6)
Creating instance 0
2023-05-19 16:43:19 [INFO] idaes.init.fs.boiler.control_volume: Initialization Complete
2023-05-19 16:43:19 [INFO] idaes.init.fs.boiler: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:19 [INFO] idaes.init.fs.turbine: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:19 [INFO] idaes.init.fs.pre_condenser.control_volume: Initialization Complete
2023-05-19 16:43:19 [INFO] idaes.init.fs.pre_condenser: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:19 [INFO] idaes.init.fs.condenser.control_volume: Initialization Complete
2023-05-19 16:43:19 [INFO] idaes.init.fs.condenser: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:19 [INFO] idaes.init.fs.bfw_pump.control_volume: Initialization Complete
2023-05-19 16:43:19 [INFO] idaes.init.fs.bfw_pump: Initialization Complete: optimal - Optimal Solution Found
2023-05-19 16:43:19 [INFO] idaes.init.fs.feed_water_heater.control_volume: Initialization Complete
2023-05-19 16:43:19 [INFO] idaes.init.fs.feed_water_heater: Initialization Complete: optimal - Optimal Solution Found
Ipopt 3.13.2: nlp_scaling_method=gradient-based
tol=1e-06
******************************************************************************
This program contains Ipopt, a library for large-scale nonlinear optimization.
Ipopt is released as open source code under the Eclipse Public License (EPL).
For more information visit http://projects.coin-or.org/Ipopt
This version of Ipopt was compiled from source code available at
https://github.com/IDAES/Ipopt as part of the Institute for the Design of
Advanced Energy Systems Process Systems Engineering Framework (IDAES PSE
Framework) Copyright (c) 2018-2019. See https://github.com/IDAES/idaes-pse.
This version of Ipopt was compiled using HSL, a collection of Fortran codes
for large-scale scientific computation. All technical papers, sales and
publicity material resulting from use of the HSL codes within IPOPT must
contain the following acknowledgement:
HSL, a collection of Fortran codes for large-scale scientific
computation. See http://www.hsl.rl.ac.uk.
******************************************************************************
This is Ipopt version 3.13.2, running with linear solver ma27.
Number of nonzeros in equality constraint Jacobian...: 107
Number of nonzeros in inequality constraint Jacobian.: 0
Number of nonzeros in Lagrangian Hessian.............: 27
Total number of variables............................: 44
variables with only lower bounds: 0
variables with lower and upper bounds: 21
variables with only upper bounds: 0
Total number of equality constraints.................: 44
Total number of inequality constraints...............: 0
inequality constraints with only lower bounds: 0
inequality constraints with lower and upper bounds: 0
inequality constraints with only upper bounds: 0
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 0.0000000e+00 1.74e+08 0.00e+00 -1.0 0.00e+00 - 0.00e+00 0.00e+00 0
1 0.0000000e+00 5.96e-08 1.83e+00 -1.0 1.74e+08 - 5.66e-01 1.00e+00h 1
Number of Iterations....: 1
(scaled) (unscaled)
Objective...............: 0.0000000000000000e+00 0.0000000000000000e+00
Dual infeasibility......: 0.0000000000000000e+00 0.0000000000000000e+00
Constraint violation....: 3.8518889586569750e-09 5.9604644775390625e-08
Complementarity.........: 0.0000000000000000e+00 0.0000000000000000e+00
Overall NLP error.......: 3.8518889586569750e-09 5.9604644775390625e-08
Number of objective function evaluations = 2
Number of objective gradient evaluations = 2
Number of equality constraint evaluations = 2
Number of inequality constraint evaluations = 0
Number of equality constraint Jacobian evaluations = 2
Number of inequality constraint Jacobian evaluations = 0
Number of Lagrangian Hessian evaluations = 1
Total CPU secs in IPOPT (w/o function evaluations) = 0.015
Total CPU secs in NLP function evaluations = 0.005
EXIT: Optimal Solution Found.
Creating instance 1
Creating instance 2
Creating instance 3
Creating instance 4
Creating instance 5
Creating instance 6
Creating instance 7
Creating instance 8
Creating instance 9
Setup Final Surrogate Inputs and Solve
We lastly fix a few surrogate inputs and solve the design problem. The solution of the design problem will determine market inputs that maximize the net plant revenue over 20 years.
[10]:
# these are representative startup costs based on startup profiles we trained on.
# it is useful to fix the startup cost to one of these values since they are technically categorical variables
startup_csts = [0., 49.66991167, 61.09068702, 101.4374234, 135.2230393]
#fix some surrogate inputs
start_cst_index=2
m.startup_cst.fix(startup_csts[start_cst_index])
m.no_load_cst.fix(1.0)
m.min_up_time.fix(4.0)
m.min_dn_multi.fix(1.0)
solver = get_solver()
solver.options = {
"tol": 1e-6
#"mu_strategy": "adaptive"
}
status = solver.solve(m, tee=True)
sol_time = status['Solver'][0]['Time']
Ipopt 3.13.2: tol=1e-06
******************************************************************************
This program contains Ipopt, a library for large-scale nonlinear optimization.
Ipopt is released as open source code under the Eclipse Public License (EPL).
For more information visit http://projects.coin-or.org/Ipopt
This version of Ipopt was compiled from source code available at
https://github.com/IDAES/Ipopt as part of the Institute for the Design of
Advanced Energy Systems Process Systems Engineering Framework (IDAES PSE
Framework) Copyright (c) 2018-2019. See https://github.com/IDAES/idaes-pse.
This version of Ipopt was compiled using HSL, a collection of Fortran codes
for large-scale scientific computation. All technical papers, sales and
publicity material resulting from use of the HSL codes within IPOPT must
contain the following acknowledgement:
HSL, a collection of Fortran codes for large-scale scientific
computation. See http://www.hsl.rl.ac.uk.
******************************************************************************
This is Ipopt version 3.13.2, running with linear solver ma27.
Number of nonzeros in equality constraint Jacobian...: 1749
Number of nonzeros in inequality constraint Jacobian.: 4
Number of nonzeros in Lagrangian Hessian.............: 713
Total number of variables............................: 641
variables with only lower bounds: 11
variables with lower and upper bounds: 317
variables with only upper bounds: 0
Total number of equality constraints.................: 638
Total number of inequality constraints...............: 2
inequality constraints with only lower bounds: 1
inequality constraints with lower and upper bounds: 0
inequality constraints with only upper bounds: 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
0 5.9764370e+02 9.25e+07 1.00e+01 -1.0 0.00e+00 - 0.00e+00 0.00e+00 0
1r 5.9764370e+02 9.25e+07 9.99e+02 8.0 0.00e+00 - 0.00e+00 7.31e-08R 2
2r 5.9765704e+02 9.25e+07 2.99e+08 8.0 1.85e+07 - 1.50e-02 2.13e-07f 1
3r 5.9781099e+02 8.66e+07 7.93e+06 3.1 9.66e+06 - 9.73e-01 6.18e-01f 1
4r 5.9865641e+02 8.48e+05 2.82e+06 3.1 2.67e+09 - 6.44e-01 3.21e-02f 1
5r 5.9865641e+02 8.48e+05 9.99e+02 5.9 0.00e+00 - 0.00e+00 7.59e-08R 2
6r 5.9867391e+02 8.48e+05 2.65e+05 5.9 7.39e+06 - 1.82e-03 2.86e-06f 1
7r 6.1547760e+02 8.45e+05 1.44e+05 5.9 7.35e+06 - 5.87e-01 3.24e-03f 1
8r 6.1825745e+02 7.75e+05 8.40e+04 3.8 4.02e+06 - 6.88e-01 7.64e-02f 1
9r 6.2787513e+02 1.72e+06 7.21e+04 3.8 1.00e+07 - 1.66e-01 7.68e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
10r 6.2787513e+02 1.72e+06 9.99e+02 4.6 0.00e+00 - 0.00e+00 7.95e-08R 2
11r 6.2807616e+02 1.69e+06 8.37e+03 4.6 3.16e+07 - 3.94e-03 1.34e-04f 1
12r 6.3353920e+02 1.19e+06 7.09e+03 4.6 1.73e+07 - 5.12e-02 3.87e-03f 1
13r 6.3353920e+02 1.19e+06 9.99e+02 4.0 0.00e+00 - 0.00e+00 8.76e-08R 2
14r 6.3577297e+02 5.67e+05 9.97e+02 4.0 6.05e+06 - 1.74e-02 1.82e-03f 1
15r 6.5429976e+02 2.99e+06 9.79e+02 4.0 2.91e+05 - 2.18e-02 1.74e-02f 1
16r 6.5429976e+02 2.99e+06 9.99e+02 3.9 0.00e+00 - 0.00e+00 5.25e-08R 2
17r 6.5542718e+02 8.75e+05 9.99e+02 3.9 7.18e+06 - 1.61e-02 1.02e-03f 1
18r 6.5542718e+02 8.75e+05 9.99e+02 3.8 0.00e+00 - 0.00e+00 5.18e-08R 2
19r 6.6031802e+02 1.91e+05 1.13e+03 3.8 3.27e+05 - 2.87e-02 4.23e-03f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
20r 6.7279643e+02 7.27e+05 1.47e+04 3.8 2.01e+05 - 9.66e-02 1.38e-02f 1
21r 6.7279643e+02 7.27e+05 9.99e+02 3.6 0.00e+00 - 0.00e+00 1.12e-07R 2
22r 6.7531771e+02 8.85e+04 4.83e+04 3.6 7.17e+05 - 3.22e-01 2.47e-03f 1
23r 6.7853728e+02 2.84e+05 1.01e+05 2.9 1.91e+05 - 2.95e-01 6.79e-02f 1
24r 6.7853728e+02 2.84e+05 9.99e+02 3.4 0.00e+00 - 0.00e+00 2.19e-07R 2
25r 6.8062861e+02 9.75e+04 4.22e+04 3.4 2.65e+05 - 2.23e-01 2.19e-03f 1
26r 6.9617345e+02 5.53e+05 6.73e+04 2.7 2.63e+05 - 1.66e-01 9.57e-02f 1
27r 6.9716796e+02 1.95e+07 6.66e+04 2.7 5.65e+05 -4.0 3.53e-03 8.67e-03f 1
28r 6.9743701e+02 1.95e+07 7.00e+04 2.7 3.93e+05 - 2.22e-01 2.13e-03f 1
29r 7.0619411e+02 1.82e+07 6.04e+04 2.7 4.51e+05 - 1.02e-01 7.49e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
30r 7.2317579e+02 1.53e+07 1.32e+04 2.7 4.23e+05 - 5.41e-01 1.63e-01f 1
31r 7.2325691e+02 1.50e+07 1.36e+04 2.7 2.93e+02 1.8 8.97e-03 1.37e-02f 1
32r 7.2330091e+02 1.48e+07 1.36e+04 2.7 2.83e+02 2.3 1.28e-02 1.43e-02f 1
33r 7.2426858e+02 1.47e+07 1.58e+04 2.7 3.80e+05 - 1.12e-01 1.12e-02f 1
34r 7.5941311e+02 8.34e+06 2.06e+04 2.7 3.67e+05 - 2.80e-01 4.29e-01f 1
35r 7.5941317e+02 8.31e+06 2.05e+04 2.7 9.48e+01 4.5 8.01e-03 4.16e-03h 1
36r 7.5941336e+02 8.27e+06 2.05e+04 2.7 8.71e+01 4.0 9.45e-03 4.07e-03h 1
37r 7.5941354e+02 8.19e+06 6.30e+04 2.7 8.71e+01 4.5 4.21e-02 1.01e-02h 1
38r 7.5941564e+02 7.87e+06 7.15e+04 2.7 8.61e+01 4.0 4.18e-02 3.91e-02h 1
39r 7.5941629e+02 7.61e+06 2.39e+05 2.7 8.28e+01 4.4 1.47e-01 3.29e-02h 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
40r 7.5942229e+02 6.87e+06 2.28e+05 2.7 8.00e+01 3.9 1.59e-01 9.77e-02h 1
41r 7.5942450e+02 6.20e+06 5.60e+05 2.7 7.22e+01 4.3 4.77e-01 9.76e-02h 1
42r 7.5945567e+02 3.48e+06 3.97e+05 2.7 6.51e+01 3.9 4.20e-01 4.38e-01h 1
43r 7.5946884e+02 1.75e+06 5.61e+05 2.7 3.66e+01 4.3 1.00e+00 4.97e-01h 1
44r 7.5950265e+02 1.62e+04 1.22e+05 2.0 1.84e+01 3.8 8.89e-01 1.00e+00h 1
45r 7.5961187e+02 1.52e+04 1.12e+03 2.0 5.08e-01 3.3 1.00e+00 1.00e+00f 1
46r 7.5984431e+02 1.51e+04 1.16e+03 1.3 1.66e+00 2.9 9.89e-01 2.99e-01f 1
47r 7.6004837e+02 1.50e+04 1.75e+03 1.3 7.60e-01 3.3 1.00e+00 4.98e-01f 1
48r 7.6036909e+02 1.49e+04 1.54e+03 1.3 1.97e+00 2.8 5.15e-01 2.02e-01f 1
49r 7.6069084e+02 1.47e+04 1.74e+03 1.3 7.36e-01 3.2 1.00e+00 5.78e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
50r 7.6077008e+02 1.47e+04 1.77e+03 1.3 3.03e+00 2.8 7.73e-02 6.96e-02f 1
51r 7.6102731e+02 1.45e+04 1.61e+03 1.3 6.87e+00 2.3 5.76e-01 1.04e-01f 1
52r 7.6128843e+02 1.38e+04 1.44e+03 1.3 2.01e+01 1.8 3.05e-01 1.10e-01f 1
53r 7.6137803e+02 1.35e+04 1.31e+03 1.3 7.07e+00 2.2 5.12e-01 1.28e-01f 1
54r 7.6138622e+02 1.35e+04 3.40e+03 1.3 2.42e+01 2.7 2.83e-02 2.83e-02f 1
55r 7.6142060e+02 1.33e+04 3.03e+03 1.3 6.73e+00 2.2 3.53e-01 7.99e-02f 1
56r 7.6138592e+02 1.30e+04 2.27e+03 1.3 2.04e+01 1.7 9.65e-02 3.60e-02f 1
57r 7.6121905e+02 1.15e+04 2.16e+03 1.3 6.12e+01 1.2 9.91e-02 6.29e-02f 1
58r 7.6096313e+02 9.38e+03 2.09e+03 1.3 1.85e+02 0.8 1.32e-01 3.01e-02f 1
59r 7.5981863e+02 4.77e+03 2.05e+03 1.3 5.57e+02 0.3 5.61e-02 2.24e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
60r 7.5929151e+02 4.11e+03 2.02e+03 1.3 2.09e+02 0.7 5.35e-03 1.53e-02f 1
61r 7.5833844e+02 4.11e+03 1.98e+03 1.3 6.26e+02 0.2 1.23e-02 1.84e-02f 1
62r 7.5787962e+02 4.11e+03 1.97e+03 1.3 1.86e+03 -0.3 3.45e-01 6.90e-03f 1
63r 7.5787192e+02 4.09e+03 1.97e+03 1.3 1.30e+07 - 3.11e-04 9.43e-05f 1
64r 7.5700780e+02 4.08e+03 1.94e+03 1.3 3.27e+03 -0.7 2.00e-02 1.09e-02f 1
65r 7.5644538e+02 4.07e+03 1.93e+03 1.3 4.91e+03 -1.2 4.03e-02 6.13e-03f 1
66r 7.5437044e+02 4.08e+03 1.89e+03 1.3 1.74e+03 -1.7 8.47e-02 1.93e-02f 1
67r 7.4796816e+02 7.73e+03 2.34e+03 1.3 3.22e+03 -2.2 1.87e-01 5.20e-02f 1
68r 7.4736690e+02 7.73e+03 2.37e+03 1.3 2.63e+03 -2.6 5.87e-02 4.62e-03f 1
69r 7.4301823e+02 7.53e+03 2.29e+03 1.3 1.36e+03 -3.1 4.41e-02 3.20e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
70r 7.4234555e+02 7.53e+03 2.28e+03 1.3 2.09e+06 - 1.38e-02 4.97e-03f 1
71r 7.1749746e+02 4.76e+04 1.86e+03 1.3 2.08e+06 - 1.52e-01 1.84e-01f 1
72r 6.8042715e+02 1.52e+05 1.44e+03 1.3 1.64e+06 - 4.27e-01 3.39e-01f 1
73r 6.7622814e+02 1.46e+05 1.31e+03 1.3 1.02e+06 - 7.64e-01 5.80e-02f 1
74r 6.7540024e+02 1.44e+05 1.27e+03 1.3 9.83e+05 - 1.00e+00 1.22e-02f 1
75r 6.1248639e+02 1.08e+04 7.91e+01 1.3 9.75e+05 - 1.00e+00 9.37e-01f 1
76r 5.7837595e+02 4.24e+03 5.18e+01 -0.1 5.84e+04 - 8.66e-01 8.82e-01f 1
77r 5.6465722e+02 4.24e+03 1.65e+02 -0.1 9.07e+06 - 2.07e-01 2.44e-02f 1
78r 5.1578952e+02 4.21e+03 2.23e+02 -0.1 8.91e+06 - 5.63e-01 8.91e-02f 1
79r 2.7178805e+02 4.10e+03 1.37e+02 -0.1 8.12e+06 - 6.14e-01 4.90e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
80r 2.4566492e+02 4.09e+03 4.75e+02 -0.1 4.14e+06 - 1.00e+00 9.73e-02f 1
81r 2.3466286e+02 4.08e+03 9.99e+02 -0.1 3.73e+06 - 5.94e-01 7.84e-02f 1
82r 2.0670237e+02 4.06e+03 2.98e+02 -0.1 3.44e+06 - 4.07e-02 2.67e-01f 1
83r 1.2649365e+02 3.99e+03 1.63e+02 -0.1 2.53e+06 - 6.15e-01 1.00e+00f 1
84r 1.0342954e+02 3.99e+03 6.50e+01 -0.1 5.79e+02 - 1.00e+00 1.00e+00f 1
85r 1.0148237e+02 3.99e+03 1.67e+00 -0.1 3.28e+01 - 1.00e+00 1.00e+00h 1
86r 8.7937703e+01 3.99e+03 2.37e+01 -1.5 3.47e+02 - 1.00e+00 8.93e-01f 1
87r-8.6514385e+01 6.13e+03 8.98e+02 -1.5 4.60e+07 - 7.42e-02 1.09e-01f 1
88r-8.7388407e+01 5.80e+03 1.21e+03 -1.5 8.03e-01 2.7 1.20e-02 5.44e-02f 1
89r-8.7388675e+01 5.62e+03 1.20e+03 -1.5 1.34e-01 4.1 6.50e-02 3.03e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
90r-8.7388552e+01 5.62e+03 2.23e+03 -1.5 1.32e-01 4.5 3.83e-02 4.89e-04f 1
91r-8.7339716e+01 5.21e+03 1.24e+03 -1.5 3.32e-01 4.0 3.77e-02 7.33e-02f 1
92r-8.7087324e+01 4.48e+03 1.76e+03 -1.5 4.54e-01 3.5 1.69e-03 1.39e-01f 1
93r-8.6964896e+01 4.48e+03 1.92e+03 -1.5 4.08e+07 - 2.48e-01 5.17e-04f 1
94r-9.2263881e+01 8.88e+03 1.60e+03 -1.5 4.08e+07 - 4.65e-03 1.20e-01f 1
95r-9.0512689e+01 1.57e+04 1.95e+03 -1.5 3.60e+07 - 1.00e+00 1.73e-01f 1
96r-6.8623710e+01 7.69e+04 7.96e+02 -1.5 2.98e+07 - 2.10e-01 5.93e-01f 1
97r-3.2989414e+01 4.13e+04 7.95e+02 -1.5 1.22e+07 - 1.00e+00 8.15e-01f 1
98r-3.2987542e+01 4.10e+04 7.88e+02 -1.5 9.16e-01 3.1 9.05e-03 8.82e-03f 1
99r-3.2984009e+01 4.07e+04 7.83e+02 -1.5 9.07e-01 2.6 1.24e-01 5.66e-03f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
100r-3.1388951e+01 5.86e+03 4.75e+02 -1.5 9.01e-01 2.1 3.43e-01 8.56e-01f 1
101r-3.1337054e+01 5.80e+03 9.62e+02 -1.5 4.30e-01 1.6 1.00e+00 9.27e-03f 1
102r-1.3614868e+01 2.74e+03 2.95e+02 -1.5 1.01e+00 1.1 9.92e-01 8.79e-01f 1
103r-1.3611671e+01 2.74e+03 3.01e+02 -1.5 1.69e-02 4.3 3.23e-01 2.17e-01f 1
104r-1.3569844e+01 2.74e+03 5.95e+02 -1.5 1.07e-01 3.8 2.05e-02 3.08e-01f 1
105r-1.3566892e+01 2.74e+03 1.02e+03 -1.5 2.93e-02 4.2 1.00e+00 2.09e-01f 1
106r-4.9865437e+00 2.74e+03 6.65e+02 -1.5 2.27e+06 - 2.31e-01 6.83e-03f 1
107r-4.2921458e+00 2.73e+03 8.35e+02 -1.5 2.26e+06 - 4.15e-01 6.70e-02f 1
108r-2.2306434e+00 2.68e+03 6.56e+02 -1.5 2.11e+06 - 9.94e-01 8.80e-01f 1
109r 2.2462867e+01 2.68e+03 6.04e+02 -1.5 2.58e+05 - 1.00e+00 2.31e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
110r 6.3028266e+01 2.68e+03 6.43e+02 -1.5 1.98e+05 - 8.41e-01 7.71e-02f 1
111r 6.7784453e+01 2.68e+03 1.32e+03 -1.5 1.83e+05 - 2.61e-01 4.75e-03f 1
112r 2.5388243e+02 2.67e+03 7.80e+02 -1.5 1.82e+05 - 2.49e-01 3.14e-01f 1
113r 2.5388338e+02 2.67e+03 7.71e+02 -1.5 1.25e+05 - 1.00e+00 1.17e-02f 1
114r 2.3081199e+02 2.67e+03 1.70e+02 -1.5 1.24e+05 - 1.00e+00 8.06e-01f 1
115r 2.0902687e+02 2.67e+03 1.73e+02 -1.5 2.41e+04 - 1.81e-01 1.00e+00f 1
116r 2.0903090e+02 2.67e+03 3.24e-05 -1.5 2.25e+00 - 1.00e+00 1.00e+00h 1
117r 2.0903090e+02 2.67e+03 1.72e+02 -3.4 3.70e+02 - 1.00e+00 1.59e-06h 20
118r 2.1611707e+02 1.38e+05 1.68e+02 -3.4 6.50e+08 - 4.93e-04 1.98e-02f 1
119r 2.1611727e+02 1.36e+05 3.57e+02 -3.4 3.09e+00 3.8 5.25e-03 1.17e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
120r 2.1612015e+02 1.36e+05 3.57e+02 -3.4 7.02e+08 - 1.82e-04 1.33e-05f 1
121r 2.1611994e+02 1.36e+05 3.57e+02 -3.4 7.03e+08 - 8.80e-05 5.07e-07f 1
122r 2.0883474e+02 3.11e+05 1.25e+03 -3.4 7.03e+08 - 1.88e-07 1.48e-02f 1
123r 2.0162667e+02 3.15e+05 1.25e+03 -3.4 7.41e+08 - 1.79e-03 1.91e-03f 1
124r 2.0162667e+02 3.15e+05 9.99e+02 3.3 0.00e+00 - 0.00e+00 8.04e-08R 2
125r 2.0171952e+02 3.14e+05 9.99e+02 3.3 1.97e+06 - 4.54e-10 4.40e-06f 1
126r 2.0158655e+02 2.62e+05 9.98e+02 1.9 1.81e+05 - 5.10e-02 6.55e-04f 1
127r 2.0227328e+02 1.09e+04 9.94e+02 1.9 8.87e+05 - 6.52e-01 5.42e-02f 1
128r 2.0235761e+02 8.94e+03 7.08e+02 1.9 7.62e+00 2.0 8.72e-01 4.74e-01f 1
129r 1.9831880e+02 7.23e+03 5.73e+02 1.9 8.38e+05 - 1.91e-01 1.47e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
130r 1.8871108e+02 4.20e+03 4.14e+02 1.2 7.15e+05 - 4.47e-01 2.49e-01f 1
131r 1.4809321e+02 1.97e+03 1.63e+03 1.2 1.63e+06 - 5.61e-01 6.45e-01f 1
132r 1.2879449e+02 1.96e+03 5.99e+02 1.2 5.81e+05 - 1.00e+00 6.33e-01f 1
133r 1.3858381e+02 1.96e+03 8.01e+02 0.5 2.13e+05 - 8.08e-01 5.34e-01f 1
134r 1.3861184e+02 1.96e+03 1.04e+03 0.5 9.02e-01 3.3 1.02e-01 1.57e-01f 1
135r 1.4897969e+02 1.92e+03 6.56e+02 0.5 2.56e+06 - 3.79e-01 3.51e-01f 1
136r 1.6482847e+02 3.68e+03 2.49e+02 0.5 1.66e+06 - 1.00e+00 7.82e-01f 1
137r 1.7046320e+02 1.86e+03 4.32e+01 0.5 3.62e+05 - 1.00e+00 1.00e+00f 1
138r 1.6746132e+02 1.86e+03 1.04e+01 -0.2 3.21e+02 - 1.00e+00 1.00e+00f 1
139r 1.3750692e+02 3.89e+04 2.50e+02 -0.9 5.57e+06 - 4.11e-01 6.89e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
140 1.6438414e+02 3.89e+04 4.19e+05 -1.0 3.39e+08 - 5.27e-06 1.50e-05h 4
141 1.9515055e+02 3.89e+04 9.88e+05 -1.0 1.07e+09 - 1.27e-06 1.56e-05h 3
142r 1.9515055e+02 3.89e+04 9.99e+02 3.2 0.00e+00 - 0.00e+00 2.79e-07R 8
143r 1.9487805e+02 3.82e+04 2.85e+03 3.2 1.98e+05 - 1.42e-06 4.74e-04f 1
144r 1.9488971e+02 3.82e+04 7.82e+04 3.2 1.35e+01 4.0 2.53e-01 3.64e-02f 1
145r 1.9852595e+02 6.35e+03 7.60e+03 3.2 1.89e+05 - 9.90e-01 2.74e-02f 1
146r 2.3796828e+02 4.50e+03 1.06e+04 3.2 1.84e+05 - 9.90e-01 3.93e-01f 1
147r 2.1543925e+02 1.69e+03 2.33e+03 2.5 1.12e+05 - 8.47e-01 1.00e+00f 1
148r 1.8555939e+02 1.69e+03 1.05e+03 1.8 2.35e+05 - 6.31e-01 6.80e-01f 1
149r 1.8556873e+02 1.69e+03 1.68e+03 1.1 4.84e+00 3.5 6.22e-02 6.16e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
150r 1.8556750e+02 1.69e+03 6.47e+03 1.1 1.12e-01 4.9 3.87e-01 1.00e+00f 1
151r 1.8556410e+02 1.69e+03 3.44e+03 1.1 7.40e-02 4.4 9.10e-01 5.65e-01f 1
152r 1.8554305e+02 1.69e+03 1.12e+03 1.1 1.13e-01 3.9 1.00e+00 7.89e-01f 1
153r 1.8547767e+02 1.69e+03 8.93e+02 1.1 2.66e-01 3.4 1.00e+00 3.63e-01f 1
154r 1.8531034e+02 1.69e+03 1.02e+03 1.1 6.07e-01 2.9 3.98e-01 4.35e-01f 1
155r 1.8529660e+02 1.69e+03 1.08e+03 1.1 5.64e-02 4.3 1.00e+00 1.00e+00f 1
156r 1.8526282e+02 1.69e+03 4.87e+02 1.1 7.44e-02 3.8 1.00e+00 8.86e-01f 1
157r 1.8515897e+02 1.69e+03 4.38e+02 0.4 2.09e-01 3.3 7.75e-01 6.38e-01f 1
158r 1.8508838e+02 1.69e+03 8.34e+02 0.4 1.36e+00 2.8 5.48e-01 1.46e-01f 1
159r 1.8492555e+02 1.69e+03 1.75e+02 0.4 1.05e-01 3.3 1.00e+00 8.94e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
160r 1.8480344e+02 1.69e+03 3.83e+02 0.4 3.42e-01 2.8 1.00e+00 2.24e-01f 1
161r 1.8131255e+02 1.69e+03 4.34e+02 0.4 4.32e+06 - 3.78e-02 4.99e-03f 1
162r 1.5746927e+02 1.69e+03 1.08e+03 0.4 4.38e+06 - 4.52e-01 4.95e-02f 1
163r 1.5661261e+02 1.69e+03 1.54e+02 0.4 2.82e-01 2.3 8.87e-01 7.63e-01f 1
164r 1.5510184e+02 1.69e+03 1.62e+02 -0.3 8.14e-01 1.8 7.89e-01 3.60e-01f 1
165r 1.5501498e+02 1.69e+03 1.77e+02 -0.3 2.76e+00 1.4 1.57e-01 1.05e-01f 1
166r 1.5508488e+02 1.69e+03 5.43e+02 -0.3 2.48e+01 0.9 9.78e-03 1.47e-02f 1
167r 1.5506405e+02 1.69e+03 2.15e+03 -0.3 2.44e+00 3.1 1.51e-01 8.88e-02f 1
168r 1.5505875e+02 1.69e+03 1.44e+03 -0.3 4.14e-02 4.4 2.67e-01 6.41e-01f 1
169r 1.5505062e+02 1.69e+03 2.34e+03 -0.3 2.53e-01 4.0 1.00e+00 4.46e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
170r 1.5504113e+02 1.69e+03 2.07e+03 -0.3 8.89e-02 3.5 9.75e-01 1.24e-01f 1
171r 1.5481416e+02 1.69e+03 1.16e+03 -0.3 1.57e-01 3.0 1.00e+00 1.00e+00f 1
172r 1.5480184e+02 1.69e+03 1.23e+03 -0.3 5.53e-02 4.3 4.51e-01 1.00e+00h 1
173r 1.5476754e+02 1.69e+03 3.80e+02 -0.3 2.71e-02 3.9 1.00e+00 1.00e+00h 1
174r 1.5466767e+02 1.69e+03 6.74e+01 -0.3 2.77e-02 3.4 1.00e+00 1.00e+00h 1
175r 1.5438143e+02 1.69e+03 3.20e+01 -0.3 3.94e-02 2.9 1.00e+00 1.00e+00f 1
176r 1.5358847e+02 1.69e+03 3.16e+01 -0.3 1.17e-01 2.4 1.00e+00 1.00e+00f 1
177r 1.5325825e+02 1.69e+03 5.91e+01 -0.3 3.51e-01 2.0 2.75e-01 2.32e-01f 1
178r 1.4565100e+02 1.68e+03 6.45e+01 -0.3 1.02e+07 - 2.62e-02 1.93e-02f 1
179r 1.4562743e+02 1.68e+03 4.63e+02 -0.3 9.18e-02 3.3 1.00e+00 2.50e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
180r 1.4854028e+02 2.59e+03 1.12e+03 -0.3 9.99e+06 - 6.46e-01 8.97e-02f 1
181r 1.5591566e+02 2.21e+04 1.01e+03 -0.3 9.09e+06 - 1.00e+00 2.78e-01f 1
182r 1.5591566e+02 2.21e+04 9.99e+02 3.2 0.00e+00 - 0.00e+00 2.84e-07R 9
183r 1.5676562e+02 1.89e+04 9.93e+02 3.2 1.99e+05 - 2.87e-03 6.23e-03f 1
184r 1.5746405e+02 1.67e+04 4.43e+03 3.2 1.98e+05 - 5.50e-02 5.18e-03f 1
185r 1.7370183e+02 1.89e+03 3.11e+03 3.2 1.97e+05 - 4.80e-01 1.02e-01f 1
186r 2.5799899e+02 1.52e+03 4.57e+03 3.2 1.76e+05 - 9.90e-01 7.65e-01f 1
187r 2.5813424e+02 1.52e+03 5.30e+02 3.2 1.49e+00 2.0 9.90e-01 1.00e+00f 1
188r 2.5706230e+02 1.52e+03 3.02e+03 2.5 1.46e+00 2.4 7.68e-01 7.35e-01f 1
189r 2.5420234e+02 1.52e+03 3.29e+03 1.8 2.07e+00 1.9 8.14e-01 7.09e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
190r 2.5421182e+02 1.52e+03 5.65e+03 1.8 2.40e-01 4.2 8.60e-01 8.50e-01f 1
191r 2.5417255e+02 1.52e+03 1.01e+03 1.8 2.49e-01 3.7 1.00e+00 1.00e+00f 1
192r 2.5417639e+02 1.52e+03 7.80e+03 1.1 1.14e+00 3.2 9.17e-01 5.16e-01f 1
193r 2.5417466e+02 1.52e+03 6.80e+03 1.1 4.25e-01 3.7 1.21e-01 1.24e-01f 1
194r 2.5417338e+02 1.52e+03 6.32e+03 1.1 8.53e-02 4.1 5.15e-01 7.32e-02f 1
195r 2.5424193e+02 1.52e+03 5.43e+03 1.1 2.10e+00 3.6 1.51e-03 1.39e-01f 1
196r 2.5418876e+02 1.52e+03 3.24e+03 1.1 2.07e+00 3.1 6.79e-01 3.77e-01f 1
197r 2.5408560e+02 1.52e+03 3.04e+03 1.1 8.49e-01 2.6 1.28e-01 6.12e-02f 1
198r 2.5367568e+02 1.52e+03 2.31e+03 1.1 1.55e+00 2.2 9.71e-01 2.27e-01f 1
199r 2.5349535e+02 1.52e+03 1.38e+03 1.1 1.44e-01 3.5 4.27e-01 8.98e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
200r 2.5326237e+02 1.52e+03 8.92e+02 1.1 2.24e-01 3.0 1.00e+00 7.21e-01f 1
201r 2.4222300e+02 1.52e+03 8.78e+02 1.1 1.99e+06 - 1.68e-01 1.65e-02f 1
202r 1.7332765e+02 1.51e+03 7.35e+02 1.1 2.03e+06 - 7.71e-01 1.61e-01f 1
203r 1.3500699e+02 1.48e+03 8.71e+02 1.1 1.71e+06 - 4.58e-01 3.71e-01f 1
204r 1.3500552e+02 1.48e+03 7.47e+02 1.1 3.32e-02 4.4 4.30e-01 6.10e-01f 1
205r 1.3498670e+02 1.48e+03 9.19e+02 1.1 5.87e-02 3.9 1.00e+00 1.00e+00f 1
206r 1.1508907e+02 1.46e+03 5.98e+02 1.1 1.08e+06 - 1.00e+00 5.79e-01f 1
207r 1.0864484e+02 1.45e+03 6.94e+02 0.4 4.53e+05 - 7.67e-01 5.99e-01f 1
208r 5.8666979e+01 2.89e+04 3.82e+02 0.4 2.96e+06 - 1.00e+00 1.00e+00f 1
209r 5.8666979e+01 2.89e+04 9.99e+02 3.1 0.00e+00 - 0.00e+00 4.68e-07R 8
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
210r 4.5915157e+01 2.62e+04 1.74e+03 3.1 2.12e+05 - 7.39e-05 2.70e-03f 1
211r 4.5927127e+01 2.61e+04 4.63e+04 3.1 6.11e+00 4.0 2.80e-01 1.22e-01f 1
212r 4.6045710e+01 2.60e+04 1.58e+04 3.1 1.22e+00 3.5 9.90e-01 4.01e-01f 1
213r 4.8808171e+01 4.29e+03 1.66e+04 3.1 2.11e+05 - 9.90e-01 3.15e-02f 1
214r 5.7159080e+01 3.92e+03 1.51e+04 2.4 2.05e+05 - 1.10e-01 8.65e-02f 1
215r 5.7233093e+01 3.75e+03 8.17e+03 2.4 1.94e+00 3.0 5.34e-01 7.52e-01f 1
216r 5.7234017e+01 3.74e+03 3.28e+03 1.7 1.43e-01 4.4 7.20e-01 1.00e+00f 1
217r 5.7236656e+01 3.73e+03 1.21e+03 1.7 1.53e-01 3.9 1.00e+00 1.00e+00f 1
218r 5.7247265e+01 3.70e+03 1.02e+03 1.7 3.85e-01 3.4 1.00e+00 9.52e-01f 1
219r 5.7263482e+01 3.60e+03 1.05e+03 1.7 1.95e+00 2.9 3.33e-01 2.56e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
220r 5.7250241e+01 3.31e+03 2.08e+03 1.7 3.37e+00 2.5 5.83e-01 2.15e-01f 1
221r 5.7253012e+01 3.28e+03 9.91e+02 1.7 3.77e+00 2.9 5.52e-02 5.77e-02f 1
222r 5.7251042e+01 3.25e+03 1.98e+03 1.0 7.67e-01 3.3 7.09e-01 1.76e-01f 1
223r 5.6572872e+01 1.32e+03 1.97e+03 1.0 2.35e+06 - 2.41e-01 3.35e-03f 1
224r 2.1999395e+01 1.29e+03 7.16e+02 1.0 2.34e+06 - 3.60e-01 2.48e-01f 1
225r 2.7181363e+01 1.42e+03 8.58e+02 1.0 1.77e+06 - 5.20e-01 2.31e-01f 1
226r 2.7218276e+01 1.27e+03 8.71e+02 1.0 3.03e-01 2.8 6.08e-01 1.41e-01f 1
227r 8.8681085e+01 2.44e+03 5.29e+02 1.0 1.36e+06 - 1.00e+00 5.52e-01f 1
228r 8.8685164e+01 1.24e+03 8.21e+02 0.3 7.31e-02 4.2 6.19e-01 7.03e-01f 1
229r 8.8689064e+01 1.24e+03 1.12e+03 0.3 2.86e-02 4.6 1.00e+00 9.86e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
230r 8.8706586e+01 1.24e+03 3.29e+02 0.3 2.48e-02 4.1 1.00e+00 9.83e-01f 1
231r 8.8738787e+01 1.24e+03 3.44e+02 0.3 9.33e-02 3.6 1.00e+00 6.15e-01f 1
232r 8.8773941e+01 1.24e+03 8.99e+02 0.3 1.94e+00 3.2 4.06e-02 1.42e-01f 1
233r 8.8958080e+01 1.24e+03 6.12e+02 0.3 4.06e-01 2.7 1.00e+00 3.79e-01f 1
234r 9.3561426e+01 1.22e+03 8.18e+02 0.3 3.58e+06 - 3.38e-01 1.33e-01f 1
235r 9.4608367e+01 1.22e+03 5.20e+02 0.3 8.75e-02 2.2 1.00e+00 6.92e-01f 1
236r 9.1204895e+01 1.28e+03 4.46e+02 0.3 3.10e+06 - 3.96e-01 1.67e-01f 1
237r 8.9433830e+01 1.20e+03 1.03e+03 0.3 2.49e-01 1.7 1.00e+00 3.54e-01f 1
238r 6.0355231e+01 2.64e+04 1.67e+02 0.3 2.58e+06 - 1.00e+00 1.00e+00f 1
239r 6.0355231e+01 2.64e+04 9.99e+02 3.0 0.00e+00 - 0.00e+00 3.09e-07R 8
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
240r 5.8433269e+01 2.17e+04 9.95e+02 3.0 2.32e+05 - 7.41e-04 4.47e-03f 1
241r 5.8055802e+01 1.93e+04 2.93e+03 3.0 2.32e+05 - 2.37e-02 2.90e-03f 1
242r 6.7118172e+01 2.51e+03 2.80e+03 3.0 2.31e+05 - 4.64e-02 4.99e-02f 1
243r 6.7170658e+01 2.45e+03 2.33e+03 3.0 2.36e-01 4.0 9.90e-01 1.00e+00f 1
244r 6.7316658e+01 2.30e+03 3.13e+03 3.0 2.49e-01 3.5 9.90e-01 1.00e+00f 1
245r 6.7357447e+01 2.05e+03 9.39e+02 2.3 8.49e-01 3.0 9.64e-01 1.00e+00f 1
246r 6.7355560e+01 2.01e+03 9.93e+02 0.9 3.36e-01 3.5 9.72e-01 1.00e+00f 1
247r 6.7373722e+01 1.90e+03 9.96e+02 0.9 1.01e+00 3.0 1.00e+00 4.75e-01f 1
248r 6.7362124e+01 1.85e+03 9.95e+02 0.9 3.77e-01 3.4 1.00e+00 4.34e-01f 1
249r 6.7330252e+01 1.69e+03 9.92e+02 0.9 1.12e+00 2.9 5.08e-01 3.67e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
250r 6.7256421e+01 1.57e+03 1.15e+03 0.9 3.37e+00 2.5 5.51e-01 8.56e-02f 1
251r 6.7425887e+01 1.08e+03 1.17e+03 0.9 2.60e+06 - 3.71e-02 3.12e-03f 1
252r 6.9106423e+01 1.07e+03 1.34e+03 0.9 2.59e+06 - 3.85e-01 4.09e-02f 1
253r 6.6844859e+01 6.14e+03 1.05e+03 0.9 2.48e+06 - 1.00e+00 4.84e-01f 1
254r 5.0945864e+01 6.76e+03 1.04e+02 0.9 1.28e+06 - 1.00e+00 9.84e-01f 1
255 4.7219269e+01 6.76e+03 2.47e+02 -1.0 9.34e+10 - 6.17e-07 1.25e-07f 1
256r 4.7219269e+01 6.76e+03 9.99e+02 3.0 0.00e+00 - 0.00e+00 2.68e-07R 8
257r 4.1631273e+01 6.08e+03 1.96e+03 3.0 2.46e+05 - 1.36e-03 8.62e-03f 1
258r 4.1633804e+01 6.08e+03 4.02e+04 3.0 7.86e+00 4.0 2.38e-01 6.24e-02f 1
259r 4.1709784e+01 5.93e+03 8.33e+03 3.0 8.46e-01 3.5 9.90e-01 5.79e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
260r 4.1715839e+01 5.92e+03 6.10e+04 3.0 6.37e-01 4.9 9.91e-01 1.00e+00h 1
261r 4.1716364e+01 5.92e+03 3.00e+04 2.3 8.17e-02 5.3 1.00e+00 1.00e+00f 1
262r 4.1716549e+01 5.92e+03 2.42e+05 2.3 6.25e-02 5.7 1.00e+00 1.00e+00f 1
263r 4.1717104e+01 5.91e+03 2.50e+04 2.3 5.82e-02 5.2 1.00e+00 1.00e+00f 1
264r 4.1718770e+01 5.91e+03 3.62e+03 2.3 3.38e-02 4.7 1.00e+00 1.00e+00f 1
265r 4.1723745e+01 5.90e+03 1.03e+03 2.3 5.50e-02 4.3 1.00e+00 1.00e+00f 1
266r 4.1727097e+01 5.88e+03 9.97e+02 1.6 1.60e-01 3.8 1.00e+00 1.00e+00f 1
267r 4.1732790e+01 5.83e+03 9.98e+02 1.6 4.80e-01 3.3 1.00e+00 6.73e-01f 1
268r 4.1736133e+01 5.67e+03 2.12e+03 1.6 2.69e+00 2.8 6.42e-01 2.33e-01f 1
269r 4.1736600e+01 5.65e+03 1.21e+03 1.6 9.10e-02 4.2 1.00e+00 6.09e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
270r 4.1736855e+01 5.64e+03 2.00e+03 1.6 5.07e-02 4.6 1.00e+00 1.00e+00f 1
271r 4.1736933e+01 5.60e+03 1.01e+03 1.6 7.66e-02 4.1 1.00e+00 1.00e+00f 1
272r 4.1736640e+01 5.58e+03 1.01e+03 1.6 2.88e-02 4.5 1.00e+00 1.00e+00f 1
273r 4.1733694e+01 5.55e+03 1.27e+03 0.9 8.81e-02 4.1 1.00e+00 9.63e-01f 1
274r 4.1723038e+01 5.49e+03 1.01e+03 0.9 2.61e-01 3.6 1.00e+00 3.63e-01f 1
275r 4.1679343e+01 5.30e+03 1.02e+03 0.9 7.83e-01 3.1 5.06e-01 4.00e-01f 1
276r 4.1494365e+01 5.16e+03 1.02e+03 0.9 2.44e+00 2.6 1.65e-01 8.94e-02f 1
277r 4.1476976e+01 5.11e+03 1.02e+03 0.9 1.10e-01 4.0 1.00e+00 8.62e-01f 1
278r 4.1428492e+01 4.97e+03 1.02e+03 0.9 3.30e-01 3.5 1.00e+00 6.67e-01f 1
279r 4.1399790e+01 9.53e+02 9.69e+02 0.9 2.76e+06 - 4.78e-02 9.49e-04f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
280r 4.0013413e+01 9.46e+02 1.19e+03 0.9 2.76e+06 - 3.34e-01 4.87e-02f 1
281r 3.8764220e+01 9.29e+02 6.59e+02 0.9 2.63e+06 - 4.14e-01 1.34e-01f 1
282r 3.8276135e+01 1.47e+04 2.94e+02 0.9 2.28e+06 - 1.00e+00 7.77e-01f 1
283r 3.8276135e+01 1.47e+04 9.99e+02 2.9 0.00e+00 - 0.00e+00 2.74e-07R 6
284r 3.1484570e+01 6.02e+03 2.99e+03 2.9 2.61e+05 - 3.98e-04 2.15e-02f 1
285r 3.1507280e+01 5.93e+03 3.24e+04 2.9 7.45e-01 4.0 9.90e-01 1.00e+00f 1
286r 3.1520470e+01 5.89e+03 9.24e+03 2.9 2.11e-01 4.4 9.90e-01 1.00e+00f 1
287r 3.1521709e+01 5.89e+03 7.69e+03 2.2 1.08e-01 4.9 9.18e-01 1.00e+00f 1
288r 3.1524540e+01 5.88e+03 9.85e+02 2.2 4.16e-02 4.4 1.00e+00 1.00e+00f 1
289r 3.1526555e+01 5.87e+03 9.98e+02 1.5 1.26e-01 3.9 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
290r 3.1530718e+01 5.82e+03 1.00e+03 1.5 3.81e-01 3.4 1.00e+00 8.70e-01f 1
291r 3.1531035e+01 5.71e+03 1.00e+03 1.5 1.25e+00 2.9 1.00e+00 2.09e-01f 1
292r 3.1508931e+01 5.18e+03 1.00e+03 1.5 3.45e+00 2.5 1.00e+00 2.48e-01f 1
293r 3.1478843e+01 4.94e+03 1.01e+03 0.8 1.30e+00 2.9 8.03e-01 3.03e-01f 1
294r 3.1347956e+01 4.13e+03 2.11e+03 0.8 3.90e+00 2.4 5.45e-01 3.33e-01f 1
295r 2.7314698e+01 8.41e+02 2.10e+03 0.8 2.92e+06 - 2.15e-02 2.09e-03f 1
296r 6.9972435e+00 8.38e+02 2.04e+03 0.8 2.92e+06 - 1.90e-01 2.38e-02f 1
297r 2.9508778e+00 1.73e+03 9.42e+02 0.8 2.85e+06 - 6.43e-01 2.07e-01f 1
298r 3.0990394e+00 1.14e+03 7.23e+02 0.8 4.35e-01 2.8 5.42e-01 3.39e-01f 1
299r 3.1099763e+00 8.24e+02 1.32e+04 0.8 3.03e-01 4.2 2.78e-01 2.78e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
300r 7.2219236e+00 8.20e+02 1.28e+04 0.8 2.26e+06 - 8.07e-02 2.72e-02f 1
301r 2.6983021e+01 1.84e+03 9.87e+03 0.8 2.20e+06 - 6.29e-02 2.26e-01f 1
302r 4.6771030e+01 4.87e+03 3.78e+03 0.8 1.70e+06 - 9.72e-01 5.28e-01f 1
303r 4.6771030e+01 4.87e+03 9.99e+02 2.9 0.00e+00 - 0.00e+00 3.89e-07R 5
304r 4.1925062e+01 3.67e+03 1.35e+03 2.9 2.79e+05 - 4.10e-04 2.28e-02f 1
305r 4.1927166e+01 3.65e+03 1.11e+04 2.9 4.31e+00 4.0 2.71e-01 2.51e-01f 1
306r 4.1975136e+01 3.47e+03 2.79e+04 2.9 1.83e+00 3.5 3.93e-01 7.44e-01f 1
307r 4.1975709e+01 3.46e+03 4.03e+04 2.9 4.57e-01 4.9 4.88e-01 1.00e+00f 1
308r 4.1985819e+01 3.42e+03 3.22e+03 2.9 1.33e-01 4.4 9.67e-01 1.00e+00f 1
309r 5.0289660e+01 7.31e+02 1.26e+04 2.9 2.72e+05 - 9.92e-01 1.36e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
310r 5.7465361e+01 7.20e+02 3.17e+03 2.2 2.35e+05 - 9.73e-01 7.16e-01f 1
311r 1.1737769e+02 7.07e+02 8.24e+02 1.5 4.12e+05 - 5.71e-01 7.45e-01f 1
312r 1.2054145e+02 8.16e+02 5.09e+02 0.8 8.78e+05 - 5.44e-01 4.11e-01f 1
313r 1.2680237e+02 1.26e+03 7.13e+02 0.8 2.26e+06 - 3.68e-01 1.43e-01f 1
314r 1.1434041e+02 5.50e+03 3.95e+02 0.8 1.93e+06 - 4.06e-01 4.92e-01f 1
315r 1.1434041e+02 5.50e+03 9.99e+02 2.8 0.00e+00 - 0.00e+00 3.93e-07R 6
316r 9.6439926e+01 2.65e+03 4.46e+03 2.8 3.02e+05 - 2.95e-04 3.61e-02f 1
317r 9.6471959e+01 2.56e+03 5.77e+04 2.8 5.55e-01 4.0 9.90e-01 1.00e+00f 1
318r 9.6475792e+01 2.52e+03 2.93e+04 2.8 2.91e-01 4.4 9.91e-01 1.00e+00f 1
319r 9.6476764e+01 2.52e+03 1.26e+04 2.1 1.67e-01 4.9 8.49e-01 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
320r 9.6478071e+01 2.51e+03 9.51e+04 2.1 7.83e-01 4.4 4.98e-01 6.45e-01f 1
321r 9.6478836e+01 2.51e+03 5.47e+04 2.1 4.21e-01 4.8 2.46e-01 1.00e+00f 1
322r 9.6479122e+01 2.51e+03 7.09e+04 2.1 1.17e-01 5.2 1.00e+00 1.00e+00f 1
323r 9.6479979e+01 2.51e+03 1.06e+04 2.1 4.90e-02 4.7 1.00e+00 1.00e+00f 1
324r 9.6482520e+01 2.49e+03 9.74e+02 2.1 5.20e-02 4.3 1.00e+00 1.00e+00f 1
325r 9.6485605e+01 2.47e+03 9.96e+02 1.4 1.59e-01 3.8 1.00e+00 1.00e+00f 1
326r 9.6487088e+01 2.41e+03 1.00e+03 1.4 4.88e-01 3.3 1.00e+00 4.53e-01f 1
327r 9.6487316e+01 2.36e+03 1.02e+03 1.4 1.86e-01 3.7 1.00e+00 6.37e-01f 1
328r 9.6494925e+01 2.27e+03 1.00e+03 1.4 9.44e-01 3.3 8.97e-01 2.77e-01f 1
329r 9.6499781e+01 2.15e+03 1.01e+03 1.4 2.05e-01 3.7 1.00e+00 9.40e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
330r 9.6502188e+01 2.12e+03 1.70e+03 0.7 2.55e+00 3.2 1.28e-01 8.06e-02f 1
331r 9.6554091e+01 1.60e+03 1.01e+03 0.7 1.85e+00 2.7 5.30e-01 4.57e-01f 1
332r 9.6726589e+01 1.45e+03 1.01e+03 0.7 3.39e+06 - 3.09e-03 5.47e-04f 1
333r 8.8753615e+01 8.92e+02 9.91e+02 0.7 3.39e+06 - 1.70e-02 3.36e-03f 1
334r 2.2942177e+01 6.18e+02 2.32e+03 0.7 3.38e+06 - 4.32e-04 1.52e-02f 1
335r 2.2947310e+01 6.18e+02 2.20e+03 0.7 6.70e-01 3.2 9.98e-01 3.39e-01f 1
336r 2.3067829e+01 6.18e+02 7.17e+03 0.7 1.83e+00 2.7 3.60e-01 2.10e-01f 1
337r 2.3254503e+01 6.18e+02 3.98e+03 0.7 3.50e+00 2.2 1.63e-01 1.27e-01f 1
338r 2.3255065e+01 6.18e+02 5.23e+03 0.7 3.46e-01 4.4 2.85e-03 8.60e-02f 1
339r 2.3256068e+01 6.18e+02 5.11e+03 0.7 1.03e+00 4.0 2.62e-01 4.13e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
340r 2.3265613e+01 6.18e+02 4.38e+03 0.7 1.22e+00 3.5 3.27e-02 1.10e-01f 1
341r 2.3323555e+01 6.18e+02 3.78e+03 0.7 1.62e+00 3.0 3.00e-01 1.95e-01f 1
342r 2.4347101e+01 6.18e+02 3.18e+03 0.7 1.67e+00 2.5 5.15e-01 1.00e+00f 1
343r 2.4388165e+01 6.18e+02 8.12e+02 0.7 1.11e-01 3.9 1.00e+00 1.00e+00f 1
344r 2.4449170e+01 6.18e+02 5.79e+02 0.7 2.35e-01 3.4 1.00e+00 4.75e-01f 1
345r 2.4542372e+01 6.18e+02 5.76e+02 0.7 7.02e-01 2.9 1.00e+00 2.70e-01f 1
346r 3.4770269e+01 6.17e+02 5.73e+02 0.7 3.32e+06 - 6.67e-03 5.14e-03f 1
347r 3.7405662e+01 6.16e+02 5.73e+02 0.7 3.32e+06 - 5.15e-04 2.79e-03f 1
348r 3.9001307e+01 6.16e+02 5.70e+02 0.7 3.30e+06 - 5.25e-03 3.80e-03f 1
349r 7.4703598e+01 6.00e+02 8.91e+02 0.7 3.28e+06 - 3.27e-01 9.25e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
350r 7.6695318e+01 2.29e+03 1.16e+03 0.7 2.98e+06 - 1.28e-01 1.94e-01f 1
351r 8.9006155e+01 1.19e+04 7.41e+02 0.7 2.40e+06 - 1.00e+00 5.74e-01f 1
352r 8.9006155e+01 1.19e+04 9.99e+02 2.7 0.00e+00 - 0.00e+00 4.02e-07R 3
353r 8.8118573e+01 3.75e+03 1.05e+03 2.7 3.32e+05 - 1.37e-04 1.93e-02f 1
354r 9.2051516e+01 2.95e+03 1.01e+03 2.7 3.25e+05 - 1.90e-01 6.69e-03f 1
355r 1.2625069e+02 5.10e+02 1.00e+03 2.7 3.24e+05 - 1.47e-01 2.36e-01f 1
356r 1.2625237e+02 5.10e+02 1.67e+04 2.7 5.52e-01 4.0 9.91e-01 1.00e+00f 1
357r 1.2622622e+02 5.10e+02 1.75e+04 2.0 7.05e-01 3.5 8.60e-01 1.00e+00f 1
358r 1.2622647e+02 5.10e+02 5.89e+03 2.0 1.52e-01 3.9 1.00e+00 1.00e+00f 1
359r 1.2622409e+02 5.10e+02 1.73e+03 2.0 2.63e-01 3.5 1.00e+00 8.23e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
360r 1.3292351e+02 5.11e+02 1.58e+03 1.3 6.60e+05 - 6.57e-01 7.35e-02f 1
361r 2.0829848e+02 5.76e+03 6.61e+02 1.3 1.53e+06 - 6.00e-01 6.40e-01f 1
362r 2.0821094e+02 1.91e+03 5.18e+03 0.6 2.69e-01 3.0 5.36e-01 6.69e-01f 1
363r 1.3927324e+02 1.82e+04 2.42e+03 0.6 2.63e+06 - 7.46e-01 6.45e-01f 1
364r 9.4716945e+01 8.12e+03 8.29e+02 0.6 9.37e+05 - 1.00e+00 7.16e-01f 1
365r 9.3313257e+01 7.68e+03 7.84e+02 0.6 2.66e+05 - 7.69e-01 5.51e-02f 1
366r 8.4410213e+01 4.86e+03 5.09e+02 0.6 2.52e+05 - 1.00e+00 3.71e-01f 1
367r 6.8833330e+01 5.73e+02 5.94e+01 0.6 1.58e+05 - 1.00e+00 1.00e+00f 1
368r 6.5458286e+01 5.73e+02 1.20e+01 -0.1 8.01e+03 - 1.00e+00 1.00e+00f 1
369r 6.3776553e+01 5.73e+02 3.52e+02 -0.8 4.74e+06 - 9.73e-01 8.02e-03f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
370r-6.2692815e+00 2.37e+04 4.81e+02 -0.8 1.59e+07 - 2.45e-01 1.02e-01f 1
371r-6.2240948e+00 2.36e+04 3.49e+02 -0.8 1.42e+07 - 8.23e-01 1.35e-03f 1
372r-4.8160388e+01 2.79e+04 4.43e+02 -0.8 1.42e+07 - 5.39e-01 5.81e-02f 1
373r-4.8202886e+01 2.79e+04 1.45e+03 -0.8 1.33e+07 - 9.03e-01 7.05e-05f 1
374r-9.1009973e+01 3.26e+04 1.34e+03 -0.8 1.32e+07 - 1.92e-01 7.62e-02f 1
375r-9.1021411e+01 3.12e+04 1.29e+03 -0.8 1.36e+00 2.5 2.10e-02 4.35e-02f 1
376r-9.1043609e+01 3.03e+04 1.25e+03 -0.8 1.71e+00 2.0 6.81e-02 2.78e-02f 1
377r-9.1049625e+01 2.97e+04 1.22e+03 -0.8 1.28e+00 2.5 9.70e-03 2.04e-02f 1
378r-9.1514873e+01 1.54e+04 6.76e+02 -0.8 8.89e-01 2.0 6.57e-03 4.81e-01f 1
379r-9.2354063e+01 1.08e+04 4.45e+02 -0.8 1.11e+00 1.5 5.33e-01 3.01e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
380r-9.2604110e+01 1.08e+04 1.46e+03 -0.8 1.21e+07 - 3.76e-01 4.67e-04f 1
381r-1.4208569e+02 1.53e+04 1.06e+03 -0.8 1.15e+07 - 5.30e-03 8.39e-02f 1
382r-1.4138165e+02 1.46e+04 1.00e+03 -0.8 7.89e-01 1.0 8.77e-02 4.83e-02f 1
383r-1.3922320e+02 1.41e+04 9.67e+02 -0.8 3.13e+00 0.6 7.62e-02 3.46e-02f 1
384r-1.3446283e+02 1.38e+04 9.46e+02 -0.8 1.59e+01 0.1 2.17e-02 2.16e-02f 1
385r-1.3459277e+02 1.37e+04 2.15e+03 -0.8 6.90e+06 - 1.00e+00 6.17e-04f 1
386r-3.0140072e+02 9.19e+04 4.68e+02 -0.8 6.95e+06 - 1.00e+00 8.79e-01f 1
387r-3.0140072e+02 9.19e+04 9.99e+02 2.6 0.00e+00 - 0.00e+00 7.19e-08R 2
388r-3.0685875e+02 1.41e+04 9.98e+02 2.6 9.20e+06 - 4.00e-04 2.14e-03f 1
389r-3.0475030e+02 1.08e+04 1.02e+03 2.6 1.05e+06 - 1.18e-01 1.59e-03f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
390r-2.9469004e+02 1.34e+03 9.58e+02 2.6 2.25e+05 - 2.05e-02 1.20e-01f 1
391r-2.9394693e+02 4.65e+02 1.96e+03 2.6 5.01e+00 2.0 6.85e-01 2.72e-01f 1
392r-2.8234772e+02 3.96e+02 1.38e+03 2.6 1.81e+05 - 5.18e-01 2.69e-01f 1
393r-2.8752702e+02 1.13e+03 2.93e+03 2.6 1.30e+05 - 4.26e-01 3.59e-01f 1
394r-2.8751742e+02 3.94e+02 1.33e+04 2.6 2.98e-01 4.2 9.90e-01 1.00e+00h 1
395r-2.8751441e+02 3.94e+02 1.30e+04 2.6 9.91e-02 4.7 9.91e-01 1.00e+00f 1
396r-2.8751618e+02 3.94e+02 9.06e+02 1.9 5.97e-02 4.2 1.00e+00 1.00e+00f 1
397r-2.8751860e+02 3.94e+02 9.11e+02 1.2 1.80e-01 3.7 1.00e+00 1.00e+00f 1
398r-2.8750642e+02 3.94e+02 8.64e+02 1.2 5.30e-01 3.2 1.00e+00 3.80e-01f 1
399r-2.8745985e+02 3.94e+02 1.35e+03 1.2 8.73e-01 2.7 1.00e+00 2.31e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
400r-2.8720824e+02 3.94e+02 1.05e+03 1.2 9.26e-01 2.3 1.00e+00 3.48e-01f 1
401r-2.8717682e+02 3.94e+02 3.68e+02 1.2 9.60e-02 3.6 1.00e+00 1.00e+00f 1
402r-2.8712360e+02 3.94e+02 2.82e+02 0.5 1.36e-01 3.1 1.00e+00 6.49e-01f 1
403r-2.8698055e+02 3.94e+02 1.85e+02 0.5 3.04e-01 2.6 1.00e+00 5.87e-01f 1
404r-2.8672919e+02 3.94e+02 3.42e+02 0.5 4.88e-01 2.2 1.00e+00 3.52e-01f 1
405r-2.2008506e+02 3.94e+02 3.39e+02 0.5 2.98e+06 - 1.56e-02 1.43e-02f 1
406r-2.1916689e+02 3.94e+02 4.21e+02 0.5 6.15e-01 1.7 1.00e+00 2.85e-01f 1
407r-2.0852386e+02 3.94e+02 6.11e+01 0.5 5.58e-01 1.2 9.54e-01 1.00e+00f 1
408r-1.4608754e+02 3.94e+02 5.12e+02 -0.2 5.20e+00 0.7 2.39e-01 6.11e-01f 1
409r-1.4606864e+02 3.94e+02 4.27e+02 -0.2 7.36e-02 3.9 1.00e+00 7.75e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
410r-1.4590666e+02 3.94e+02 7.19e+02 -0.2 1.76e-01 3.4 6.18e-02 1.00e+00f 1
411r-1.4583401e+02 3.94e+02 4.83e+02 -0.2 8.22e-01 2.9 1.00e+00 2.15e-01f 1
412r-1.4218386e+02 3.93e+02 5.91e+02 -0.2 5.64e+06 - 2.59e-01 6.22e-03f 1
413r-1.5480102e+02 3.47e+03 7.24e+02 -0.2 5.11e+06 - 6.00e-01 2.33e-01f 1
414r-1.5470254e+02 3.25e+03 6.80e+02 -0.2 8.01e-02 2.4 6.26e-02 6.49e-02h 1
415r-1.5582097e+02 3.21e+03 1.02e+03 -0.2 3.91e+06 - 1.00e+00 1.56e-02f 1
416r-1.9939207e+02 2.15e+04 3.98e+02 -0.2 3.86e+06 - 1.00e+00 7.58e-01f 1
417r-1.9939207e+02 2.15e+04 9.99e+02 2.4 0.00e+00 - 0.00e+00 3.08e-07R 4
418r-1.9862009e+02 3.05e+03 9.94e+02 2.4 8.37e+05 - 8.32e-02 5.34e-03f 1
419r-1.9354186e+02 2.28e+03 1.02e+03 2.4 2.45e+05 - 1.76e-02 3.32e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
420r-1.7944941e+02 4.10e+02 1.24e+03 2.4 2.36e+05 - 1.31e-01 1.33e-01f 1
421r-1.7934667e+02 3.49e+02 4.26e+03 2.4 2.03e+01 2.0 5.06e-02 7.33e-02f 1
422r-1.7934407e+02 3.38e+02 5.29e+03 2.4 8.06e-01 4.2 2.49e-01 3.12e-01h 1
423r-1.7934171e+02 3.26e+02 1.28e+05 2.4 1.03e+00 4.7 9.90e-01 7.49e-01f 1
424r-1.7934120e+02 3.24e+02 2.10e+05 2.4 1.52e+00 5.1 1.27e-01 4.31e-01f 1
425r-1.7934116e+02 3.24e+02 7.69e+05 2.4 3.42e-01 6.4 9.90e-01 7.33e-01h 1
426r-1.7934114e+02 3.24e+02 5.43e+05 2.4 6.94e-02 6.8 8.36e-01 6.22e-01h 1
427r-1.7934108e+02 3.23e+02 2.01e+05 2.4 4.00e-02 6.4 1.00e+00 1.00e+00f 1
428r-1.7934090e+02 3.23e+02 8.75e+03 2.4 1.14e-02 5.9 1.00e+00 1.00e+00f 1
429r-1.7934090e+02 3.22e+02 1.45e+03 1.7 5.69e-03 5.4 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
430r-1.7934067e+02 3.20e+02 1.16e+03 1.7 1.37e-02 4.9 1.00e+00 1.00e+00f 1
431r-1.7933997e+02 3.14e+02 1.07e+03 1.7 3.76e-02 4.5 1.00e+00 1.00e+00f 1
432r-1.7933784e+02 2.88e+02 1.01e+03 1.7 1.06e-01 4.0 1.00e+00 1.00e+00f 1
433r-1.7933168e+02 2.77e+02 1.27e+03 1.7 3.18e-01 3.5 1.00e+00 8.71e-01f 1
434r-1.8379638e+02 2.76e+02 1.22e+03 1.0 5.26e+05 - 1.21e-01 7.76e-02f 1
435r-1.8379595e+02 2.76e+02 9.46e+02 1.0 6.27e-01 3.0 7.87e-01 3.02e-01f 1
436r-1.8379482e+02 2.76e+02 7.86e+02 1.0 3.75e-01 3.4 1.00e+00 8.11e-01f 1
437r-1.8379043e+02 2.76e+02 1.21e+03 1.0 1.63e+00 3.0 1.54e-01 1.92e-01f 1
438r-1.8608440e+02 4.63e+02 7.71e+02 1.0 1.18e+06 - 6.08e-01 3.73e-01f 1
439r-1.8643541e+02 1.31e+03 1.21e+03 1.0 7.42e+05 - 4.42e-01 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
440r-1.8643541e+02 1.31e+03 9.99e+02 2.4 0.00e+00 - 0.00e+00 3.18e-07R 8
441r-1.9199738e+02 8.44e+02 1.35e+03 2.4 2.63e+05 - 5.72e-04 2.71e-02f 1
442r-1.9661848e+02 7.90e+02 1.34e+03 2.4 2.55e+05 - 1.91e-02 5.00e-03f 1
443r-1.7784312e+02 2.43e+02 1.20e+03 2.4 2.54e+05 - 9.63e-03 1.02e-01f 1
444r-1.7783739e+02 2.43e+02 3.48e+04 2.4 7.28e-01 4.0 6.59e-01 4.28e-01f 1
445r-1.7780178e+02 2.42e+02 1.60e+03 2.4 3.71e-01 3.5 9.90e-01 1.00e+00f 1
446r-1.7780031e+02 2.42e+02 6.09e+04 2.4 1.65e-01 4.9 9.90e-01 1.00e+00h 1
447r-1.7779972e+02 2.42e+02 1.40e+04 2.4 7.40e-02 5.3 1.00e+00 1.00e+00f 1
448r-1.7779793e+02 2.42e+02 4.06e+03 2.4 5.30e-02 4.8 1.00e+00 1.00e+00f 1
449r-1.8099730e+02 2.42e+02 3.59e+03 1.7 2.28e+05 - 7.55e-01 1.25e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
450r-1.7664600e+02 3.00e+03 4.47e+03 1.7 5.36e+05 - 1.00e+00 1.00e+00f 1
451r-1.6756663e+02 8.96e+02 2.56e+03 1.7 3.01e+03 - 6.78e-01 7.01e-01f 1
452r-1.6756662e+02 7.83e+02 2.40e+03 1.0 3.39e-01 4.3 3.15e-02 1.26e-01f 1
453r-1.6756570e+02 5.51e+02 2.08e+03 1.0 9.10e-01 3.8 2.87e-01 2.96e-01f 1
454r-1.6756544e+02 4.34e+02 1.28e+03 1.0 2.45e-01 4.3 1.25e-01 2.13e-01f 1
455r-1.6756529e+02 2.43e+02 6.49e+03 1.0 3.20e-01 4.7 7.21e-01 4.39e-01f 1
456r-1.6762391e+02 2.40e+02 6.16e+03 1.0 7.47e+05 - 4.52e-01 1.52e-02f 1
457r-1.7188854e+02 2.18e+02 3.92e+03 1.0 7.33e+05 - 1.00e+00 3.80e-01f 1
458r-1.7188854e+02 2.18e+02 9.99e+02 2.3 0.00e+00 - 0.00e+00 3.06e-07R 9
459r-1.7040756e+02 2.18e+02 2.37e+03 2.3 2.77e+05 - 1.43e-03 1.51e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
460r-1.7040541e+02 2.17e+02 2.82e+04 2.3 1.03e+00 4.0 4.17e-01 2.56e-01h 1
461r-1.7039171e+02 2.17e+02 3.22e+03 2.3 1.79e+00 3.5 3.85e-01 4.93e-01f 1
462r-1.7039153e+02 2.17e+02 8.97e+03 2.3 6.75e-01 4.9 7.44e-02 1.44e-01h 1
463r-1.7039070e+02 2.17e+02 1.13e+04 2.3 8.41e-01 4.4 5.94e-01 2.14e-01h 1
464r-1.7038266e+02 2.17e+02 1.58e+04 2.3 2.25e+00 3.9 5.32e-01 6.88e-01h 1
465r-1.7037860e+02 2.17e+02 2.18e+04 2.3 4.83e-01 4.3 9.63e-01 9.25e-01f 1
466r-1.7037695e+02 2.17e+02 5.69e+03 2.3 1.01e-01 4.7 9.90e-01 1.00e+00f 1
467r-1.7037597e+02 2.17e+02 1.12e+03 1.6 5.97e-02 4.3 9.94e-01 1.00e+00f 1
468r-1.7037275e+02 2.17e+02 1.05e+03 1.6 1.69e-01 3.8 1.00e+00 1.00e+00f 1
469r-1.7036732e+02 2.17e+02 1.07e+03 1.6 5.34e-01 3.3 1.00e+00 4.56e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
470r-1.7036475e+02 2.17e+02 1.16e+03 0.9 2.09e-01 3.7 1.00e+00 9.75e-01f 1
471r-1.7035408e+02 4.05e+02 2.51e+04 0.9 6.54e-01 3.3 4.21e-01 5.23e-01f 1
472r-1.7034453e+02 2.17e+02 1.21e+03 0.9 2.45e-01 3.7 1.00e+00 1.00e+00f 1
473r-1.7034113e+02 2.17e+02 1.73e+03 0.9 9.17e-02 4.1 1.00e+00 1.00e+00f 1
474r-1.7033990e+02 2.17e+02 1.31e+03 0.9 3.75e-02 4.5 1.00e+00 1.00e+00f 1
475r-1.7033652e+02 2.18e+02 1.84e+03 0.9 1.03e-01 4.1 1.00e+00 1.00e+00f 1
476r-1.7033530e+02 2.18e+02 1.20e+03 0.9 3.84e-02 4.5 1.00e+00 1.00e+00f 1
477r-1.7033485e+02 2.18e+02 2.31e+05 0.9 2.01e-02 4.9 1.00e+00 1.00e+00f 1
478r-1.7033355e+02 2.18e+02 1.18e+03 0.9 4.25e-02 4.4 1.00e+00 1.00e+00f 1
479r-1.7033002e+02 2.18e+02 1.17e+03 0.9 1.26e-01 4.0 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
480r-1.7032196e+02 2.18e+02 1.14e+03 0.9 3.71e-01 3.5 1.00e+00 1.00e+00f 1
481r-1.7031778e+02 2.18e+02 1.08e+03 0.9 9.81e-01 3.0 3.71e-01 4.89e-01f 1
482r-1.7033831e+02 2.18e+02 3.40e+03 0.9 7.66e-01 2.5 4.92e-01 6.96e-01f 1
483r-1.7033825e+02 2.18e+02 1.41e+03 0.9 2.20e-01 3.9 1.00e+00 8.97e-01f 1
484r-1.7033837e+02 2.18e+02 1.10e+03 0.9 2.67e-01 3.4 8.51e-01 2.66e-01f 1
485r-1.7033849e+02 2.18e+02 1.21e+03 0.9 6.79e-02 3.8 1.00e+00 1.00e+00h 1
486r-1.7034004e+02 2.18e+02 5.31e+01 0.9 2.45e-02 3.3 1.00e+00 1.00e+00h 1
487r-1.7035259e+02 2.18e+02 7.19e+01 -0.5 9.99e-02 2.9 9.64e-01 8.44e-01f 1
488r-1.7041323e+02 2.18e+02 2.82e+02 -0.5 1.98e-01 2.4 1.00e+00 4.05e-01f 1
489r-1.7105560e+02 2.17e+02 4.57e+02 -0.5 6.99e+06 - 6.01e-01 4.94e-03f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
490r-2.1357866e+02 2.52e+04 3.42e+02 -0.5 6.78e+06 - 1.00e+00 4.79e-01f 1
491r-2.1357866e+02 2.52e+04 9.99e+02 2.1 0.00e+00 - 0.00e+00 1.48e-07R 2
492r-2.1344993e+02 4.00e+03 9.98e+02 2.1 2.67e+06 - 7.83e-05 2.15e-03f 1
493r-2.1319997e+02 3.48e+03 1.04e+03 2.1 4.29e+05 - 1.46e-01 2.37e-03f 1
494r-2.0075024e+02 4.30e+02 4.09e+03 2.1 3.57e+05 - 9.86e-01 1.63e-01f 1
495r-1.8116445e+02 2.57e+02 2.02e+03 2.1 2.94e+05 - 8.05e-01 5.09e-01f 1
496r-1.8079784e+02 1.19e+02 3.12e+03 2.1 9.87e-01 2.0 5.73e-01 9.81e-01f 1
497r-1.8079560e+02 1.19e+02 6.57e+03 2.1 2.72e-01 4.2 9.90e-01 1.00e+00f 1
498r-1.8079584e+02 1.19e+02 8.12e+03 1.4 1.68e-01 4.7 8.82e-01 1.00e+00f 1
499r-1.8079548e+02 1.19e+02 4.71e+03 1.4 1.19e-01 4.2 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
500r-1.8027341e+02 1.17e+02 3.52e+03 1.4 5.99e+05 - 2.81e-01 1.46e-01f 1
501r-1.8233523e+02 3.19e+02 1.17e+03 1.4 5.09e+05 - 1.00e+00 7.26e-01f 1
502r-1.8233523e+02 3.19e+02 9.99e+02 2.0 0.00e+00 - 0.00e+00 4.19e-07R 9
503r-1.7683966e+02 1.42e+02 1.53e+03 2.0 3.87e+05 - 9.34e-04 7.88e-02f 1
504r-1.7737496e+02 1.33e+02 1.62e+03 2.0 3.56e+05 - 7.08e-01 9.10e-03f 1
505r-1.6339288e+02 1.09e+02 3.70e+03 2.0 3.52e+05 - 9.89e-01 2.79e-01f 1
506r-1.2813317e+02 2.89e+02 1.46e+03 2.0 2.52e+05 - 9.90e-01 1.00e+00f 1
507r-1.2813317e+02 2.89e+02 9.99e+02 2.0 0.00e+00 - 0.00e+00 4.49e-07R 5
508r-1.2806155e+02 2.84e+02 9.96e+02 2.0 4.30e+05 - 3.56e-01 3.44e-03f 1
509r-9.8143126e+01 1.20e+02 8.87e+02 2.0 4.28e+05 - 9.89e-01 5.30e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
510r-9.8142852e+01 8.64e+01 1.47e+03 2.0 4.63e-01 4.0 2.68e-01 2.79e-01h 1
511r-9.8142729e+01 8.52e+01 7.50e+03 2.0 4.36e-01 4.4 9.88e-01 3.36e-01h 1
512r-9.8142591e+01 8.52e+01 2.48e+04 2.0 3.49e-01 4.9 9.71e-01 1.00e+00h 1
513r-9.8142178e+01 8.52e+01 4.02e+03 2.0 4.81e-02 4.4 1.00e+00 1.00e+00f 1
514r-9.8142983e+01 8.52e+01 7.51e+02 1.3 9.51e-02 3.9 1.00e+00 1.00e+00f 1
515r-1.0255535e+02 3.42e+02 8.52e+02 1.3 7.42e+05 - 1.00e+00 2.57e-01f 1
516r-1.0255681e+02 9.01e+01 2.12e+02 0.6 5.84e-02 3.4 1.00e+00 9.85e-01f 1
517r-1.0255760e+02 9.02e+01 2.02e+02 0.6 2.87e-02 3.8 1.00e+00 1.00e+00f 1
518r-1.0256625e+02 9.02e+01 2.07e+02 -0.1 8.27e-02 3.4 7.10e-01 6.56e-01f 1
519r-1.0259901e+02 9.02e+01 2.13e+02 -0.1 1.73e-01 2.9 1.00e+00 4.17e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
520r-1.0266822e+02 9.02e+01 3.14e+02 -0.1 3.24e-01 2.4 1.00e+00 2.61e-01f 1
521r-1.0276331e+02 9.02e+01 1.31e+02 -0.1 1.23e-01 2.8 1.00e+00 1.00e+00f 1
522r-1.0278604e+02 9.02e+01 1.37e+02 -0.1 4.44e-02 3.3 1.00e+00 8.12e-01f 1
523r-1.0279526e+02 9.02e+01 1.36e+02 -0.1 2.48e-02 3.7 1.00e+00 9.37e-01f 1
524r-1.0279887e+02 9.02e+01 1.52e+02 -0.1 1.13e-02 4.1 1.00e+00 1.00e+00f 1
525r-1.0280730e+02 9.02e+01 1.19e+03 -0.1 1.58e-01 3.6 6.16e-01 8.38e-01f 1
526r-1.0283145e+02 9.02e+01 3.32e+03 -0.1 1.32e-01 3.2 5.65e-01 1.00e+00f 1
527r-1.0286025e+02 9.02e+01 6.45e+03 -0.1 3.28e+00 2.7 2.34e-01 5.19e-01f 1
528r-1.0554001e+02 8.86e+01 6.35e+03 -0.1 4.36e+06 - 1.59e-01 1.54e-02f 1
529r-1.3346031e+02 1.31e+03 5.09e+03 -0.1 4.20e+06 - 1.00e+00 1.78e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
530r-1.3346031e+02 1.31e+03 9.99e+02 1.8 0.00e+00 - 0.00e+00 1.15e-07R 2
531r-1.3284115e+02 2.29e+02 2.30e+03 1.8 4.95e+05 - 5.51e-05 1.96e-02f 1
532r-1.3300032e+02 2.26e+02 2.41e+03 1.8 4.77e+05 - 4.22e-01 1.88e-03f 1
533r-1.0631396e+02 6.24e+02 1.19e+04 1.8 4.76e+05 - 5.67e-01 5.24e-01f 1
534 -1.0416558e+02 6.22e+02 6.84e+01 -1.0 7.64e+09 - 2.19e-05 4.87e-06h 4
535r-1.0416558e+02 6.22e+02 9.99e+02 1.8 0.00e+00 - 0.00e+00 3.34e-07R 4
536r-1.0430418e+02 5.21e+02 9.98e+02 1.8 5.15e+05 - 1.20e-01 7.97e-04f 1
537r-9.6937904e+01 8.96e+01 8.83e+02 1.8 5.14e+05 - 4.15e-01 1.16e-01f 1
538r-8.9119544e+01 1.96e+02 8.64e+03 1.8 4.54e+05 - 4.68e-01 1.95e-01f 1
539r-8.4132491e+01 1.64e+02 7.22e+03 1.8 3.65e+05 - 1.06e-02 1.65e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
540r-8.4132194e+01 1.01e+02 4.70e+03 1.8 3.54e-01 4.0 8.47e-02 3.61e-01h 1
541r-8.4132131e+01 7.77e+01 1.87e+04 1.8 6.72e-01 4.4 1.58e-01 2.08e-01h 1
542r-8.4131653e+01 6.09e+01 5.98e+03 1.8 4.16e-01 3.9 5.12e-01 4.98e-01h 1
543r-8.4131643e+01 6.09e+01 7.13e+03 1.8 1.29e-01 5.3 9.90e-01 2.24e-01h 1
544r-8.4131559e+01 6.09e+01 6.98e+03 1.8 1.72e-01 4.8 5.31e-01 6.30e-01h 1
545r-7.7038613e+01 6.09e+01 4.87e+03 1.8 3.05e+05 - 9.91e-01 2.95e-01f 1
546r-6.3047680e+01 1.04e+02 2.30e+04 1.8 2.15e+05 - 1.00e+00 8.04e-01f 1
547r-6.0492990e+01 6.09e+01 5.40e+03 1.8 4.21e+04 - 1.00e+00 7.52e-01f 1
548r-5.9588530e+01 6.09e+01 3.87e+02 1.8 1.04e+04 - 1.00e+00 1.00e+00f 1
549r-9.3439891e+01 6.44e+01 2.21e+02 0.4 3.12e+04 - 8.30e-01 7.99e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
550r-9.4883661e+01 6.20e+01 7.59e+02 0.4 2.04e+06 - 9.36e-01 4.60e-02f 1
551r-9.4883466e+01 6.09e+01 4.87e+02 0.4 2.31e-02 4.3 1.00e+00 1.00e+00f 1
552r-9.5765669e+01 6.09e+01 5.59e+02 0.4 1.95e+06 - 1.00e+00 1.78e-02f 1
553r-9.8803815e+01 6.09e+01 6.26e+02 0.4 1.91e+06 - 1.00e+00 6.82e-02f 1
554r-1.0298482e+02 1.15e+02 5.91e+02 0.4 1.78e+06 - 1.00e+00 1.03e-01f 1
555r-1.1056500e+02 3.51e+02 4.70e+02 0.4 1.60e+06 - 1.00e+00 2.10e-01f 1
556r-1.3205168e+02 2.18e+03 1.76e+02 0.4 1.26e+06 - 1.00e+00 7.56e-01f 1
557r-1.3236543e+02 2.08e+03 1.68e+02 0.4 2.52e+05 - 1.00e+00 4.62e-02f 1
558r-1.3438735e+02 1.42e+03 1.20e+02 0.4 4.57e+04 - 1.00e+00 3.18e-01f 1
559r-1.5171007e+02 6.08e+01 4.40e+01 -0.3 3.38e+03 - 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
560r-1.6144004e+02 6.08e+01 3.72e+02 -0.3 3.21e+04 - 1.00e+00 2.43e-01f 1
561r-1.8909523e+02 6.08e+01 4.45e+00 -0.3 2.43e+04 - 1.00e+00 1.00e+00f 1
562r-1.9709267e+02 6.08e+01 2.54e+01 -1.7 5.87e+01 - 1.00e+00 8.92e-01f 1
563r-2.1826498e+02 6.08e+01 4.40e+02 -1.7 2.32e+05 - 1.00e+00 1.16e-01f 1
564r-2.2546789e+02 6.08e+01 5.05e+02 -1.7 2.05e+05 - 1.00e+00 2.59e-01f 1
565r-2.2572113e+02 6.07e+01 2.17e+02 -1.7 1.52e+05 - 1.00e+00 6.15e-01f 1
566r-2.2550412e+02 6.06e+01 2.33e-01 -1.7 5.86e+04 - 1.00e+00 1.00e+00f 1
567r-2.2561121e+02 6.06e+01 8.48e-02 -2.6 5.55e+00 - 1.00e+00 1.00e+00f 1
568r-2.2447833e+02 1.74e+02 3.64e+01 -3.8 4.86e+05 - 1.00e+00 8.27e-01f 1
569r-2.1752933e+02 6.47e+03 1.25e+01 -3.8 2.71e+06 - 1.00e+00 9.03e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
570r-2.1677685e+02 7.17e+01 1.48e-03 -3.8 2.63e+05 - 1.00e+00 1.00e+00f 1
571r-2.1677672e+02 5.83e+01 4.64e-08 -3.8 5.39e-01 - 1.00e+00 1.00e+00h 1
572r-2.1677712e+02 5.83e+01 1.03e-02 -5.8 5.40e-01 - 1.00e+00 1.00e+00f 1
573r-1.7730743e+02 1.76e+05 2.19e+02 -5.8 2.84e+07 - 1.00e+00 4.50e-01f 1
574r-1.6136095e+02 1.48e+05 7.88e+01 -5.8 1.56e+07 - 1.00e+00 2.37e-01f 1
575r-1.6136092e+02 1.48e+05 8.84e+01 -5.8 1.19e+07 - 1.00e+00 5.64e-07h 1
576r-1.5512337e+02 1.35e+05 2.36e+02 -5.8 1.19e+07 - 1.00e+00 9.49e-02f 1
577r-1.2386941e+02 9.56e+04 1.30e+02 -5.8 1.07e+07 - 1.00e+00 5.69e-01f 1
578r-1.2386929e+02 9.56e+04 1.03e+02 -5.8 4.62e+06 - 1.00e+00 2.34e-06h 1
579r-7.4273439e+01 2.08e+04 5.98e+00 -5.8 4.62e+06 - 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
580r-7.4268141e+01 2.06e+04 1.23e+02 -5.8 1.83e+00 3.8 9.94e-03 9.59e-03h 1
581r-7.4251342e+01 1.99e+04 8.99e+02 -5.8 1.81e+00 3.4 1.15e-02 3.07e-02h 1
582r-7.4120744e+01 1.50e+04 3.72e+03 -5.8 1.75e+00 2.9 4.30e-04 2.46e-01h 1
583r-7.3987860e+01 1.00e+04 2.95e+03 -5.8 1.31e+00 2.4 1.41e-04 3.33e-01h 1
584 -7.2336516e+01 1.00e+04 1.56e+02 -1.0 2.35e+08 - 3.76e-08 1.89e-06f 6
585r-7.2336516e+01 1.00e+04 9.99e+02 1.6 0.00e+00 - 0.00e+00 5.12e-08R 2
586r-7.2339594e+01 9.89e+03 1.88e+03 1.6 6.01e+05 - 5.41e-08 3.01e-05f 1
587r-7.1515348e+01 3.04e+03 1.78e+03 1.6 1.09e+06 - 3.69e-05 1.75e-03f 1
588r-7.1521801e+01 1.81e+03 1.79e+03 1.6 4.76e+05 - 5.04e-02 2.03e-03f 1
589r-7.1801318e+01 1.45e+03 1.78e+03 1.6 1.66e+04 - 2.26e-02 1.12e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
590r-6.8865381e+01 1.41e+02 1.73e+03 1.6 8.96e+03 - 3.14e-01 7.60e-02f 1
591r-6.6321005e+01 1.07e+02 3.29e+03 1.6 7.14e+03 - 9.90e-01 9.91e-02f 1
592r-4.5024248e+01 4.09e+01 3.38e+02 1.6 6.43e+03 - 9.90e-01 1.00e+00f 1
593r-4.4080018e+01 4.09e+01 1.44e+02 1.6 3.92e+02 - 1.00e+00 1.00e+00f 1
594r-6.3413622e+01 4.09e+01 6.96e+01 0.2 3.45e+02 - 8.57e-01 8.27e-01f 1
595r-6.3412285e+01 4.09e+01 1.87e+02 -0.5 1.70e-02 4.0 1.00e+00 7.38e-01f 1
596r-6.5201099e+01 4.09e+01 4.99e+02 -0.5 8.18e+04 - 1.00e+00 2.09e-01f 1
597r-6.9799490e+01 4.09e+01 8.59e+01 -0.5 6.47e+04 - 1.00e+00 8.35e-01f 1
598r-7.0449737e+01 4.09e+01 4.56e+00 -0.5 1.07e+04 - 1.00e+00 1.00e+00f 1
599r-7.1346044e+01 4.09e+01 3.62e+00 -1.2 3.44e+01 - 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
600r-7.0099748e+01 4.09e+01 2.10e+00 -1.2 1.11e+05 - 1.00e+00 1.00e+00f 1
601r-7.0059800e+01 4.09e+01 7.30e-01 -1.9 1.64e+01 - 1.00e+00 1.00e+00f 1
602r-6.7388145e+01 6.08e+01 5.49e-02 -1.9 2.48e+05 - 1.00e+00 1.00e+00f 1
603r-6.7378045e+01 4.09e+01 2.70e-04 -1.9 2.96e-01 - 1.00e+00 1.00e+00h 1
604r-6.7372992e+01 4.09e+01 7.29e-02 -4.3 1.09e+01 - 9.67e-01 9.65e-01f 1
605r-4.6258787e+01 2.33e+02 2.78e+02 -4.3 6.79e+06 - 3.65e-01 7.16e-02f 1
606r-4.6258773e+01 2.32e+02 2.67e+02 -4.3 2.05e-02 3.5 3.88e-02 2.43e-03h 1
607r-4.6258562e+01 2.25e+02 6.85e+02 -4.3 2.04e-02 3.0 1.00e+00 3.37e-02f 1
608r-4.6256872e+01 1.71e+02 5.12e+02 -4.3 1.97e-02 2.6 1.00e+00 2.40e-01h 1
609r-4.6253262e+01 7.46e+01 8.96e+01 -4.3 1.50e-02 2.1 1.00e+00 5.63e-01h 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
610r-4.6251317e+01 4.08e+01 2.70e-01 -4.3 6.55e-03 1.6 1.00e+00 1.00e+00h 1
611r-4.6251262e+01 4.08e+01 2.20e-03 -4.3 1.60e-04 1.1 1.00e+00 1.00e+00h 1
612r-4.6251119e+01 4.08e+01 6.90e+00 -6.5 6.23e-04 0.7 1.00e+00 9.29e-01f 1
613r-4.6250553e+01 4.08e+01 2.35e+01 -6.5 1.92e-03 0.2 9.44e-01 9.13e-01f 1
614r-4.6248692e+01 4.08e+01 2.93e-03 -6.5 5.76e-03 -0.3 1.00e+00 1.00e+00f 1
615r-4.6243109e+01 4.08e+01 2.93e-03 -6.5 1.73e-02 -0.8 1.00e+00 1.00e+00f 1
616r-4.6226355e+01 4.08e+01 2.93e-03 -6.5 5.19e-02 -1.2 1.00e+00 1.00e+00f 1
617r-4.6176048e+01 4.08e+01 2.93e-03 -6.5 1.56e-01 -1.7 1.00e+00 1.00e+00f 1
618r-4.6074536e+01 4.08e+01 1.27e+02 -6.5 4.68e-01 -2.2 1.00e+00 6.71e-01f 1
619r-4.5928671e+01 4.08e+01 2.33e+02 -6.5 1.41e+00 -2.7 1.00e+00 3.19e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
620r-4.4530864e+01 4.08e+01 2.01e-02 -6.5 4.31e+00 -3.2 1.00e+00 1.00e+00f 1
621r-3.9981187e+01 4.08e+01 2.62e-01 -6.5 1.39e+01 -3.6 1.00e+00 1.00e+00f 1
622r-3.8872315e+01 4.08e+01 1.15e+02 -6.5 5.77e+01 -4.1 1.00e+00 5.64e-02f 1
623r-3.4742289e+01 4.08e+01 1.86e+02 -6.5 1.74e+03 -4.6 3.12e-02 7.21e-03f 1
624r-2.5429876e+01 4.08e+01 3.28e+02 -6.5 1.15e+02 -4.2 1.00e+00 2.53e-01f 1
625r-2.3932487e+01 4.08e+01 1.48e+02 -6.5 4.66e+01 -3.7 8.10e-01 9.12e-02f 1
626r-2.3075895e+01 4.08e+01 1.12e+02 -6.5 1.11e+01 -3.3 1.80e-04 2.44e-01f 1
627r-2.3059718e+01 4.08e+01 1.11e+02 -6.5 1.52e+02 -3.8 1.20e-07 3.12e-04f 1
628r-2.3060073e+01 4.08e+01 1.11e+02 -6.5 1.00e+01 -3.4 8.29e-02 1.91e-04f 1
629r-2.1680878e+01 4.08e+01 1.03e+02 -6.5 9.69e+01 -3.8 6.52e-05 4.07e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
630r-2.1676418e+01 4.08e+01 1.03e+02 -6.5 1.84e+01 -3.4 2.40e-01 6.94e-04f 1
631r-2.1728085e+01 4.08e+01 1.30e+02 -6.5 1.21e+03 -3.9 4.99e-02 1.25e-04h 1
632r-1.6251835e+01 4.08e+01 3.36e+01 -6.5 2.45e+01 -3.5 1.00e+00 6.41e-01f 1
633r-1.2189075e+01 4.08e+01 1.93e+02 -6.5 1.63e+01 -3.0 1.00e+00 7.07e-01f 1
634r-1.1964986e+01 4.08e+01 1.28e+02 -6.5 3.28e+00 -1.7 1.14e-01 1.90e-01f 1
635r-1.1964986e+01 4.08e+01 1.28e+02 -6.5 4.08e-01 3.2 5.20e-06 2.89e-04f 1
636r-1.1964986e+01 4.08e+01 2.96e+02 -6.5 2.97e-02 3.7 1.00e+00 1.94e-06f 2
637r-1.1964986e+01 4.08e+01 2.06e+02 -6.5 1.13e-03 4.1 1.00e+00 1.40e-02f 1
638r-1.1964983e+01 4.08e+01 1.76e+02 -6.5 4.67e-03 3.6 1.00e+00 5.37e-01f 1
639r-1.0438751e+01 7.68e+02 1.89e+02 -6.5 9.17e+07 - 3.86e-01 9.59e-03f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
640r-1.0438409e+01 7.52e+02 9.00e+02 -6.5 6.73e-02 4.0 1.72e-01 1.99e-02h 1
641r-1.0433752e+01 5.44e+02 8.84e+02 -6.5 6.60e-02 3.6 2.09e-02 2.77e-01h 1
642r-1.0425314e+01 1.66e+02 2.86e+02 -6.5 4.77e-02 3.1 1.00e+00 6.95e-01h 1
643r-1.0424057e+01 1.09e+02 5.64e+02 -6.5 1.46e-02 2.6 5.83e-01 3.44e-01h 1
644r-1.0421691e+01 4.06e+01 8.86e+03 -6.5 9.56e-03 2.1 1.00e+00 1.00e+00h 1
645r-1.0425132e+01 4.06e+01 6.33e+02 -6.5 9.56e-03 1.6 4.70e-01 1.00e+00h 1
646r-1.0426751e+01 4.06e+01 2.01e+02 -6.5 4.48e-03 1.2 1.00e+00 1.00e+00h 1
647r-1.0431214e+01 4.06e+01 7.69e+01 -6.5 1.23e-02 0.7 1.00e+00 1.00e+00h 1
648r-1.0444069e+01 4.06e+01 2.33e+01 -6.5 3.55e-02 0.2 1.00e+00 1.00e+00h 1
649r-1.0487217e+01 4.06e+01 2.41e+00 -6.5 1.19e-01 -0.3 4.55e-01 1.00e+00h 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
650r-1.0487537e+01 4.06e+01 2.41e+00 -6.5 1.70e+00 -0.7 1.26e-01 5.19e-04f 1
651r-1.0558263e+01 4.06e+01 1.47e+01 -6.5 1.96e-01 -0.3 2.82e-01 1.00e+00f 1
652r-1.0558281e+01 4.06e+01 1.47e+01 -6.5 1.95e+00 0.1 2.28e-04 2.56e-05h 1
653r-1.0558365e+01 4.06e+01 1.47e+01 -6.5 7.36e-02 0.5 6.54e-01 3.16e-03f 1
654r-1.0559047e+01 4.06e+01 1.33e+01 -6.5 2.13e-02 1.0 7.60e-01 8.84e-02f 1
655r-1.0573516e+01 4.06e+01 1.28e+01 -6.5 1.04e-01 0.5 1.00e+00 3.83e-01f 1
656r-1.0577668e+01 4.06e+01 2.75e+01 -6.5 1.15e-02 1.8 5.86e-01 1.00e+00h 1
657r-1.0595056e+01 4.06e+01 1.72e+01 -6.5 9.57e-02 2.2 2.38e-01 1.00e+00f 1
658r-1.0595250e+01 4.06e+01 1.74e+01 -6.5 5.64e-02 2.7 1.00e+00 1.94e-02f 1
659r-1.0598792e+01 4.06e+01 2.47e+02 -6.5 1.71e+00 2.2 1.00e+00 1.05e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
660r-1.0598792e+01 4.06e+01 1.53e+03 -6.5 5.90e-02 1.7 1.00e+00 3.98e-05h 1
661r-1.0584740e+01 4.06e+01 1.14e+01 -6.5 3.89e-02 1.2 1.00e+00 1.00e+00f 1
662r-1.0589850e+01 4.06e+01 3.87e+02 -6.5 1.63e+00 0.8 6.32e-01 4.17e-01f 1
663r-1.0589851e+01 4.06e+01 2.25e+02 -6.5 1.37e+00 2.1 1.00e+00 2.17e-06f 1
664r-1.0588841e+01 4.06e+01 6.35e+02 -6.5 2.16e-01 1.6 1.00e+00 1.91e-01f 1
665r-1.0583975e+01 4.06e+01 1.13e+02 -6.5 6.59e-01 1.1 1.00e+00 8.63e-01f 1
666r-1.0588263e+01 4.06e+01 1.19e+02 -6.5 2.30e+00 0.7 1.00e+00 1.00e+00f 1
667r-1.0582064e+01 4.06e+01 2.83e+02 -6.5 5.63e+00 0.2 1.00e+00 8.19e-01f 1
668r-1.0583755e+01 4.06e+01 7.63e+02 -6.5 1.51e+01 -0.3 1.00e+00 7.60e-01f 1
669r-1.0574272e+01 4.06e+01 1.14e+03 -6.5 3.73e+01 -0.8 1.00e+00 1.79e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
670r-1.0574272e+01 4.06e+01 1.23e+03 -6.5 1.13e+02 -1.2 1.00e+00 2.31e-08f 1
671r-1.0219618e+01 2.97e+03 1.29e+03 -6.5 7.12e+01 -1.7 3.68e-01 5.51e-01f 1
672r-1.0219618e+01 2.97e+03 1.29e+03 -6.5 5.26e+02 -1.3 8.41e-02 3.04e-08h 1
673r-1.0219408e+01 2.40e+03 1.27e+03 -6.5 9.12e+01 -0.9 6.84e-01 2.63e-01f 1
674r-1.0219408e+01 2.40e+03 9.26e+02 -6.5 1.54e+02 -1.4 4.75e-01 5.89e-08h 1
675r-1.0161965e+01 2.37e+03 8.70e+02 -6.5 1.13e+02 -1.8 1.55e-07 3.40e-02f 1
676r-1.0161965e+01 2.37e+03 8.70e+02 -6.5 3.41e-01 3.1 4.34e-01 9.47e-07h 1
677r-1.0161960e+01 2.37e+03 8.70e+02 -6.5 2.06e+00 2.6 2.56e-01 1.10e-04f 1
678r-1.0105352e+01 4.06e+01 6.45e+02 -6.5 2.48e-01 2.2 4.70e-01 1.00e+00h 1
679r-1.0102490e+01 4.06e+01 3.60e+02 -6.5 1.70e-01 1.7 7.08e-06 4.62e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
680r-1.0100146e+01 4.06e+01 1.66e+02 -6.5 5.10e-01 1.2 1.00e+00 5.40e-01f 1
681r-1.0097229e+01 4.06e+01 8.15e+00 -6.5 1.53e+00 0.7 1.00e+00 1.00e+00f 1
682r-1.0094211e+01 4.06e+01 1.49e+01 -6.5 4.55e+00 0.3 1.00e+00 7.90e-01f 1
683r-1.0094211e+01 4.06e+01 2.32e+02 -6.5 1.24e+01 -0.2 1.00e+00 9.96e-08h 1
684r-1.0093722e+01 4.06e+01 1.92e+03 -6.5 1.31e+01 -0.7 9.83e-01 4.02e-02f 1
685r-1.0093721e+01 4.06e+01 1.75e+03 -6.5 1.35e+01 -1.2 1.00e+00 4.37e-06h 1
686r-9.8543013e+00 6.54e+02 7.57e+02 -6.5 4.04e+01 -1.7 1.00e+00 1.00e+00f 1
687r-9.5857754e+00 3.76e+03 5.71e+02 -6.5 1.90e+02 -2.1 4.87e-01 4.73e-01f 1
688r-9.5857754e+00 3.76e+03 5.71e+02 -6.5 1.39e+04 -2.6 3.14e-09 1.95e-09h 1
689r-9.5857754e+00 3.76e+03 1.29e+03 -6.5 1.77e+02 -2.2 1.00e+00 5.01e-08h 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
690r-9.7410705e+00 1.20e+04 1.49e+03 -6.5 8.25e+02 -2.7 5.55e-01 1.98e-01f 1
691r-9.7410705e+00 1.20e+04 1.49e+03 -6.5 2.85e+02 -2.2 3.34e-06 1.14e-07h 1
692r-9.7599203e+00 1.15e+04 1.24e+03 -6.5 5.72e+02 -2.7 5.40e-01 3.23e-02f 1
693r-9.7756841e+00 9.16e+03 1.94e+03 -6.5 1.36e+02 -2.3 3.64e-03 1.83e-01f 1
694r-9.7756842e+00 9.16e+03 1.94e+03 -6.5 5.74e+02 -2.8 1.76e-02 9.75e-08h 1
695r-9.7579617e+00 6.47e+03 1.18e+03 -6.5 1.51e+02 -2.3 9.91e-05 4.18e-01f 1
696r-9.7975565e+00 6.47e+03 1.17e+03 -6.5 1.16e+03 -2.8 4.82e-02 9.01e-03f 1
697r-9.9343736e+00 7.56e+03 7.59e+02 -6.5 2.49e+02 -2.4 1.00e+00 3.54e-01f 1
698r-9.9343736e+00 7.56e+03 7.59e+02 -6.5 1.21e+02 -2.0 3.37e-01 2.61e-07h 1
699r-1.0684358e+01 8.76e+04 5.93e+03 -6.5 2.27e+03 -2.4 1.30e-02 2.05e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
700r-1.2031126e+01 7.88e+04 5.83e+03 -6.5 2.33e+03 -2.9 9.78e-08 5.87e-02f 1
701r-1.1405079e+01 5.21e+04 4.78e+03 -6.5 5.57e+02 -2.5 6.29e-02 3.75e-01h 1
702r-1.1405079e+01 5.21e+04 4.71e+03 -6.5 7.31e+03 -1.2 4.37e-04 2.30e-09h 1
703r-1.1439635e+01 4.89e+04 5.87e+03 -6.5 2.14e+02 -0.7 1.01e-07 7.77e-02f 1
704r-1.1439453e+01 4.89e+04 5.87e+03 -6.5 3.09e+02 0.6 2.47e-03 1.40e-04h 1
705r-1.1438275e+01 4.88e+04 5.87e+03 -6.5 3.01e+02 0.1 7.74e-09 1.30e-03h 1
706r-1.1437459e+01 4.76e+04 5.77e+03 -6.5 2.99e+01 0.5 4.20e-03 2.45e-02h 1
707r-1.1437733e+01 4.75e+04 5.75e+03 -6.5 1.21e+02 0.1 7.10e-09 3.06e-03h 1
708r-1.1440114e+01 4.49e+04 5.30e+03 -6.5 2.02e+01 0.5 2.60e-03 5.40e-02h 1
709r-1.1440114e+01 4.49e+04 5.30e+03 -6.5 5.10e+00 4.5 0.00e+00 2.99e-07R 3
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
710r-1.1440114e+01 4.49e+04 5.30e+03 -6.5 1.12e-02 5.0 9.76e-01 3.12e-07f 1
711r-1.1431615e+01 4.49e+04 1.00e+03 -6.5 3.35e-02 4.5 5.60e-01 8.77e-01F 1
712r-1.1431615e+01 4.49e+04 1.00e+03 -6.5 1.01e-01 4.0 9.82e-01 1.74e-06f 1
713r-1.1429060e+01 4.49e+04 1.00e+03 -6.5 3.02e-01 3.5 9.86e-01 3.35e-02f 4
714r-1.1424058e+01 4.49e+04 1.00e+03 -6.5 1.13e-01 3.9 1.00e+00 1.77e-01f 3
715r-1.1423548e+01 4.49e+04 1.00e+03 -6.5 3.40e-01 3.5 8.67e-01 8.41e-03f 5
716r-1.1422298e+01 4.49e+04 1.00e+03 -6.5 1.27e-01 3.9 1.00e+00 5.45e-02f 4
717r-1.1421887e+01 4.49e+04 1.00e+03 -6.5 3.82e-01 3.4 4.87e-01 1.03e-02f 4
718r-1.1416008e+01 4.49e+04 1.00e+03 -6.5 1.43e-01 3.8 1.00e+00 2.90e-01f 1
719r-1.1416007e+01 4.49e+04 1.00e+03 -6.5 4.30e-01 3.4 7.04e-01 8.50e-05f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
720r-1.1399159e+01 4.49e+04 1.54e+03 -6.5 2.03e+01 2.9 1.25e-02 1.55e-02f 1
721r-1.1399159e+01 4.49e+04 2.40e+03 -6.5 4.83e-01 3.3 3.51e-01 3.42e-07f 1
722r-1.1409397e+01 4.49e+04 2.95e+03 -6.5 3.54e+00 2.8 1.54e-03 2.29e-02f 1
723r-1.1409397e+01 4.49e+04 3.48e+03 -6.5 5.44e-01 3.3 2.47e-01 9.03e-09f 1
724r-1.1413180e+01 4.49e+04 3.00e+03 -6.5 1.63e+00 2.8 2.46e-01 2.67e-02f 1
725r-1.1413529e+01 4.49e+04 2.22e+03 -6.5 4.90e+00 2.3 3.36e-01 5.63e-03f 1
726r-1.1471285e+01 4.47e+04 2.34e+03 -6.5 1.47e+01 1.8 1.73e-07 3.27e-01f 1
727r-1.1469249e+01 4.47e+04 1.63e+03 -6.5 6.88e-01 3.2 6.72e-01 1.39e-01f 1
728r-1.1469249e+01 4.47e+04 1.86e+03 -6.5 2.07e+00 2.7 5.32e-01 1.76e-07f 1
729r-1.1506594e+01 4.45e+04 1.00e+03 -6.5 6.20e+00 2.2 7.40e-01 5.56e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
730r-1.1506595e+01 4.45e+04 1.00e+03 -6.5 1.86e+01 1.7 1.00e+00 5.05e-07f 1
731r-1.1671896e+01 4.39e+04 1.00e+03 -6.5 5.58e+01 1.3 1.00e+00 2.75e-01f 1
732r-1.1671896e+01 4.39e+04 1.00e+03 -6.5 1.67e+02 0.8 5.73e-01 2.86e-08f 1
733r-1.1681907e+01 4.38e+04 1.00e+03 -6.5 5.02e+02 0.3 5.43e-01 2.30e-03f 1
734r-1.1681938e+01 4.38e+04 1.00e+03 -6.5 1.50e+03 -0.2 4.65e-08 3.14e-06f 1
735r-1.1698677e+01 4.37e+04 1.00e+03 -6.5 4.52e+03 -0.7 1.60e-03 5.88e-04f 1
736r-1.1698682e+01 4.37e+04 1.00e+03 -6.5 1.69e+03 -0.2 4.51e-03 7.71e-07f 1
737r-1.1827485e+01 4.24e+04 1.00e+03 -6.5 5.09e+03 -0.7 4.19e-03 6.13e-03f 1
738r-1.1899068e+01 4.16e+04 1.00e+03 -6.5 1.90e+03 -0.3 9.94e-02 8.93e-03f 1
739r-1.1899068e+01 4.16e+04 1.00e+03 -6.5 5.76e+03 -0.8 4.31e-03 4.20e-09f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
740r-1.2003296e+01 3.94e+04 1.91e+03 -6.5 2.14e+03 -0.3 4.01e-01 2.36e-02f 1
741r-1.2003296e+01 3.94e+04 1.68e+03 -6.5 8.03e+02 0.1 3.74e-01 2.28e-08f 1
742r-1.2046277e+01 3.90e+04 1.94e+03 -6.5 2.41e+03 -0.4 3.95e-02 4.37e-03f 1
743r-1.2046278e+01 3.90e+04 1.66e+03 -6.5 9.03e+02 0.0 1.00e+00 5.20e-07f 1
744r-1.2136287e+01 3.70e+04 1.80e+03 -6.5 2.71e+03 -0.4 3.64e-02 1.71e-02f 1
745r-1.2136288e+01 3.70e+04 1.95e+03 -6.5 1.02e+03 -0.0 1.19e-01 2.77e-07f 1
746r-1.2210891e+01 3.60e+04 1.36e+03 -6.5 3.06e+03 -0.5 6.37e-10 7.74e-03f 1
747r-1.2210891e+01 3.60e+04 1.36e+03 -6.5 2.79e-01 3.6 1.00e+00 2.65e-05f 1
748r-1.2212699e+01 3.60e+04 1.20e+03 -6.5 8.37e-01 3.1 7.62e-01 1.15e-01f 1
749r-1.2212699e+01 3.60e+04 1.00e+03 -6.5 2.51e+00 2.6 2.83e-01 3.55e-07f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
750r-1.2217809e+01 3.59e+04 1.00e+03 -6.5 7.53e+00 2.1 1.00e+00 2.88e-01f 1
751r-1.2265888e+01 3.49e+04 1.00e+03 -6.5 2.26e+01 1.6 1.00e+00 1.00e+00f 1
752r-1.2403850e+01 3.20e+04 1.00e+03 -6.5 6.78e+01 1.2 1.00e+00 1.00e+00f 1
753r-1.2600136e+01 2.78e+04 1.00e+03 -6.5 2.03e+02 0.7 1.00e+00 4.73e-01f 1
754r-1.2600136e+01 2.78e+04 1.16e+03 -6.5 6.10e+02 0.2 6.86e-01 4.90e-08f 1
755r-1.2599207e+01 2.51e+04 1.00e+03 -6.5 1.83e+03 -0.3 1.51e-01 3.39e-02f 1
756r-1.2599207e+01 2.51e+04 1.00e+03 -6.5 5.49e+03 -0.7 1.00e+00 3.69e-09f 1
757r-1.1368220e+01 4.56e+01 1.00e+03 -6.5 1.65e+04 -1.2 2.57e-01 3.54e-02f 1
758r-1.1368220e+01 4.56e+01 1.00e+03 -6.5 4.82e+04 -1.7 1.46e-08 3.74e-09f 1
759r-1.1355116e+01 4.06e+01 6.89e+02 -6.5 3.90e-01 3.2 1.00e+00 1.00e+00h 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
760r-1.1353827e+01 4.06e+01 1.20e+02 -6.5 2.03e-01 2.8 1.00e+00 9.42e-01f 1
761r-1.1352949e+01 4.06e+01 4.23e+01 -6.5 2.15e-01 2.3 1.00e+00 9.19e-01f 1
762r-1.1352512e+01 4.06e+01 2.05e+00 -6.5 3.17e-02 1.8 1.00e+00 9.63e-01f 1
763r-1.1351644e+01 4.06e+01 2.08e+00 -6.5 9.50e-02 1.3 1.00e+00 1.00e+00f 1
764r-1.1349188e+01 4.06e+01 2.07e+00 -6.5 2.85e-01 0.9 1.00e+00 1.00e+00f 1
765r-1.1342433e+01 4.06e+01 4.98e+00 -6.5 8.50e-01 0.4 1.00e+00 1.00e+00f 1
766r-1.1326798e+01 4.06e+01 1.35e+01 -6.5 2.52e+00 -0.1 1.00e+00 1.00e+00f 1
767r-1.1291394e+01 4.06e+01 3.00e+00 -6.5 7.50e+00 -0.6 1.00e+00 1.00e+00f 1
768r-1.1197570e+01 2.16e+02 2.11e+01 -6.5 2.27e+01 -1.0 1.00e+00 9.51e-01f 1
769r-1.1165334e+01 2.76e+02 3.32e+02 -6.5 7.22e+01 -1.5 3.10e-06 1.46e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
770r-1.1165334e+01 2.76e+02 1.32e+03 -6.5 2.58e+01 -1.1 7.02e-01 1.06e-07h 1
771r-1.1165334e+01 2.76e+02 6.81e+02 -6.5 7.87e+01 -1.6 3.45e-01 1.16e-07h 1
772r-1.1160900e+01 2.75e+02 6.31e+02 -6.5 2.39e+02 -2.1 5.34e-08 4.06e-03f 1
773r-1.1160897e+01 2.75e+02 6.31e+02 -6.5 1.66e+00 0.2 6.39e-01 7.29e-04h 1
774r-1.1148228e+01 4.06e+01 1.25e+01 -6.5 3.77e+00 -0.3 1.00e+00 1.00e+00f 1
775r-1.1128012e+01 4.06e+01 4.34e+01 -6.5 1.04e+01 -0.8 4.53e-01 5.09e-01f 1
776r-1.1118199e+01 4.06e+01 1.70e+02 -6.5 4.73e+00 -0.3 1.29e-03 6.41e-01f 1
777r-1.1116660e+01 1.28e+02 5.63e+02 -6.5 1.63e+02 -0.8 7.98e-03 4.22e-02f 1
778r-1.1116660e+01 1.28e+02 5.63e+02 -6.5 7.52e+00 -0.4 0.00e+00 2.50e-11R 12
779r-1.1115325e+01 1.25e+02 5.51e+02 -6.5 1.39e+01 -0.9 1.00e+00 2.13e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
780r-1.1115325e+01 1.25e+02 1.78e+03 -6.5 3.94e+01 -1.4 1.00e+00 2.57e-06f 1
781r-1.1081277e+01 1.11e+02 1.83e+03 -6.5 1.27e+02 -1.8 9.76e-01 5.65e-02f 5
782r-1.1064056e+01 1.03e+02 2.42e+03 -6.5 5.15e+02 -2.3 9.27e-01 1.05e-02f 7
783r-1.1035672e+01 1.00e+02 2.05e+03 -6.5 1.47e+02 -1.9 5.02e-01 4.27e-02f 5
784r-1.1021970e+01 1.10e+02 2.32e+03 -6.5 6.49e+02 -2.4 5.00e-01 7.90e-03f 7
785r-1.1021840e+01 1.10e+02 2.31e+03 -6.5 1.76e+02 -1.9 1.00e+00 1.49e-04f 1
786r-1.1016303e+01 1.12e+02 2.31e+03 -6.5 8.27e+02 -2.4 9.68e-04 3.00e-03f 8
787r-1.1016303e+01 1.12e+02 2.09e+03 -6.5 1.89e+02 -2.0 1.00e+00 8.33e-09h 1
788r-1.1014139e+01 1.13e+02 2.20e+03 -6.5 1.09e+03 -2.5 8.08e-02 1.12e-03f 9
789r-1.0138221e+01 2.39e+04 7.80e+02 -6.5 2.28e+02 -2.0 1.00e+00 9.64e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
790r-1.0123778e+01 2.02e+04 1.33e+03 -6.5 1.05e+02 -1.6 3.99e-01 1.61e-01f 1
791r-1.0123778e+01 2.02e+04 1.33e+03 -6.5 2.62e+00 0.6 4.84e-02 1.37e-07h 1
792r-1.0123764e+01 2.01e+04 1.57e+03 -6.5 2.66e+00 0.1 1.00e+00 5.85e-03f 1
793r-1.0123734e+01 2.01e+04 1.99e+03 -6.5 5.36e+00 -0.3 2.92e-01 8.82e-04h 1
794r-1.0086367e+01 5.23e+03 3.21e+03 -6.5 1.92e+01 -0.8 6.14e-08 7.62e-01f 1
795r-1.0086356e+01 5.23e+03 3.21e+03 -6.5 2.42e+00 0.5 4.13e-01 9.49e-04h 1
796r-1.0081866e+01 4.48e+03 2.66e+03 -6.5 3.92e+00 0.0 7.62e-01 1.43e-01f 1
797r-1.0081866e+01 4.48e+03 2.66e+03 -6.5 1.90e+00 0.5 0.00e+00 9.74e-10R 7
798r-1.0081866e+01 4.48e+03 2.66e+03 -6.5 3.11e-02 4.5 5.90e-01 2.39e-07f 1
799r-1.0081521e+01 4.48e+03 2.67e+03 -6.5 9.32e-02 4.0 2.41e-01 8.55e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
800r-1.0081596e+01 4.48e+03 2.65e+03 -6.5 2.79e-01 3.6 2.10e-07 1.47e-02f 1
801r-1.0082824e+01 4.47e+03 2.12e+03 -6.5 8.38e-01 3.1 1.97e-05 2.10e-01f 1
802r-1.0082426e+01 4.46e+03 1.83e+03 -6.5 2.51e+00 2.6 1.70e-01 1.34e-01f 1
803r-1.0079462e+01 4.30e+03 2.12e+03 -6.5 7.54e+00 2.1 7.57e-08 4.88e-01f 1
804r-1.0079284e+01 4.30e+03 1.54e+03 -6.5 4.42e-02 4.4 1.00e+00 2.80e-01f 1
805r-1.0076745e+01 4.30e+03 1.01e+03 -6.5 1.33e-01 3.9 1.00e+00 3.18e-01f 1
806r-1.0076745e+01 4.30e+03 1.91e+03 -6.5 3.98e-01 3.4 3.60e-01 1.00e-07f 1
807r-1.0077618e+01 4.28e+03 1.19e+03 -6.5 1.19e+00 2.9 3.27e-07 2.06e-01f 1
808r-1.0077912e+01 4.28e+03 1.70e+03 -6.5 5.60e-02 4.3 3.29e-01 1.16e-01f 1
809r-1.0078784e+01 4.28e+03 1.00e+03 -6.5 1.68e-01 3.8 4.06e-07 9.35e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
810r-1.0078636e+01 4.28e+03 1.00e+03 -6.5 5.04e-01 3.3 5.15e-02 3.02e-02f 1
811r-1.0078027e+01 4.26e+03 1.00e+03 -6.5 1.51e+00 2.8 1.00e+00 2.87e-01f 1
812r-1.0077740e+01 4.21e+03 1.00e+03 -6.5 4.53e+00 2.3 1.00e+00 2.46e-01f 1
813r-1.0077895e+01 3.71e+03 1.00e+03 -6.5 1.36e+01 1.9 1.00e+00 8.52e-01f 1
814r-1.0077937e+01 3.57e+03 1.00e+03 -6.5 4.08e+01 1.4 1.00e+00 8.03e-02f 1
815r-1.0077937e+01 3.57e+03 1.00e+03 -6.5 1.22e+02 0.9 1.00e+00 1.08e-07f 1
816r-1.0064953e+01 3.21e+03 1.00e+03 -6.5 3.67e+02 0.4 5.92e-01 2.27e-02f 1
817r-1.0064950e+01 3.21e+03 1.00e+03 -6.5 1.10e+03 -0.0 1.00e+00 1.01e-06f 1
818r-1.0043898e+01 2.66e+03 2.02e+03 -6.5 3.30e+03 -0.5 2.96e-01 3.82e-03f 1
819r-1.0041187e+01 2.59e+03 2.02e+03 -6.5 9.91e+03 -1.0 9.85e-08 1.83e-04f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
820r-1.0040450e+01 2.54e+03 2.02e+03 -6.5 2.97e+04 -1.5 2.16e-02 3.34e-05f 1
821r-1.0038827e+01 2.49e+03 1.99e+03 -6.5 8.92e+04 -2.0 2.15e-02 1.26e-05f 1
822r-1.0057488e+01 5.98e+02 1.69e+03 -6.5 3.34e+04 -1.5 4.09e-01 1.32e-03f 1
823r-1.0080589e+01 1.49e+02 1.68e+03 -6.5 1.59e+04 -2.0 6.97e-02 8.73e-04f 1
824r-1.0111349e+01 1.65e+02 1.67e+03 -6.5 7.78e+03 -2.5 2.10e-03 4.51e-04f 7
825r-1.0130232e+01 1.97e+02 1.65e+03 -6.5 1.93e+02 -2.1 1.70e-04 3.07e-02f 2
826r-1.0130253e+01 1.97e+02 1.10e+03 -6.5 6.84e+01 -1.6 1.00e+00 2.58e-04f 1
827r-1.0956235e+01 5.98e+04 1.09e+03 -6.5 2.85e+02 -2.1 1.50e-01 7.81e-01f 1
828r-1.0951852e+01 5.39e+04 9.61e+02 -6.5 5.37e+01 -0.8 3.11e-01 9.81e-02f 1
829r-1.0921291e+01 4.44e+04 1.01e+03 -6.5 4.58e+01 -0.3 1.90e-07 1.76e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
830r-1.0846912e+01 2.12e+04 2.50e+03 -6.5 7.27e+01 -0.8 1.28e-01 4.99e-01f 1
831r-1.0846912e+01 2.12e+04 2.50e+03 -6.5 7.65e+01 -1.3 0.00e+00 1.09e-09R 5
832r-1.0846902e+01 2.12e+04 2.50e+03 -6.5 3.40e+04 -1.8 2.85e-01 6.79e-07f 1
833r-1.0784314e+01 1.59e+04 2.49e+03 -6.5 5.44e+04 -2.3 8.38e-02 2.26e-03f 3
834r-1.0799724e+01 1.51e+04 2.49e+03 -6.5 4.06e+04 -2.7 1.54e-02 4.49e-04f 4
835r-1.0782314e+01 1.32e+04 2.49e+03 -6.5 3.31e+04 -2.3 2.20e-08 1.32e-03f 4
836r-1.0754140e+01 9.90e+03 2.48e+03 -6.5 2.38e+04 -1.9 8.36e-01 3.22e-03f 3
837r-1.0754148e+01 9.90e+03 2.48e+03 -6.5 1.59e+04 -2.4 7.20e-04 1.28e-06f 1
838r-1.0738844e+01 7.42e+03 2.47e+03 -6.5 1.42e+04 -1.9 4.63e-01 4.04e-03f 3
839r-1.0745389e+01 7.19e+03 2.47e+03 -6.5 9.25e+03 -2.4 1.86e-08 5.54e-04f 6
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
840r-1.0742700e+01 6.29e+03 2.46e+03 -6.5 8.48e+03 -2.0 4.95e-02 2.46e-03f 4
841r-1.0747031e+01 6.28e+03 2.46e+03 -6.5 9.36e+03 -2.5 6.19e-05 4.70e-05f 9
842r-1.0746585e+01 5.88e+03 2.46e+03 -6.5 6.47e+03 -2.0 1.00e+00 1.41e-03f 5
843r-1.0744959e+01 1.35e+03 2.39e+03 -6.5 5.26e+03 -1.6 1.20e-02 2.60e-02f 1
844r-1.0745177e+01 1.32e+03 2.39e+03 -6.5 5.79e+02 -2.1 1.00e+00 1.27e-03f 1
845r-1.0592214e+01 6.86e+03 1.98e+03 -6.5 7.34e+02 -2.6 7.09e-05 1.03e-01f 1
846r-1.0443926e+01 7.01e+03 1.73e+03 -6.5 2.26e+02 -2.1 6.10e-02 1.03e-01f 1
847r-1.0435037e+01 6.91e+03 1.69e+03 -6.5 1.16e+02 -1.7 5.47e-02 1.45e-02f 1
848r-1.0435031e+01 6.91e+03 1.68e+03 -6.5 2.22e+02 -2.2 1.00e+00 3.56e-06f 1
849r-1.0391564e+01 7.54e+03 1.58e+03 -6.5 8.79e+02 -2.7 5.06e-02 5.55e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
850r-1.0391537e+01 7.54e+03 1.40e+03 -6.5 1.84e+03 -3.1 4.76e-01 3.07e-06f 1
851r-1.0265881e+01 7.44e+03 1.39e+03 -6.5 4.91e+03 -3.6 7.29e-01 6.01e-03f 1
852r-1.0265889e+01 7.44e+03 1.39e+03 -6.5 1.70e+03 -4.1 1.28e-01 3.25e-07h 1
853r-1.3116502e+01 8.84e+03 1.33e+03 -6.5 5.92e+03 -4.6 1.23e-01 2.41e-02f 1
854r-1.9227388e+01 2.44e+04 1.06e+03 -6.5 2.42e+03 -4.1 3.59e-01 1.40e-01f 1
855r-1.9227389e+01 2.44e+04 1.23e+03 -6.5 9.24e+02 -3.7 1.00e+00 8.01e-08h 1
856r-1.9223341e+01 2.44e+04 1.29e+03 -6.5 9.08e+07 - 5.38e-02 9.18e-06f 1
857r-2.1765718e+00 6.72e+04 1.54e+03 -6.5 9.08e+07 - 8.59e-02 3.85e-02f 1
858r-1.8987327e+00 1.82e+04 1.28e+03 -6.5 7.86e+01 -1.5 1.23e-07 7.29e-01f 1
859r-1.8916213e+00 1.72e+04 1.21e+03 -6.5 3.99e+00 0.7 1.47e-01 5.48e-02f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
860r-1.8663649e+00 1.39e+04 1.41e+03 -6.5 9.91e+00 0.3 4.88e-07 1.95e-01f 1
861r-1.8619247e+00 1.33e+04 1.36e+03 -6.5 1.08e+01 -0.2 2.44e-01 3.83e-02f 1
862r-1.8619247e+00 1.33e+04 1.36e+03 -6.5 3.29e+01 -0.7 7.61e-06 2.21e-08h 1
863r-1.8619894e+00 1.33e+04 2.84e+03 -6.5 8.70e+00 -0.3 6.14e-02 2.78e-03f 1
864r-1.8568665e+00 1.22e+04 1.24e+03 -6.5 3.80e+01 -0.7 2.70e-08 8.15e-02f 1
865r-1.8568541e+00 1.22e+04 1.24e+03 -6.5 9.16e+00 -0.3 4.61e-02 1.64e-04f 1
866r-1.8902346e+00 8.16e+03 1.24e+03 -6.5 7.09e+00 0.1 2.08e-04 3.32e-01f 1
867r-1.9216910e+00 3.92e+03 2.01e+03 -6.5 5.16e+00 0.5 1.99e-04 5.20e-01f 1
868r-1.9206165e+00 3.68e+03 1.88e+03 -6.5 3.16e-01 1.9 6.20e-01 6.07e-02f 1
869r-1.9179301e+00 1.75e+03 1.31e+03 -6.5 4.29e-01 1.4 1.70e-07 5.24e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
870r-1.9243325e+00 1.51e+03 1.13e+03 -6.5 5.62e-01 0.9 8.97e-03 1.36e-01f 1
871r-1.9401167e+00 4.58e+02 3.42e+02 -6.5 1.65e+00 0.4 9.92e-01 6.97e-01f 1
872 -1.3725059e+00 3.85e+03 1.32e+06 -1.0 2.31e+08 - 1.84e-07 4.21e-03f 6
873r-1.3725059e+00 3.85e+03 9.99e+02 2.0 0.00e+00 - 0.00e+00 3.91e-07R 7
874r-8.7215869e-01 8.69e+02 9.99e+02 2.0 2.45e+05 - 1.24e-05 9.91e-04f 1
875 -6.4260566e-01 3.65e+03 1.15e+02 -1.0 1.43e+08 - 4.42e-03 4.42e-03s 15
876r-6.4260566e-01 3.65e+03 9.99e+02 1.9 0.00e+00 - 0.00e+00 0.00e+00R 1
877r 3.0718000e-02 9.22e+02 9.99e+02 1.9 3.19e+05 - 6.69e-05 9.91e-04f 1
878r 6.6419154e-02 5.66e+02 1.08e+03 1.9 2.55e+04 - 1.40e-01 1.17e-03f 1
879r 6.5122646e+00 2.08e+02 8.54e+02 1.9 5.29e+03 - 9.44e-01 1.41e-01f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
880 6.5292457e+00 2.10e+02 1.53e+03 -1.0 1.50e+08 - 5.73e-02 3.60e-04h 1
881 6.5712919e+00 2.12e+02 1.81e+04 -1.0 2.30e+08 - 5.46e-02 3.12e-04h 5
882 6.5878718e+00 2.13e+02 3.35e+04 -1.0 2.30e+08 - 5.73e-02 1.23e-04h 10
883 6.5960879e+00 2.13e+02 1.26e+05 -1.0 2.30e+08 - 3.92e-01 6.07e-05h 12
884 6.6023608e+00 2.13e+02 1.27e+05 -1.0 2.30e+08 - 6.72e-02 4.62e-05h 13
885 3.8711414e+01 6.54e+06 1.03e+05 -1.0 2.30e+08 - 1.89e-01 1.89e-01s 20
886 4.5742106e+01 6.20e+06 9.81e+04 -1.0 2.12e+08 - 4.88e-02 4.88e-02s 20
887 6.1790486e+01 5.45e+06 8.70e+04 -1.0 2.07e+08 - 1.12e-01 1.12e-01s 20
888r 6.1790486e+01 5.45e+06 1.00e+03 5.1 0.00e+00 -4.0 0.00e+00 0.00e+00R 1
889r 5.9095753e+01 5.10e+06 1.05e+04 5.1 7.66e+07 - 4.13e-03 1.05e-04f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
890r 2.9038017e+01 8.94e+05 2.99e+04 5.1 2.92e+07 - 6.71e-02 4.01e-03f 1
891r 1.1569213e+02 4.96e+05 6.71e+04 5.1 1.98e+05 - 6.48e-01 6.76e-02f 1
892r 3.6863587e+02 1.18e+05 3.11e+03 5.1 2.52e+03 - 9.91e-01 7.16e-01f 1
893r 2.2217222e+02 4.71e+04 2.71e+03 3.0 3.95e+02 - 8.08e-01 8.53e-01f 1
894r 1.9858400e+02 4.29e+04 5.81e+03 2.3 1.50e+03 - 8.24e-01 7.99e-01f 1
895r 1.4232498e+02 3.27e+04 2.57e+02 2.3 2.73e+03 - 5.72e-01 6.35e-01f 1
896r 1.3176095e+02 2.92e+04 1.10e+04 1.6 2.91e+03 - 7.75e-01 1.04e-01f 1
897r 5.3816329e+01 2.13e+04 7.14e+02 1.6 8.14e+03 - 6.93e-01 8.04e-01f 1
898r 5.3821532e+01 6.03e+03 8.98e+03 1.6 1.25e+00 4.0 5.27e-01 7.17e-01f 1
899r 5.3825745e+01 2.71e+01 1.09e+03 1.6 3.28e-01 3.5 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
900 5.4269705e+01 2.71e+01 1.02e+03 -1.0 2.07e+08 - 6.76e-01 6.56e-04h 1
901 5.4403375e+01 2.71e+01 5.22e+04 -1.0 1.93e+08 - 7.36e-01 2.27e-04h 7
902 5.4476720e+01 2.71e+01 2.40e+04 -1.0 1.93e+08 - 7.38e-01 1.25e-04h 13
903 5.4551462e+01 2.71e+01 7.50e+03 -1.0 1.93e+08 - 7.87e-01 1.27e-04h 13
904 5.4627179e+01 2.71e+01 1.56e+03 -1.0 1.93e+08 - 8.51e-01 1.29e-04h 13
905 5.4703245e+01 2.71e+01 1.91e+02 -1.0 1.93e+08 - 8.58e-01 1.29e-04h 13
906 5.4779352e+01 2.71e+01 3.71e+01 -1.0 1.93e+08 - 5.51e-01 1.30e-04h 13
907 5.4855409e+01 2.71e+01 3.74e+01 -1.0 1.93e+08 - 8.47e-03 1.30e-04h 13
908 5.5005769e+01 2.71e+01 3.74e+01 -1.0 2.09e+08 - 2.73e-04 2.65e-04f 12
909 5.5080157e+01 2.71e+01 4.12e+01 -1.0 2.22e+08 - 1.32e-01 1.35e-04h 13
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
910 5.5156007e+01 2.71e+01 4.25e+01 -1.0 1.93e+08 - 3.45e-02 1.30e-04h 13
911 3.6408532e+02 4.73e+06 2.11e+04 -1.0 2.20e+08 - 3.71e-04 5.51e-01w 1
912 3.4333206e+02 4.63e+06 1.80e+04 -1.0 7.86e+08 - 2.13e-01 1.85e-02w 1
913 3.4301479e+02 4.62e+06 1.80e+04 -1.0 4.29e+07 - 5.68e-02 1.23e-03w 1
914 5.5230432e+01 2.71e+01 4.25e+01 -1.0 1.26e+08 - 3.71e-04 1.35e-04h 12
915 5.5247199e+01 2.70e+01 4.25e+01 -1.0 5.97e+08 - 2.96e-04 1.35e-04f 11
916 5.5254012e+01 2.70e+01 4.24e+01 -1.0 6.13e+08 - 9.22e-03 6.47e-05h 12
917 5.5269186e+01 2.70e+01 4.24e+01 -1.0 4.56e+08 - 3.64e-04 5.35e-05h 13
918 5.5269980e+01 2.70e+01 4.22e+01 -1.0 6.99e+08 - 1.28e-02 1.07e-04h 11
919 5.5274801e+01 2.70e+01 4.22e+01 -1.0 6.35e+08 - 3.85e-04 6.13e-05h 12
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
920 5.5274877e+01 2.70e+01 3.99e+01 -1.0 7.02e+08 - 1.50e-01 5.29e-05h 12
921 5.5309940e+01 2.70e+01 4.00e+01 -1.0 3.60e+08 - 1.29e-03 8.91e-05h 13
922 5.5309987e+01 2.70e+01 3.99e+01 -1.0 7.02e+08 - 6.47e-04 5.29e-05h 12
923 5.5309380e+01 2.70e+01 3.98e+01 -1.0 7.12e+08 - 8.31e-03 5.17e-05h 12
924 5.7629043e+01 6.89e+05 1.28e+03 -1.0 6.81e+08 - 6.18e-04 1.13e-01w 1
925 5.7736878e+01 6.86e+05 1.72e+03 -1.0 4.41e+08 - 4.16e-01 4.96e-03w 1
926 7.3974358e+01 6.86e+05 2.57e+04 -1.0 3.04e+08 - 7.94e-01 3.54e-02w 1
927 5.5312003e+01 2.70e+01 3.98e+01 -1.0 2.80e+08 - 6.18e-04 1.10e-04h 10
928 5.5310892e+01 2.70e+01 3.71e+01 -1.0 7.19e+08 - 1.71e-01 5.10e-05h 12
929 5.5324450e+01 2.70e+01 3.71e+01 -1.0 4.69e+08 - 1.32e-03 5.07e-05h 13
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
930 5.5323528e+01 2.70e+01 3.71e+01 -1.0 7.15e+08 - 1.32e-03 5.14e-05h 12
931 5.5322243e+01 2.70e+01 3.70e+01 -1.0 7.21e+08 - 8.59e-03 5.08e-05h 12
932 5.5321683e+01 2.70e+01 3.70e+01 -1.0 7.07e+08 - 9.18e-04 5.21e-05h 12
933 5.5319971e+01 2.70e+01 4.18e+01 -1.0 7.27e+08 - 2.71e-01 5.01e-05h 12
934 5.5334487e+01 2.70e+01 4.20e+01 -1.0 4.60e+08 - 1.81e-03 5.25e-05h 13
935 5.5333499e+01 2.70e+01 4.26e+01 -1.0 7.13e+08 - 1.50e-03 5.15e-05h 12
936 5.5332071e+01 2.70e+01 4.57e+01 -1.0 7.20e+08 - 8.62e-03 5.07e-05h 12
937 5.2350226e+01 6.21e+05 7.99e+02 -1.0 7.17e+08 - 1.19e-03 1.05e-01w 1
938 5.2814646e+01 6.19e+05 5.00e+03 -1.0 4.27e+08 - 5.30e-01 5.83e-03w 1
939 2.1073021e+02 2.95e+06 1.91e+04 -1.0 3.09e+08 - 9.55e-01 3.56e-01w 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
940 5.5330807e+01 2.70e+01 4.61e+01 -1.0 2.25e+08 - 1.19e-03 5.11e-05h 11
941 5.5328745e+01 2.70e+01 1.69e+02 -1.0 7.30e+08 - 3.01e-01 4.97e-05h 12
942 5.5340973e+01 2.70e+01 1.70e+02 -1.0 4.80e+08 - 2.11e-03 4.85e-05h 13
943 5.5340471e+01 2.70e+01 1.73e+02 -1.0 7.02e+08 - 2.06e-03 5.25e-05h 12
944 5.5338605e+01 2.70e+01 1.86e+02 -1.0 7.09e+08 - 9.13e-03 1.04e-04h 11
945 5.5337689e+01 2.70e+01 1.88e+02 -1.0 7.07e+08 - 1.57e-03 5.19e-05h 12
946 5.5336149e+01 2.70e+01 8.84e+02 -1.0 7.17e+08 - 4.33e-01 5.09e-05h 12
947 5.5352840e+01 2.70e+01 8.92e+02 -1.0 4.43e+08 - 3.01e-03 5.65e-05h 13
948 5.5357306e+01 2.70e+01 9.10e+02 -1.0 6.30e+08 - 2.97e-03 6.16e-05h 12
949 5.5361256e+01 2.70e+01 9.79e+02 -1.0 6.36e+08 - 1.09e-02 6.07e-05h 12
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
950 6.3469526e+01 7.80e+05 7.22e+02 -1.0 6.34e+08 - 2.22e-03 1.25e-01w 1
951 9.1494386e+01 8.09e+05 1.43e+04 -1.0 3.09e+08 - 9.90e-01 6.25e-02w 1
952 4.8206184e+02 7.33e+06 2.27e+04 -1.0 2.62e+08 - 8.33e-01 6.78e-01w 1
953 5.5365363e+01 2.70e+01 9.93e+02 -1.0 5.28e+08 - 2.22e-03 6.11e-05h 11
954 5.5368932e+01 2.70e+01 3.56e+03 -1.0 6.40e+08 - 3.66e-01 6.01e-05h 12
955 5.5379603e+01 2.70e+01 3.60e+03 -1.0 4.14e+08 - 3.96e-03 3.24e-05h 14
956 5.5390965e+01 2.70e+01 3.85e+03 -1.0 4.87e+08 - 1.72e-02 4.71e-05h 13
957 5.5402543e+01 2.70e+01 1.95e+04 -1.0 4.85e+08 - 9.90e-01 4.75e-05h 13
958 5.5437349e+01 2.70e+01 1.95e+04 -1.0 2.93e+08 - 9.00e-03 7.34e-05h 14
959 5.5464865e+01 2.70e+01 2.01e+04 -1.0 3.20e+08 - 3.08e-02 6.22e-05h 14
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
960 5.5491907e+01 2.70e+01 2.10e+04 -1.0 3.21e+08 - 4.81e-02 6.14e-05h 14
961 5.5519808e+01 2.70e+01 2.12e+04 -1.0 3.18e+08 - 8.67e-03 6.29e-05h 14
962 5.5547281e+01 2.70e+01 4.11e+04 -1.0 3.20e+08 - 9.93e-01 6.21e-05h 14
963 3.4768279e+02 7.75e+06 1.66e+04 -1.0 2.90e+08 - 9.67e-03 6.02e-01w 1
964 2.6985413e+02 6.71e+06 6.38e+04 -1.0 5.12e+08 - 6.36e-01 1.08e-01w 1
965 2.6899993e+02 6.70e+06 4.47e+04 -1.0 2.95e+08 - 1.54e-01 1.75e-03w 1
966 5.5582339e+01 2.70e+01 4.13e+04 -1.0 1.15e+08 - 9.67e-03 7.35e-05h 13
967 5.5616784e+01 2.70e+01 4.24e+04 -1.0 3.02e+08 - 4.62e-02 7.45e-05h 14
968 5.5650801e+01 2.70e+01 6.96e+04 -1.0 3.03e+08 - 1.00e+00 7.37e-05h 14
969 5.5685845e+01 2.70e+01 6.99e+04 -1.0 2.98e+08 - 9.37e-03 7.49e-05h 14
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
970 5.5720873e+01 2.70e+01 7.19e+04 -1.0 3.01e+08 - 4.68e-02 7.54e-05h 14
971 5.5755831e+01 2.70e+01 7.24e+04 -1.0 3.01e+08 - 1.27e-02 7.54e-05h 14
972 5.5790676e+01 2.70e+01 1.16e+05 -1.0 3.01e+08 - 1.00e+00 7.52e-05h 14
973 5.5825767e+01 2.70e+01 1.17e+05 -1.0 2.96e+08 - 9.68e-03 7.48e-05h 14
974 5.5860830e+01 2.70e+01 1.20e+05 -1.0 3.00e+08 - 4.86e-02 7.55e-05h 14
975 5.5895826e+01 2.70e+01 1.93e+05 -1.0 3.00e+08 - 1.00e+00 7.54e-05h 14
976 3.4840961e+02 8.43e+06 1.23e+05 -1.0 2.99e+08 - 9.98e-03 6.18e-01w 1
977 2.9019863e+02 7.45e+06 9.41e+05 -1.0 3.94e+08 - 6.57e-01 1.01e-01w 1
978 8.4928662e+02 1.68e+07 1.15e+08 -1.0 6.38e+08 - 1.68e-02 1.00e+00w 1
979 5.5930908e+01 2.70e+01 1.94e+05 -1.0 1.80e+07 -4.5 9.98e-03 7.54e-05h 13
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
980 5.5965995e+01 2.71e+01 2.00e+05 -1.0 3.00e+08 - 4.77e-02 7.56e-05h 14
981 5.6001071e+01 2.71e+01 2.01e+05 -1.0 3.00e+08 - 1.27e-02 7.56e-05h 14
982 5.6036096e+01 2.71e+01 3.22e+05 -1.0 3.00e+08 - 1.00e+00 7.55e-05h 14
983 5.6071204e+01 2.72e+01 3.24e+05 -1.0 2.98e+08 - 1.00e-02 7.53e-05h 14
984 5.6106310e+01 2.72e+01 3.34e+05 -1.0 2.99e+08 - 4.96e-02 7.56e-05h 14
985 5.6141379e+01 2.73e+01 3.57e+05 -1.0 2.99e+08 - 1.18e-01 7.55e-05h 14
986 5.6176452e+01 2.73e+01 3.74e+05 -1.0 2.99e+08 - 7.75e-02 7.55e-05h 14
987 5.6211560e+01 2.73e+01 3.76e+05 -1.0 2.99e+08 - 1.17e-02 7.56e-05h 14
988 5.6246620e+01 2.74e+01 6.03e+05 -1.0 2.99e+08 - 1.00e+00 7.55e-05h 14
989 3.4907080e+02 8.39e+06 3.81e+05 -1.0 2.98e+08 - 1.05e-02 6.18e-01w 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
990 2.8667207e+02 7.35e+06 2.71e+06 -1.0 4.01e+08 - 6.49e-01 1.06e-01w 1
991 8.6964576e+02 1.73e+07 4.34e+08 -1.0 6.59e+08 - 1.34e-02 1.00e+00w 1
992 5.6281741e+01 2.74e+01 6.06e+05 -1.0 1.89e+07 -5.0 1.05e-02 7.55e-05h 13
993 5.6316867e+01 2.74e+01 6.24e+05 -1.0 2.99e+08 - 4.98e-02 7.56e-05h 14
994 5.6351969e+01 2.75e+01 6.38e+05 -1.0 2.98e+08 - 3.71e-02 7.56e-05h 14
995 5.6387064e+01 2.75e+01 1.02e+06 -1.0 2.98e+08 - 1.00e+00 7.55e-05h 14
996 5.6422196e+01 2.75e+01 1.03e+06 -1.0 2.98e+08 - 1.08e-02 7.55e-05h 14
997 5.6457323e+01 2.76e+01 1.07e+06 -1.0 2.98e+08 - 6.31e-02 7.56e-05h 14
998 5.6492442e+01 2.76e+01 1.09e+06 -1.0 2.98e+08 - 4.02e-02 7.56e-05h 14
999 5.6527556e+01 2.76e+01 1.75e+06 -1.0 2.98e+08 - 1.00e+00 7.56e-05h 14
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1000 5.6562698e+01 2.77e+01 1.76e+06 -1.0 2.97e+08 - 1.11e-02 7.56e-05h 14
1001 5.6597837e+01 2.77e+01 1.83e+06 -1.0 2.98e+08 - 6.87e-02 7.56e-05h 14
1002 3.4952577e+02 8.39e+06 1.19e+06 -1.0 2.97e+08 - 3.18e-02 6.19e-01w 1
1003 3.1036180e+02 7.78e+06 1.50e+07 -1.0 4.07e+08 - 6.42e-01 6.58e-02w 1
1004 7.3018011e+02 1.41e+07 1.48e+09 -1.0 5.26e+08 - 1.71e-02 9.95e-01w 1
1005 5.6632971e+01 2.77e+01 1.86e+06 -1.0 1.61e+07 -4.5 3.18e-02 7.56e-05h 13
1006 5.6668100e+01 2.78e+01 2.98e+06 -1.0 2.97e+08 - 1.00e+00 7.56e-05h 14
1007 5.6703252e+01 2.78e+01 3.00e+06 -1.0 2.97e+08 - 1.15e-02 7.56e-05h 14
1008 5.6738405e+01 2.78e+01 3.10e+06 -1.0 2.97e+08 - 5.70e-02 7.56e-05h 14
1009 5.6773552e+01 2.79e+01 3.17e+06 -1.0 2.97e+08 - 3.32e-02 7.56e-05h 14
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1010 5.6808693e+01 2.79e+01 5.06e+06 -1.0 2.97e+08 - 1.00e+00 7.56e-05h 14
1011 5.6843855e+01 2.79e+01 5.10e+06 -1.0 2.96e+08 - 1.18e-02 7.56e-05h 14
1012 5.6879018e+01 2.80e+01 5.28e+06 -1.0 2.96e+08 - 5.89e-02 7.56e-05h 14
1013 5.6914176e+01 2.80e+01 5.37e+06 -1.0 2.96e+08 - 3.06e-02 7.56e-05h 14
1014 5.6949329e+01 2.80e+01 8.59e+06 -1.0 2.96e+08 - 1.00e+00 7.56e-05h 14
1015 3.5015675e+02 8.35e+06 5.57e+06 -1.0 2.96e+08 - 1.22e-02 6.19e-01w 1
1016 3.0365672e+02 7.63e+06 5.73e+07 -1.0 4.16e+08 - 6.37e-01 7.65e-02w 1
1017 7.6030302e+02 1.50e+07 7.71e+09 -1.0 5.63e+08 - 1.35e-02 9.97e-01w 1
1018 5.6984499e+01 2.81e+01 8.65e+06 -1.0 1.83e+07 -4.1 1.22e-02 7.56e-05h 13
1019 5.7019672e+01 2.81e+01 8.95e+06 -1.0 2.96e+08 - 5.77e-02 7.56e-05h 14
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1020 5.7054840e+01 2.81e+01 9.11e+06 -1.0 2.96e+08 - 2.99e-02 7.56e-05h 14
1021 5.7090003e+01 2.82e+01 1.46e+07 -1.0 2.96e+08 - 1.00e+00 7.56e-05h 14
1022 5.7125181e+01 2.82e+01 1.47e+07 -1.0 2.95e+08 - 1.25e-02 7.56e-05h 14
1023 5.7160364e+01 2.82e+01 1.52e+07 -1.0 2.95e+08 - 5.81e-02 7.56e-05h 14
1024 5.7195541e+01 2.83e+01 1.54e+07 -1.0 2.95e+08 - 2.92e-02 7.56e-05h 14
1025 5.7230713e+01 2.83e+01 2.47e+07 -1.0 2.95e+08 - 1.00e+00 7.56e-05h 14
1026 5.7265899e+01 2.83e+01 2.49e+07 -1.0 2.95e+08 - 1.29e-02 7.56e-05h 14
1027 5.7301090e+01 2.84e+01 2.58e+07 -1.0 2.95e+08 - 5.84e-02 7.56e-05h 14
1028 3.5062173e+02 8.30e+06 1.68e+07 -1.0 2.95e+08 - 2.88e-02 6.20e-01w 1
1029 3.0177884e+02 7.56e+06 1.66e+08 -1.0 4.24e+08 - 6.33e-01 7.90e-02w 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1030 7.4707356e+02 1.44e+07 2.42e+10 -1.0 5.43e+08 - 1.24e-02 9.97e-01w 1
1031 5.7336275e+01 2.84e+01 2.62e+07 -1.0 1.75e+07 -3.7 2.88e-02 7.56e-05h 13
1032 5.7371456e+01 2.84e+01 4.19e+07 -1.0 2.95e+08 - 1.00e+00 7.56e-05h 14
1033 5.7406649e+01 2.85e+01 4.22e+07 -1.0 2.94e+08 - 1.33e-02 7.56e-05h 14
1034 5.7441848e+01 2.85e+01 4.37e+07 -1.0 2.94e+08 - 5.89e-02 7.56e-05h 14
1035 5.7477041e+01 2.85e+01 4.44e+07 -1.0 2.94e+08 - 2.86e-02 7.56e-05h 14
1036 5.7512229e+01 2.86e+01 7.10e+07 -1.0 2.94e+08 - 1.00e+00 7.56e-05h 14
1037 5.7547430e+01 2.86e+01 7.16e+07 -1.0 2.94e+08 - 1.36e-02 7.56e-05h 14
1038 5.7582635e+01 2.86e+01 7.41e+07 -1.0 2.94e+08 - 5.94e-02 7.56e-05h 14
1039 5.7617836e+01 2.87e+01 7.54e+07 -1.0 2.94e+08 - 2.85e-02 7.56e-05h 14
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1040 5.7653031e+01 2.87e+01 1.21e+08 -1.0 2.94e+08 - 1.00e+00 7.56e-05h 14
1041 3.5113055e+02 8.26e+06 7.86e+07 -1.0 2.93e+08 - 1.40e-02 6.20e-01w 1
1042 2.9375033e+02 7.39e+06 6.42e+08 -1.0 4.33e+08 - 6.32e-01 9.09e-02w 1
1043 4.5802362e+02 9.96e+06 9.58e+10 -1.0 3.86e+08 - 1.16e-02 9.98e-01w 1
1044 3.3595036e+02 2.50e+03 1.53e+09 -1.0 3.94e+07 - 4.71e-02 1.00e+00h 1
1045 3.3583554e+02 3.11e-01 8.97e+06 -1.0 4.09e+04 -4.2 9.94e-01 1.00e+00h 1
1046 3.3583394e+02 1.01e-04 1.13e+04 -1.0 4.39e+01 -4.6 9.89e-01 1.00e+00h 1
1047 3.3583393e+02 3.58e-07 3.10e-01 -1.0 1.61e-02 -5.1 9.90e-01 1.00e+00h 1
1048 3.3284606e+02 3.31e-01 2.64e+01 -1.0 6.45e+06 - 9.80e-01 1.00e+00f 1
1049 3.3284564e+02 1.26e+02 5.31e+00 -1.0 4.61e+01 -5.6 1.00e+00 1.00e+00H 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1050 3.3287267e+02 1.22e+02 5.16e+00 -1.0 2.24e+01 -6.1 1.00e+00 3.12e-02h 6
1051 3.3287427e+02 1.22e+02 5.15e+00 -1.0 2.11e+01 -6.5 1.00e+00 1.95e-03h 10
1052 3.3287467e+02 1.22e+02 5.15e+00 -1.0 3.17e+01 -7.0 1.00e+00 4.88e-04h 12
1053 3.3287486e+02 1.22e+02 5.15e+00 -1.0 9.52e+01 -7.5 1.00e+00 2.44e-04h 13
1054 3.3287491e+02 1.22e+02 5.15e+00 -1.0 2.86e+02 -8.0 1.00e+00 6.10e-05h 15
1055 3.3367782e+02 4.41e+01 1.42e+01 -1.0 8.66e+02 -8.4 1.00e+00 1.00e+00h 1
1056 3.3346366e+02 2.29e+00 1.27e+00 -1.0 1.95e+03 -8.9 1.00e+00 1.00e+00h 1
1057 3.3335169e+02 5.02e-02 3.36e-02 -1.0 5.77e+03 -9.4 1.00e+00 1.00e+00h 1
1058 3.3300763e+02 3.66e-01 8.35e-01 -2.5 1.79e+04 -9.9 9.23e-01 1.00e+00f 1
1059 3.3184539e+02 7.81e+00 9.12e-01 -2.5 5.79e+04 -10.4 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1060 3.2820788e+02 8.12e+01 7.28e+00 -2.5 1.73e+05 -10.8 1.00e+00 1.00e+00f 1
1061 3.0911770e+02 3.65e+03 1.36e+02 -2.5 1.83e+07 - 1.00e+00 1.00e+00f 1
1062 3.2815512e+02 3.55e+02 1.46e+01 -2.5 3.46e+05 -11.3 1.00e+00 1.00e+00h 1
1063 3.1276588e+02 3.11e+03 5.54e+01 -2.5 5.80e+06 - 2.89e-01 4.45e-01f 1
1064 3.3365584e+02 6.37e+02 1.97e+01 -2.5 1.71e+05 -10.9 1.00e+00 1.00e+00h 1
1065 3.1447219e+02 1.64e+03 1.61e+02 -2.5 8.86e+05 -11.4 5.39e-01 1.00e+00f 1
1066 3.2139869e+02 4.38e+02 1.99e+02 -2.5 1.08e+06 -11.8 1.00e+00 1.00e+00h 1
1067 3.2507697e+02 2.56e+02 4.18e+01 -2.5 6.04e+05 -11.4 1.00e+00 1.00e+00h 1
1068 3.2320339e+02 9.40e+01 1.51e+01 -2.5 1.74e+05 -11.9 1.00e+00 1.00e+00h 1
1069 3.1263875e+02 1.49e+00 8.74e+02 -2.5 3.68e+07 - 1.00e+00 1.00e+00f 1
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1070 3.0458683e+02 1.52e+00 1.35e+02 -2.5 1.78e+07 - 1.00e+00 1.00e+00f 1
1071 3.2100236e+02 2.88e-01 2.11e+01 -2.5 4.90e+06 - 1.00e+00 1.00e+00h 1
1072 3.2169820e+02 1.23e-02 9.34e-01 -2.5 1.11e+06 - 1.00e+00 1.00e+00h 1
1073 3.2172657e+02 2.59e-05 2.12e-03 -2.5 5.38e+04 - 1.00e+00 1.00e+00h 1
1074 3.2172663e+02 2.26e-05 3.39e-06 -3.8 1.65e+03 - 1.00e+00 1.00e+00h 1
1075 3.2172663e+02 2.38e-07 9.83e-09 -5.7 8.47e+01 - 1.00e+00 1.00e+00h 1
1076 3.2172663e+02 4.77e-07 3.02e-12 -7.0 5.40e+01 - 1.00e+00 1.00e+00h 1
Number of Iterations....: 1076
(scaled) (unscaled)
Objective...............: 3.2172662923377101e+02 3.2172662923377101e+02
Dual infeasibility......: 3.0248858973180859e-12 3.0248858973180859e-12
Constraint violation....: 2.9802322387695312e-08 4.7683715820312500e-07
Complementarity.........: 9.0909154426932121e-08 9.0909154426932121e-08
Overall NLP error.......: 9.0909154426932121e-08 4.7683715820312500e-07
Number of objective function evaluations = 3129
Number of objective gradient evaluations = 255
Number of equality constraint evaluations = 3133
Number of inequality constraint evaluations = 3133
Number of equality constraint Jacobian evaluations = 1112
Number of inequality constraint Jacobian evaluations = 1112
Number of Lagrangian Hessian evaluations = 1076
Total CPU secs in IPOPT (w/o function evaluations) = 3.193
Total CPU secs in NLP function evaluations = 370.594
EXIT: Optimal Solution Found.
Display Solution
[11]:
#get value of surrogate inputs
x = [value(m.pmax),
value(m.pmin_multi),
value(m.ramp_multi),
value(m.min_up_time),
value(m.min_dn_multi),
value(m.marg_cst),
value(m.no_load_cst),
value(m.startup_cst)
]
print("Market Inputs:", x)
#calculate revenues, costs, etc...
optimal_objective = -value(m.obj)
zone_hours = [value(m.zone_off.zone_hours)]
scaled_zone_hours = [value(m.zone_off.scaled_zone_hours)]
op_cost = [value(m.zone_off.fs.operating_cost)]
op_expr = value(m.zone_off.fs.operating_cost) # in dollars [$]
for zone in m.op_zones:
zone_hours.append(value(zone.zone_hours))
scaled_zone_hours.append(value(zone.scaled_zone_hours))
op_cost.append(value(zone.fs.operating_cost))
op_expr += value(zone.scaled_zone_hours)*value(zone.fs.operating_cost)
#more calculations of revenue and cost
revenue_per_year = value(m.revenue)
cap_expr = value(m.cap_fs.fs.capital_cost)/capital_payment_years
startup_expr = value(m.startup_expr)
total_cost = plant_lifetime*op_expr/1e6 + capital_payment_years*cap_expr
total_revenue = plant_lifetime*revenue_per_year
print("Annual Revenue [MM$]: ", revenue_per_year)
print("Capital Cost [MM$]: ",cap_expr)
print("Annual Startup Cost [MM$]: ", startup_expr)
print("Total Cost (20 years): ", total_cost)
print("Net Revenue (20 years)", optimal_objective)
Market Inputs: [259.91910248189976, 0.4423848332093146, 0.6392249722616824, 4.0, 1.0, 17.03405649932845, 1.0, 61.09068702]
Annual Revenue [MM$]: 42.350552833052404
Capital Cost [MM$]: 109.78197558234294
Annual Startup Cost [MM$]: 2.4502773280683162e-11
Total Cost (20 years): 1166.2990079394763
Net Revenue (20 years) -321.72662923377106