167 lines
4.7 KiB
Perl
167 lines
4.7 KiB
Perl
#!/usr/local/bin/perl
|
|
# ********************************************************************
|
|
# * COPYRIGHT:
|
|
# * Copyright (c) 2002, International Business Machines Corporation and
|
|
# * others. All Rights Reserved.
|
|
# ********************************************************************
|
|
|
|
my $PLUS_MINUS = "±";
|
|
|
|
#|#---------------------------------------------------------------------
|
|
#|# Format a confidence interval, as given by a Dataset. Output is as
|
|
#|# as follows:
|
|
#|# 241.23 - 241.98 => 241.5 +/- 0.3
|
|
#|# 241.2 - 243.8 => 242 +/- 1
|
|
#|# 211.0 - 241.0 => 226 +/- 15 or? 230 +/- 20
|
|
#|# 220.3 - 234.3 => 227 +/- 7
|
|
#|# 220.3 - 300.3 => 260 +/- 40
|
|
#|# 220.3 - 1000 => 610 +/- 390 or? 600 +/- 400
|
|
#|# 0.022 - 0.024 => 0.023 +/- 0.001
|
|
#|# 0.022 - 0.032 => 0.027 +/- 0.005
|
|
#|# 0.022 - 1.000 => 0.5 +/- 0.5
|
|
#|# In other words, take one significant digit of the error value and
|
|
#|# display the mean to the same precision.
|
|
#|sub formatDataset {
|
|
#| my $ds = shift;
|
|
#| my $lower = $ds->getMean() - $ds->getError();
|
|
#| my $upper = $ds->getMean() + $ds->getError();
|
|
#| my $scale = 0;
|
|
#| # Find how many initial digits are the same
|
|
#| while ($lower < 1 ||
|
|
#| int($lower) == int($upper)) {
|
|
#| $lower *= 10;
|
|
#| $upper *= 10;
|
|
#| $scale++;
|
|
#| }
|
|
#| while ($lower >= 10 &&
|
|
#| int($lower) == int($upper)) {
|
|
#| $lower /= 10;
|
|
#| $upper /= 10;
|
|
#| $scale--;
|
|
#| }
|
|
#|}
|
|
|
|
#---------------------------------------------------------------------
|
|
# Format a number, optionally with a +/- delta, to n significant
|
|
# digits.
|
|
#
|
|
# @param significant digit, a value >= 1
|
|
# @param multiplier
|
|
# @param time in seconds to be formatted
|
|
# @optional delta in seconds
|
|
#
|
|
# @return string of the form "23" or "23 +/- 10".
|
|
#
|
|
sub formatNumber {
|
|
my $sigdig = shift;
|
|
my $mult = shift;
|
|
my $a = shift;
|
|
my $delta = shift; # may be undef
|
|
|
|
my $result = formatSigDig($sigdig, $a*$mult);
|
|
if (defined($delta)) {
|
|
my $d = formatSigDig($sigdig, $delta*$mult);
|
|
# restrict PRECISION of delta to that of main number
|
|
if ($result =~ /\.(\d+)/) {
|
|
# TODO make this work for values with all significant
|
|
# digits to the left of the decimal, e.g., 1234000.
|
|
|
|
# TODO the other thing wrong with this is that it
|
|
# isn't rounding the $delta properly. Have to put
|
|
# this logic into formatSigDig().
|
|
my $x = length($1);
|
|
$d =~ s/\.(\d{$x})\d+/.$1/;
|
|
}
|
|
$result .= " $PLUS_MINUS " . $d;
|
|
}
|
|
$result;
|
|
}
|
|
|
|
#---------------------------------------------------------------------
|
|
# Format a time, optionally with a +/- delta, to n significant
|
|
# digits.
|
|
#
|
|
# @param significant digit, a value >= 1
|
|
# @param time in seconds to be formatted
|
|
# @optional delta in seconds
|
|
#
|
|
# @return string of the form "23 ms" or "23 +/- 10 ms".
|
|
#
|
|
sub formatSeconds {
|
|
my $sigdig = shift;
|
|
my $a = shift;
|
|
my $delta = shift; # may be undef
|
|
|
|
my @MULT = (1 , 1e3, 1e6, 1e9);
|
|
my @SUFF = ('s' , 'ms', 'us', 'ns');
|
|
|
|
# Determine our scale
|
|
my $i = 0;
|
|
#always do seconds if the following line is commented out
|
|
++$i while ($a*$MULT[$i] < 1 && $i < @MULT);
|
|
|
|
formatNumber($sigdig, $MULT[$i], $a, $delta) . ' ' . $SUFF[$i];
|
|
}
|
|
|
|
#---------------------------------------------------------------------
|
|
# Format a percentage, optionally with a +/- delta, to n significant
|
|
# digits.
|
|
#
|
|
# @param significant digit, a value >= 1
|
|
# @param value to be formatted, as a fraction, e.g. 0.5 for 50%
|
|
# @optional delta, as a fraction
|
|
#
|
|
# @return string of the form "23 %" or "23 +/- 10 %".
|
|
#
|
|
sub formatPercent {
|
|
my $sigdig = shift;
|
|
my $a = shift;
|
|
my $delta = shift; # may be undef
|
|
|
|
formatNumber($sigdig, 100, $a, $delta) . '%';
|
|
}
|
|
|
|
#---------------------------------------------------------------------
|
|
# Format a number to n significant digits without using exponential
|
|
# notation.
|
|
#
|
|
# @param significant digit, a value >= 1
|
|
# @param number to be formatted
|
|
#
|
|
# @return string of the form "1234" "12.34" or "0.001234". If
|
|
# number was negative, prefixed by '-'.
|
|
#
|
|
sub formatSigDig {
|
|
my $n = shift() - 1;
|
|
my $a = shift;
|
|
|
|
local $_ = sprintf("%.${n}e", $a);
|
|
my $sign = (s/^-//) ? '-' : '';
|
|
|
|
my $a_e;
|
|
my $result;
|
|
if (/^(\d)\.(\d+)e([-+]\d+)$/) {
|
|
my ($d, $dn, $e) = ($1, $2, $3);
|
|
$a_e = $e;
|
|
$d .= $dn;
|
|
$e++;
|
|
$d .= '0' while ($e > length($d));
|
|
while ($e < 1) {
|
|
$e++;
|
|
$d = '0' . $d;
|
|
}
|
|
if ($e == length($d)) {
|
|
$result = $sign . $d;
|
|
} else {
|
|
$result = $sign . substr($d, 0, $e) . '.' . substr($d, $e);
|
|
}
|
|
} else {
|
|
die "Can't parse $_";
|
|
}
|
|
$result;
|
|
}
|
|
|
|
1;
|
|
|
|
#eof
|