44 #ifndef ROL_ABSOLUTEVALUE_HPP
45 #define ROL_ABSOLUTEVALUE_HPP
49 #include "ROL_ParameterList.hpp"
74 Real param = parlist.get(
"Smoothing Parameter",1.);
75 param_ = 1./((param > 0.) ? param : 1.);
76 std::string type = parlist.get(
"Absolute Value Approximation",
"true");
92 ROL_TEST_FOR_EXCEPTION(
true, std::invalid_argument,
93 ">>> ERROR (ROL::AbsoluteValue): Absolute value approximation not defined!");
100 Real output = 0.0, e = 0.0;
101 if ( std::abs(
param_) > ROL_EPSILON<Real>() ) { e = 0.5/
param_; }
104 if ( input < -e ) { region = -1; }
105 else if ( input > e ) { region = 1; }
107 if ( deriv == 0 ) { output = std::abs(input); }
108 else if ( deriv == 1 ) { output = (input < 0.0 ? -1.0 : 1.0); }
109 else if ( deriv == 2 ) {
110 if ( region == -1 ) { output = 0.0; }
111 else if ( region == 0 ) { output = e*std::exp( -1.0/(1.0 - std::pow(std::abs(input/e),2.0)) ); }
112 else if ( region == 1 ) { output = 0.0; }
119 if ( deriv == 0 ) { output = std::sqrt(input*input + 1.0/
param_); }
120 else if ( deriv == 1 ) { output = input/std::sqrt(input*input+1.0/
param_); }
121 else if ( deriv == 2 ) { output = (1.0/
param_)/std::pow(input*input+1.0/
param_,1.5); }
127 if ( deriv == 0 ) { output = input*input/std::sqrt(input*input + 1.0/
param_); }
128 else if ( deriv == 1 ) { output = (2.0/
param_*input+std::pow(input,3.0)) /
129 std::pow(input*input+1.0/
param_,1.5); }
130 else if ( deriv == 2 ) { output = ((2.0/
param_-input*input)/
param_) /
131 std::pow(input*input+1.0/
param_,2.5); }
136 Real output = 0.0, e = 1.0;
137 if ( std::abs(
param_) > ROL_EPSILON<Real>() ) { e = 0.5/
param_; }
140 if ( input < -e ) { region = -1; }
141 else if ( input > e ) { region = 1; }
144 if ( std::abs(region) == 1 ) { output = std::abs(input); }
145 else if ( region == 0 ) { output = 1.875*std::pow(input*e,2.0) -
146 1.25 *std::pow(input*e,4.0) +
147 0.375*std::pow(input*e,6.0); }
149 else if ( deriv == 1 ) {
150 if ( std::abs(region) == 1 ) { output = (input < 0.0 ? -1.0 : 1.0); }
151 else if ( region == 0 ) { output = e*2.0*1.875*input*e -
152 e*4.0*1.25 *std::pow(input*e,3.0) +
153 e*6.0*0.375*std::pow(input*e,5.0); }
155 else if ( deriv == 2 ) {
156 if ( std::abs(region) == 1 ) { output = 0.0; }
157 else if ( region == 0 ) { output = e* 2.0*1.875 -
158 e*e*12.0*1.25 *std::pow(input*e,2.0) +
159 e*e*30.0*0.375*std::pow(input*e,4.0); }
Real sqrtd_absolute_value(Real input, int deriv)
Contains definitions of custom data types in ROL.
Real c2_absolute_value(Real input, int deriv)
AbsoluteValue(Real param=1., EAbsoluteValue eav=ABSOLUTEVALUE_TRUE)
Real sqrt_absolute_value(Real input, int deriv)
AbsoluteValue(ROL::ParameterList &parlist)
Real evaluate(Real input, int deriv)
Real true_absolute_value(Real input, int deriv)