我正在尝试实现一些C++代码,以使用单纯形算法找到函数的最大值。不幸的是,我对C++的经验为零。
我遇到了这个错误,似乎无法从类似问题的答案中找到解决方案。
也有与链接文件“simplex615.h”有关的警告。
simplex.cpp
#include <vector>
#include <cmath>
#include <iostream>
#include "simplex615.h"
#define ZEPS 1e-10
// function object used as an argument
class arbitraryFunc {
public:
double operator() (std::vector<double>& x) {
// f(x0,x1) = 100*(x1-x0^2)^2 + (1-x0)^2
return 100*(x[1]-x[0]*x[0])*(x[1]-x[0]*x[0])+(1-x[0])*(1-x[0]);
}
};
int main(int main, char** argv) {
double point[2] = {-1.2, 1};
arbitraryFunc foo;
// initial point to start
// WILL BE DISCUSSED LATER
simplex615 <arbitraryFunc> simplex(point, 2); // create a simplex
simplex.amoeba(foo, 1e-7); // optimize for a function
// print outputs
std::cout << "Minimum = " << simplex.ymin() << ", at ("
<< simplex.xmin()[0] << ", " << simplex.xmin()[1]
<< ")" << std::endl;
return 0;
}
simplex615.h
template <class F> // F is a function object
void simplex615 <F>::amoeba(optFunc& foo, double tol) {
evaluateFunction(foo);
while(true) {
evaluateExtremes();
prepareUpdate();
if ( check_tol(Y[idxHi],Y[idxLo],tol) ) break;
updateSimplex(foo, -1.0); // reflection
if ( Y[idxHi] < Y[idxLo] ) {
updateSimplex(foo, -2.0); // expansion
}
else if ( Y[idxHi] >= Y[idxNextHi] ) {
if ( !updateSimplex(foo, 0.5) ) {
contractSimplex(foo);
}
}
}
}
simplex615.h
class optFunc {
public:
virtual double operator() (std::vector<double>& x) = 0;
};
链接到完整文件simplex.cpp和simplex.h:Source code
任何帮助将不胜感激。谢谢。
最佳答案
在我看来,您在simplex615.h
中忘记了在amoeba
方法中使用“F类”。只需将optFunc
替换为F
,它就可以解决此问题。
template <class F> // F is a function object
void simplex615 <F>::amoeba(F& foo, double tol) {
...
}
C++中的模板类参数定义了一种通用类型,可以在使用模板时将其替换。
同样从此示例中,您可以从头文件中删除optFunc的声明。
关于c++ - C++错误:没有匹配的函数可以调用‘simplex615 <arbitraryFunc>::amoeba,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/54362898/