#ifndef VALARRAY_LINEAR_INTERPOLATION_04_07_2006
#define VALARRAY_LINEAR_INTERPOLATION_04_07_2006

/**
      @file ValInterp.h
      @author Brian Magill
      @datecreated 4/07/2006

	  $Date:$
      $Revision:$

      @copyright
      (©) Copyright 2006 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.
 
      @brief Class for linear interpolation     (double precision version)
*/
//----------------------------------------------------------------------
//

#include <valarray>

    int Vlocate(std::valarray<double> xin, double x);

    double linear_interpolate(double x, double xa, double xb, 
                                                  double ya, double yb);

    double Interpol(std::valarray<double> const &yin, 
                                std::valarray<double> const &xin, double xout);

    std::valarray<double> VInterpol(std::valarray<double> const &yin, 
                 std::valarray<double> const &xin, std::valarray<double> const &xout);


#endif