#!/bin/sh
#
# multest
#
# Run multiply and square timing tests, to compute a chart for the
# current processor and compiler combination.

# ***** BEGIN LICENSE BLOCK *****
# Version: MPL 1.1/GPL 2.0/LGPL 2.1
#
# The contents of this file are subject to the Mozilla Public License Version
# 1.1 (the "License"); you may not use this file except in compliance with
# the License. You may obtain a copy of the License at
# http://www.mozilla.org/MPL/
#
# Software distributed under the License is distributed on an "AS IS" basis,
# WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
# for the specific language governing rights and limitations under the
# License.
#
# The Original Code is the MPI Arbitrary Precision Integer Arithmetic library.
#
# The Initial Developer of the Original Code is
# Michael J. Fromberger <sting@linguist.dartmouth.edu>.
# Portions created by the Initial Developer are Copyright (C) 2000
# the Initial Developer. All Rights Reserved.
#
# Contributor(s):
#
# Alternatively, the contents of this file may be used under the terms of
# either the GNU General Public License Version 2 or later (the "GPL"), or
# the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
# in which case the provisions of the GPL or the LGPL are applicable instead
# of those above. If you wish to allow use of your version of this file only
# under the terms of either the GPL or the LGPL, and not to allow others to
# use your version of this file under the terms of the MPL, indicate your
# decision by deleting the provisions above and replace them with the notice
# and other provisions required by the GPL or the LGPL. If you do not delete
# the provisions above, a recipient may use your version of this file under
# the terms of any one of the MPL, the GPL or the LGPL.
#
# ***** END LICENSE BLOCK *****

# $Id: multest,v 1.3 2005/02/02 22:28:22 gerv%gerv.net Exp $
#

ECHO=/bin/echo
MAKE=gmake

$ECHO "\n** Running multiply and square timing tests\n"

$ECHO "Bringing 'mulsqr' up to date ... "
if $MAKE mulsqr ; then
    :
else
    $ECHO "\nMake failed to build mulsqr.\n"
    exit 1
fi

if [ ! -x ./mulsqr ] ; then
    $ECHO "\nCannot find 'mulsqr' program, testing cannot continue.\n"
    exit 1
fi

sizes='64 128 192 256 320 384 448 512 640 768 896 1024 1536 2048'
ntests=500000

$ECHO "Running timing tests, please wait ... "

trap 'echo "oop!";rm -f tt*.tmp;exit 0' INT HUP

touch tt$$.tmp
$ECHO $ntests tests >> tt$$.tmp
for size in $sizes ; do
    $ECHO "$size bits ... \c"
    set -A res `./mulsqr $ntests $size|head -3|tr -d '%'|awk '{print $2}'`
    $ECHO $size"\t"${res[0]}"\t"${res[1]}"\t"${res[2]} >> tt$$.tmp
    $ECHO "(done)"
done
mv tt$$.tmp mulsqr-results.txt
rm -f tt$$.tmp

$ECHO "\n** Running Karatsuba-Ofman multiplication tests\n"

$ECHO "Brining 'karatsuba' up to date ... "
if $MAKE karatsuba ; then
    :
else
    $ECHO "\nMake failed to build karatsuba.\n"
    exit 1
fi

if [ ! -x ./karatsuba ] ; then
    $ECHO "\nCannot find 'karatsuba' program, testing cannot continue.\n"
    exit 1
fi

ntests=100000

trap 'echo "oop!";rm -f tt*.tmp;exit 0' INT HUP

touch tt$$.tmp
for size in $sizes ; do
    $ECHO "$size bits ... "
    ./karatsuba $ntests $size >> tt$$.tmp
    tail -2 tt$$.tmp
done
mv tt$$.tmp karatsuba-results.txt
rm -f tt$$.tmp

exit 0
