我今天才开始使用boost,发现this post非常有用。我正在尝试使用boost::bisect解决一系列值的参数方程。如果我想求解值为0.8,则以下工作:
#include <boost/math/tools/roots.hpp>
//declare constants
const double Y_r = 0.2126;
const double Y_g = 0.7152;
const double Y_b = 0.0722;
struct FunctionToApproximate {
double operator() (double x) {
return (pow (x, 2.4) * Y_r) + (pow ((x*0.6)+0.4, 2.4) * Y_g) + (1 * Y_b) - 0.8;
}
};
struct TerminationCondition {
bool operator() (double min, double max) {
return fabs(min - max) <= t_c;
}
};
using boost::math::tools::bisect;
std::pair<double, double> result = bisect(FunctionToApproximate(), 0.0, 1.0, TerminationCondition());
double root = (result.first + result.second) / 2;
我想将其包装在一个循环中,以便可以求解除0.8以外的值。我将如何处理?
非常感谢!
最佳答案
您可以将状态/数据成员添加到FunctionToApproximate
中,以保存要在每次调用中减去的值:
struct FunctionToApproximate {
FunctionToApproximate(double val) :val(val){}
double operator() (double x) {
return (pow (x, 2.4) * Y_r) + (pow ((x*0.6)+0.4, 2.4) * Y_g) + (1 * Y_b) - val;
}
double val;
};
然后将计算结果包装在一个循环中应该很简单。
for(double i = 0.1; i <= 1.0; i+=1.0) {
std::pair<double, double> result = bisect(FunctionToApproximate(i), 0.0, 1.0, TerminationCondition());
double root = (result.first + result.second) / 2;
}
关于c++ - 在一系列方程上使用boost::bisect,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/50195549/