FunctionQuotient.cpp

//==========================================================================//
// Imported from the QAT library (qat.pitt.edu) by the copyright holder     //
// Joe Boudreau.  Software licensed under the terms and conditions of the   //
// GNU Lesser Public License v3.                                            //
//==========================================================================//

#include "GeoGenericFunctions/FunctionQuotient.h"
#include <stdexcept>

namespace GeoGenfun {
FUNCTION_OBJECT_IMP(FunctionQuotient)

FunctionQuotient::FunctionQuotient(const AbsFunction *arg1, const AbsFunction *arg2):
  _arg1(arg1->clone()),
  _arg2(arg2->clone())
{
  if (arg1->dimensionality()!=arg2->dimensionality()) {
    throw std::runtime_error ("FunctionQuotient:  dimension mismatch");
  }
}

FunctionQuotient::FunctionQuotient( const FunctionQuotient & right) :
AbsFunction(right),
_arg1(right._arg1->clone()),
_arg2(right._arg2->clone())
{}
FunctionQuotient::~FunctionQuotient()
{
  delete _arg1;
  delete _arg2;
}

unsigned int FunctionQuotient::dimensionality() const {
  return _arg1->dimensionality();
}


double FunctionQuotient::operator ()(double x) const
{
  return (*_arg1)(x)/(*_arg2)(x);
}

double FunctionQuotient::operator ()(const Argument & x) const
{
  return (*_arg1)(x)/(*_arg2)(x);
}


Derivative FunctionQuotient::partial(unsigned int index) const {
  const Derivative & d1=_arg1->partial(index);
  const Derivative & d2=_arg2->partial(index);
  
  const AbsFunction & fPrime  = 
    ((d1)*(*_arg2)-(*_arg1)*d2) / 
    (*_arg2)/ (*_arg2);
  return Derivative(&fPrime);
}

} // namespace GeoGenfun