floraachy
/
OpenBLAS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Restore the remaining utests

This commit is contained in:
Martin Kroeker 2018-02-15 23:54:59 +01:00 committed by GitHub
parent 6940c59a88
commit df06c59ee4
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
8 changed files with 368 additions and 97 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
utest/CMakeLists.txt
View File

@ -8,9 +8,28 @@ else ()
utest_main.c
test_amax.c
test_rotmg.c
test_rot.c
test_axpy.c
test_dotu.c
test_dsdot.c
test_swap.c
)
endif ()
if (NOT USE_OPENMP)
if (NOT OS_WINDOWS)
set(OpenBLAS_utest_src
${OpenBLAS_utest_src}
test_fork.c
)
elseif (OS_CYGWIN_NT)
set(OpenBLAS_utest_src
${OpenBLAS_utest_src}
test_fork.c
)
endif()
endif()
if (NOT NO_LAPACK)
set(OpenBLAS_utest_src
${OpenBLAS_utest_src}

4
utest/Makefile
View File

@ -8,13 +8,14 @@ UTESTBIN=openblas_utest
include $(TOPDIR)/Makefile.system
OBJS=utest_main.o test_amax.o test_rotmg.o
OBJS=utest_main.o test_amax.o test_rotmg.o test_axpy.o test_dotu.o test_dsdot.o test_swap.o test_rot.o
#test_rot.o test_swap.o test_axpy.o test_dotu.o test_dsdot.o test_fork.o
ifneq ($(NO_LAPACK), 1)
#OBJS += test_potrs.o
endif
ifndef USE_OPENMP
ifndef OS_WINDOWS
OBJS += test_fork.o
else
@ -22,6 +23,7 @@ ifdef OS_CYGWIN_NT
OBJS += test_fork.o
endif
endif
endif
all : run_test

43
utest/test_axpy.c
View File

@ -31,30 +31,29 @@ USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
**********************************************************************************/
#include "common_utest.h"
#include "openblas_utest.h"
void test_daxpy_inc_0(void)
CTEST(axpy,daxpy_inc_0)
{
int i;
int 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[]={2.0,4.0,6.0,8.0,2.0,4.0,6.0,8.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);
//reference
BLASFUNC_REF(daxpy)(&N,&a,x2,&incX,y2,&incY);
for(i=0; i<N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
void test_zaxpy_inc_0(void)
CTEST(axpy,zaxpy_inc_0)
{
int i;
int N=4,incX=0,incY=0;
@ -62,20 +61,18 @@ void test_zaxpy_inc_0(void)
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};
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);
//reference
BLASFUNC_REF(zaxpy)(&N,a,x2,&incX,y2,&incY);
for(i=0; i<2*N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
void test_saxpy_inc_0(void)
CTEST(axpy,saxpy_inc_0)
{
int i;
int N=8,incX=0,incY=0;
@ -83,20 +80,18 @@ void test_saxpy_inc_0(void)
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};
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);
//reference
BLASFUNC_REF(saxpy)(&N,&a,x2,&incX,y2,&incY);
for(i=0; i<N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
void test_caxpy_inc_0(void)
CTEST(axpy,caxpy_inc_0)
{
int i;
int N=4,incX=0,incY=0;
@ -104,15 +99,13 @@ void test_caxpy_inc_0(void)
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};
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);
//reference
BLASFUNC_REF(caxpy)(&N,a,x2,&incX,y2,&incY);
for(i=0; i<2*N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}

22
utest/test_dotu.c
View File

@ -31,10 +31,10 @@ USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
**********************************************************************************/
#include "common_utest.h"
#include "openblas_utest.h"
#include <complex.h>
void test_zdotu_n_1(void)
CTEST( zdotu,zdotu_n_1)
{
int N=1,incX=1,incY=1;
double x1[]={1.0,1.0};
@ -42,19 +42,17 @@ void test_zdotu_n_1(void)
double x2[]={1.0,1.0};
double y2[]={1.0,2.0};
double _Complex result1=0.0;
double _Complex result2=0.0;
double _Complex result2={-1.0000+3.0000*I};
//OpenBLAS
result1=BLASFUNC(zdotu)(&N,x1,&incX,y1,&incY);
//reference
result2=BLASFUNC_REF(zdotu)(&N,x2,&incX,y2,&incY);
CU_ASSERT_DOUBLE_EQUAL(creal(result1), creal(result2), CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(cimag(result1), cimag(result2), CHECK_EPS);
ASSERT_DBL_NEAR_TOL(creal(result1), creal(result2), DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(cimag(result1), cimag(result2), DOUBLE_EPS);
// printf("\%lf,%lf\n",creal(result1),cimag(result1));
}
void test_zdotu_offset_1(void)
CTEST(zdotu, zdotu_offset_1)
{
int N=1,incX=1,incY=1;
double x1[]={1.0,2.0,3.0,4.0};
@ -62,14 +60,12 @@ void test_zdotu_offset_1(void)
double x2[]={1.0,2.0,3.0,4.0};
double y2[]={5.0,6.0,7.0,8.0};
double _Complex result1=0.0;
double _Complex result2=0.0;
double _Complex result2={-9.0000+32.0000*I};
//OpenBLAS
result1=BLASFUNC(zdotu)(&N,x1+1,&incX,y1+1,&incY);
//reference
result2=BLASFUNC_REF(zdotu)(&N,x2+1,&incX,y2+1,&incY);
CU_ASSERT_DOUBLE_EQUAL(creal(result1), creal(result2), CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(cimag(result1), cimag(result2), CHECK_EPS);
ASSERT_DBL_NEAR_TOL(creal(result1), creal(result2), DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(cimag(result1), cimag(result2), DOUBLE_EPS);
// printf("\%lf,%lf\n",creal(result1),cimag(result1));
}

10
utest/test_dsdot.c
View File

@ -31,9 +31,9 @@ USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
**********************************************************************************/
#include "common_utest.h"
#include "openblas_utest.h"
void test_dsdot_n_1()
CTEST(dsdot,dsdot_n_1)
{
float x= 0.172555164;
float y= -0.0138700781;
@ -41,11 +41,9 @@ void test_dsdot_n_1()
int incy=1;
int n=1;
double res1=0.0f, res2=0.0f;
double res1=0.0f, res2=-0.00239335360107;
res1=BLASFUNC(dsdot)(&n, &x, &incx, &y, &incy);
res2=BLASFUNC_REF(dsdot)(&n, &x, &incx, &y, &incy);
CU_ASSERT_DOUBLE_EQUAL(res1, res2, CHECK_EPS);
ASSERT_DBL_NEAR_TOL(res1, res2, DOUBLE_EPS);
}

54
utest/test_rot.c
View File

@ -31,88 +31,82 @@ USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
**********************************************************************************/
#include "common_utest.h"
#include "openblas_utest.h"
void test_drot_inc_0(void)
CTEST(rot,drot_inc_0)
{
int i=0;
int 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[]={1.0,3.0,5.0,7.0};
double y2[]={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);
//reference
BLASFUNC_REF(drot)(&N,x2,&incX,y2,&incY,&c,&s);
for(i=0; i<N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
void test_zdrot_inc_0(void)
CTEST(rot,zdrot_inc_0)
{
int i=0;
int 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[]={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};
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);
//reference
BLASFUNC_REF(zdrot)(&N,x2,&incX,y2,&incY,&c,&s);
for(i=0; i<2*N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
void test_srot_inc_0(void)
CTEST(rot,srot_inc_0)
{
int i=0;
int 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[]={1.0,3.0,5.0,7.0};
float y2[]={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);
//reference
BLASFUNC_REF(srot)(&N,x2,&incX,y2,&incY,&c,&s);
for(i=0; i<N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], SINGLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], SINGLE_EPS);
}
}
void test_csrot_inc_0(void)
CTEST(rot, csrot_inc_0)
{
int i=0;
int 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[]={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};
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);
//reference
BLASFUNC_REF(csrot)(&N,x2,&incX,y2,&incY,&c,&s);
for(i=0; i<2*N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], SINGLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], SINGLE_EPS);
}
}

36
utest/test_swap.c
View File

@ -31,9 +31,9 @@ USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
**********************************************************************************/
#include "common_utest.h"
#include "openblas_utest.h"
void test_dswap_inc_0(void)
CTEST(swap,dswap_inc_0)
{
int i=0;
int N=4,incX=0,incY=0;
@ -41,19 +41,17 @@ void test_dswap_inc_0(void)
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);
//reference
BLASFUNC_REF(dswap)(&N,x2,&incX,y2,&incY);
for(i=0; i<N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
void test_zswap_inc_0(void)
CTEST(swap,zswap_inc_0)
{
int i=0;
int N=4,incX=0,incY=0;
@ -64,16 +62,14 @@ void test_zswap_inc_0(void)
//OpenBLAS
BLASFUNC(zswap)(&N,x1,&incX,y1,&incY);
//reference
BLASFUNC_REF(zswap)(&N,x2,&incX,y2,&incY);
for(i=0; i<2*N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
void test_sswap_inc_0(void)
CTEST(swap,sswap_inc_0)
{
int i=0;
int N=4,incX=0,incY=0;
@ -84,16 +80,14 @@ void test_sswap_inc_0(void)
//OpenBLAS
BLASFUNC(sswap)(&N,x1,&incX,y1,&incY);
//reference
BLASFUNC_REF(sswap)(&N,x2,&incX,y2,&incY);
for(i=0; i<N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], SINGLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], SINGLE_EPS);
}
}
void test_cswap_inc_0(void)
CTEST(swap,cswap_inc_0)
{
int i=0;
int N=4,incX=0,incY=0;
@ -104,11 +98,9 @@ void test_cswap_inc_0(void)
//OpenBLAS
BLASFUNC(cswap)(&N,x1,&incX,y1,&incY);
//reference
BLASFUNC_REF(cswap)(&N,x2,&incX,y2,&incY);
for(i=0; i<2*N; i++){
CU_ASSERT_DOUBLE_EQUAL(x1[i], x2[i], CHECK_EPS);
CU_ASSERT_DOUBLE_EQUAL(y1[i], y2[i], CHECK_EPS);
ASSERT_DBL_NEAR_TOL(x1[i], x2[i], SINGLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], SINGLE_EPS);
}
}

277
utest/utest_main2.c
View File

@ -32,6 +32,7 @@ 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
@ -183,6 +184,282 @@ CTEST(drotmg, rotmg_D1eqD2_X1eqX2){
}
}
CTEST(axpy,daxpy_inc_0)
{
int i;
int 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(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
CTEST(axpy,zaxpy_inc_0)
{
int i;
int 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(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
CTEST(axpy,saxpy_inc_0)
{
int i;
int 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(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
CTEST(axpy,caxpy_inc_0)
{
int i;
int 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(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
CTEST( zdotu,zdotu_n_1)
{
int N=1,incX=1,incY=1;
double x1[]={1.0,1.0};
double y1[]={1.0,2.0};
double x2[]={1.0,1.0};
double y2[]={1.0,2.0};
double _Complex result1=0.0;
double _Complex result2={-1.0000+3.0000*I};
//OpenBLAS
result1=BLASFUNC(zdotu)(&N,x1,&incX,y1,&incY);
ASSERT_DBL_NEAR_TOL(creal(result1), creal(result2), DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(cimag(result1), cimag(result2), DOUBLE_EPS);
}
CTEST(zdotu, zdotu_offset_1)
{
int 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};
double x2[]={1.0,2.0,3.0,4.0};
double y2[]={5.0,6.0,7.0,8.0};
double _Complex result1=0.0;
double _Complex result2={-9.0000+32.0000*I};
//OpenBLAS
result1=BLASFUNC(zdotu)(&N,x1+1,&incX,y1+1,&incY);
ASSERT_DBL_NEAR_TOL(creal(result1), creal(result2), DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(cimag(result1), cimag(result2), DOUBLE_EPS);
}
CTEST(dsdot,dsdot_n_1)
{
float x= 0.172555164;
float y= -0.0138700781;
int incx=1;
int incy=1;
int n=1;
double res1=0.0f, res2=-0.00239335360107;
res1=BLASFUNC(dsdot)(&n, &x, &incx, &y, &incy);
ASSERT_DBL_NEAR_TOL(res1, res2, DOUBLE_EPS);
}
CTEST(rot,drot_inc_0)
{
int i=0;
int 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(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
CTEST(rot,zdrot_inc_0)
{
int i=0;
int 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(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
CTEST(rot,srot_inc_0)
{
int i=0;
int 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(x1[i], x2[i], SINGLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], SINGLE_EPS);
}
}
CTEST(rot, csrot_inc_0)
{
int i=0;
int 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(x1[i], x2[i], SINGLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], SINGLE_EPS);
}
}
CTEST(swap,dswap_inc_0)
{
int i=0;
int 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(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
CTEST(swap,zswap_inc_0)
{
int i=0;
int 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(x1[i], x2[i], DOUBLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], DOUBLE_EPS);
}
}
CTEST(swap,sswap_inc_0)
{
int i=0;
int 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(x1[i], x2[i], SINGLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], SINGLE_EPS);
}
}
CTEST(swap,cswap_inc_0)
{
int i=0;
int 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(x1[i], x2[i], SINGLE_EPS);
ASSERT_DBL_NEAR_TOL(y1[i], y2[i], SINGLE_EPS);
}
}
int main(int argc, const char ** argv){
CTEST_ADD(amax, samax);