OLD | NEW |
| (Empty) |
1 #!/usr/local/bin/perl | |
2 # ******************************************************************** | |
3 # * COPYRIGHT: | |
4 # * Copyright (c) 2002, International Business Machines Corporation and | |
5 # * others. All Rights Reserved. | |
6 # ******************************************************************** | |
7 | |
8 my $PLUS_MINUS = "±"; | |
9 | |
10 #|#--------------------------------------------------------------------- | |
11 #|# Format a confidence interval, as given by a Dataset. Output is as | |
12 #|# as follows: | |
13 #|# 241.23 - 241.98 => 241.5 +/- 0.3 | |
14 #|# 241.2 - 243.8 => 242 +/- 1 | |
15 #|# 211.0 - 241.0 => 226 +/- 15 or? 230 +/- 20 | |
16 #|# 220.3 - 234.3 => 227 +/- 7 | |
17 #|# 220.3 - 300.3 => 260 +/- 40 | |
18 #|# 220.3 - 1000 => 610 +/- 390 or? 600 +/- 400 | |
19 #|# 0.022 - 0.024 => 0.023 +/- 0.001 | |
20 #|# 0.022 - 0.032 => 0.027 +/- 0.005 | |
21 #|# 0.022 - 1.000 => 0.5 +/- 0.5 | |
22 #|# In other words, take one significant digit of the error value and | |
23 #|# display the mean to the same precision. | |
24 #|sub formatDataset { | |
25 #| my $ds = shift; | |
26 #| my $lower = $ds->getMean() - $ds->getError(); | |
27 #| my $upper = $ds->getMean() + $ds->getError(); | |
28 #| my $scale = 0; | |
29 #| # Find how many initial digits are the same | |
30 #| while ($lower < 1 || | |
31 #| int($lower) == int($upper)) { | |
32 #| $lower *= 10; | |
33 #| $upper *= 10; | |
34 #| $scale++; | |
35 #| } | |
36 #| while ($lower >= 10 && | |
37 #| int($lower) == int($upper)) { | |
38 #| $lower /= 10; | |
39 #| $upper /= 10; | |
40 #| $scale--; | |
41 #| } | |
42 #|} | |
43 | |
44 #--------------------------------------------------------------------- | |
45 # Format a number, optionally with a +/- delta, to n significant | |
46 # digits. | |
47 # | |
48 # @param significant digit, a value >= 1 | |
49 # @param multiplier | |
50 # @param time in seconds to be formatted | |
51 # @optional delta in seconds | |
52 # | |
53 # @return string of the form "23" or "23 +/- 10". | |
54 # | |
55 sub formatNumber { | |
56 my $sigdig = shift; | |
57 my $mult = shift; | |
58 my $a = shift; | |
59 my $delta = shift; # may be undef | |
60 | |
61 my $result = formatSigDig($sigdig, $a*$mult); | |
62 if (defined($delta)) { | |
63 my $d = formatSigDig($sigdig, $delta*$mult); | |
64 # restrict PRECISION of delta to that of main number | |
65 if ($result =~ /\.(\d+)/) { | |
66 # TODO make this work for values with all significant | |
67 # digits to the left of the decimal, e.g., 1234000. | |
68 | |
69 # TODO the other thing wrong with this is that it | |
70 # isn't rounding the $delta properly. Have to put | |
71 # this logic into formatSigDig(). | |
72 my $x = length($1); | |
73 $d =~ s/\.(\d{$x})\d+/.$1/; | |
74 } | |
75 $result .= " $PLUS_MINUS " . $d; | |
76 } | |
77 $result; | |
78 } | |
79 | |
80 #--------------------------------------------------------------------- | |
81 # Format a time, optionally with a +/- delta, to n significant | |
82 # digits. | |
83 # | |
84 # @param significant digit, a value >= 1 | |
85 # @param time in seconds to be formatted | |
86 # @optional delta in seconds | |
87 # | |
88 # @return string of the form "23 ms" or "23 +/- 10 ms". | |
89 # | |
90 sub formatSeconds { | |
91 my $sigdig = shift; | |
92 my $a = shift; | |
93 my $delta = shift; # may be undef | |
94 | |
95 my @MULT = (1 , 1e3, 1e6, 1e9); | |
96 my @SUFF = ('s' , 'ms', 'us', 'ns'); | |
97 | |
98 # Determine our scale | |
99 my $i = 0; | |
100 #always do seconds if the following line is commented out | |
101 ++$i while ($a*$MULT[$i] < 1 && $i < @MULT); | |
102 | |
103 formatNumber($sigdig, $MULT[$i], $a, $delta) . ' ' . $SUFF[$i]; | |
104 } | |
105 | |
106 #--------------------------------------------------------------------- | |
107 # Format a percentage, optionally with a +/- delta, to n significant | |
108 # digits. | |
109 # | |
110 # @param significant digit, a value >= 1 | |
111 # @param value to be formatted, as a fraction, e.g. 0.5 for 50% | |
112 # @optional delta, as a fraction | |
113 # | |
114 # @return string of the form "23 %" or "23 +/- 10 %". | |
115 # | |
116 sub formatPercent { | |
117 my $sigdig = shift; | |
118 my $a = shift; | |
119 my $delta = shift; # may be undef | |
120 | |
121 formatNumber($sigdig, 100, $a, $delta) . '%'; | |
122 } | |
123 | |
124 #--------------------------------------------------------------------- | |
125 # Format a number to n significant digits without using exponential | |
126 # notation. | |
127 # | |
128 # @param significant digit, a value >= 1 | |
129 # @param number to be formatted | |
130 # | |
131 # @return string of the form "1234" "12.34" or "0.001234". If | |
132 # number was negative, prefixed by '-'. | |
133 # | |
134 sub formatSigDig { | |
135 my $n = shift() - 1; | |
136 my $a = shift; | |
137 | |
138 local $_ = sprintf("%.${n}e", $a); | |
139 my $sign = (s/^-//) ? '-' : ''; | |
140 | |
141 my $a_e; | |
142 my $result; | |
143 if (/^(\d)\.(\d+)e([-+]\d+)$/) { | |
144 my ($d, $dn, $e) = ($1, $2, $3); | |
145 $a_e = $e; | |
146 $d .= $dn; | |
147 $e++; | |
148 $d .= '0' while ($e > length($d)); | |
149 while ($e < 1) { | |
150 $e++; | |
151 $d = '0' . $d; | |
152 } | |
153 if ($e == length($d)) { | |
154 $result = $sign . $d; | |
155 } else { | |
156 $result = $sign . substr($d, 0, $e) . '.' . substr($d, $e); | |
157 } | |
158 } else { | |
159 die "Can't parse $_"; | |
160 } | |
161 $result; | |
162 } | |
163 | |
164 1; | |
165 | |
166 #eof | |
OLD | NEW |