C++ Определение перегруженного оператора виртуального базового класса

Я пытаюсь написать набор общих математических служебных классов (искатели корней, интеграторы и т. д.), которые при построении принимают указатель на базовый класс, определяющий функцию, с которой я хочу, чтобы работал конкретный алгоритм. Базовый класс должен определять только общедоступный виртуальный интерфейс (абстрактный или с тривиальной функциональностью по умолчанию) type operator()(type inputArg), который может быть реализован пользователем по мере необходимости. Это позволит пользователю просто реализовать необходимые функторы и выполнить необходимые вычисления. Мой mwe ниже:

Этот первый заголовок определяет абстрактный класс интерфейса.

// BaseFunctor.h

#ifndef _BASE_FUNCTOR_H_
#define _BASE_FUNCTOR_H_

class BaseFunctor
{
public:
   virtual double operator() (double x) = 0;
};
#endif

Это класс для одного из методов решателя

// NewtonsMethod.h

#ifndef _NEWTONS_METHOD_H_
#define _NEWTONS_METHOD_H_

class BaseFunctor;

class NewtonsMethod
{
public: 

   NewtonsMethod(BaseFunctor *rhsIn,
                 BaseFunctor *rhsPrimeIn,
                 double       x0In);

   ~NewtonsMethod();

   bool DetermineRoot(double &root);

 private:

   double       x0;
   BaseFunctor *rhs;
   BaseFunctor *rhsPrime;

   static const double       EPSILON;
   static const unsigned int MAX_ITER;

};
#endif

Это реализация решателя: // NewtonsMethod.cpp

#include "NewtonsMethod.h"
#include "BaseFunctor.h"
#include <cmath>

const double       NewtonsMethod::EPSILON  = 1e-9;
const unsigned int NewtonsMethod::MAX_ITER = 30;

NewtonsMethod::NewtonsMethod(BaseFunctor *rhsIn,
                             BaseFunctor *rhsPrimeIn,
                             double       x0In) :
   rhs     (rhsIn),
   rhsPrime(rhsPrimeIn),
   x0      (x0In)
{ }

NewtonsMethod::~NewtonsMethod() { }

bool NewtonsMethod::DetermineRoot(double &root)
{
   // This is obviously not implemented
   root = rhs(1.0) / rhsPrime(2.0);
   return false;
}

И основная функция, в которой я делаю реализации производного класса: // Main.cpp

#include "BaseFunctor.h"
#include "NewtonsMethod.h"
#include <iostream>
#include <iomanip>

class fOfX : public BaseFunctor
{
   double operator() (double x)
   {
      return x * x - 2.0;
   }
};

class fPrimeOfX : public BaseFunctor
{
   double operator() (double x)
   {
      return 2.0 * x;
   }
};


int main()
{
   double x0 = 2.0;

   BaseFunctor *f      = new fOfX();
   BaseFunctor *fPrime = new fPrimeOfX(); 

   NewtonsMethod newton(f, fPrime, x0);

   double root      = 0.0;
   bool   converged = newton.DetermineRoot(root);

   if (converged)
   {
      std::cout << "SUCCESS: root == " << std::setprecision(16) << root << std::endl;
   }
   else
   {
      std::cout << "FAILED: root == " << std::setprecision(16) << root << std::endl;
   }
   delete f;
   delete fPrime;
}

Я постарался сделать это как можно короче, так что извините, если получилось слишком длинно. В основном я получаю ошибку:

g++ Main.cpp NewtonsMethod.cpp -o main
NewtonsMethod.cpp: In member function ‘bool NewtonsMethod::DetermineRoot(double&)’: 
NewtonsMethod.cpp:29: error: ‘((NewtonsMethod*)this)->NewtonsMethod::rhs’ cannot be used    as a function
NewtonsMethod.cpp:29: error: ‘((NewtonsMethod*)this)->NewtonsMethod::rhsPrime’ cannot be   used as a function

Как я могу решить эту проблему, сохранив желаемую функциональность или создав класс для различных необходимых функций?

Спасибо