#======================================================================

#                    N R U T I L . P L 

#                    doc: Wed Feb 24 17:44:49 1999

#                    dlm: Sun Jul  2 00:47:04 2006

#                    (c) 1999 A.M. Thurnherr

#                    uE-Info: 45 31 NIL 0 0 72 2 2 4 NIL ofnI

#======================================================================

# extract from nrutil.c/nrutil.h (Numerical Recipes) adapted for ANTS

# HISTORY:

#	Feb 24, 1999: -	created from c-source

#	Aug 01, 1999: - added macros from nrutil.h

#	Sep 26, 1999: - added &dumpMatrix()

#   Jul  1, 2006: - Version 3.3 [HISTORY]

# NOTES:

#	- allocation routines &vector, &matrix needed to make sure

#	  right number of elts is allocated (for $# op)

#	- array elts are initialized with nan

#	- array indices MUST start with 1 (in the spirit of FORTRAN IV, bless)

#	- instead of pointer return, we use refs

sub vector($$$)

{

	my($vR,$nl,$nh) = @_;

	my($i);

	croak("vector must be 1-relative")

		unless ($nl == 1);

	for ($i=1; $i<=$nh; $i++) {

		$vR->[$i] = nan;

	}

}

sub matrix($$$$$)

{

	my($mR,$nrl,$nrh,$ncl,$nch) = @_;

	my($i,$j);

	croak("matrix must be 1-relative")

		unless ($nrl == 1 && $ncl == 1);

	for ($i=1; $i<=$nrh; $i++) {

		for ($j=1; $j<=$nch; $j++) {

			$mR->[$i][$j] = nan;

		}

	}

}

#----------------------------------------------------------------------

sub dumpMatrix($$)

{

	my($msg,$mR) = @_;

	my($rows) = $#{$mR};

	my($cols) = $#{$mR->[1]};

	my($r,$c);

	print(STDERR "$msg: $rows x $cols (rows x cols)\n");

	for ($r=1; $r<=$rows; $r++) {

		for ($c=1; $c<=$cols; $c++) {

			printf(STDERR "%.3e\t",$mR->[$r][$c]);

		}

		print(STDERR "\n");

	}

}

#----------------------------------------------------------------------

# Macros

#----------------------------------------------------------------------

sub SQR($) { return $_[0] * $_[0]; }						# D?SQR

sub MAX($$) { return ($_[0] > $_[1]) ? $_[0] : $_[1]; }		# [DF]MAX

sub MIN($$) { return ($_[0] < $_[1]) ? $_[0] : $_[1]; }		# [DF]MIN

sub SIGN($$) { return ($_[1] >= 0) ? $_[0] : -$_[0]; }		# SIGN

1;