Real TaylorApproximation::value(const Variables& vars)
{
short bdo = sharedDataRep->buildDataOrder;
if (bdo == 1)
return approxData.anchor_function();
else { // build up approx value from constant and derivative terms
Real approx_val = (bdo & 1) ? approxData.anchor_function() : 0.;
if (bdo & 6) {
const RealVector& x = vars.continuous_variables();
const RealVector& x0 = approxData.anchor_continuous_variables();
const RealVector& grad = approxData.anchor_gradient();
const RealSymMatrix& hess = approxData.anchor_hessian();
size_t num_v = sharedDataRep->numVars;
for (size_t i=0; i<num_v; i++) {
Real dist_i = x[i] - x0[i];
if (bdo & 2) // include gradient terms
approx_val += grad[i] * dist_i;
if (bdo & 4) // include Hessian terms
for (size_t j=0; j<num_v; j++)
approx_val += dist_i * hess(i,j) * (x[j] - x0[j])/2.;
}
}
return approx_val;
}
}
void Approximation::
add(const Variables& vars, bool anchor_flag, bool deep_copy)
{
if (approxRep)
approxRep->add(vars, anchor_flag, deep_copy);
else { // not virtual: all derived classes use following definition
// Approximation does not know about view mappings; therefore, take the
// simple approach of matching up active or all counts with numVars.
size_t num_v = sharedDataRep->numVars;
if (vars.cv() + vars.div() + vars.drv() == num_v)
add(vars.continuous_variables(), vars.discrete_int_variables(),
vars.discrete_real_variables(), anchor_flag, deep_copy);
else if (vars.acv() + vars.adiv() + vars.adrv() == num_v)
add(vars.all_continuous_variables(), vars.all_discrete_int_variables(),
vars.all_discrete_real_variables(), anchor_flag, deep_copy);
/*
else if (vars.cv() == num_v) { // compactMode does not affect vars
IntVector empty_iv; RealVector empty_rv;
add(vars.continuous_variables(), empty_iv, empty_rv,
anchor_flag, deep_copy);
}
else if (vars.acv() == num_v) { // potential conflict with cv/div/drv
IntVector empty_iv; RealVector empty_rv;
add(vars.all_continuous_variables(), empty_iv, empty_rv,
anchor_flag, deep_copy);
}
*/
else {
Cerr << "Error: variable size mismatch in Approximation::add()."
<< std::endl;
abort_handler(-1);
}
}
}
void SharedSurfpackApproxData::
vars_to_realarray(const Variables& vars, RealArray& ra)
{
// check incoming vars for correct length (active or all views)
if (vars.cv() + vars.div() + vars.drv() == numVars)
merge_variable_arrays(vars.continuous_variables(),
vars.discrete_int_variables(),
vars.discrete_real_variables(), ra);
else if (vars.acv() + vars.adiv() + vars.adrv() == numVars)
merge_variable_arrays(vars.all_continuous_variables(),
vars.all_discrete_int_variables(),
vars.all_discrete_real_variables(), ra);
else {
Cerr << "Error: bad parameter set length in SharedSurfpackApproxData::"
<< "vars_to_realarray()." << std::endl;
abort_handler(-1);
}
}
