707 lines
18 KiB
C
707 lines
18 KiB
C
/*****************************************************************************
|
|
Copyright (c) 2011-2016, The OpenBLAS Project
|
|
All rights reserved.
|
|
|
|
Redistribution and use in source and binary forms, with or without
|
|
modification, are permitted provided that the following conditions are
|
|
met:
|
|
|
|
1. Redistributions of source code must retain the above copyright
|
|
notice, this list of conditions and the following disclaimer.
|
|
|
|
2. Redistributions in binary form must reproduce the above copyright
|
|
notice, this list of conditions and the following disclaimer in
|
|
the documentation and/or other materials provided with the
|
|
distribution.
|
|
3. Neither the name of the OpenBLAS project nor the names of
|
|
its contributors may be used to endorse or promote products
|
|
derived from this software without specific prior written
|
|
permission.
|
|
|
|
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
|
|
AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
|
|
IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
|
|
ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
|
|
LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
|
|
DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
|
|
SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
|
|
CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
|
|
OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE
|
|
USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|
|
|
**********************************************************************************/
|
|
|
|
#include <stdio.h>
|
|
#include <complex.h>
|
|
|
|
#define CTEST_MAIN
|
|
#define CTEST_SEGFAULT
|
|
#define CTEST_ADD_TESTS_MANUALLY
|
|
|
|
#include "cblas.h"
|
|
#include "openblas_utest.h"
|
|
|
|
CTEST(amax, samax){
|
|
blasint N=3, inc=1;
|
|
float te_max=0.0, tr_max=0.0;
|
|
float x[]={-1.1, 2.2, -3.3};
|
|
te_max=BLASFUNC(samax)(&N, x, &inc);
|
|
tr_max=3.3;
|
|
|
|
ASSERT_DBL_NEAR_TOL((double)(tr_max), (double)(te_max), SINGLE_EPS);
|
|
}
|
|
|
|
CTEST(amax, damax){
|
|
blasint N=3, inc=1;
|
|
double te_max=0.0, tr_max=0.0;
|
|
double x[]={-1.1, 2.2, -3.3};
|
|
|
|
te_max=BLASFUNC(damax)(&N, x, &inc);
|
|
tr_max=3.3;
|
|
|
|
ASSERT_DBL_NEAR_TOL((double)(tr_max), (double)(te_max), DOUBLE_EPS);
|
|
}
|
|
|
|
CTEST (drotmg,rotmg)
|
|
{
|
|
double te_d1, tr_d1;
|
|
double te_d2, tr_d2;
|
|
double te_x1, tr_x1;
|
|
double te_y1, tr_y1;
|
|
double te_param[5];
|
|
double tr_param[5];
|
|
int i=0;
|
|
// original test case for libGoto bug fixed by feb2014 rewrite
|
|
te_d1= 0.21149573940783739;
|
|
te_d2= 0.046892057172954082;
|
|
te_x1= -0.42272687517106533;
|
|
te_y1= 0.42211309121921659;
|
|
|
|
|
|
for(i=0; i<5; i++){
|
|
te_param[i]=tr_param[i]=0.0;
|
|
}
|
|
|
|
//reference values as calculated by netlib blas
|
|
|
|
tr_d1= 0.1732048;
|
|
tr_d2= 0.03840234;
|
|
tr_x1= -0.516180;
|
|
tr_y1= 0.422113;
|
|
tr_d1= 0.17320483687975;
|
|
tr_d2= 0.03840233915037;
|
|
tr_x1= -0.51618034832329;
|
|
tr_y1= 0.42211309121922;
|
|
|
|
tr_param[0]= 0.0;
|
|
tr_param[1]= 0.0;
|
|
tr_param[2]= 0.99854803659786;
|
|
tr_param[3]= -0.22139439665872;
|
|
tr_param[4]= 0.0;
|
|
|
|
BLASFUNC(drotmg)(&te_d1, &te_d2, &te_x1, &te_y1, te_param);
|
|
ASSERT_DBL_NEAR_TOL(tr_d1, te_d1, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_d2, te_d2, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_x1, te_x1, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_y1, te_y1, DOUBLE_EPS);
|
|
|
|
for(i=0; i<5; i++){
|
|
ASSERT_DBL_NEAR_TOL(tr_param[i], te_param[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST (drotmg,rotmg_issue1452)
|
|
{
|
|
double te_d1, tr_d1;
|
|
double te_d2, tr_d2;
|
|
double te_x1, tr_x1;
|
|
double te_y1, tr_y1;
|
|
double te_param[5];
|
|
double tr_param[5];
|
|
int i=0;
|
|
|
|
// from issue #1452
|
|
te_d1 = 5.9e-8;
|
|
te_d2 = 5.960464e-8;
|
|
te_x1 = 1.0;
|
|
te_y1 = 150.0;
|
|
|
|
for(i=0; i<5; i++){
|
|
te_param[i]=tr_param[i]=0.0;
|
|
}
|
|
te_param[3]=1./4096.;
|
|
//reference values as calculated by gonum blas with rotmg rewritten to Hopkins' algorithm
|
|
tr_d1= 0.99995592822897;
|
|
tr_d2= 0.98981219860583;
|
|
tr_x1= 0.03662270484346;
|
|
tr_y1= 150.000000000000;
|
|
|
|
tr_param[0]= -1.0;
|
|
tr_param[1]= 0.00000161109346;
|
|
tr_param[2]= -0.00024414062500;
|
|
tr_param[3]= 0.00024414062500;
|
|
tr_param[4]= 0.00000162760417;
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(drotmg)(&te_d1, &te_d2, &te_x1, &te_y1, te_param);
|
|
|
|
ASSERT_DBL_NEAR_TOL(tr_d1, te_d1, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_d2, te_d2, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_x1, te_x1, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_y1, te_y1, DOUBLE_EPS);
|
|
|
|
for(i=0; i<5; i++){
|
|
ASSERT_DBL_NEAR_TOL(tr_param[i], te_param[i], DOUBLE_EPS);
|
|
}
|
|
|
|
}
|
|
|
|
CTEST(drotmg, rotmg_D1eqD2_X1eqX2)
|
|
{
|
|
double te_d1, tr_d1;
|
|
double te_d2, tr_d2;
|
|
double te_x1, tr_x1;
|
|
double te_y1, tr_y1;
|
|
double te_param[5];
|
|
double tr_param[5];
|
|
int i=0;
|
|
te_d1= tr_d1=2.;
|
|
te_d2= tr_d2=2.;
|
|
te_x1= tr_x1=8.;
|
|
te_y1= tr_y1=8.;
|
|
|
|
for(i=0; i<5; i++){
|
|
te_param[i]=tr_param[i]=0.0;
|
|
}
|
|
|
|
//reference values as calculated by netlib blas
|
|
tr_d1= 1.0;
|
|
tr_d2= 1.0;
|
|
tr_x1= 16.0;
|
|
tr_y1= 8.0;
|
|
|
|
tr_param[0]=1.0;
|
|
tr_param[1]=1.0;
|
|
tr_param[2]=0.0;
|
|
tr_param[3]=0.0;
|
|
tr_param[4]=1.0;
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(drotmg)(&te_d1, &te_d2, &te_x1, &te_y1, te_param);
|
|
|
|
ASSERT_DBL_NEAR_TOL(tr_d1, te_d1, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_d2, te_d2, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_x1, te_x1, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_y1, te_y1, DOUBLE_EPS);
|
|
|
|
for(i=0; i<5; i++){
|
|
ASSERT_DBL_NEAR_TOL(tr_param[i], te_param[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(drotmg, drotmg_D1_big_D2_big_flag_zero)
|
|
{
|
|
double te_d1, tr_d1;
|
|
double te_d2, tr_d2;
|
|
double te_x1, tr_x1;
|
|
double te_y1, tr_y1;
|
|
double te_param[5]={1.,4096.,-4096.,1.,4096.};
|
|
double tr_param[5]={-1.,4096.,-3584.,1792.,4096.};
|
|
int i=0;
|
|
te_d1= tr_d1=1600000000.;
|
|
te_d2= tr_d2=800000000.;
|
|
te_x1= tr_x1=8.;
|
|
te_y1= tr_y1=7.;
|
|
|
|
|
|
//reference values as calculated by gonum
|
|
tr_d1= 68.96627824858757;
|
|
tr_d2= 34.483139124293785;
|
|
tr_x1= 45312.;
|
|
tr_y1= 7.0;
|
|
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(drotmg)(&te_d1, &te_d2, &te_x1, &te_y1, te_param);
|
|
|
|
ASSERT_DBL_NEAR_TOL(tr_d1, te_d1, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_d2, te_d2, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_x1, te_x1, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(tr_y1, te_y1, DOUBLE_EPS);
|
|
|
|
for(i=0; i<5; i++){
|
|
ASSERT_DBL_NEAR_TOL(tr_param[i], te_param[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(axpy,daxpy_inc_0)
|
|
{
|
|
blasint i;
|
|
blasint N=8,incX=0,incY=0;
|
|
double a=0.25;
|
|
double x1[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
double y1[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
double x2[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
double y2[]={4.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(daxpy)(&N,&a,x1,&incX,y1,&incY);
|
|
|
|
for(i=0; i<N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(axpy,zaxpy_inc_0)
|
|
{
|
|
blasint i;
|
|
blasint N=4,incX=0,incY=0;
|
|
double a[2]={0.25,0.5};
|
|
double x1[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
double y1[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
double x2[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
double y2[]={-3.0,9.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(zaxpy)(&N,a,x1,&incX,y1,&incY);
|
|
|
|
for(i=0; i<2*N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(axpy,saxpy_inc_0)
|
|
{
|
|
blasint i;
|
|
blasint N=8,incX=0,incY=0;
|
|
float a=0.25;
|
|
float x1[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
float y1[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
float x2[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
float y2[]={4.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(saxpy)(&N,&a,x1,&incX,y1,&incY);
|
|
|
|
for(i=0; i<N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(axpy,caxpy_inc_0)
|
|
{
|
|
blasint i;
|
|
blasint N=4,incX=0,incY=0;
|
|
float a[2]={0.25,0.5};
|
|
float x1[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
float y1[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
float x2[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
float y2[]={-3.0,9.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(caxpy)(&N,a,x1,&incX,y1,&incY);
|
|
|
|
for(i=0; i<2*N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST( zdotu,zdotu_n_1)
|
|
{
|
|
blasint N=1,incX=1,incY=1;
|
|
double x1[]={1.0,1.0};
|
|
double y1[]={1.0,2.0};
|
|
openblas_complex_double result1=openblas_make_complex_double(0.0,0.0);
|
|
openblas_complex_double result2=openblas_make_complex_double(-1.0,3.0);
|
|
#ifdef RETURN_BY_STACK
|
|
BLASFUNC(zdotu)(&result1,&N,x1,&incX,y1,&incY);
|
|
#else
|
|
result1=BLASFUNC(zdotu)(&N,x1,&incX,y1,&incY);
|
|
#endif
|
|
|
|
#ifdef OPENBLAS_COMPLEX_STRUCT
|
|
ASSERT_DBL_NEAR_TOL(result2.real, result1.real, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(result2.imag, result1.imag, DOUBLE_EPS);
|
|
#else
|
|
ASSERT_DBL_NEAR_TOL(creal(result2), creal(result1), DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(cimag(result2), cimag(result1), DOUBLE_EPS);
|
|
#endif
|
|
}
|
|
|
|
CTEST(zdotu, zdotu_offset_1)
|
|
{
|
|
blasint N=1,incX=1,incY=1;
|
|
double x1[]={1.0,2.0,3.0,4.0};
|
|
double y1[]={5.0,6.0,7.0,8.0};
|
|
openblas_complex_double result1=openblas_make_complex_double(0.0,0.0);
|
|
openblas_complex_double result2=openblas_make_complex_double(-9.0,32.0);
|
|
#ifdef RETURN_BY_STACK
|
|
BLASFUNC(zdotu)(&result1,&N,x1+1,&incX,y1+1,&incY);
|
|
#else
|
|
result1=BLASFUNC(zdotu)(&N,x1+1,&incX,y1+1,&incY);
|
|
#endif
|
|
|
|
#ifdef OPENBLAS_COMPLEX_STRUCT
|
|
ASSERT_DBL_NEAR_TOL(result2.real, result1.real, DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(result2.imag, result1.imag, DOUBLE_EPS);
|
|
#else
|
|
ASSERT_DBL_NEAR_TOL(creal(result2), creal(result1), DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(cimag(result2), cimag(result1), DOUBLE_EPS);
|
|
#endif
|
|
}
|
|
|
|
CTEST(dsdot,dsdot_n_1)
|
|
{
|
|
float x= 0.172555164F;
|
|
float y= -0.0138700781F;
|
|
blasint incx=1;
|
|
blasint incy=1;
|
|
blasint n=1;
|
|
|
|
double res1=0.0f, res2=-0.00239335360107;
|
|
|
|
res1=BLASFUNC(dsdot)(&n, &x, &incx, &y, &incy);
|
|
ASSERT_DBL_NEAR_TOL(res2, res1, DOUBLE_EPS);
|
|
|
|
}
|
|
|
|
#if defined(BUILD_DOUBLE)
|
|
CTEST(dnrm2,dnrm2_inf)
|
|
{
|
|
#ifndef INFINITY
|
|
#define INFINITY HUGE_VAL
|
|
#endif
|
|
int i;
|
|
double x[29];
|
|
blasint incx=1;
|
|
blasint n=28;
|
|
double res1=0.0f, res2=INFINITY;
|
|
|
|
for (i=0;i<n;i++)x[i]=0.0f;
|
|
x[10]=-INFINITY;
|
|
res1=BLASFUNC(dnrm2)(&n, x, &incx);
|
|
ASSERT_DBL_NEAR_TOL(res2, res1, DOUBLE_EPS);
|
|
|
|
}
|
|
CTEST(dnrm2,dnrm2_tiny)
|
|
{
|
|
int i;
|
|
double x[29];
|
|
blasint incx=1;
|
|
blasint n=28;
|
|
double res1=0.0f, res2=0.0f;
|
|
|
|
for (i=0;i<n;i++)x[i]=7.457008414e-310;
|
|
res1=BLASFUNC(dnrm2)(&n, x, &incx);
|
|
ASSERT_DBL_NEAR_TOL(res2, res1, DOUBLE_EPS);
|
|
}
|
|
#endif
|
|
|
|
CTEST(rot,drot_inc_0)
|
|
{
|
|
blasint i=0;
|
|
blasint N=4,incX=0,incY=0;
|
|
double c=0.25,s=0.5;
|
|
double x1[]={1.0,3.0,5.0,7.0};
|
|
double y1[]={2.0,4.0,6.0,8.0};
|
|
double x2[]={-0.21484375000000,3.0,5.0,7.0};
|
|
double y2[]={ 0.03906250000000,4.0,6.0,8.0};
|
|
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(drot)(&N,x1,&incX,y1,&incY,&c,&s);
|
|
|
|
for(i=0; i<N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(rot,zdrot_inc_0)
|
|
{
|
|
blasint i=0;
|
|
blasint N=4,incX=0,incY=0;
|
|
double c=0.25,s=0.5;
|
|
double x1[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
double y1[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
double x2[]={-0.21484375000000,-0.45703125000000 ,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
double y2[]={ 0.03906250000000, 0.17187500000000 ,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(zdrot)(&N,x1,&incX,y1,&incY,&c,&s);
|
|
|
|
for(i=0; i<2*N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(rot,srot_inc_0)
|
|
{
|
|
blasint i=0;
|
|
blasint N=4,incX=0,incY=0;
|
|
float c=0.25,s=0.5;
|
|
float x1[]={1.0,3.0,5.0,7.0};
|
|
float y1[]={2.0,4.0,6.0,8.0};
|
|
float x2[]={-0.21484375000000,3.0,5.0,7.0};
|
|
float y2[]={ 0.03906250000000,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(srot)(&N,x1,&incX,y1,&incY,&c,&s);
|
|
|
|
for(i=0; i<N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], SINGLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], SINGLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(rot, csrot_inc_0)
|
|
{
|
|
blasint i=0;
|
|
blasint N=4,incX=0,incY=0;
|
|
float c=0.25,s=0.5;
|
|
float x1[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
float y1[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
float x2[]={-0.21484375000000,-0.45703125000000 ,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
float y2[]={ 0.03906250000000, 0.17187500000000 ,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(csrot)(&N,x1,&incX,y1,&incY,&c,&s);
|
|
|
|
for(i=0; i<2*N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], SINGLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], SINGLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(swap,dswap_inc_0)
|
|
{
|
|
blasint i=0;
|
|
blasint N=4,incX=0,incY=0;
|
|
double x1[]={1.0,3.0,5.0,7.0};
|
|
double y1[]={2.0,4.0,6.0,8.0};
|
|
double x2[]={1.0,3.0,5.0,7.0};
|
|
double y2[]={2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(dswap)(&N,x1,&incX,y1,&incY);
|
|
|
|
for(i=0; i<N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(swap,zswap_inc_0)
|
|
{
|
|
blasint i=0;
|
|
blasint N=4,incX=0,incY=0;
|
|
double x1[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
double y1[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
double x2[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
double y2[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(zswap)(&N,x1,&incX,y1,&incY);
|
|
|
|
for(i=0; i<2*N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], DOUBLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], DOUBLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(swap,sswap_inc_0)
|
|
{
|
|
blasint i=0;
|
|
blasint N=4,incX=0,incY=0;
|
|
float x1[]={1.0,3.0,5.0,7.0};
|
|
float y1[]={2.0,4.0,6.0,8.0};
|
|
float x2[]={1.0,3.0,5.0,7.0};
|
|
float y2[]={2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(sswap)(&N,x1,&incX,y1,&incY);
|
|
|
|
for(i=0; i<N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], SINGLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], SINGLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(swap,cswap_inc_0)
|
|
{
|
|
blasint i=0;
|
|
blasint N=4,incX=0,incY=0;
|
|
float x1[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
float y1[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
float x2[]={1.0,3.0,5.0,7.0,1.0,3.0,5.0,7.0};
|
|
float y2[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.0};
|
|
|
|
//OpenBLAS
|
|
BLASFUNC(cswap)(&N,x1,&incX,y1,&incY);
|
|
|
|
for(i=0; i<2*N; i++){
|
|
ASSERT_DBL_NEAR_TOL(x2[i], x1[i], SINGLE_EPS);
|
|
ASSERT_DBL_NEAR_TOL(y2[i], y1[i], SINGLE_EPS);
|
|
}
|
|
}
|
|
|
|
CTEST(min, smin_negative){
|
|
blasint N=3, inc=1;
|
|
float te_min=0.0, tr_min=0.0;
|
|
float x[]={-1.1, -2.2, -3.3};
|
|
|
|
te_min=BLASFUNC(smin)(&N, x, &inc);
|
|
tr_min=-3.3;
|
|
|
|
ASSERT_DBL_NEAR_TOL((double)(tr_min), (double)(te_min), SINGLE_EPS);
|
|
}
|
|
|
|
CTEST(min, dmin_positive){
|
|
blasint N=3, inc=1;
|
|
double te_min=0.0, tr_min=0.0;
|
|
double x[]={1.1, 0.0, 3.3};
|
|
|
|
te_min=BLASFUNC(dmin)(&N, x, &inc);
|
|
tr_min=0.0;
|
|
|
|
ASSERT_DBL_NEAR_TOL((double)(tr_min), (double)(te_min), DOUBLE_EPS);
|
|
}
|
|
|
|
CTEST(min, smin_zero){
|
|
blasint N=3, inc=1;
|
|
float te_min=0.0, tr_min=0.0;
|
|
float x[]={1.1, 2.2, 0.0};
|
|
|
|
te_min=BLASFUNC(smin)(&N, x, &inc);
|
|
tr_min=0.0;
|
|
|
|
ASSERT_DBL_NEAR_TOL((double)(tr_min), (double)(te_min), SINGLE_EPS);
|
|
}
|
|
|
|
CTEST(max, smax_negative){
|
|
blasint N=3, inc=1;
|
|
float te_max=0.0, tr_max=0.0;
|
|
float x[]={-1.1, -2.2, -3.3};
|
|
|
|
te_max=BLASFUNC(smax)(&N, x, &inc);
|
|
tr_max=-1.1;
|
|
|
|
ASSERT_DBL_NEAR_TOL((double)(tr_max), (double)(te_max), SINGLE_EPS);
|
|
}
|
|
|
|
CTEST(max, dmax_positive){
|
|
blasint N=3, inc=1;
|
|
double te_max=0.0, tr_max=0.0;
|
|
double x[]={1.1, 0.0, 3.3};
|
|
|
|
te_max=BLASFUNC(dmax)(&N, x, &inc);
|
|
tr_max=3.3;
|
|
|
|
ASSERT_DBL_NEAR_TOL((double)(tr_max), (double)(te_max), DOUBLE_EPS);
|
|
}
|
|
|
|
CTEST(max, smax_zero){
|
|
blasint N=3, inc=1;
|
|
float te_max=0.0, tr_max=0.0;
|
|
float x[]={-1.1, -2.2, 0.0};
|
|
|
|
te_max=BLASFUNC(smax)(&N, x, &inc);
|
|
tr_max=0.0;
|
|
|
|
ASSERT_DBL_NEAR_TOL((double)(tr_max), (double)(te_max), SINGLE_EPS);
|
|
}
|
|
|
|
|
|
CTEST(zscal, i_nan)
|
|
{
|
|
double i[] = {0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1 };
|
|
double nan[] = {NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0};
|
|
cblas_zscal(9, i, &nan, 1);
|
|
ASSERT_TRUE(isnan(nan[0]));
|
|
ASSERT_TRUE(isnan(nan[1]));
|
|
ASSERT_TRUE(isnan(nan[16]));
|
|
ASSERT_TRUE(isnan(nan[17]));
|
|
}
|
|
|
|
CTEST(zscal, nan_i)
|
|
{
|
|
double i[] = {0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1 };
|
|
double nan[] = {NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0, NAN,0};
|
|
cblas_zscal(9, &nan, &i, 1);
|
|
ASSERT_TRUE(isnan(i[0]));
|
|
ASSERT_TRUE(isnan(i[1]));
|
|
ASSERT_TRUE(isnan(i[16]));
|
|
ASSERT_TRUE(isnan(i[17]));
|
|
}
|
|
|
|
CTEST(zscal, i_inf)
|
|
{
|
|
double i[] = {0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1 };
|
|
double inf[] = {INFINITY, 0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0};
|
|
cblas_zscal(9, i, &inf, 1);
|
|
ASSERT_TRUE(isnan(inf[0]));
|
|
ASSERT_TRUE(isinf(inf[1]));
|
|
ASSERT_TRUE(isnan(inf[16]));
|
|
ASSERT_TRUE(isinf(inf[17]));
|
|
}
|
|
|
|
CTEST(zscal, inf_i)
|
|
{
|
|
double i[] = {0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1, 0,1 };
|
|
double inf[] = {INFINITY, 0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0, INFINITY,0};
|
|
cblas_zscal(9, &inf, &i, 1);
|
|
ASSERT_TRUE(isnan(i[0]));
|
|
ASSERT_TRUE(isinf(i[1]));
|
|
ASSERT_TRUE(isnan(i[16]));
|
|
ASSERT_TRUE(isinf(i[17]));
|
|
}
|
|
|
|
int main(int argc, const char ** argv){
|
|
|
|
CTEST_ADD (amax, samax);
|
|
CTEST_ADD (amax, damax);
|
|
CTEST_ADD (min, smin_negative);
|
|
CTEST_ADD (min, dmin_positive);
|
|
CTEST_ADD (min, smin_zero);
|
|
CTEST_ADD (max, smax_negative);
|
|
CTEST_ADD (max, dmax_positive);
|
|
CTEST_ADD (max, smax_zero);
|
|
CTEST_ADD (drotmg,rotmg);
|
|
CTEST_ADD (drotmg,rotmg_issue1452);
|
|
CTEST_ADD (drotmg,rotmg_D1eqD2_X1eqX2);
|
|
CTEST_ADD (drotmg,drotmg_D1_big_D2_big_flag_zero);
|
|
CTEST_ADD (axpy,daxpy_inc_0);
|
|
CTEST_ADD (axpy,zaxpy_inc_0);
|
|
CTEST_ADD (axpy,saxpy_inc_0);
|
|
CTEST_ADD (axpy,caxpy_inc_0);
|
|
CTEST_ADD (zdotu,zdotu_n_1);
|
|
CTEST_ADD (zdotu,zdotu_offset_1);
|
|
CTEST_ADD (dsdot,dsdot_n_1);
|
|
CTEST_ADD (dnrm2,dnrm2_inf);
|
|
CTEST_ADD (dnrm2,dnrm2_tiny);
|
|
CTEST_ADD (rot,drot_inc_0);
|
|
CTEST_ADD (rot,zdrot_inc_0);
|
|
CTEST_ADD (rot,srot_inc_0);
|
|
CTEST_ADD (rot,csrot_inc_0);
|
|
CTEST_ADD (swap,dswap_inc_0);
|
|
CTEST_ADD (swap,zswap_inc_0);
|
|
CTEST_ADD (swap,sswap_inc_0);
|
|
CTEST_ADD (swap,cswap_inc_0);
|
|
CTEST_ADD (zscal, i_nan);
|
|
CTEST_ADD (zscal, nan_i);
|
|
CTEST_ADD (zscal, i_inf);
|
|
CTEST_ADD (zscal, inf_i);
|
|
int num_fail=0;
|
|
|
|
num_fail=ctest_main(argc, argv);
|
|
|
|
return num_fail;
|
|
}
|