"be consistent" PerlMonks

### Re^13: Algorithm for cancelling common factors between two lists of multiplicands

by tmoertel (Chaplain)
 on Aug 09, 2005 at 06:16 UTC Need Help??

How do I get that output reduced to a real in scientific format?
If you want arbitrary precision, the safest bet is to extract the digits yourself:
```module ScientificFormat (sciformat) where

import Ratio

-- return prec+1 digits because I don't round

sciformat prec r =
sign ++ digits ++ "e" ++ show exp
where
(ds, exp) = tosci (abs numer) (denominator r) 0
numer     = numerator r
sign      = if signum numer < 0 then "-" else "+"
digits    = show (head ds') ++ "." ++ concatMap show (tail ds')
ds'       = take (prec+1) (ds ++ repeat 0)

tosci 0     _     _   = ([0], 0)
tosci numer denom exp =
case numer `div` denom of
0           -> tosci (10 * numer) denom (exp-1)
x | x >= 10 -> tosci numer (10*denom) (exp+1)
_           -> (digits numer denom, exp)

digits n d =
let (q,r) = n `quotRem` d in q : digits (10*r) d
And:
```module Main (main) where
import FishersExactTest
import ScientificFormat
main = interact \$ sciformat 60 . pCutoff . map (map read . words) . li
+nes
Example:
```[thor@arinmir fishers-exact-test]\$ time ./fet < ex1.dat
+8.070604647867604097576877675243109941729476950993498765160880e-7030
real    0m0.839s
user    0m0.819s
sys     0m0.016s

Cheers,
Tom

Create A New User
Domain Nodelet?
Node Status?
node history
Node Type: note [id://482119]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others avoiding work at the Monastery: (6)
As of 2024-03-01 14:15 GMT
Voting Booth?
My favourite way to spend a leap day ...

Results (28 votes). Check out past polls.