Browse Source

Added zdotu with x & y offset=1 test case.

tags/v0.1alpha1
Xianyi Zhang 14 years ago
parent
commit
44acb7503e
3 changed files with 22 additions and 0 deletions
  1. +1
    -0
      utest/common_utest.h
  2. +2
    -0
      utest/main.c
  3. +19
    -0
      utest/test_dotu.c

+ 1
- 0
utest/common_utest.h View File

@@ -55,5 +55,6 @@ void test_saxpy_inc_0(void);
void test_caxpy_inc_0(void);

void test_zdotu_n_1(void);
void test_zdotu_offset_1(void);

#endif

+ 2
- 0
utest/main.c View File

@@ -53,6 +53,8 @@ CU_TestInfo test_level1[]={
{"Testing zaxpy with incx || incy == 0",test_zaxpy_inc_0},

{"Testing zdotu with n == 1",test_zdotu_n_1},
{"Testing zdotu with input x & y offset == 1",test_zdotu_offset_1},
CU_TEST_INFO_NULL,
};



+ 19
- 0
utest/test_dotu.c View File

@@ -53,4 +53,23 @@ void test_zdotu_n_1(void)
}

void test_zdotu_offset_1(void)
{
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=0.0;
//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);
// printf("\%lf,%lf\n",creal(result1),cimag(result1));
}


Loading…
Cancel
Save