2016-01-11 17:19:48 +00:00
|
|
|
|
/*
|
|
|
|
|
* rational numbers
|
|
|
|
|
* Copyright (c) 2003 Michael Niedermayer <michaelni@gmx.at>
|
|
|
|
|
*
|
|
|
|
|
* 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
|
|
|
|
|
*/
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* @file
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* @ingroup lavu_math_rational
|
|
|
|
|
* Utilties for rational number calculation.
|
2016-01-11 17:19:48 +00:00
|
|
|
|
* @author Michael Niedermayer <michaelni@gmx.at>
|
|
|
|
|
*/
|
|
|
|
|
|
|
|
|
|
#ifndef AVUTIL_RATIONAL_H
|
|
|
|
|
#define AVUTIL_RATIONAL_H
|
|
|
|
|
|
|
|
|
|
#include <stdint.h>
|
|
|
|
|
#include <limits.h>
|
|
|
|
|
#include "attributes.h"
|
|
|
|
|
|
|
|
|
|
/**
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* @defgroup lavu_math_rational AVRational
|
|
|
|
|
* @ingroup lavu_math
|
|
|
|
|
* Rational number calculation.
|
|
|
|
|
*
|
|
|
|
|
* While rational numbers can be expressed as floating-point numbers, the
|
|
|
|
|
* conversion process is a lossy one, so are floating-point operations. On the
|
|
|
|
|
* other hand, the nature of FFmpeg demands highly accurate calculation of
|
|
|
|
|
* timestamps. This set of rational number utilities serves as a generic
|
|
|
|
|
* interface for manipulating rational numbers as pairs of numerators and
|
|
|
|
|
* denominators.
|
|
|
|
|
*
|
|
|
|
|
* Many of the functions that operate on AVRational's have the suffix `_q`, in
|
|
|
|
|
* reference to the mathematical symbol "ℚ" (Q) which denotes the set of all
|
|
|
|
|
* rational numbers.
|
|
|
|
|
*
|
2016-01-11 17:19:48 +00:00
|
|
|
|
* @{
|
|
|
|
|
*/
|
|
|
|
|
|
|
|
|
|
/**
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* Rational number (pair of numerator and denominator).
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
typedef struct AVRational{
|
2018-07-30 02:07:02 +00:00
|
|
|
|
int num; ///< Numerator
|
|
|
|
|
int den; ///< Denominator
|
2016-01-11 17:19:48 +00:00
|
|
|
|
} AVRational;
|
|
|
|
|
|
|
|
|
|
/**
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* Create an AVRational.
|
|
|
|
|
*
|
2016-01-11 17:19:48 +00:00
|
|
|
|
* Useful for compilers that do not support compound literals.
|
2018-07-30 02:07:02 +00:00
|
|
|
|
*
|
|
|
|
|
* @note The return value is not reduced.
|
|
|
|
|
* @see av_reduce()
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
static inline AVRational av_make_q(int num, int den)
|
|
|
|
|
{
|
|
|
|
|
AVRational r = { num, den };
|
|
|
|
|
return r;
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* Compare two rationals.
|
2018-07-30 02:07:02 +00:00
|
|
|
|
*
|
|
|
|
|
* @param a First rational
|
|
|
|
|
* @param b Second rational
|
|
|
|
|
*
|
|
|
|
|
* @return One of the following values:
|
|
|
|
|
* - 0 if `a == b`
|
|
|
|
|
* - 1 if `a > b`
|
|
|
|
|
* - -1 if `a < b`
|
|
|
|
|
* - `INT_MIN` if one of the values is of the form `0 / 0`
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
static inline int av_cmp_q(AVRational a, AVRational b){
|
|
|
|
|
const int64_t tmp= a.num * (int64_t)b.den - b.num * (int64_t)a.den;
|
|
|
|
|
|
|
|
|
|
if(tmp) return (int)((tmp ^ a.den ^ b.den)>>63)|1;
|
|
|
|
|
else if(b.den && a.den) return 0;
|
|
|
|
|
else if(a.num && b.num) return (a.num>>31) - (b.num>>31);
|
|
|
|
|
else return INT_MIN;
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
/**
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* Convert an AVRational to a `double`.
|
|
|
|
|
* @param a AVRational to convert
|
|
|
|
|
* @return `a` in floating-point form
|
|
|
|
|
* @see av_d2q()
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
static inline double av_q2d(AVRational a){
|
|
|
|
|
return a.num / (double) a.den;
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* Reduce a fraction.
|
2018-07-30 02:07:02 +00:00
|
|
|
|
*
|
2016-01-11 17:19:48 +00:00
|
|
|
|
* This is useful for framerate calculations.
|
2018-07-30 02:07:02 +00:00
|
|
|
|
*
|
|
|
|
|
* @param[out] dst_num Destination numerator
|
|
|
|
|
* @param[out] dst_den Destination denominator
|
|
|
|
|
* @param[in] num Source numerator
|
|
|
|
|
* @param[in] den Source denominator
|
|
|
|
|
* @param[in] max Maximum allowed values for `dst_num` & `dst_den`
|
|
|
|
|
* @return 1 if the operation is exact, 0 otherwise
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
int av_reduce(int *dst_num, int *dst_den, int64_t num, int64_t den, int64_t max);
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* Multiply two rationals.
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* @param b First rational
|
|
|
|
|
* @param c Second rational
|
2016-01-11 17:19:48 +00:00
|
|
|
|
* @return b*c
|
|
|
|
|
*/
|
|
|
|
|
AVRational av_mul_q(AVRational b, AVRational c) av_const;
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* Divide one rational by another.
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* @param b First rational
|
|
|
|
|
* @param c Second rational
|
2016-01-11 17:19:48 +00:00
|
|
|
|
* @return b/c
|
|
|
|
|
*/
|
|
|
|
|
AVRational av_div_q(AVRational b, AVRational c) av_const;
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* Add two rationals.
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* @param b First rational
|
|
|
|
|
* @param c Second rational
|
2016-01-11 17:19:48 +00:00
|
|
|
|
* @return b+c
|
|
|
|
|
*/
|
|
|
|
|
AVRational av_add_q(AVRational b, AVRational c) av_const;
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* Subtract one rational from another.
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* @param b First rational
|
|
|
|
|
* @param c Second rational
|
2016-01-11 17:19:48 +00:00
|
|
|
|
* @return b-c
|
|
|
|
|
*/
|
|
|
|
|
AVRational av_sub_q(AVRational b, AVRational c) av_const;
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* Invert a rational.
|
|
|
|
|
* @param q value
|
|
|
|
|
* @return 1 / q
|
|
|
|
|
*/
|
|
|
|
|
static av_always_inline AVRational av_inv_q(AVRational q)
|
|
|
|
|
{
|
|
|
|
|
AVRational r = { q.den, q.num };
|
|
|
|
|
return r;
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* Convert a double precision floating point number to a rational.
|
|
|
|
|
*
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* In case of infinity, the returned value is expressed as `{1, 0}` or
|
|
|
|
|
* `{-1, 0}` depending on the sign.
|
|
|
|
|
*
|
|
|
|
|
* @param d `double` to convert
|
|
|
|
|
* @param max Maximum allowed numerator and denominator
|
|
|
|
|
* @return `d` in AVRational form
|
|
|
|
|
* @see av_q2d()
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
AVRational av_d2q(double d, int max) av_const;
|
|
|
|
|
|
|
|
|
|
/**
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* Find which of the two rationals is closer to another rational.
|
|
|
|
|
*
|
|
|
|
|
* @param q Rational to be compared against
|
|
|
|
|
* @param q1,q2 Rationals to be tested
|
|
|
|
|
* @return One of the following values:
|
|
|
|
|
* - 1 if `q1` is nearer to `q` than `q2`
|
|
|
|
|
* - -1 if `q2` is nearer to `q` than `q1`
|
|
|
|
|
* - 0 if they have the same distance
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
int av_nearer_q(AVRational q, AVRational q1, AVRational q2);
|
|
|
|
|
|
|
|
|
|
/**
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* Find the value in a list of rationals nearest a given reference rational.
|
|
|
|
|
*
|
|
|
|
|
* @param q Reference rational
|
|
|
|
|
* @param q_list Array of rationals terminated by `{0, 0}`
|
|
|
|
|
* @return Index of the nearest value found in the array
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
int av_find_nearest_q_idx(AVRational q, const AVRational* q_list);
|
|
|
|
|
|
|
|
|
|
/**
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* Convert an AVRational to a IEEE 32-bit `float` expressed in fixed-point
|
|
|
|
|
* format.
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*
|
2018-07-30 02:07:02 +00:00
|
|
|
|
* @param q Rational to be converted
|
|
|
|
|
* @return Equivalent floating-point value, expressed as an unsigned 32-bit
|
|
|
|
|
* integer.
|
|
|
|
|
* @note The returned value is platform-indepedant.
|
2016-01-11 17:19:48 +00:00
|
|
|
|
*/
|
|
|
|
|
uint32_t av_q2intfloat(AVRational q);
|
|
|
|
|
|
|
|
|
|
/**
|
|
|
|
|
* @}
|
|
|
|
|
*/
|
|
|
|
|
|
|
|
|
|
#endif /* AVUTIL_RATIONAL_H */
|