46 double (*
const *
funcs1)(
void *,
double a);
48 double (*
const *
funcs2)(
void *,
double a,
double b);
96 if(numstr[0]==
'0' && (numstr[1]|0x20)==
'x') {
97 d = strtoul(numstr, &next, 16);
102 if (next[0] ==
'd' && next[1] ==
'B') {
106 }
else if (*next >=
'E' && *next <=
'z') {
109 if (next[1] ==
'i') {
131 #define IS_IDENTIFIER_CHAR(c) ((c) - '0' <= 9U || (c) - 'a' <= 25U || (c) - 'A' <= 25U || (c) == '_')
136 for (i=0; prefix[i]; i++) {
137 if (prefix[i] != s[i])
return 0;
158 double (*
func2)(
void *, double, double);
206 av_log(p, level,
"%f\n", x);
212 r= r*1664525+1013904223;
214 return e->
value * (r * (1.0/UINT64_MAX));
223 double t = 1, d = 0,
v;
227 double var0 = p->
var[
id];
228 for(i=0; i<1000; i++) {
242 double low = -1, high = -1,
v, low_v = -DBL_MAX, high_v = DBL_MAX;
243 double var0 = p->
var[0];
245 for(i=-1; i<1024; i++) {
249 p->
var[0] = x_max*pow(0.9, i-255);
250 if (i&1) p->
var[0] *= -1;
251 if (i&2) p->
var[0] += low;
252 else p->
var[0] += high;
255 if (v<=0 && v>low_v) {
259 if (
v>=0 &&
v<high_v) {
263 if (low>=0 && high>=0){
264 for (j=0; j<1000; j++) {
265 p->
var[0] = (low+high)*0.5;
266 if (low == p->
var[0] || high == p->
var[0])
269 if (
v<=0) low = p->
var[0];
270 if (
v>=0) high= p->
var[0];
280 return -low_v<high_v ? low : high;
286 case e_mod:
return e->
value * (d - floor((!CONFIG_FTRAPV || d2) ? d / d2 : d *
INFINITY) * d2);
288 case e_max:
return e->
value * (d > d2 ? d : d2);
289 case e_min:
return e->
value * (d < d2 ? d : d2);
290 case e_eq:
return e->
value * (d == d2 ? 1.0 : 0.0);
291 case e_gt:
return e->
value * (d > d2 ? 1.0 : 0.0);
292 case e_gte:
return e->
value * (d >= d2 ? 1.0 : 0.0);
293 case e_lt:
return e->
value * (d < d2 ? 1.0 : 0.0);
294 case e_lte:
return e->
value * (d <= d2 ? 1.0 : 0.0);
325 char *next = p->
s, *
s0 = p->
s;
361 p->
s= strchr(p->
s,
'(');
373 if (p->
s[0] !=
')') {
394 if (p->
s[0] !=
')') {
490 *sign= (*p->
s ==
'+') - (*p->
s ==
'-');
502 if (next != p->
s && next[0] ==
'd' && next[1] ==
'B') {
512 int sign, sign2,
ret;
514 if ((ret =
parse_dB(&e0, p, &sign)) < 0)
519 if ((ret =
parse_dB(&e2, p, &sign2)) < 0) {
531 if (e0) e0->
value *= (sign|1);
543 while (p->
s[0]==
'*' || p->
s[0]==
'/') {
567 while (*p->
s ==
'+' || *p->
s ==
'-') {
595 while (*p->
s ==
';') {
653 const char *
const *const_names,
655 const char *
const *func2_names,
double (*
const *funcs2)(
void *,
double,
double),
656 int log_offset,
void *log_ctx)
718 const char *
const *const_names,
const double *const_values,
719 const char *
const *func1_names,
double (*
const *
funcs1)(
void *,
double),
720 const char *
const *func2_names,
double (*
const *funcs2)(
void *,
double,
double),
721 void *opaque,
int log_offset,
void *log_ctx)
724 int ret =
av_expr_parse(&e, s, const_names, func1_names,
funcs1, func2_names, funcs2, log_offset, log_ctx);
738 static const double const_values[] = {
744 static const char *
const const_names[] = {
750 int main(
int argc,
char **argv)
754 const char *
const *expr;
755 static const char *
const exprs[] = {
761 "1+(5-2)^(3-1)+1/2+sin(PI)-max(-2.2,-3.1)",
782 "bar + PI + E + 100f*2 + foo",
783 "13k + 12f - foo(1, 2)",
797 "st(0, 1); while(lte(ld(0), 100), st(1, ld(1)+ld(0));st(0, ld(0)+1)); ld(1)",
799 "st(1, 1); st(2, 2); st(0, 1); while(lte(ld(0),10), st(3, ld(1)+ld(2)); st(1, ld(2)); st(2, ld(3)); st(0, ld(0)+1)); ld(3)",
801 "st(0, 1); while(lte(ld(0),100), st(1, ld(1)+ld(0)); st(0, ld(0)+1))",
830 "ifnot(1, NaN) + if(0, 1)",
834 "taylor(eq(mod(ld(1),4),1)-eq(mod(ld(1),4),3), PI/2, 1)",
835 "root(sin(ld(0))-1, 2)",
836 "root(sin(ld(0))+6+sin(ld(0)/12)-log(ld(0)), 100)",
837 "7000000B*random(0)",
841 "gcd(30,55)*print(min(9,1))",
845 "between(10, -3, 10)",
846 "between(-4, -2, -1)",
854 for (expr = exprs; *expr; expr++) {
855 printf(
"Evaluating '%s'\n", *expr);
857 const_names, const_values,
860 printf(
"'%s' -> nan\n\n", *expr);
862 printf(
"'%s' -> %f\n\n", *expr, d);
866 const_names, const_values,
868 printf(
"%f == 12.7\n", d);
870 const_names, const_values,
872 printf(
"%f == 0.931322575\n", d);
874 if (argc > 1 && !strcmp(argv[1],
"-t")) {
875 for (i = 0; i < 1050; i++) {
878 const_names, const_values,
static int verify_expr(AVExpr *e)
const char *const * func1_names
static const char *const func1_names[]
#define LIBAVUTIL_VERSION_INT
#define FF_ARRAY_ELEMS(a)
int av_expr_parse(AVExpr **expr, const char *s, const char *const *const_names, const char *const *func1_names, double(*const *funcs1)(void *, double), const char *const *func2_names, double(*const *funcs2)(void *, double, double), int log_offset, void *log_ctx)
Parse an expression.
Macro definitions for various function/variable attributes.
static int parse_dB(AVExpr **e, Parser *p, int *sign)
int av_isspace(int c)
Locale-independent conversion of ASCII isspace.
static av_always_inline av_const int isnan(float x)
static av_cold int end(AVCodecContext *avctx)
static int parse_expr(AVExpr **e, Parser *p)
static double(*const funcs1[])(void *, double)
double strtod(const char *, char **)
const char *const * const_names
double(* func1)(void *, double)
static int parse_pow(AVExpr **e, Parser *p, int *sign)
static int parse_factor(AVExpr **e, Parser *p)
high precision timer, useful to profile code
static int parse_subexpr(AVExpr **e, Parser *p)
double(* func2)(void *, double, double)
int av_expr_parse_and_eval(double *d, const char *s, const char *const *const_names, const double *const_values, const char *const *func1_names, double(*const *funcs1)(void *, double), const char *const *func2_names, double(*const *funcs2)(void *, double, double), void *opaque, int log_offset, void *log_ctx)
Parse and evaluate an expression.
const double * const_values
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
#define IS_IDENTIFIER_CHAR(c)
static int parse_term(AVExpr **e, Parser *p)
int64_t av_gcd(int64_t a, int64_t b)
Return the greatest common divisor of a and b.
static double etime(double v)
static const int8_t si_prefixes['z'- 'E'+1]
GLsizei GLboolean const GLfloat * value
static av_always_inline av_const double trunc(double x)
static int strmatch(const char *s, const char *prefix)
static int parse_primary(AVExpr **e, Parser *p)
static const struct @177 constants[]
int64_t av_gettime(void)
Get the current time in microseconds.
double av_strtod(const char *numstr, char **tail)
Parse the string in numstr and return its value as a double.
#define AV_LOG_INFO
Standard information.
static double eval_expr(Parser *p, AVExpr *e)
void av_expr_free(AVExpr *e)
Free a parsed expression previously created with av_expr_parse().
Describe the class of an AVClass context structure.
static const AVClass eval_class
double(*const funcs1)(void *, double a)
const char *const * func2_names
attribute_deprecated const uint8_t av_reverse[256]
Reverse the order of the bits of an 8-bits unsigned integer.
common internal and external API header
double(*const funcs2)(void *, double a, double b)
double av_expr_eval(AVExpr *e, const double *const_values, void *opaque)
Evaluate a previously parsed expression.
#define FF_QP2LAMBDA
factor to convert from H.263 QP to lambda
static AVExpr * make_eval_expr(int type, int value, AVExpr *p0, AVExpr *p1)
static av_always_inline av_const int isinf(float x)
int main(int argc, char **argv)
void * av_mallocz(size_t size)
Allocate a block of size bytes with alignment suitable for all memory accesses (including vectors if ...
simple arithmetic expression evaluator