/*
  * This file is part of FFmpeg.
  *
  * FFmpeg is free software; you can redistribute it and/or
  * modify it under the terms of the GNU Lesser General Public
  * License as published by the Free Software Foundation; either
  * version 2.1 of the License, or (at your option) any later version.
  *
  * FFmpeg 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 GNU
  * Lesser General Public License for more details.
  *
  * You should have received a copy of the GNU Lesser General Public
  * License along with FFmpeg; if not, write to the Free Software
  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  */
 

 #ifndef AVFILTER_HERMITE_H
 #define AVFILTER_HERMITE_H
 

 static inline double hermite_interpolation(double x, double x0, double x1,

                                     double p0, double p1,
                                     double m0, double m1)
 {
     double width = x1 - x0;
     double t = (x - x0) / width;
     double t2, t3;
     double ct0, ct1, ct2, ct3;
 
     m0 *= width;
     m1 *= width;
 
     t2 = t*t;
     t3 = t2*t;
     ct0 = p0;
     ct1 = m0;
 
     ct2 = -3 * p0 - 2 * m0 + 3 * p1 - m1;
     ct3 = 2 * p0 + m0  - 2 * p1 + m1;
 
     return ct3 * t3 + ct2 * t2 + ct1 * t + ct0;
 }

 
 #endif /* AVFILTER_HERMITE_H */