/*****************************************************************/
/*    NAME: Michael Benjamin                                     */
/*    ORGN: Dept of Mechanical Eng / CSAIL, MIT Cambridge MA     */
/*    FILE: AOF_Quadratic.cpp                                    */
/*    DATE: Feb 2rd, 2008                                        */
/*                                                               */
/* This file is part of IvP Helm Core Libs                       */
/*                                                               */
/* IvP Helm Core Libs is free software: you can redistribute it  */
/* and/or modify it under the terms of the Lesser GNU General    */
/* Public License as published by the Free Software Foundation,  */
/* either version 3 of the License, or (at your option) any      */
/* later version.                                                */
/*                                                               */
/* IvP Helm Core Libs is distributed in the hope that it will    */
/* be useful but WITHOUT ANY WARRANTY; without even the implied  */
/* warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR       */
/* PURPOSE. See the Lesser GNU General Public License for more   */
/* details.                                                      */
/*                                                               */
/* You should have received a copy of the Lesser GNU General     */
/* Public License along with MOOS-IvP.  If not, see              */
/* <http://www.gnu.org/licenses/>.                               */
/*****************************************************************/

#include "AOF_Quadratic.h"

using namespace std;

//----------------------------------------------------------
// Procedure: Constructor

AOF_Quadratic::AOF_Quadratic(IvPDomain domain) : AOF(domain)
{
  m_coeff  = 0;
  n_coeff  = 0;
  x_center = 0;
  y_center = 0;
  x_index  = domain.getIndex("x");
  y_index  = domain.getIndex("y");
}

//----------------------------------------------------------------
// Procedure: setParam

bool AOF_Quadratic::setParam(const std::string& param, double val) 
{
  if(param == "mcoeff")
    m_coeff = val;
  else if(param == "ncoeff")
    n_coeff = val;
  else if(param == "x_center")
    x_center = val;
  else if(param == "y_center")
    y_center = val;
  else
    return(false);
  return(true);
}

//----------------------------------------------------------------
// Procedure: initialize

bool AOF_Quadratic::initialize()
{
  if((x_index != -1) && (y_index != -1))
    return(true);
  return(false);
}

//----------------------------------------------------------------
// Procedure: evalBox

double AOF_Quadratic::evalBox(const IvPBox *ptbox) const
{
  double x_val = m_domain.getVal(x_index, ptbox->pt(x_index));
  double y_val = m_domain.getVal(y_index, ptbox->pt(y_index));
  
  return((m_coeff * (x_val - x_center) * (x_val - x_center)) + 
	 (n_coeff * (y_val - y_center) * (y_val - y_center)));
}