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 "common.h"
10 #include "cblas_test.h"
11
F77_dasum(const int * N,double * X,const int * incX)12 double F77_dasum(const int *N, double *X, const int *incX)
13 {
14 return cblas_dasum(*N, X, *incX);
15 }
16
F77_daxpy(const int * N,const double * alpha,const double * X,const int * incX,double * Y,const int * incY)17 void F77_daxpy(const int *N, const double *alpha, const double *X,
18 const int *incX, double *Y, const int *incY)
19 {
20 cblas_daxpy(*N, *alpha, X, *incX, Y, *incY);
21 return;
22 }
23
F77_dcopy(const int * N,double * X,const int * incX,double * Y,const int * incY)24 void F77_dcopy(const int *N, double *X, const int *incX,
25 double *Y, const int *incY)
26 {
27 cblas_dcopy(*N, X, *incX, Y, *incY);
28 return;
29 }
30
F77_ddot(const int * N,const double * X,const int * incX,const double * Y,const int * incY)31 double F77_ddot(const int *N, const double *X, const int *incX,
32 const double *Y, const int *incY)
33 {
34 return cblas_ddot(*N, X, *incX, Y, *incY);
35 }
36
F77_dnrm2(const int * N,const double * X,const int * incX)37 double F77_dnrm2(const int *N, const double *X, const int *incX)
38 {
39 return cblas_dnrm2(*N, X, *incX);
40 }
41
F77_drotg(double * a,double * b,double * c,double * s)42 void F77_drotg( double *a, double *b, double *c, double *s)
43 {
44 cblas_drotg(a,b,c,s);
45 return;
46 }
47
F77_drot(const int * N,double * X,const int * incX,double * Y,const int * incY,const double * c,const double * s)48 void F77_drot( const int *N, double *X, const int *incX, double *Y,
49 const int *incY, const double *c, const double *s)
50 {
51
52 cblas_drot(*N,X,*incX,Y,*incY,*c,*s);
53 return;
54 }
55
F77_dscal(const int * N,const double * alpha,double * X,const int * incX)56 void F77_dscal(const int *N, const double *alpha, double *X,
57 const int *incX)
58 {
59 cblas_dscal(*N, *alpha, X, *incX);
60 return;
61 }
62
F77_dswap(const int * N,double * X,const int * incX,double * Y,const int * incY)63 void F77_dswap( const int *N, double *X, const int *incX,
64 double *Y, const int *incY)
65 {
66 cblas_dswap(*N,X,*incX,Y,*incY);
67 return;
68 }
69
F77_dzasum(const int * N,void * X,const int * incX)70 double F77_dzasum(const int *N, void *X, const int *incX)
71 {
72 return cblas_dzasum(*N, X, *incX);
73 }
74
F77_dznrm2(const int * N,const void * X,const int * incX)75 double F77_dznrm2(const int *N, const void *X, const int *incX)
76 {
77 return cblas_dznrm2(*N, X, *incX);
78 }
79
F77_idamax(const int * N,const double * X,const int * incX)80 int F77_idamax(const int *N, const double *X, const int *incX)
81 {
82 if (*N < 1 || *incX < 1) return(0);
83 return (cblas_idamax(*N, X, *incX)+1);
84 }
85