Laboratoire de l'Informatique et du Parallélisme » MEPlib

#! /usr/bin/perl

# printTCPointsList.pl - ST - 2005-07-18

# COMMAND LINE SYNTAX:

# printPointsList.pl lowerbBound upperBound numPoints

# NOTE(S):

# 

##############################################################################

# PROGRAM DECLARATIONS

##############################################################################

BEGIN

{

  # Modify @INC here if necessary

  push(@INC,"../utils");

  push(@INC,"../utils/tchebychev");

  $ENV{PATH} .= ":../utils/tchebychev";

};

# core modules

use strict; # must be first to force strict in all others

use Carp;

use Getopt::Long;

use Math::BigFloat;

use Math::BigRat;

# custom modules

# Uncomment if you want to use standard exit error codes.

# the stdExitErrors.pm file must be available around.

use stdExitErrors;

use Tchebychev qw(:BASIC);

#use Matrix qw(:BASIC);

##############################################################################

# PROGRAM CONSTANT DECLARATIONS

##############################################################################

my @test = (1,2,10);

my $maxErrContFrac = Math::BigFloat->new("0.001");

##############################################################################

# PROGRAM VARIABLE DECLARATIONS

##############################################################################

my $lowerBoundAsBigRat;

my $upperBoundAsBigRat;

my $numPoints;

my $refPoints;

##############################################################################

# TEST ENVIRONMENT AND CONFIGURATION VARIABLES

##############################################################################

##############################################################################

# INITIALIZE DEBUG FILE

##############################################################################

##############################################################################

# INITIALIZE LOG FILE

##############################################################################

##############################################################################

# GET COMMAND LINE ARGUMENTS

##############################################################################

##############################################################################

# PROGRAM MAIN

##############################################################################

&usage();

# Use $numPoints -2 because we add the lower and the upper bounds.

$refPoints = &tchebychevPointsForInterval($lowerBoundAsBigRat,

                                          $upperBoundAsBigRat,

                                          $numPoints - 2,

                                          $maxErrContFrac);

# Add the lower bound (at the tail) and the upper bound (at the front)

# of the list

unshift(@$refPoints, $upperBoundAsBigRat);

push(@$refPoints, $lowerBoundAsBigRat);

for (my $i = (@$refPoints - 1) ; $i >= 0 ; $i--)

{

  print STDOUT "$$refPoints[$i] ";

}

print STDOUT "\n";

# Uncomment this line and remove the next one if using stdExitErrors.

exit($EX_OK);

##############################################################################

# PROGRAM SUBS

##############################################################################

##############################################################################

# sub mySub

##############################################################################

# Performed task:

#

# input   :

# output  :

# globals :

#

sub mySub {

  my ($someParameter) = @_;

  if (! defined($someParameter))

  {

    my @call_info = caller(0);

    print "\n\n", $call_info[3], " missing parameter. Aborting program!\n\n";

    # Uncomment this line and remove the next one if using stdExitErrors.

    #exit($EX_SOFTWARE);

    exit(1);

  }

} # End mySub

##############################################################################

# sub $retv = _privateSubName($argv)

##############################################################################

sub _privateSubName

{

	my $argv = shift();

	my $retv = undef;

	return($retv);

};

#

##############################################################################

# sub usage

##############################################################################

# Performed task: checks the command line parameters

#

# input   : none

# output  : none

# globals : $ARGV,others...

#

sub usage {

  if (! defined($ARGV[2]))

  {

    my $scriptName = `basename $0`;

    chomp $scriptName;

    print STDERR "\n\nUsage : $scriptName lowerBound upperBound num_points\n\n";

    print STDERR "  - lowerBound: the lower bound of the interval;\n";

    print STDERR "  - upperBound: the upper bound of the interval;\n";

    print STDERR "  - num_points: the number of points.\n\n";

    # Uncomment this line and remove the next one if using stdExitErrors.

    exit($EX_USAGE);

  } # End if

  # Some initializations

  $ARGV[0] =~ s/"//g;

  $lowerBoundAsBigRat = Math::BigRat->new($ARGV[0]);

  if ($lowerBoundAsBigRat->is_nan())

  {

    print STDERR "\n\n";

    print STDERR "The lower bound \"$ARGV[0]\" could not be converted into a ";

    print STDERR "Math::BigRat. Aborting the program!\n\n";

    exit ($EX_DATAERR);

  }

  $ARGV[1] =~ s/"//g;

  $upperBoundAsBigRat = Math::BigRat->new($ARGV[1]);

  if ($upperBoundAsBigRat->is_nan())

  {

    print STDERR "\n\n";

    print STDERR "The upper bound \"$ARGV[1]\" could not be converted into a ";

    print STDERR "Math::BigRat. Aborting the program!\n\n";

    exit ($EX_DATAERR);

  }

  $ARGV[2] =~ s/"//g;

  $numPoints = $ARGV[2];

  if (Math::BigInt->new($numPoints)->is_nan())

  {

    print STDERR "\n\n";

    print STDERR "The number of points \"$ARGV[3]\" is not a valid integer. ";

    print STDERR "Aborting the program!\n\n";

    exit ($EX_DATAERR);

  }

} # End usage

__END__