Add DNRM2 regression test for issues 2998 and 3654
authorMartin Kroeker <martin@ruby.chemie.uni-freiburg.de>
Sun, 3 Jul 2022 15:56:49 +0000 (17:56 +0200)
committerGitHub <noreply@github.com>
Sun, 3 Jul 2022 15:56:49 +0000 (17:56 +0200)
utest/CMakeLists.txt
utest/Makefile
utest/test_dnrm2.c [new file with mode: 0644]
utest/utest_main2.c

index 0c99e0d..5511015 100644 (file)
@@ -13,6 +13,7 @@ else ()
     test_rot.c
     test_axpy.c
     test_dsdot.c
+    test_dnrm2.c
     test_swap.c
   )
 endif ()
index fad3607..f990354 100644 (file)
@@ -11,7 +11,7 @@ UTESTBIN=openblas_utest
 
 include $(TOPDIR)/Makefile.system
 
-OBJS=utest_main.o test_min.o test_amax.o test_ismin.o test_rotmg.o test_axpy.o test_dotu.o test_dsdot.o test_swap.o test_rot.o
+OBJS=utest_main.o test_min.o test_amax.o test_ismin.o test_rotmg.o test_axpy.o test_dotu.o test_dsdot.o test_swap.o test_rot.o test_dnrm2.o
 #test_rot.o test_swap.o test_axpy.o test_dotu.o test_dsdot.o test_fork.o
 
 ifneq ($(NO_LAPACK), 1)
diff --git a/utest/test_dnrm2.c b/utest/test_dnrm2.c
new file mode 100644 (file)
index 0000000..f710431
--- /dev/null
@@ -0,0 +1,60 @@
+/*****************************************************************************
+Copyright (c) 2011-2022, 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 "openblas_utest.h"
+#if defined(BUILD_DOUBLE)
+CTEST(dnrm2,dnrm2_inf)
+{
+       double x[29];
+       blasint incx=1;
+       blasint n=28;
+       double res1=0.0f, res2=INFINITY;
+
+       for (int 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)
+{
+       double x[29];
+       blasint incx=1;
+       blasint n=28;
+       double res1=0.0f, res2=0.0f;
+
+       for (int 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
index 6b25286..e382237 100644 (file)
@@ -606,6 +606,7 @@ int main(int argc, const char ** argv){
   CTEST_ADD (zdotu,zdotu_n_1);
   CTEST_ADD (zdotu,zdotu_offset_1);
   CTEST_ADD (dsdot,dsdot_n_1);
+  CTEST_ADD (dnrm2,dnrm2_tiny);
   CTEST_ADD (rot,drot_inc_0);
   CTEST_ADD (rot,zdrot_inc_0);
   CTEST_ADD (rot,srot_inc_0);