v0.9.0
cblas_sspmv.c
Go to the documentation of this file.
1 /*
2  *
3  * cblas_sspmv.c
4  * This program is a C interface to sspmv.
5  * Written by Keita Teranishi
6  * 4/6/1998
7  *
8  */
9 #include "cblas.h"
10 #include "cblas_f77.h"
11 void cblas_sspmv(const enum CBLAS_ORDER order,
12  const enum CBLAS_UPLO Uplo, const int N,
13  const float alpha, const float *AP,
14  const float *X, const int incX, const float beta,
15  float *Y, const int incY)
16 {
17  char UL;
18 #ifdef F77_CHAR
19  F77_CHAR F77_UL;
20 #else
21  #define F77_UL &UL
22 #endif
23 #ifdef F77_INT
24  F77_INT F77_N=N, F77_incX=incX, F77_incY=incY;
25 #else
26  #define F77_N N
27  #define F77_incX incX
28  #define F77_incY incY
29 #endif
30  extern int CBLAS_CallFromC;
31  extern int RowMajorStrg;
32  RowMajorStrg = 0;
33 
34  CBLAS_CallFromC = 1;
35  if (order == CblasColMajor)
36  {
37  if (Uplo == CblasUpper) UL = 'U';
38  else if (Uplo == CblasLower) UL = 'L';
39  else
40  {
41  cblas_xerbla(2, "cblas_sspmv","Illegal Uplo setting, %d\n",Uplo );
42  CBLAS_CallFromC = 0;
43  RowMajorStrg = 0;
44  return;
45  }
46  #ifdef F77_CHAR
47  F77_UL = C2F_CHAR(&UL);
48  #endif
49  F77_sspmv(F77_UL, &F77_N, &alpha, AP, X,
50  &F77_incX, &beta, Y, &F77_incY);
51  }
52  else if (order == CblasRowMajor)
53  {
54  RowMajorStrg = 1;
55  if (Uplo == CblasUpper) UL = 'L';
56  else if (Uplo == CblasLower) UL = 'U';
57  else
58  {
59  cblas_xerbla(2, "cblas_sspmv","Illegal Uplo setting, %d\n", Uplo);
60  CBLAS_CallFromC = 0;
61  RowMajorStrg = 0;
62  return;
63  }
64  #ifdef F77_CHAR
65  F77_UL = C2F_CHAR(&UL);
66  #endif
67  F77_sspmv(F77_UL, &F77_N, &alpha,
68  AP, X,&F77_incX, &beta, Y, &F77_incY);
69  }
70  else cblas_xerbla(1, "cblas_sspmv", "Illegal Order setting, %d\n", order);
71  CBLAS_CallFromC = 0;
72  RowMajorStrg = 0;
73 }
#define F77_incX
CBLAS_UPLO
Definition: cblas.h:12
int CBLAS_CallFromC
Definition: cblas_globals.c:1
void cblas_sspmv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const int N, const float alpha, const float *AP, const float *X, const int incX, const float beta, float *Y, const int incY)
Definition: cblas_sspmv.c:11
CBLAS_ORDER
Definition: cblas.h:10
void cblas_xerbla(int p, const char *rout, const char *form,...)
Definition: cblas_xerbla.c:8
#define F77_incY
#define F77_N
const int N
Definition: speed_test.cpp:3
void F77_sspmv(FCHAR, FINT, const float *, const float *, const float *, FINT, const float *, float *, FINT)
#define F77_UL
int RowMajorStrg
Definition: cblas_globals.c:2