#! /usr/bin/python -t
# _*_ coding: iso-8859-1 _*_
# Last edited on 2009-05-02 19:03:06 by stolfilocal

PROG_NAME = "hrn_test"
PROG_DESC = "Tests {hrn.py}"
PROG_VERS = "1.0"

PROG_COPYRIGHT = "Copyright © 2009 State University of Campinas"

PROG_INFO = "!!! PROG_INFO to be written"

import sys;
import os;
sys.path[1:0] = [ sys.path[0] + '/../lib', os.path.expandvars('${STOLFIHOME}/lib'), '.' ] 
import hrn;
import rmxn;
import rn;

def data_error(msg) :
  "Prints the error message {msg} about the current input line, and aborts."
  sys.stderr.write("%s:%d: %s\n" % (in_file_name, in_line_number, msg));
  sys.exit(1)
 
def arg_error(msg):
  "Prints the error message {msg} about the command line arguments, and aborts."
  sys.stderr.write("%s\n" % msg);
  sys.stderr.write("usage: %s\n" % PROG_HELP);
  sys.exit(1)

# Testing {rn}:

s = rn.norm_sqr([1,2,3]);
sys.stderr.write("rn.norm_sqr([1,2,3]) = %s\n" % s);

s = rn.norm([1,2,3]);
sys.stderr.write("rn.norm([1,2,3]) = %s\n" % s);

d = rn.dir([1,2,3]);
sys.stderr.write("rn.dir([1,2,3]) = %s\n" % d);

s = rn.dot([1,2,3],[4,5,6]);
sys.stderr.write("rn.dot([1,2,3],[4,5,6]) = %s\n" % s);

w = rn.cross2([1,2,3],[4,5,6]);
sys.stderr.write("rn.cross2([1,2,3],[4,5,6]) = %s\n" % w);

# Testing {rmxn}:

M = rmxn.zero_matrix(2,3);
sys.stderr.write("rmxn.zero_matrix(2,3) = %s\n" % M);

M = rmxn.ident_matrix(3,4);
sys.stderr.write("rmxn.ident_matrix(3,4) = %s\n" % M);

M = rmxn.diag_matrix([2,3,4]);
sys.stderr.write("rmxn.diag_matrix([2,3,4]) = %s\n" % M);

A = [[1 2 3], [4,5,6]];
sys.stderr.write("A = %s\n" % A);

B = [[1 2], [3, 4], [5, 6]];
sys.stderr.write("B = %s\n" % B);

M = rmxn.mul(A,B);
sys.stderr.write("rmxn.mul(A,B) = %s\n" % M);

v = rmxn.map_row([10,100],A);
sys.stderr.write("rmxn.map_row([10,100],A) = %s\n" % v);

v = rmxn.map_col(B,[10,100]);
sys.stderr.write("rmxn.map_col(B,[10,100]) = %s\n" % v);

# Testing {hrn}:

q = hrn.pt_scale([2, 3, 4, 5], [10,100,1000]);
sys.stderr.write("hrn.pt_scale([2, 3, 4, 5], [10,100,1000]) = %s\n" % q);

M = hrn.trans_matrix([2, 3, 4, 5]);
sys.stderr.write("hrn.trans_matrix([2, 3, 4, 5]) = %s\n" % M);

M = hrn.scale_matrix([3, 4, 5]);
sys.stderr.write("hrn.scale_matrix([3, 4, 5]) = %s\n" % M);

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