//
//-----------------------------------------------------------------------
/// @copyright
///            (c) Copyright 2008 by GATS, Inc.,
///     11864 Canon Blvd, Suite 101, Newport News VA 23606
///
///  All Rights Reserved. No part of this software or publication may be
///  reproduced, stored in a retrieval system, or transmitted, in any form
///  or by any means, electronic, mechanical, photocopying, recording, or
///  otherwise without the prior written permission of GATS, Inc.
///
//-----------------------------------------------------------------------
///
/// @file	EdgeLocation.cpp
///
/// @author	John Burton	<j.c.burton@gats-inc.com>
///
/// @date  	Fri May 16 09:28:34 2008
///
//-----------------------------------------------------------------------
//
//-----------------------------------------------------------------------
// Include Files:
//-----------------------------------------------------------------------
//
#include "EventVar.h"


//
//-----------------------------------------------------------------------
// Defines and Macros:
//-----------------------------------------------------------------------
//



//
//-----------------------------------------------------------------------
// Global Variables:
//-----------------------------------------------------------------------
//



//
//-----------------------------------------------------------------------
// Class Methods:
//-----------------------------------------------------------------------
//



static double linterp(double &x, double &xa, double &xb, double &ya, double &yb)
{
  double dx, y;
  double a;

  dx = xb - xa;
  if(dx != 0)
  {
    a = (x - xa) / dx;
    y = ya + a * (yb - ya);
  }else
    y = (ya * yb) / 2.0;

  return y;
}

double findLowXEdge(double &edge, EventVar<double> &X, EventVar<double> &Y)
{
  int len = X.size();
  int i,j;
  double x,xa,xb,y,ya,yb;

  i = 0;
  while((i<len-1)&&(X[i]<edge)&&(X[i]<=X[i+1]))
    i++;
  
  ya = Y[i-1];
  yb = Y[i];
  xa = X[i-1];
  xb = X[i];
  x  = edge;
  y  = linterp(x,xa,xb,ya,yb);
//  std::cerr << "findLowXEdge:isamp, x,xa,xb,min,max,max2,y,ya,yb: ";
//  std::cerr << isamp << ", " << x << ", " << xa << ", " << xb << ", ";
//  std::cerr << min << ", " << max << ", " << sample.max() << ", ";
//  std::cerr << y << ", " << ya << ", " << yb << "\n";
  return y;
}


double findHighXEdge(double &edge, EventVar<double> &X, EventVar<double> &Y)
{
  int len = X.size();
  int i,j;
  double x,xa,xb,y,ya,yb;

  i = len-2;
  while((i>0)&&(X[i]<edge)&&(X[i]<=X[i-1]))
    i--;
  
  ya = Y[i];
  yb = Y[i+1];
  xa = X[i];
  xb = X[i+1];
  x  = edge;
  y  = linterp(x,xa,xb,ya,yb);
//  std::cerr << "findHighXEdge:isamp, x,xa,xb,min,max,maxndx,y,ya,yb: ";
//  std::cerr << isamp << ", " << x << ", " << xa << ", " << xb << ", ";
//  std::cerr << min << ", " << max << ", " << maxndx << ", ";
//  std::cerr << y << ", " << ya << ", " << yb << "\n";
  return y;
}