OptimizationReporter

Computes objective function, gradient and contains reporters for communicating between optimizeSolve and subapps

Overview

This optimization reporter provides a basic interface for inverse optimization with measurement data where all data is provided in the input file or from csv file. The objective function is defined as:

var element = document.getElementById("moose-equation-dd6a790e-a702-4487-a93a-275bdd521cc0");katex.render("f(\\mathbf{u}, \\mathbf{p}) = \\frac{1}{2}\\sum^N_i \\left(u_i - \\tilde u_i\\right)^2 +R(\\mathbf{u}, \\mathbf{p}),", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});

where $var element = document.getElementById("moose-equation-502eae8a-fdb2-455d-9d86-04bbcc56f755");katex.render("\\mathbf{u}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ represents the simulation solution with $var element = document.getElementById("moose-equation-0d4525a1-6efe-46ea-a1ba-787142babce9");katex.render("u_i", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ being the solution value at measurement location $var element = document.getElementById("moose-equation-97ae0d53-62a8-430d-b030-8184d9cc117b");katex.render("i", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ . $var element = document.getElementById("moose-equation-9488f322-d358-4eb1-8f81-5e33b43016b5");katex.render("\\tilde u_i", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the measurement value at location $var element = document.getElementById("moose-equation-af658084-80b1-482a-a312-4b613d98429f");katex.render("i", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ . $var element = document.getElementById("moose-equation-ee8db4bc-861d-4ce0-81de-dfd149cd4c4a");katex.render("\\mathbf{p}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the vector of parameters being optimized that the simulation solution depends on. $var element = document.getElementById("moose-equation-ce1a1414-6eba-4036-a050-679ab1ad02bd");katex.render("R(\\mathbf{u}, \\mathbf{p})", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ provides regularization of the objective function. We currently only support the Tikhonov regularization given by

var element = document.getElementById("moose-equation-f1e5650e-10c3-42c6-8a3c-393edb08cafb");katex.render("R(\\mathbf{u}, \\mathbf{p}) = \\frac{\\alpha}{2}\\|\\mathbf{p}\\|^2,", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});

where $var element = document.getElementById("moose-equation-21bdd623-a961-419b-a7f0-6198128eac40");katex.render("\\alpha", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the Tikhonov coefficient specified in the input file by "tikhonov_coeff".

Measurement Data

The measurement locations and values are defined either by a CSV file or are specified at input.

The CSV file is specified with the "measurement_file". This file must be at least four columns: x-coordinates indicated by "file_xcoord", y-coordinates indicated by "file_ycoord", z-coordinates inidicated by "file_zcoord", and values indicated by "file_value". A column for time coordinates can also be included and indicated by "file_time".

Additionally, locations and values can be specified at input using "measurement_points" for x-y-z coordinates, "measurement_times" for time coordinates, and "measurement_values" for the values.

Optimization Parameters

OptimizationReporter is also responsible for creating parameter vector(s) for optimization, setting the initial condition for the optimization, and setting parameter bounds. Although the Optimize executioner holds a single vector for parameter values, this vector can be split into groups of parameters. This is done by specifying a name for each group with "parameter_names" and the number of parameters in each group with "num_values". The total number of parameters is ultimately defined by the sum of "num_values". The initial condition for the optimization can then be defined with "initial_condition", where a vector of data must defined for each group. This vector an be a single value in which case all parameters in that group are set to that value or a value can be set for every parameter in that group. The lower and upper bounds for the parameters can then specified by "lower_bounds" and "upper_bounds", respectively. The bounds follow the same input format rules as the initial_condtion. If no initial conditions are provided, the parameters are initialized with 0. Default values for upper_bounds and lower_bounds are std::numeric<Real>::max() and std::numeric<Real>::lower(), respectively. These bounds are only applied if a bounded optimization algorithm is used.

Declared Data

OptimizationReporter declares a number of vector reporter values that can be read by other objects and/or transferred to sub-applications with MultiAppReporterTransfer. Table 1 lists each of these vectors. These vectors can be supplied or transferred by specifying OptimizationReporter/<Vector Name> at input.

Table 1: List of vectors declared by OptimizationReporter. $var element = document.getElementById("moose-equation-494bc24d-09d7-43fb-936f-3e17dc1aa0da");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the number of measurement points.

Description	Vector Name	Size
$var element = document.getElementById("moose-equation-994196d0-8c89-4973-8b1f-a685abc1951e");katex.render("x", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ -coordinate	`measurement_xcoord`	$var element = document.getElementById("moose-equation-a85ec705-9b11-4b4f-b8e5-b717a3f14895");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$
$var element = document.getElementById("moose-equation-5f843572-6a2e-41b9-8789-181fbe600e23");katex.render("y", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ -coordinate	`measurement_ycoord`	$var element = document.getElementById("moose-equation-3d78bdf9-a367-4661-81f6-f0a8eabcf07f");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$
$var element = document.getElementById("moose-equation-f80fa588-61c2-4380-aeb8-cd40996f5a6f");katex.render("z", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ -coordinate	`measurement_zcoord`	$var element = document.getElementById("moose-equation-b0f0decf-ffe3-4214-9413-16da3a669f01");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$
Time coordinate	`measurement_time`	$var element = document.getElementById("moose-equation-114df022-2aca-4b2f-be02-9cb96b0f8d54");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$
Measured values	`measurement_values`	$var element = document.getElementById("moose-equation-94ea8311-ef9a-4c99-b30d-82f6f7b4ed99");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$
Simulation values	`simulation_values`	$var element = document.getElementById("moose-equation-796fdb23-3f66-401e-b3f1-fe82c598307c");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$
$var element = document.getElementById("moose-equation-6bfef8e1-de87-4c41-b8a4-afa148666de8");katex.render("u_i - \\tilde{u}_i", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	`misfit_values`	$var element = document.getElementById("moose-equation-b6464ec1-8898-484b-9c88-f5d69d52b87f");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$
Values of parameter group $var element = document.getElementById("moose-equation-3b7a3e14-0b01-4a3c-9a42-a1b1c415616c");katex.render("g", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	"parameter_names" $var element = document.getElementById("moose-equation-58f8aec0-ffe3-4956-8cf1-8795b6c76d7d");katex.render("_g", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	"num_values" $var element = document.getElementById("moose-equation-2481412c-d82c-4f33-b994-67b5ca2b400a");katex.render("_g", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$
Parameter Gradient of parameter group $var element = document.getElementById("moose-equation-071c83ca-3696-4d90-bcb2-8ee50f1d0a39");katex.render("g", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	`grad_`"parameter_names" $var element = document.getElementById("moose-equation-baa4d212-ccda-4722-9cfe-1c2e1204e936");katex.render("_g", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	"num_values" $var element = document.getElementById("moose-equation-e941631a-053d-44f0-b15f-b2a72f6e8d18");katex.render("_g", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$

Example Input Syntax

The following example creates three groups of parameters named p1, p2, and p3 with two, four, and six parameters, respectively. There are several ways to define the initial conditions and upper and lower bounds. Starting with the initial_condition, all parameters are specified for the p1 parameter group, only one parameter is defined for the p2 group making all parameter values in this group set to this value, and all parameters are defined for the p3 group. A similar mix of defining one or all of the lower and upper bounds are demonstrated in this input file.

It sets an initial condition and upper and lower bounds for every parameter. The measurement data is taken from a CSV file where the x-coordinate, y-coordinate, z-coordinate, and measurement value columns are specified with the header names of the file.

[OptimizationReporter]
  type = GeneralOptimization
  objective_name = objective_value
  parameter_names = 'p1 p2 p3'
  num_values = '2 4 6'
  initial_condition = '1 2; 3; 7 8 9 10 11 12'
  upper_bounds = '101; 103; 107 108 109 110 111 112'
  lower_bounds = '-1 -2; -3 -4 -5 -6; -7'

  outputs = out
[]

measured_value,coordx,y,z
11,0.1,0.2,0.3
12,0.4,0.5,0.6
13,0.7,0.8,0.9
14,1.0,1.1,1.2

The measurement data can also be given directly, like in the following example.

[OptimizationReporter]
  type = GeneralOptimization
  objective_name = objective_value
  parameter_names = 'parameter_results'
  num_values = '3'
[]

Input Parameters

num_valuesNumber of parameter values associated with each parameter group in 'parameter_names'.
C++ Type:std::vector<unsigned long>
Unit:(no unit assumed)
Controllable:No
Description:Number of parameter values associated with each parameter group in 'parameter_names'.
parameter_namesList of parameter names, one for each group of parameters.
C++ Type:std::vector<ReporterValueName>
Unit:(no unit assumed)
Controllable:No
Description:List of parameter names, one for each group of parameters.

Required Parameters

equality_namesList of equality names.
C++ Type:std::vector<ReporterValueName>
Unit:(no unit assumed)
Controllable:No
Description:List of equality names.
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Unit:(no unit assumed)
Options:XFEM_MARK, FORWARD, ADJOINT, HOMOGENEOUS_FORWARD, ADJOINT_TIMESTEP_BEGIN, ADJOINT_TIMESTEP_END, NONE, INITIAL, LINEAR, NONLINEAR_CONVERGENCE, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM
Controllable:No
Description:The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
inequality_namesList of inequality names.
C++ Type:std::vector<ReporterValueName>
Unit:(no unit assumed)
Controllable:No
Description:List of inequality names.
initial_conditionInitial conditions for each parameter. A vector is given for each parameter group. A single value can be given for each group and all parameters in that group will be set to that value. The default value is 0.
C++ Type:std::vector<std::vector<double>>
Unit:(no unit assumed)
Controllable:No
Description:Initial conditions for each parameter. A vector is given for each parameter group. A single value can be given for each group and all parameters in that group will be set to that value. The default value is 0.
lower_boundsLower bound for each parameter. A vector is given for each parameter group. A single value can be given for each group and all parameters in that group will be set to that value
C++ Type:std::vector<std::vector<double>>
Unit:(no unit assumed)
Controllable:No
Description:Lower bound for each parameter. A vector is given for each parameter group. A single value can be given for each group and all parameters in that group will be set to that value
objective_nameName of reporter value defining the objective.
C++ Type:ReporterValueName
Unit:(no unit assumed)
Controllable:No
Description:Name of reporter value defining the objective.
prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
tikhonov_coeff0Coefficient for Tikhonov Regularization.
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Coefficient for Tikhonov Regularization.
upper_boundsUpper bound for each parameter. A vector is given for each parameter group. A single value can be given for each group and all parameters in that group will be set to that value
C++ Type:std::vector<std::vector<double>>
Unit:(no unit assumed)
Controllable:No
Description:Upper bound for each parameter. A vector is given for each parameter group. A single value can be given for each group and all parameters in that group will be set to that value
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
variableVector of variable names to sample at measurement points.
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:Vector of variable names to sample at measurement points.
variable_weight_namesVector of weight reporter names that will create a reporter to transfer weights into. The ordering of these weight reporter names corresponds to the ordering used in variable.
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
Controllable:No
Description:Vector of weight reporter names that will create a reporter to transfer weights into. The ordering of these weight reporter names corresponds to the ordering used in variable.

Optional Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Unit:(no unit assumed)
Controllable:Yes
Description:Set the enabled status of the MooseObject.
execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Unit:(no unit assumed)
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup
outputsVector of output names where you would like to restrict the output of variables(s) associated with this object
C++ Type:std::vector<OutputName>
Unit:(no unit assumed)
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

file_timetimetime column name from csv file being read in.
Default:time
C++ Type:std::string
Unit:(no unit assumed)
Controllable:No
Description:time column name from csv file being read in.
file_valuevaluemeasurement value column name from csv file being read in.
Default:value
C++ Type:std::string
Unit:(no unit assumed)
Controllable:No
Description:measurement value column name from csv file being read in.
file_variable_weightsvariable weight column names from csv file being read in.
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
Controllable:No
Description:variable weight column names from csv file being read in.
file_xcoordxx coordinate column name from measurement_file csv being read in.
Default:x
C++ Type:std::string
Unit:(no unit assumed)
Controllable:No
Description:x coordinate column name from measurement_file csv being read in.
file_ycoordyy coordinate column name from csv file being read in.
Default:y
C++ Type:std::string
Unit:(no unit assumed)
Controllable:No
Description:y coordinate column name from csv file being read in.
file_zcoordzz coordinate column name from csv file being read in.
Default:z
C++ Type:std::string
Unit:(no unit assumed)
Controllable:No
Description:z coordinate column name from csv file being read in.
measurement_fileCSV file with measurement value and coordinates (value, x, y, z).
C++ Type:FileName
Unit:(no unit assumed)
Controllable:No
Description:CSV file with measurement value and coordinates (value, x, y, z).

File Measurement Data Parameters

measurement_pointsPoint locations corresponding to each measurement value
C++ Type:std::vector<libMesh::Point>
Unit:(no unit assumed)
Controllable:No
Description:Point locations corresponding to each measurement value
measurement_timesTimes corresponding to each measurement value
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Times corresponding to each measurement value
measurement_valuesMeasurement values collected from locations given by measurement_points
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Measurement values collected from locations given by measurement_points

Input Measurement Data Parameters

(modules/optimization/test/tests/optimizationreporter/optimization_reporter_base/optRep_fromCsv_mixBounds.i)

[Mesh]
  type = GeneratedMesh
  dim = 1
[]

[Variables]
  [u]
  []
[]

[Kernels]
  [null]
    type = NullKernel
    variable = u
  []
[]

[OptimizationReporter]
  type = GeneralOptimization
  objective_name = objective_value
  parameter_names = 'p1 p2 p3'
  num_values = '2 4 6'
  initial_condition = '1 2; 3; 7 8 9 10 11 12'
  upper_bounds = '101; 103; 107 108 109 110 111 112'
  lower_bounds = '-1 -2; -3 -4 -5 -6; -7'

  outputs = out
[]

[UserObjects]
  [optReporterTester]
    type = OptimizationReporterTest
    values_to_set_parameters_to = '10 20 30 40 50 60 70 80 90 100 110 120'
    expected_lower_bounds = '-1 -2 -3 -4 -5 -6 -7 -7 -7 -7 -7 -7'
    expected_upper_bounds = '101 101 103 103 103 103 107 108 109 110 111 112'
  []
[]

[Executioner]
  type = Steady
[]

[Outputs]
  [out]
    type = JSON
    execute_system_information_on = none
  []
[]

(modules/optimization/test/tests/optimizationreporter/point_loads/main.i)

# DO NOT CHANGE THIS TEST
# this test is documented as an example in forceInv_pointLoads.md
# if this test is changed, the figures will need to be updated.
measurement_points = '0.5 0.28 0
   0.5 0.6 0
   0.5 0.8 0
   0.5 1.1 0'
measurement_values = '293 304 315 320'
[Optimization]
[]

[OptimizationReporter]
  type = GeneralOptimization
  objective_name = objective_value
  parameter_names = 'parameter_results'
  num_values = '3'
[]

[Reporters]
  [main]
    type = OptimizationData
    measurement_points = ${measurement_points}
    measurement_values = ${measurement_values}
  []
[]

[Executioner]
  type = Optimize
  tao_solver = taonls
  petsc_options_iname = '-tao_gttol -tao_max_it -tao_nls_pc_type -tao_nls_ksp_type'
  petsc_options_value = '1e-5 10 none cg'
  verbose = true
[]

[MultiApps]
  [forward]
    type = FullSolveMultiApp
    input_files = forward.i
    execute_on = "FORWARD"
  []
  [adjoint]
    type = FullSolveMultiApp
    input_files = adjoint.i
    execute_on = "ADJOINT"
  []
  [homogeneousForward]
    type = FullSolveMultiApp
    input_files = forward_homogeneous.i
    execute_on = "HOMOGENEOUS_FORWARD"
  []
[]

[Transfers]
  # FORWARD transfers
  [toForward_measument]
    type = MultiAppReporterTransfer
    to_multi_app = forward
    from_reporters = 'main/measurement_xcoord
                      main/measurement_ycoord
                      main/measurement_zcoord
                      main/measurement_time
                      main/measurement_values
                      OptimizationReporter/parameter_results'
    to_reporters = 'measure_data/measurement_xcoord
                    measure_data/measurement_ycoord
                    measure_data/measurement_zcoord
                    measure_data/measurement_time
                    measure_data/measurement_values
                    point_source/value'
  []
  [fromForward]
    type = MultiAppReporterTransfer
    from_multi_app = forward
    # Note: We are transferring the misfit values into main misfit
    from_reporters = 'measure_data/objective_value measure_data/misfit_values'
    to_reporters = 'OptimizationReporter/objective_value main/misfit_values'
  []

  # ADJOINT transfers
  #NOTE:  the adjoint variable we are transferring is actually the gradient
  [toAdjoint]
    type = MultiAppReporterTransfer
    to_multi_app = adjoint
    from_reporters = 'main/measurement_xcoord
                      main/measurement_ycoord
                      main/measurement_zcoord
                      main/measurement_time
                      main/misfit_values'
    to_reporters = 'misfit/measurement_xcoord
                    misfit/measurement_ycoord
                    misfit/measurement_zcoord
                    misfit/measurement_time
                    misfit/misfit_values'
  []
  [fromAdjoint]
    type = MultiAppReporterTransfer
    from_multi_app = adjoint
    from_reporters = 'gradient/adjoint'
    to_reporters = 'OptimizationReporter/grad_parameter_results'
  []

  # HESSIAN transfers.  Same as forward.
  [toHomoForward]
    type = MultiAppReporterTransfer
    multi_app = homogeneousForward
    direction = to_multiapp
    from_reporters = 'main/measurement_xcoord
                      main/measurement_ycoord
                      main/measurement_zcoord
                      main/measurement_time
                      main/measurement_values
                      OptimizationReporter/parameter_results'
    to_reporters = 'measure_data/measurement_xcoord
                    measure_data/measurement_ycoord
                    measure_data/measurement_zcoord
                    measure_data/measurement_time
                    measure_data/measurement_values
                    point_source/value'
  []
  [fromHomoForward]
    type = MultiAppReporterTransfer
    multi_app = homogeneousForward
    direction = from_multiapp
    # Note: We are transferring the simulation values into misfit
    # this has to be done when using general opt and homogenous forward.
    from_reporters = 'measure_data/simulation_values'
    to_reporters = 'main/misfit_values'
  []
[]

[Reporters]
  [optInfo]
    type = OptimizationInfo
  []
[]

[Outputs]
  csv = true
[]

Overview
Measurement Data
Optimization Parameters
Declared Data
Example Input Syntax
Input Parameters

OptimizationReporter

Overview

Measurement Data

Optimization Parameters

Declared Data

Example Input Syntax

Input Parameters

Required Parameters

Optional Parameters

Advanced Parameters

File Measurement Data Parameters

Input Measurement Data Parameters

tikhonov_coeff

measurement_file

file_xcoord

file_ycoord

file_zcoord

file_value

file_time

measurement_points

measurement_times

measurement_values

parameter_names

num_values

num_values

initial_condition

lower_bounds

upper_bounds

parameter_names

num_values

parameter_names

num_values

(modules/optimization/test/tests/optimizationreporter/optimization_reporter_base/optRep_fromCsv_mixBounds.i)

(modules/optimization/test/tests/optimizationreporter/optimization_reporter_base/measurementData.csv)

(modules/optimization/test/tests/optimizationreporter/point_loads/main.i)