mirror of
git://sourceware.org/git/glibc.git
synced 2024-12-15 04:20:28 +08:00
e6b6a85705
As recently discussed <https://sourceware.org/ml/libc-alpha/2014-02/msg00670.html>, it doesn't seem particularly useful for libm-test-ulps files to contain huge amounts of data on ulps for individual tests; just the global maximum observed ulps for each function, together with the verification of exceptions, errno and special results such as infinities and NaNs for each test, suffices to verify that a function's behavior on the given test inputs is within the expected accuracy. Removing this data reduces source tree churn caused by updates to these files when libm tests are added, and reduces the frequency with which testsuite additions actually need libm-test-ulps changes at all. Accordingly, this patch removes that data, so that individual tests get checked against the global bounds for the given function and only generate an error if those are exceeded. Tested x86_64 (including verifying that if an ulps value is artificially reduced, the tests do indeed fail as they should and "make regen-ulps" generates the expected changes). * math/libm-test.inc (struct ulp_data): Don't refer to ulps for individual tests in comment. (libm-test-ulps.h): Don't refer to test_ulps in #include comment. (prev_max_error): New variable. (prev_real_max_error): Likewise. (prev_imag_max_error): Likewise. (compare_ulp_data): Don't refer to test names in comment. (find_test_ulps): Remove function. (find_function_ulps): Likewise. (find_complex_function_ulps): Likewise. (init_max_error): Take function name as argument. Look up ulps for that function. (print_ulps): Remove function. (print_max_error): Use prev_max_error instead of calling find_function_ulps. (print_complex_max_error): Use prev_real_max_error and prev_imag_max_error instead of calling find_complex_function_ulps. (check_float_internal): Take max_ulp parameter instead of calling find_test_ulps. Don't call print_ulps. (check_float): Update call to check_float_internal. (check_complex): Update calls to check_float_internal. (START): Pass argument to init_max_error. * math/gen-libm-test.pl (%results): Don't include "kind" information. (parse_ulps): Don't handle ulps of individual tests. (print_ulps_file): Likewise. (output_ulps): Likewise. * math/README.libm-test: Update. * manual/libm-err-tab.pl (parse_ulps): Don't handle ulps of individual tests. * sysdeps/aarch64/libm-test-ulps: Remove individual test ulps. * sysdeps/alpha/fpu/libm-test-ulps: Likewise. * sysdeps/arm/libm-test-ulps: Likewise. * sysdeps/i386/fpu/libm-test-ulps: Likewise. * sysdeps/ia64/fpu/libm-test-ulps: Likewise. * sysdeps/m68k/coldfire/fpu/libm-test-ulps: Likewise. * sysdeps/m68k/m680x0/fpu/libm-test-ulps: Likewise. * sysdeps/microblaze/libm-test-ulps: Likewise. * sysdeps/mips/mips32/libm-test-ulps: Likewise. * sysdeps/mips/mips64/libm-test-ulps: Likewise. * sysdeps/powerpc/fpu/libm-test-ulps: Likewise. * sysdeps/powerpc/nofpu/libm-test-ulps: Likewise. * sysdeps/s390/fpu/libm-test-ulps: Likewise. * sysdeps/sh/libm-test-ulps: Likewise. * sysdeps/sparc/fpu/libm-test-ulps: Likewise. * sysdeps/tile/libm-test-ulps: Likewise. * sysdeps/x86_64/fpu/libm-test-ulps: Likewise. * sysdeps/hppa/fpu/libm-test-ulps: Remove individual test ulps.
234 lines
6.6 KiB
Perl
Executable File
234 lines
6.6 KiB
Perl
Executable File
#!/usr/bin/perl -w
|
|
# Copyright (C) 1999-2014 Free Software Foundation, Inc.
|
|
# This file is part of the GNU C Library.
|
|
# Contributed by Andreas Jaeger <aj@suse.de>, 1999.
|
|
|
|
# The GNU C Library 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.
|
|
|
|
# The GNU C Library 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 the GNU C Library; if not, see
|
|
# <http://www.gnu.org/licenses/>.
|
|
|
|
# Information about tests are stored in: %results
|
|
# $results{$test}{"type"} is the result type, e.g. normal or complex.
|
|
# In the following description $platform, $type and $float are:
|
|
# - $platform is the used platform
|
|
# - $type is either "normal", "real" (for the real part of a complex number)
|
|
# or "imag" (for the imaginary part # of a complex number).
|
|
# - $float is either of float, ifloat, double, idouble, ldouble, ildouble;
|
|
# It represents the underlying floating point type (float, double or long
|
|
# double) and if inline functions (the leading i stands for inline)
|
|
# are used.
|
|
# $results{$test}{$platform}{$type}{$float} is defined and has a delta
|
|
# or 'fail' as value.
|
|
|
|
use File::Find;
|
|
|
|
use strict;
|
|
|
|
use vars qw ($sources @platforms %pplatforms);
|
|
use vars qw (%results @all_floats %suffices @all_functions);
|
|
|
|
|
|
# all_floats is in output order and contains all recognised float types that
|
|
# we're going to output
|
|
@all_floats = ('float', 'double', 'ldouble');
|
|
%suffices =
|
|
( 'float' => 'f',
|
|
'double' => '',
|
|
'ldouble' => 'l'
|
|
);
|
|
|
|
# Pretty description of platform
|
|
%pplatforms =
|
|
( "i386/fpu" => "ix86",
|
|
"generic" => "Generic",
|
|
"alpha/fpu" => "Alpha",
|
|
"ia64/fpu" => "IA64",
|
|
"m68k/fpu" => "M68k",
|
|
"mips/fpu" => "MIPS",
|
|
"powerpc/fpu" => "PowerPC",
|
|
"sparc/sparc32/fpu" => "Sparc 32-bit",
|
|
"sparc/sparc64/fpu" => "Sparc 64-bit",
|
|
"sh/sh4/fpu" => "SH4",
|
|
"s390/fpu" => "S/390",
|
|
"arm" => "ARM"
|
|
);
|
|
|
|
@all_functions =
|
|
( "acos", "acosh", "asin", "asinh", "atan", "atanh",
|
|
"atan2", "cabs", "cacos", "cacosh", "carg", "casin", "casinh",
|
|
"catan", "catanh", "cbrt", "ccos", "ccosh", "ceil", "cexp", "cimag",
|
|
"clog", "clog10", "conj", "copysign", "cos", "cosh", "cpow", "cproj",
|
|
"creal", "csin", "csinh", "csqrt", "ctan", "ctanh", "erf", "erfc",
|
|
"exp", "exp10", "exp2", "expm1", "fabs", "fdim", "floor", "fma",
|
|
"fmax", "fmin", "fmod", "frexp", "gamma", "hypot",
|
|
"ilogb", "j0", "j1", "jn", "lgamma", "lrint",
|
|
"llrint", "log", "log10", "log1p", "log2", "logb", "lround",
|
|
"llround", "modf", "nearbyint", "nextafter", "nexttoward", "pow",
|
|
"remainder", "remquo", "rint", "round", "scalb", "scalbn", "scalbln",
|
|
"sin", "sincos", "sinh", "sqrt", "tan", "tanh", "tgamma",
|
|
"trunc", "y0", "y1", "yn" );
|
|
# fpclassify, isnormal, isfinite, isinf, isnan, issignaling, signbit,
|
|
# isgreater, isgreaterequal, isless, islessequal, islessgreater, isunordered
|
|
# are not tabulated.
|
|
|
|
if ($#ARGV == 0) {
|
|
$sources = $ARGV[0];
|
|
} else {
|
|
$sources = '/usr/src/cvs/libc';
|
|
}
|
|
|
|
find (\&find_files, $sources);
|
|
|
|
@platforms = sort by_platforms @platforms;
|
|
|
|
&print_all;
|
|
|
|
sub find_files {
|
|
if ($_ eq 'libm-test-ulps') {
|
|
# print "Parsing $File::Find::name\n";
|
|
push @platforms, $File::Find::dir;
|
|
&parse_ulps ($File::Find::name, $File::Find::dir);
|
|
}
|
|
}
|
|
|
|
# Parse ulps file
|
|
sub parse_ulps {
|
|
my ($file, $platform) = @_;
|
|
my ($test, $type, $float, $eps);
|
|
|
|
# $type has the following values:
|
|
# "normal": No complex variable
|
|
# "real": Real part of complex result
|
|
# "imag": Imaginary part of complex result
|
|
open ULP, $file or die ("Can't open $file: $!");
|
|
while (<ULP>) {
|
|
chop;
|
|
# ignore comments and empty lines
|
|
next if /^#/;
|
|
next if /^\s*$/;
|
|
if (/^Function: /) {
|
|
if (/Real part of/) {
|
|
s/Real part of //;
|
|
$type = 'real';
|
|
} elsif (/Imaginary part of/) {
|
|
s/Imaginary part of //;
|
|
$type = 'imag';
|
|
} else {
|
|
$type = 'normal';
|
|
}
|
|
($test) = ($_ =~ /^Function:\s*\"([a-zA-Z0-9_]+)\"/);
|
|
next;
|
|
}
|
|
if (/^i?(float|double|ldouble):/) {
|
|
($float, $eps) = split /\s*:\s*/,$_,2;
|
|
if ($eps eq 'fail') {
|
|
$results{$test}{$platform}{$type}{$float} = 'fail';
|
|
} elsif ($eps eq "0") {
|
|
# ignore
|
|
next;
|
|
} elsif (!exists $results{$test}{$platform}{$type}{$float}
|
|
|| $results{$test}{$platform}{$type}{$float} ne 'fail') {
|
|
$results{$test}{$platform}{$type}{$float} = $eps;
|
|
}
|
|
if ($type =~ /^real|imag$/) {
|
|
$results{$test}{'type'} = 'complex';
|
|
} elsif ($type eq 'normal') {
|
|
$results{$test}{'type'} = 'normal';
|
|
}
|
|
next;
|
|
}
|
|
print "Skipping unknown entry: `$_'\n";
|
|
}
|
|
close ULP;
|
|
}
|
|
|
|
sub get_value {
|
|
my ($fct, $platform, $type, $float) = @_;
|
|
|
|
return (exists $results{$fct}{$platform}{$type}{$float}
|
|
? $results{$fct}{$platform}{$type}{$float} : "0");
|
|
}
|
|
|
|
sub canonicalize_platform {
|
|
my ($platform) = @_;
|
|
|
|
$platform =~ s|^(.*/sysdeps/)||;
|
|
|
|
|
|
return exists $pplatforms{$platform} ? $pplatforms{$platform} : $platform;
|
|
}
|
|
|
|
sub print_platforms {
|
|
my (@p) = @_;
|
|
my ($fct, $platform, $float, $first, $i, $platform_no, $platform_total);
|
|
|
|
print '@multitable {nexttowardf} ';
|
|
foreach (@p) {
|
|
print ' {1000 + i 1000}';
|
|
}
|
|
print "\n";
|
|
|
|
print '@item Function ';
|
|
foreach (@p) {
|
|
print ' @tab ';
|
|
print &canonicalize_platform ($_);
|
|
}
|
|
print "\n";
|
|
|
|
|
|
foreach $fct (@all_functions) {
|
|
foreach $float (@all_floats) {
|
|
print "\@item $fct$suffices{$float} ";
|
|
foreach $platform (@p) {
|
|
print ' @tab ';
|
|
if (exists $results{$fct}{$platform}{'normal'}{$float}
|
|
|| exists $results{$fct}{$platform}{'real'}{$float}
|
|
|| exists $results{$fct}{$platform}{'imag'}{$float}) {
|
|
if ($results{$fct}{'type'} eq 'complex') {
|
|
print &get_value ($fct, $platform, 'real', $float),
|
|
' + i ', &get_value ($fct, $platform, 'imag', $float);
|
|
} else {
|
|
print $results{$fct}{$platform}{'normal'}{$float};
|
|
}
|
|
} else {
|
|
print '-';
|
|
}
|
|
}
|
|
print "\n";
|
|
}
|
|
}
|
|
|
|
print "\@end multitable\n";
|
|
}
|
|
|
|
sub print_all {
|
|
my ($i, $max);
|
|
|
|
my ($columns) = 5;
|
|
|
|
# Print only 5 platforms at a time.
|
|
for ($i=0; $i < $#platforms; $i+=$columns) {
|
|
$max = $i+$columns-1 > $#platforms ? $#platforms : $i+$columns-1;
|
|
print_platforms (@platforms[$i .. $max]);
|
|
}
|
|
}
|
|
|
|
sub by_platforms {
|
|
my ($pa, $pb);
|
|
|
|
$pa = $pplatforms{$a} ? $pplatforms{$a} : $a;
|
|
$pb = $pplatforms{$b} ? $pplatforms{$b} : $b;
|
|
|
|
return $pa cmp $pb;
|
|
}
|