• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1 /*
2  * c_dblas1.c
3  *
4  * The program is a C wrapper for dcblat1.
5  *
6  * Written by Keita Teranishi.  2/11/1998
7  *
8  */
9 #include "cblas_test.h"
10 #include "cblas.h"
F77_dasum(const int * N,double * X,const int * incX)11 double F77_dasum(const int *N, double *X, const int *incX)
12 {
13    return cblas_dasum(*N, X, *incX);
14 }
15 
F77_daxpy(const int * N,const double * alpha,const double * X,const int * incX,double * Y,const int * incY)16 void F77_daxpy(const int *N, const double *alpha, const double *X,
17                     const int *incX, double *Y, const int *incY)
18 {
19    cblas_daxpy(*N, *alpha, X, *incX, Y, *incY);
20    return;
21 }
22 
F77_dcopy(const int * N,double * X,const int * incX,double * Y,const int * incY)23 void F77_dcopy(const int *N, double *X, const int *incX,
24                     double *Y, const int *incY)
25 {
26    cblas_dcopy(*N, X, *incX, Y, *incY);
27    return;
28 }
29 
F77_ddot(const int * N,const double * X,const int * incX,const double * Y,const int * incY)30 double F77_ddot(const int *N, const double *X, const int *incX,
31                 const double *Y, const int *incY)
32 {
33    return cblas_ddot(*N, X, *incX, Y, *incY);
34 }
35 
F77_dnrm2(const int * N,const double * X,const int * incX)36 double F77_dnrm2(const int *N, const double *X, const int *incX)
37 {
38    return cblas_dnrm2(*N, X, *incX);
39 }
40 
F77_drotg(double * a,double * b,double * c,double * s)41 void F77_drotg( double *a, double *b, double *c, double *s)
42 {
43    cblas_drotg(a,b,c,s);
44    return;
45 }
46 
F77_drot(const int * N,double * X,const int * incX,double * Y,const int * incY,const double * c,const double * s)47 void F77_drot( const int *N, double *X, const int *incX, double *Y,
48        const int *incY, const double *c, const double *s)
49 {
50 
51    cblas_drot(*N,X,*incX,Y,*incY,*c,*s);
52    return;
53 }
54 
F77_dscal(const int * N,const double * alpha,double * X,const int * incX)55 void F77_dscal(const int *N, const double *alpha, double *X,
56                          const int *incX)
57 {
58    cblas_dscal(*N, *alpha, X, *incX);
59    return;
60 }
61 
F77_dswap(const int * N,double * X,const int * incX,double * Y,const int * incY)62 void F77_dswap( const int *N, double *X, const int *incX,
63                           double *Y, const int *incY)
64 {
65    cblas_dswap(*N,X,*incX,Y,*incY);
66    return;
67 }
68 
F77_dzasum(const int * N,void * X,const int * incX)69 double F77_dzasum(const int *N, void *X, const int *incX)
70 {
71    return cblas_dzasum(*N, X, *incX);
72 }
73 
F77_dznrm2(const int * N,const void * X,const int * incX)74 double F77_dznrm2(const int *N, const void *X, const int *incX)
75 {
76    return cblas_dznrm2(*N, X, *incX);
77 }
78 
F77_idamax(const int * N,const double * X,const int * incX)79 int F77_idamax(const int *N, const double *X, const int *incX)
80 {
81    if (*N < 1 || *incX < 1) return(0);
82    return (cblas_idamax(*N, X, *incX)+1);
83 }
84