76 Real f1 = std::pow((*xp)[0],2.0) + std::pow((*xp)[1],4.0);
77 Real f2 = std::pow(2.0-(*xp)[0],2.0) + std::pow(2.0-(*xp)[1],2.0);
78 Real f3 = 2.0*std::exp(-(*xp)[0] + (*xp)[1]);
79 return std::max(f1,std::max(f2,f3));
85 Real f1 = std::pow((*xp)[0],2.0) + std::pow((*xp)[1],4.0);
86 Real f2 = std::pow(2.0-(*xp)[0],2.0) + std::pow(2.0-(*xp)[1],2.0);
87 Real f3 = 2.0*std::exp(-(*xp)[0] + (*xp)[1]);
90 if( f1 >= std::max(f2,f3) ) {
91 (*gp)[0] = 2.0*(*xp)[0];
92 (*gp)[1] = 4.0*std::pow((*xp)[1],3.0);
94 else if ( f2 >= std::max(f1,f3) ) {
95 (*gp)[0] = 2.0*(*xp)[0]-4.0;
96 (*gp)[1] = 2.0*(*xp)[1]-4.0;
98 else if ( f3 >= std::max(f1,f2) ) {
99 (*gp)[0] = -2.0*std::exp(-(*xp)[0]+(*xp)[1]);
100 (*gp)[1] = 2.0*std::exp(-(*xp)[0]+(*xp)[1]);