6821269217

Committed 10 Nov 2023 05:48AM UTC coverage: 88.456% (-0.004%) from 88.46%

Build # 6821269217

Build Type

push

github

Committed by

igormunkin

Commit Message

Fix FOLD rule for x-0.

Reported by XmiliaH.

(cherry-picked from commit 7b994e0ee)

Fold optimization x - (-0) ==> x is INVALID for x = -0 in FP arithmetic.
Its result is -0 instead of +0. This patch allows only x - (+0) ==> x
optimization.

Sergey Kaplun:
* added the description and the test for the problem

Part of tarantool/tarantool#9145

Run Details

5359 of 5975 branches covered (0.0%)

Branch coverage included in aggregate %.

1 of 1 new or added line in 1 file covered. (100.0%)

3 existing lines in 3 files now uncovered.

20540 of 23304 relevant lines covered (88.14%)

2753453.46 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

92.11

/src/lj_vmmath.c

/*
** Math helper functions for assembler VM.
** Copyright (C) 2005-2017 Mike Pall. See Copyright Notice in luajit.h
*/

#define lj_vmmath_c
#define LUA_CORE

#include <errno.h>
#include <math.h>

#include "lj_obj.h"
#include "lj_ir.h"
#include "lj_vm.h"

/* -- Wrapper functions --------------------------------------------------- */

#if LJ_TARGET_X86 && __ELF__ && __PIC__
/* Wrapper functions to deal with the ELF/x86 PIC disaster. */
LJ_FUNCA double lj_wrap_log(double x) { return log(x); }
LJ_FUNCA double lj_wrap_log10(double x) { return log10(x); }
LJ_FUNCA double lj_wrap_exp(double x) { return exp(x); }
LJ_FUNCA double lj_wrap_sin(double x) { return sin(x); }
LJ_FUNCA double lj_wrap_cos(double x) { return cos(x); }
LJ_FUNCA double lj_wrap_tan(double x) { return tan(x); }
LJ_FUNCA double lj_wrap_asin(double x) { return asin(x); }
LJ_FUNCA double lj_wrap_acos(double x) { return acos(x); }
LJ_FUNCA double lj_wrap_atan(double x) { return atan(x); }
LJ_FUNCA double lj_wrap_sinh(double x) { return sinh(x); }
LJ_FUNCA double lj_wrap_cosh(double x) { return cosh(x); }
LJ_FUNCA double lj_wrap_tanh(double x) { return tanh(x); }
LJ_FUNCA double lj_wrap_atan2(double x, double y) { return atan2(x, y); }
LJ_FUNCA double lj_wrap_pow(double x, double y) { return pow(x, y); }
LJ_FUNCA double lj_wrap_fmod(double x, double y) { return fmod(x, y); }
#endif

/* -- Helper functions ---------------------------------------------------- */

double lj_vm_foldarith(double x, double y, int op)
{
  switch (op) {
  case IR_ADD - IR_ADD: return x+y; break;
  case IR_SUB - IR_ADD: return x-y; break;
  case IR_MUL - IR_ADD: return x*y; break;
  case IR_DIV - IR_ADD: return x/y; break;
  case IR_MOD - IR_ADD: return x-lj_vm_floor(x/y)*y; break;
  case IR_POW - IR_ADD: return pow(x, y); break;
  case IR_NEG - IR_ADD: return -x; break;
  case IR_ABS - IR_ADD: return fabs(x); break;
#if LJ_HASJIT
  case IR_LDEXP - IR_ADD: return ldexp(x, (int)y); break;
  case IR_MIN - IR_ADD: return x < y ? x : y; break;
  case IR_MAX - IR_ADD: return x > y ? x : y; break;
#endif
  default: return x;
  }
}

/* -- Helper functions for generated machine code ------------------------- */

#if (LJ_HASJIT && !(LJ_TARGET_ARM || LJ_TARGET_ARM64 || LJ_TARGET_PPC)) || LJ_TARGET_MIPS
int32_t LJ_FASTCALL lj_vm_modi(int32_t a, int32_t b)
{
  uint32_t y, ua, ub;
  /* This must be checked before using this function. */
  lj_assertX(b != 0, "modulo with zero divisor");
  ua = a < 0 ? (uint32_t)-a : (uint32_t)a;
  ub = b < 0 ? (uint32_t)-b : (uint32_t)b;
  y = ua % ub;
  if (y != 0 && (a^b) < 0) y = y - ub;
  if (((int32_t)y^b) < 0) y = (uint32_t)-(int32_t)y;
  return (int32_t)y;
}
#endif

#if LJ_HASJIT

#ifdef LUAJIT_NO_LOG2
double lj_vm_log2(double a)
{
  return log(a) * 1.4426950408889634074;
}
#endif

/* Computes fpm(x) for extended math functions. */
double lj_vm_foldfpm(double x, int fpm)
{
  switch (fpm) {
  case IRFPM_FLOOR: return lj_vm_floor(x);
  case IRFPM_CEIL: return lj_vm_ceil(x);
  case IRFPM_TRUNC: return lj_vm_trunc(x);
  case IRFPM_SQRT: return sqrt(x);
  case IRFPM_LOG: return log(x);
  case IRFPM_LOG2: return lj_vm_log2(x);
  default: lj_assertX(0, "bad fpm %d", fpm);
  }
  return 0;
}

#if LJ_HASFFI
int lj_vm_errno(void)
{
  return errno;
}
#endif

