Informatyka UJ forum

Robson

/******************* Klasa Newton **********************/
/******** (c) Marian Mrozek, Kraków 2006 **************/

template<typename Function, typename Scalar>
class Newton{
private:
const Function& f;
const Function& derf;
const Scalar epsilon;
const int limit;
Scalar findRoot(const Scalar x,const int depth);
public:
Newton(const Function& A_f, const Function& A_derf, Scalar A_epsilon, int A_limit):
f(A_f),derf(A_derf),epsilon(A_epsilon),limit(A_limit){}
Scalar findRoot(const Scalar x){
return findRoot(x,0);
}
};

template<typename Function, typename Scalar>
Scalar Newton<Function,Scalar>::findRoot(const Scalar x,const int depth){
const Scalar y=f(x);
cout << "Step " << depth << "\n";
cout << " x =" << x << "\n";
cout << " f(x)=" << y << "\n";
if(y<epsilon && y>-epsilon || depth>=limit){
return x;
}else{
return findRoot(x-y/derf(x),depth+1);
}
}

/**************** Metoda Newtona (sin) *****************/
/******** (c) Marian Mrozek, Kraków 2006 **************/

#include <iostream>
using namespace std;

#include <cmath>

typedef double doubleFdouble(double);

int main(){
const Newton<doubleFdouble,double> n(sin,cos,0.0001,5);
cout << n.findRoot(3.0) << std::endl;
}

kompilator (gcc z dev 4.9.costam ) rzuca sie o:
main.cpp: In function `int main()':
main.cpp:54: passing `const Newton<double ()(double), double>' as `this'

argument of `Scalar Newton<Function, Scalar>::findRoot(Scalar) [with
Function = double ()(double), Scalar = double]' discards qualifiers

mogłby ktoś przystępnie wyjaśnic skąd to sie bierze... bo nie moge zrozumiec tej informacji...

rafal

Funkcje findRoot() zarówno ta jedno, jak i dwuargumentowa muszą być const, skoro podajesz jako Function stałą. Albo dopisz const po funkcjach findRoot(..) albo usuń const, które jest przed Newton<> w main..

rafal

Sorry, nie jako Function (to tutaj nie ma nic do rzeczy ;)) tylko po prostu n jest stałą, więc może wywoływać tylko funkcje, które jej nie zmodyfikują czyli findRoot(..) musi być const.