psblas3/test/computational_routines/geaxpby

Introduction

This is a directory developed by Luca Pepè Sciarria and Simone Staccone froma Tor Vergata University to start to create some unit tests for PSBLAS 3.9, in particular for psb_geaxpby routine.

Getting started

Steps to reproduce the tests:

• Compile the code using make (Optional)
• Launch the script ./autotest.sh or with source ./autotest.sh if you want to add modules to the .bashrc file permenently.
• Check the output log file psb_geaxpby_test.log to collect results

NOTE: If the code is changed and a new compilation is needed to show the changes, the autotest.sh script isn't aware of this scenario, therefore it is necessary to manually recompile the code.

Test Suite

Overall Analysys

The psb_geaxpby. The signature of the function is:

call psb_geaxpby(alpha, x, beta, y, desc_a, info)

The strategy to validate the correctness of the computation is to compare single precision result and double precision result in the test cases in which the test should not give an error. In this way it is possible to have a correctness check of the computation comparing the two results considering a number of significant digits which is tuned on the single precision computation.

Parameters Values

x vectors are located in the vectors/ directory. They are generated randomly using the same seed and then saved on different files based on their characteristics. The size of the vector is choosen accordingly to the size of the matrix column space considered for the single test instance.

Vector	File Name	Coefficients	Coefficients Description
x_1	x1.txt	x_i> 0, \forall i	Positive coefficients
x_2	x2.txt	x_i < 0, \forall i	Negative coefficients
x_3	x3.txt	x_i \ne 0, \forall i	Random coefficients
x_4	x4.txt	x_i = 0, \forall i	Null coefficients

y vectors are located in the vectors/ directory. They are generated randomly using the same seed and then saved on different files based on their characteristics. The size of the vector is choosen accordingly to the size of the matrix rows space considered for the single test instance.

Vector	File Name	Coefficients	Coefficients Description
y_1	y1.txt	y_i> 0, \forall i	Positive coefficients
y_2	y2.txt	y_i < 0, \forall i	Negative coefficients
y_3	y3.txt	y_i \ne 0, \forall i	Random coefficients
y_4	y4.txt	y_i = 0, \forall i	Null coefficients

\alpha real coefficient multiplied by vector x

\alpha	Value	Coefficients Description
\alpha_1	1.0	Positive value
\alpha_2	-1.0	Negative value
\alpha_3	0.0	Null value

\beta real coefficient multiplied by vector y

\alpha	Value	Coefficients Description
\beta_1	1.0	Positive value
\beta_2	-1.0	Negative value
\beta_3	0.0	Null value

TODO

• Add computation with broken descriptor and catch the errore result
• Test using complex data
• Try multiple distributions
• Try using a matrix instead of a vector