#!/bin/bash
# ULP error check script.
#
# Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
# See https://llvm.org/LICENSE.txt for license information.
# SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
#set -x
set -eu
# cd to bin directory.
cd "${0%/*}"
rmodes='n u d z'
#rmodes=n
flags="${ULPFLAGS:--q}"
emu="$@"
FAIL=0
PASS=0
t() {
[ $r = "n" ] && Lt=$L || Lt=$Ldir
$emu ./ulp -r $r -e $Lt $flags "$@" && PASS=$((PASS+1)) || FAIL=$((FAIL+1))
}
check() {
$emu ./ulp -f -q "$@" >/dev/null
}
Ldir=0.5
for r in $rmodes
do
L=0.01
t exp 0 0xffff000000000000 10000
t exp 0x1p-6 0x1p6 40000
t exp -0x1p-6 -0x1p6 40000
t exp 633.3 733.3 10000
t exp -633.3 -777.3 10000
L=0.01
t exp2 0 0xffff000000000000 10000
t exp2 0x1p-6 0x1p6 40000
t exp2 -0x1p-6 -0x1p6 40000
t exp2 633.3 733.3 10000
t exp2 -633.3 -777.3 10000
L=0.02
t log 0 0xffff000000000000 10000
t log 0x1p-4 0x1p4 40000
t log 0 inf 40000
L=0.05
t log2 0 0xffff000000000000 10000
t log2 0x1p-4 0x1p4 40000
t log2 0 inf 40000
L=0.05
t pow 0.5 2.0 x 0 inf 20000
t pow -0.5 -2.0 x 0 inf 20000
t pow 0.5 2.0 x -0 -inf 20000
t pow -0.5 -2.0 x -0 -inf 20000
t pow 0.5 2.0 x 0x1p-10 0x1p10 40000
t pow 0.5 2.0 x -0x1p-10 -0x1p10 40000
t pow 0 inf x 0.5 2.0 80000
t pow 0 inf x -0.5 -2.0 80000
t pow 0x1.fp-1 0x1.08p0 x 0x1p8 0x1p17 80000
t pow 0x1.fp-1 0x1.08p0 x -0x1p8 -0x1p17 80000
t pow 0 0x1p-1000 x 0 1.0 50000
t pow 0x1p1000 inf x 0 1.0 50000
t pow 0x1.ffffffffffff0p-1 0x1.0000000000008p0 x 0x1p60 0x1p68 50000
t pow 0x1.ffffffffff000p-1 0x1p0 x 0x1p50 0x1p52 50000
t pow -0x1.ffffffffff000p-1 -0x1p0 x 0x1p50 0x1p52 50000
L=0.01
t expf 0 0xffff0000 10000
t expf 0x1p-14 0x1p8 50000
t expf -0x1p-14 -0x1p8 50000
L=0.01
t exp2f 0 0xffff0000 10000
t exp2f 0x1p-14 0x1p8 50000
t exp2f -0x1p-14 -0x1p8 50000
L=0.32
t logf 0 0xffff0000 10000
t logf 0x1p-4 0x1p4 50000
t logf 0 inf 50000
L=0.26
t log2f 0 0xffff0000 10000
t log2f 0x1p-4 0x1p4 50000
t log2f 0 inf 50000
#L=0.06
#t sinf 0 0xffff0000 10000
#t sinf 0x1p-14 0x1p54 50000
#t sinf -0x1p-14 -0x1p54 50000
#
#L=0.06
#t cosf 0 0xffff0000 10000
#t cosf 0x1p-14 0x1p54 50000
#t cosf -0x1p-14 -0x1p54 50000
#
#L=0.06
#t sincosf_sinf 0 0xffff0000 10000
#t sincosf_sinf 0x1p-14 0x1p54 50000
#t sincosf_sinf -0x1p-14 -0x1p54 50000
#
#L=0.06
#t sincosf_cosf 0 0xffff0000 10000
#t sincosf_cosf 0x1p-14 0x1p54 50000
#t sincosf_cosf -0x1p-14 -0x1p54 50000
L=0.4
t powf 0x1p-1 0x1p1 x 0x1p-7 0x1p7 50000
t powf 0x1p-1 0x1p1 x -0x1p-7 -0x1p7 50000
t powf 0x1p-70 0x1p70 x 0x1p-1 0x1p1 50000
t powf 0x1p-70 0x1p70 x -0x1p-1 -0x1p1 50000
t powf 0x1.ep-1 0x1.1p0 x 0x1p8 0x1p14 50000
t powf 0x1.ep-1 0x1.1p0 x -0x1p8 -0x1p14 50000
done
# vector functions
Ldir=0.5
r='n'
flags="${ULPFLAGS:--q} -f"
runs=
check __s_exp 1 && runs=1
runv=
check __v_exp 1 && runv=1
runvn=
check __vn_exp 1 && runvn=1
range_exp='
0 0xffff000000000000 10000
0x1p-6 0x1p6 400000
-0x1p-6 -0x1p6 400000
633.3 733.3 10000
-633.3 -777.3 10000
'
range_log='
0 0xffff000000000000 10000
0x1p-4 0x1p4 400000
0 inf 400000
'
range_pow='
0x1p-1 0x1p1 x 0x1p-10 0x1p10 50000
0x1p-1 0x1p1 x -0x1p-10 -0x1p10 50000
0x1p-500 0x1p500 x 0x1p-1 0x1p1 50000
0x1p-500 0x1p500 x -0x1p-1 -0x1p1 50000
0x1.ep-1 0x1.1p0 x 0x1p8 0x1p16 50000
0x1.ep-1 0x1.1p0 x -0x1p8 -0x1p16 50000
'
range_sin='
0 0xffff000000000000 10000
0x1p-4 0x1p4 400000
-0x1p-23 0x1p23 400000
'
range_cos="$range_sin"
range_expf='
0 0xffff0000 10000
0x1p-14 0x1p8 500000
-0x1p-14 -0x1p8 500000
'
range_expf_1u="$range_expf"
range_exp2f="$range_expf"
range_exp2f_1u="$range_expf"
range_logf='
0 0xffff0000 10000
0x1p-4 0x1p4 500000
'
range_sinf='
0 0xffff0000 10000
0x1p-4 0x1p4 300000
-0x1p-9 -0x1p9 300000
'
range_cosf="$range_sinf"
range_powf='
0x1p-1 0x1p1 x 0x1p-7 0x1p7 50000
0x1p-1 0x1p1 x -0x1p-7 -0x1p7 50000
0x1p-70 0x1p70 x 0x1p-1 0x1p1 50000
0x1p-70 0x1p70 x -0x1p-1 -0x1p1 50000
0x1.ep-1 0x1.1p0 x 0x1p8 0x1p14 50000
0x1.ep-1 0x1.1p0 x -0x1p8 -0x1p14 50000
'
# error limits
L_exp=1.9
L_log=1.2
L_pow=0.05
L_sin=3.0
L_cos=3.0
L_expf=1.49
L_expf_1u=0.4
L_exp2f=1.49
L_exp2f_1u=0.4
L_logf=2.9
L_sinf=1.4
L_cosf=1.4
L_powf=2.1
while read G F R
do
[ "$R" = 1 ] || continue
case "$G" in \#*) continue ;; esac
eval range="\${range_$G}"
eval L="\${L_$G}"
while read X
do
[ -n "$X" ] || continue
case "$X" in \#*) continue ;; esac
t $F $X
done << EOF
$range
EOF
done << EOF
# group symbol run
exp __s_exp $runs
exp __v_exp $runv
exp __vn_exp $runvn
exp _ZGVnN2v_exp $runvn
log __s_log $runs
log __v_log $runv
log __vn_log $runvn
log _ZGVnN2v_log $runvn
pow __s_pow $runs
pow __v_pow $runv
pow __vn_pow $runvn
pow _ZGVnN2vv_pow $runvn
sin __s_sin $runs
sin __v_sin $runv
sin __vn_sin $runvn
sin _ZGVnN2v_sin $runvn
cos __s_cos $runs
cos __v_cos $runv
cos __vn_cos $runvn
cos _ZGVnN2v_cos $runvn
expf __s_expf $runs
expf __v_expf $runv
expf __vn_expf $runvn
expf _ZGVnN4v_expf $runvn
expf_1u __s_expf_1u $runs
expf_1u __v_expf_1u $runv
expf_1u __vn_expf_1u $runvn
exp2f __s_exp2f $runs
exp2f __v_exp2f $runv
exp2f __vn_exp2f $runvn
exp2f _ZGVnN4v_exp2f $runvn
exp2f_1u __s_exp2f_1u $runs
exp2f_1u __v_exp2f_1u $runv
exp2f_1u __vn_exp2f_1u $runvn
logf __s_logf $runs
logf __v_logf $runv
logf __vn_logf $runvn
logf _ZGVnN4v_logf $runvn
sinf __s_sinf $runs
sinf __v_sinf $runv
sinf __vn_sinf $runvn
sinf _ZGVnN4v_sinf $runvn
cosf __s_cosf $runs
cosf __v_cosf $runv
cosf __vn_cosf $runvn
cosf _ZGVnN4v_cosf $runvn
powf __s_powf $runs
powf __v_powf $runv
powf __vn_powf $runvn
powf _ZGVnN4vv_powf $runvn
EOF
[ 0 -eq $FAIL ] || {
echo "FAILED $FAIL PASSED $PASS"
exit 1
}
Codebase Browser
llvm