66 wts_[0] = 0.1527533871307258;
pts_[0] = -0.0765265211334973;
67 wts_[1] = 0.1527533871307258;
pts_[1] = 0.0765265211334973;
68 wts_[2] = 0.1491729864726037;
pts_[2] = -0.2277858511416451;
69 wts_[3] = 0.1491729864726037;
pts_[3] = 0.2277858511416451;
70 wts_[4] = 0.1420961093183820;
pts_[4] = -0.3737060887154195;
71 wts_[5] = 0.1420961093183820;
pts_[5] = 0.3737060887154195;
72 wts_[6] = 0.1316886384491766;
pts_[6] = -0.5108670019508271;
73 wts_[7] = 0.1316886384491766;
pts_[7] = 0.5108670019508271;
74 wts_[8] = 0.1181945319615184;
pts_[8] = -0.6360536807265150;
75 wts_[9] = 0.1181945319615184;
pts_[9] = 0.6360536807265150;
76 wts_[10] = 0.1019301198172404;
pts_[10] = -0.7463319064601508;
77 wts_[11] = 0.1019301198172404;
pts_[11] = 0.7463319064601508;
78 wts_[12] = 0.0832767415767048;
pts_[12] = -0.8391169718222188;
79 wts_[13] = 0.0832767415767048;
pts_[13] = 0.8391169718222188;
80 wts_[14] = 0.0626720483341091;
pts_[14] = -0.9122344282513259;
81 wts_[15] = 0.0626720483341091;
pts_[15] = 0.9122344282513259;
82 wts_[16] = 0.0406014298003869;
pts_[16] = -0.9639719272779138;
83 wts_[17] = 0.0406014298003869;
pts_[17] = 0.9639719272779138;
84 wts_[18] = 0.0176140071391521;
pts_[18] = -0.9931285991850949;
85 wts_[19] = 0.0176140071391521;
pts_[19] = 0.9931285991850949;
86 for (
size_t i = 0; i < 20; i++) {
133 if ( input <= ROL_EPSILON<Real>() ) {
136 if ( input >= 1.-ROL_EPSILON<Real>() ) {
139 Real a = ROL_EPSILON<Real>(), b = 1.-ROL_EPSILON<Real>(), c = 0.;
142 Real sa = ((fa < 0.) ? -1. : ((fa > 0.) ? 1. : 0.));
144 for (
size_t i = 0; i < 100; i++) {
147 sc = ((fc < 0.) ? -1. : ((fc > 0.) ? 1. : 0.));
148 if ( fc == 0. || (b-a)*0.5 < ROL_EPSILON<Real>() ) {
151 if ( sc == sa ) { a = c; fa = fc; sa = sc; }
179 void test(std::ostream &outStream = std::cout )
const {
181 std::vector<Real> X(size,0.);
182 std::vector<int> T(size,0);
183 X[0] = -4.*(Real)rand()/(Real)RAND_MAX;
187 X[2] = 0.5*(Real)rand()/(Real)RAND_MAX;
191 X[4] = 1.+4.*(Real)rand()/(Real)RAND_MAX;