01: /*
02:  *
03:  * cblas_dspr.c
04:  * This program is a C interface to dspr.
05:  * Written by Keita Teranishi
06:  * 4/6/1998
07:  *
08:  */
09: 
10: #include "cblas.h"
11: #include "cblas_f77.h"
12: void cblas_dspr(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo,
13:                 const int N, const double alpha, const double *X,
14:                 const int incX, double *Ap)
15: {
16:    char UL;
17: #ifdef F77_CHAR
18:    F77_CHAR F77_UL;
19: #else
20:    #define F77_UL &UL
21: #endif
22: 
23: #ifdef F77_INT
24:    F77_INT F77_N=N, F77_incX=incX;
25: #else
26:    #define F77_N N
27:    #define F77_incX incX
28: #endif
29:    extern int CBLAS_CallFromC;
30:    extern int RowMajorStrg;
31:    RowMajorStrg = 0;
32:    CBLAS_CallFromC = 1;
33:    if (order == CblasColMajor)
34:    {
35:       if (Uplo == CblasLower) UL = 'L';
36:       else if (Uplo == CblasUpper) UL = 'U';
37:       else
38:       {
39:          cblas_xerbla(2, "cblas_dspr","Illegal Uplo setting, %d\n",Uplo );
40:          CBLAS_CallFromC = 0;
41:          RowMajorStrg = 0;
42:          return;
43:       }
44:       #ifdef F77_CHAR
45:          F77_UL = C2F_CHAR(&UL);
46:       #endif
47: 
48:       F77_dspr(F77_UL, &F77_N, &alpha, X, &F77_incX, Ap);
49: 
50:    }  else if (order == CblasRowMajor)
51:    {
52:       RowMajorStrg = 1;
53:       if (Uplo == CblasLower) UL = 'U';
54:       else if (Uplo == CblasUpper) UL = 'L';
55:       else
56:       {
57:          cblas_xerbla(2, "cblas_dspr","Illegal Uplo setting, %d\n",Uplo );
58:          CBLAS_CallFromC = 0;
59:          RowMajorStrg = 0;
60:          return;
61:       }
62:       #ifdef F77_CHAR
63:          F77_UL = C2F_CHAR(&UL);
64:       #endif
65:       F77_dspr(F77_UL, &F77_N, &alpha, X, &F77_incX, Ap);
66:    } else cblas_xerbla(1, "cblas_dspr", "Illegal Order setting, %d\n", order);
67:    CBLAS_CallFromC = 0;
68:    RowMajorStrg = 0;
69:    return;
70: }
71: