Gecko [ armaneco.co.uk ]

Name	Size	Permission	Date
App	[ DIR ]	drwxr-xr-x	2021-08-24 16:58
Archive	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
B	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Bundle	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Business	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
CGI	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
CPAN	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
CPANPLUS	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Carp	[ DIR ]	drwxr-xr-x	2021-08-24 16:57
DBIx	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Date	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Devel	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Digest	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Encode	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Error	[ DIR ]	drwxr-xr-x	2021-08-24 17:01
Exporter	[ DIR ]	drwxr-xr-x	2021-08-24 16:57
ExtUtils	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
File	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Getopt	[ DIR ]	drwxr-xr-x	2021-08-24 16:57
Git	[ DIR ]	drwxr-xr-x	2026-02-03 06:51
HTML	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
HTTP	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
IO	[ DIR ]	drwxr-xr-x	2018-04-11 04:22
IPC	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
JSON	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
LWP	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Locale	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Log	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Module	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Mozilla	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Net	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
POD2	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Params	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Parse	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Perl	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Perl4	[ DIR ]	drwxr-xr-x	2021-08-25 16:57
Pod	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
RPC	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
TAP	[ DIR ]	drwxr-xr-x	2021-08-24 16:58
Term	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Test	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Text	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Thread	[ DIR ]	drwxr-xr-x	2021-08-24 16:58
Time	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
Try	[ DIR ]	drwxr-xr-x	2022-04-06 12:48
Types	[ DIR ]	drwxr-xr-x	2022-04-06 12:48
Version	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
WWW	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
autodie	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
common	[ DIR ]	drwxr-xr-x	2022-04-06 12:48
inc	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
lib	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
local	[ DIR ]	drwxr-xr-x	2021-08-24 17:00
CGI.pm	255.24 KB	-rw-r--r--	2012-11-15 03:12
CPANPLUS.pm	7.05 KB	-rw-r--r--	2013-05-21 02:16
Carp.pm	19.93 KB	-rw-r--r--	2012-06-18 16:02
Digest.pm	10.35 KB	-rw-r--r--	2011-10-02 13:36
Env.pm	5.39 KB	-rw-r--r--	2013-03-02 20:40
Error.pm	24.27 KB	-rw-r--r--	2014-06-10 00:49
Expect.pm	53.44 KB	-rw-r--r--	2007-07-19 15:53
Expect.pod	42.71 KB	-rw-r--r--	2007-07-19 16:04
Exporter.pm	18.31 KB	-rw-r--r--	2013-03-27 20:04
Fatal.pm	45.16 KB	-rw-r--r--	2013-02-22 18:19
Git.pm	42.55 KB	-r--r--r--	2013-06-11 00:31
abbrev.pl	819 B	-rw-r--r--	2017-06-14 04:08
assert.pl	1.27 KB	-rw-r--r--	2017-06-14 04:08
autodie.pm	11.83 KB	-rw-r--r--	2013-02-22 18:19
bigfloat.pl	7.14 KB	-rw-r--r--	2017-06-14 04:08
bigint.pl	8.71 KB	-rw-r--r--	2017-06-14 04:08
bigrat.pl	4.35 KB	-rw-r--r--	2017-06-14 04:08
cacheout.pl	1.1 KB	-rw-r--r--	2017-06-14 04:08
chat2.pl	9.73 KB	-rw-r--r--	2017-06-14 04:08
complete.pl	3.12 KB	-rw-r--r--	2017-06-14 04:08
constant.pm	13.04 KB	-rw-r--r--	2014-06-10 07:21
ctime.pl	1.93 KB	-rw-r--r--	2017-06-14 04:08
dotsh.pl	2.12 KB	-rw-r--r--	2017-06-14 04:08
exceptions.pl	1.7 KB	-rw-r--r--	2017-06-14 04:08
fastcwd.pl	1019 B	-rw-r--r--	2017-06-14 04:08
find.pl	1.16 KB	-rw-r--r--	2017-06-14 04:08
finddepth.pl	1.1 KB	-rw-r--r--	2017-06-14 04:08
flush.pl	642 B	-rw-r--r--	2017-06-14 04:08
ftp.pl	23.53 KB	-rw-r--r--	2017-06-14 04:08
getcwd.pl	1.37 KB	-rw-r--r--	2017-06-14 04:08
getopt.pl	1.27 KB	-rw-r--r--	2017-06-14 04:08
getopts.pl	1.35 KB	-rw-r--r--	2017-06-14 04:08
hostname.pl	727 B	-rw-r--r--	2017-06-14 04:08
importenv.pl	283 B	-rw-r--r--	2017-06-14 04:08
look.pl	1.23 KB	-rw-r--r--	2017-06-14 04:08
newgetopt.pl	2.15 KB	-rw-r--r--	2010-07-09 16:56
open2.pl	185 B	-rw-r--r--	2017-06-14 04:08
open3.pl	185 B	-rw-r--r--	2017-06-14 04:08
parent.pm	2.83 KB	-rw-r--r--	2011-03-08 20:56
perldoc.pod	8.25 KB	-rw-r--r--	2013-01-28 07:02
pwd.pl	1.44 KB	-rw-r--r--	2017-06-14 04:08
shellwords.pl	280 B	-rw-r--r--	2017-06-14 04:08
stat.pl	525 B	-rw-r--r--	2017-06-14 04:08
syslog.pl	4.69 KB	-rw-r--r--	2017-06-14 04:08
tainted.pl	164 B	-rw-r--r--	2017-06-14 04:08
termcap.pl	4.02 KB	-rw-r--r--	2017-06-14 04:08
timelocal.pl	690 B	-rw-r--r--	2017-06-14 04:08
validate.pl	3.64 KB	-rw-r--r--	2017-06-14 04:08

Rename

package bigrat;
require "bigint.pl";
#
# This library is no longer being maintained, and is included for backward
# compatibility with Perl 4 programs which may require it.
#
# In particular, this should not be used as an example of modern Perl
# programming techniques.
#
# Arbitrary size rational math package
#
# by Mark Biggar
#
# Input values to these routines consist of strings of the form 
#   m|^\s*[+-]?[\d\s]+(/[\d\s]+)?$|.
# Examples:
#   "+0/1"                          canonical zero value
#   "3"                             canonical value "+3/1"
#   "   -123/123 123"               canonical value "-1/1001"
#   "123 456/7890"                  canonical value "+20576/1315"
# Output values always include a sign and no leading zeros or
#   white space.
# This package makes use of the bigint package.
# The string 'NaN' is used to represent the result when input arguments 
#   that are not numbers, as well as the result of dividing by zero and
#       the sqrt of a negative number.
# Extreamly naive algorthims are used.
#
# Routines provided are:
#
#   rneg(RAT) return RAT                negation
#   rabs(RAT) return RAT                absolute value
#   rcmp(RAT,RAT) return CODE           compare numbers (undef,<0,=0,>0)
#   radd(RAT,RAT) return RAT            addition
#   rsub(RAT,RAT) return RAT            subtraction
#   rmul(RAT,RAT) return RAT            multiplication
#   rdiv(RAT,RAT) return RAT            division
#   rmod(RAT) return (RAT,RAT)          integer and fractional parts
#   rnorm(RAT) return RAT               normalization
#   rsqrt(RAT, cycles) return RAT       square root

# Convert a number to the canonical string form m|^[+-]\d+/\d+|.
sub main'rnorm { #(string) return rat_num
    local($_) = @_;
    s/\s+//g;
    if (m#^([+-]?\d+)(/(\d*[1-9]0*))?$#) {
	&norm($1, $3 ? $3 : '+1');
    } else {
	'NaN';
    }
}

# Normalize by reducing to lowest terms
sub norm { #(bint, bint) return rat_num
    local($num,$dom) = @_;
    if ($num eq 'NaN') {
	'NaN';
    } elsif ($dom eq 'NaN') {
	'NaN';
    } elsif ($dom =~ /^[+-]?0+$/) {
	'NaN';
    } else {
	local($gcd) = &'bgcd($num,$dom);
	$gcd =~ s/^-/+/;
	if ($gcd ne '+1') { 
	    $num = &'bdiv($num,$gcd);
	    $dom = &'bdiv($dom,$gcd);
	} else {
	    $num = &'bnorm($num);
	    $dom = &'bnorm($dom);
	}
	substr($dom,0,1) = '';
	"$num/$dom";
    }
}

# negation
sub main'rneg { #(rat_num) return rat_num
    local($_) = &'rnorm(@_);
    tr/-+/+-/ if ($_ ne '+0/1');
    $_;
}

# absolute value
sub main'rabs { #(rat_num) return $rat_num
    local($_) = &'rnorm(@_);
    substr($_,0,1) = '+' unless $_ eq 'NaN';
    $_;
}

# multipication
sub main'rmul { #(rat_num, rat_num) return rat_num
    local($xn,$xd) = split('/',&'rnorm($_[0]));
    local($yn,$yd) = split('/',&'rnorm($_[1]));
    &norm(&'bmul($xn,$yn),&'bmul($xd,$yd));
}

# division
sub main'rdiv { #(rat_num, rat_num) return rat_num
    local($xn,$xd) = split('/',&'rnorm($_[0]));
    local($yn,$yd) = split('/',&'rnorm($_[1]));
    &norm(&'bmul($xn,$yd),&'bmul($xd,$yn));
}

# addition
sub main'radd { #(rat_num, rat_num) return rat_num
    local($xn,$xd) = split('/',&'rnorm($_[0]));
    local($yn,$yd) = split('/',&'rnorm($_[1]));
    &norm(&'badd(&'bmul($xn,$yd),&'bmul($yn,$xd)),&'bmul($xd,$yd));
}

# subtraction
sub main'rsub { #(rat_num, rat_num) return rat_num
    local($xn,$xd) = split('/',&'rnorm($_[0]));
    local($yn,$yd) = split('/',&'rnorm($_[1]));
    &norm(&'bsub(&'bmul($xn,$yd),&'bmul($yn,$xd)),&'bmul($xd,$yd));
}

# comparison
sub main'rcmp { #(rat_num, rat_num) return cond_code
    local($xn,$xd) = split('/',&'rnorm($_[0]));
    local($yn,$yd) = split('/',&'rnorm($_[1]));
    &bigint'cmp(&'bmul($xn,$yd),&'bmul($yn,$xd));
}

# int and frac parts
sub main'rmod { #(rat_num) return (rat_num,rat_num)
    local($xn,$xd) = split('/',&'rnorm(@_));
    local($i,$f) = &'bdiv($xn,$xd);
    if (wantarray) {
	("$i/1", "$f/$xd");
    } else {
	"$i/1";
    }   
}

# square root by Newtons method.
#   cycles specifies the number of iterations default: 5
sub main'rsqrt { #(fnum_str[, cycles]) return fnum_str
    local($x, $scale) = (&'rnorm($_[0]), $_[1]);
    if ($x eq 'NaN') {
	'NaN';
    } elsif ($x =~ /^-/) {
	'NaN';
    } else {
	local($gscale, $guess) = (0, '+1/1');
	$scale = 5 if (!$scale);
	while ($gscale++ < $scale) {
	    $guess = &'rmul(&'radd($guess,&'rdiv($x,$guess)),"+1/2");
	}
	"$guess";          # quotes necessary due to perl bug
    }
}