#endif

1	/*
2	** Math helper functions for assembler VM.
3	** Copyright (C) 2005-2017 Mike Pall. See Copyright Notice in luajit.h
4	*/
5
6	#define lj_vmmath_c
7	#define LUA_CORE
8
9	#include <errno.h>
10	#include <math.h>
11
12	#include "lj_obj.h"
13	#include "lj_ir.h"
14	#include "lj_vm.h"
15
16	/* -- Wrapper functions --------------------------------------------------- */
17
18	#if LJ_TARGET_X86 && __ELF__ && __PIC__
19	/* Wrapper functions to deal with the ELF/x86 PIC disaster. */
20	LJ_FUNCA double lj_wrap_log(double x) { return log(x); }
21	LJ_FUNCA double lj_wrap_log10(double x) { return log10(x); }
22	LJ_FUNCA double lj_wrap_exp(double x) { return exp(x); }
23	LJ_FUNCA double lj_wrap_sin(double x) { return sin(x); }
24	LJ_FUNCA double lj_wrap_cos(double x) { return cos(x); }
25	LJ_FUNCA double lj_wrap_tan(double x) { return tan(x); }
26	LJ_FUNCA double lj_wrap_asin(double x) { return asin(x); }
27	LJ_FUNCA double lj_wrap_acos(double x) { return acos(x); }
28	LJ_FUNCA double lj_wrap_atan(double x) { return atan(x); }
29	LJ_FUNCA double lj_wrap_sinh(double x) { return sinh(x); }
30	LJ_FUNCA double lj_wrap_cosh(double x) { return cosh(x); }
31	LJ_FUNCA double lj_wrap_tanh(double x) { return tanh(x); }
32	LJ_FUNCA double lj_wrap_atan2(double x, double y) { return atan2(x, y); }
33	LJ_FUNCA double lj_wrap_pow(double x, double y) { return pow(x, y); }
34	LJ_FUNCA double lj_wrap_fmod(double x, double y) { return fmod(x, y); }
35	#endif
36
37	/* -- Helper functions ---------------------------------------------------- */
38
39	double lj_vm_foldarith(double x, double y, int op)	1,770✔
40	{
41	switch (op) {	1,770✔
42	case IR_ADD - IR_ADD: return x+y; break;	554✔
43	case IR_SUB - IR_ADD: return x-y; break;	143✔
44	case IR_MUL - IR_ADD: return x*y; break;	161✔
45	case IR_DIV - IR_ADD: return x/y; break;	54✔
46	case IR_MOD - IR_ADD: return x-lj_vm_floor(x/y)*y; break;	19✔
47	case IR_POW - IR_ADD: return pow(x, y); break;	744✔
48	case IR_NEG - IR_ADD: return -x; break;	13✔
49	case IR_ABS - IR_ADD: return fabs(x); break;	5✔
50	#if LJ_HASJIT
51	case IR_LDEXP - IR_ADD: return ldexp(x, (int)y); break;	41✔
52	case IR_MIN - IR_ADD: return x < y ? x : y; break;	4✔
53	case IR_MAX - IR_ADD: return x > y ? x : y; break;	32✔
54	#endif
55	default: return x;
56	}
57	}
58
59	/* -- Helper functions for generated machine code ------------------------- */
60
61	#if (LJ_HASJIT && !(LJ_TARGET_ARM \|\| LJ_TARGET_ARM64 \|\| LJ_TARGET_PPC)) \|\| LJ_TARGET_MIPS
62	int32_t LJ_FASTCALL lj_vm_modi(int32_t a, int32_t b)	6,892✔
63	{
64	uint32_t y, ua, ub;	6,892✔
65	/* This must be checked before using this function. */
66	lj_assertX(b != 0, "modulo with zero divisor");	6,892✔
67	ua = a < 0 ? (uint32_t)-a : (uint32_t)a;	6,892✔
68	ub = b < 0 ? (uint32_t)-b : (uint32_t)b;	6,892✔
69	y = ua % ub;	6,892✔
70	if (y != 0 && (a^b) < 0) y = y - ub;	6,892✔
71	if (((int32_t)y^b) < 0) y = (uint32_t)-(int32_t)y;	6,892✔
72	return (int32_t)y;	6,892✔
73	}
74	#endif
75
76	#if LJ_HASJIT
77
78	#ifdef LUAJIT_NO_LOG2
79	double lj_vm_log2(double a)
80	{
81	return log(a) * 1.4426950408889634074;
82	}
83	#endif
84
85	/* Computes fpm(x) for extended math functions. */
86	double lj_vm_foldfpm(double x, int fpm)	8✔
87	{
88	switch (fpm) {	8✔
89	case IRFPM_FLOOR: return lj_vm_floor(x);	4✔
90	case IRFPM_CEIL: return lj_vm_ceil(x);	3✔
91	case IRFPM_TRUNC: return lj_vm_trunc(x);	×
UNCOV 92	case IRFPM_SQRT: return sqrt(x);	×
93	case IRFPM_LOG: return log(x);	1✔
94	case IRFPM_LOG2: return lj_vm_log2(x);	×
95	default: lj_assertX(0, "bad fpm %d", fpm);
96	}
97	return 0;
98	}
99
100	#if LJ_HASFFI
101	int lj_vm_errno(void)	42✔
102	{
103	return errno;	42✔
104	}
105	#endif
106
107	#endif

tarantool / luajit / 6821269217

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous