Fix undefined reference to cblas_[sd]cabs1

cblas.h defines

   double cblas_dcabs1(const void  *z);
   float  cblas_scabs1(const void  *c);

but does not provide an implementation. This commit adds
the missing implementation in the common CBLAS style and
eliminates the linker error when calling the above function.

Author: angsch

CBLAS/include/cblas_f77.h

@@ -97,6 +97,8 @@
 #define F77_csrot_base		F77_GLOBAL(csrot, CSROT)
 #define F77_zrotg_base		F77_GLOBAL(zrotg, ZROTG)
 #define F77_zdrot_base		F77_GLOBAL(zdrot, ZDROT)
+#define F77_scabs1_sub_base F77_GLOBAL(scabs1sub, SCABS1SUB)
+#define F77_dcabs1_sub_base F77_GLOBAL(dcabs1sub, DCABS1SUB)
 
 /*
  * Level 2 BLAS
@@ -221,6 +223,7 @@
    #define F77_snrm2_sub(...) F77_snrm2_sub_base(__VA_ARGS__)
    #define F77_sasum_sub(...) F77_sasum_sub_base(__VA_ARGS__)
    #define F77_isamax_sub(...) F77_isamax_sub_base(__VA_ARGS__)
+   #define F77_scabs1_sub(...) F77_scabs1_sub_base(__VA_ARGS__)
 
 /* Double Precision */
 
@@ -238,6 +241,7 @@
    #define F77_dnrm2_sub(...) F77_dnrm2_sub_base(__VA_ARGS__)
    #define F77_dasum_sub(...) F77_dasum_sub_base(__VA_ARGS__)
    #define F77_idamax_sub(...) F77_idamax_sub_base(__VA_ARGS__)
+   #define F77_dcabs1_sub(...) F77_dcabs1_sub_base(__VA_ARGS__)
 
 /* Single Complex Precision */
 
@@ -602,6 +606,7 @@ void F77_xerbla_base(FCHAR, void *
    void F77_csscal_base( FINT, const float *, void *, FINT);
    void F77_scnrm2_sub_base( FINT, const void *, FINT, float *);
    void F77_scasum_sub_base( FINT, const void *, FINT, float *);
+   void F77_scabs1_sub_base( const void *, float *);
 
 /* Double Complex Precision */
 
@@ -618,6 +623,7 @@ void F77_xerbla_base(FCHAR, void *
    void F77_dznrm2_sub_base( FINT, const void *, FINT, double *);
    void F77_dzasum_sub_base( FINT, const void *, FINT, double *);
    void F77_izamax_sub_base( FINT, const void *, FINT, FINT2);
+   void F77_dcabs1_sub_base( const void *, double *);
 
 /*
  * Level 2 Fortran Prototypes

CBLAS/src/CMakeLists.txt

@@ -29,14 +29,16 @@ set(DLEV1 cblas_drotg.c cblas_drotmg.c cblas_drot.c cblas_drotm.c
 set(CLEV1 cblas_crotg.c cblas_csrot.c
          cblas_cswap.c cblas_cscal.c cblas_csscal.c cblas_ccopy.c
          cblas_caxpy.c cblas_cdotu_sub.c cblas_cdotc_sub.c
-         cblas_icamax.c cdotcsub.f cdotusub.f icamaxsub.f)
+         cblas_icamax.c cdotcsub.f cdotusub.f icamaxsub.f
+         cblas_scabs1.c scabs1sub.f )
 
 # Files for level 1 double precision complex
 set(ZLEV1 cblas_zrotg.c cblas_zdrot.c
          cblas_zswap.c cblas_zscal.c cblas_zdscal.c cblas_zcopy.c
          cblas_zaxpy.c cblas_zdotu_sub.c cblas_zdotc_sub.c cblas_dznrm2.c
          cblas_dzasum.c cblas_izamax.c zdotcsub.f zdotusub.f
-         dzasumsub.f dznrm2sub.f izamaxsub.f)
+         dzasumsub.f dznrm2sub.f izamaxsub.f
+         cblas_dcabs1.c dcabs1sub.f)
 
 # Common files for level 1 single precision
 set(SCLEV1 cblas_scasum.c scasumsub.f cblas_scnrm2.c scnrm2sub.f)

CBLAS/src/Makefile

@@ -39,14 +39,16 @@ dlev1 = cblas_drotg.o cblas_drotmg.o cblas_drot.o  cblas_drotm.o \
 clev1 = cblas_crotg.o cblas_csrot.o \
         cblas_cswap.o cblas_cscal.o cblas_csscal.o cblas_ccopy.o \
         cblas_caxpy.o cblas_cdotu_sub.o cblas_cdotc_sub.o \
-        cblas_icamax.o cdotcsub.o cdotusub.o icamaxsub.o
+        cblas_icamax.o cdotcsub.o cdotusub.o icamaxsub.o \
+        cblas_scabs1.o scabs1sub.o
 
 # Files for level 1 double precision complex
 zlev1 = cblas_zrotg.o cblas_zdrot.o \
         cblas_zswap.o cblas_zscal.o cblas_zdscal.o cblas_zcopy.o \
         cblas_zaxpy.o cblas_zdotu_sub.o cblas_zdotc_sub.o cblas_dznrm2.o \
         cblas_dzasum.o cblas_izamax.o zdotcsub.o zdotusub.o \
-        dzasumsub.o dznrm2sub.o izamaxsub.o
+        dzasumsub.o dznrm2sub.o izamaxsub.o \
+        cblas_dcabs1.o dcabs1sub.o
 
 # Common files for level 1 single precision
 sclev1 = cblas_scasum.o scasumsub.o cblas_scnrm2.o scnrm2sub.o

CBLAS/src/cblas_dcabs1.c

@@ -0,0 +1,15 @@
+/*
+ * cblas_scabs1.c
+ *
+ * The program is a C interface to scabs1.
+ *
+ */
+#include "cblas.h"
+#include "cblas_f77.h"
+double cblas_dcabs1(const void *c)
+{
+   double cabs1 = 0.0;
+   F77_dcabs1_sub(c, &cabs1);
+   return cabs1; 
+}
+

CBLAS/src/cblas_scabs1.c

@@ -0,0 +1,15 @@
+/*
+ * cblas_scabs1.c
+ *
+ * The program is a C interface to scabs1.
+ *
+ */
+#include "cblas.h"
+#include "cblas_f77.h"
+float cblas_scabs1(const void *c)
+{
+   float cabs1 = 0.0;
+   F77_scabs1_sub(c, &cabs1);
+   return cabs1; 
+}
+

CBLAS/src/dcabs1sub.f

@@ -0,0 +1,13 @@
+c     dcabs1.f
+c
+c     The program is a fortran wrapper for dcabs1.
+c
+      subroutine dcabs1sub(z, cabs1)
+c
+      external dcabs1
+      double complex z
+      double precision dcabs1, cabs1
+c
+      cabs1=dcabs1(z)
+      return
+      end

CBLAS/src/scabs1sub.f

@@ -0,0 +1,13 @@
+c     scabs1.f
+c
+c     The program is a fortran wrapper for scabs1.
+c
+      subroutine scabs1sub(z, cabs1)
+c
+      external scabs1
+      complex z
+      real scabs1, cabs1
+c
+      cabs1=scabs1(z)
+      return
+      end