1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
|
/*
* c_sblas1.c
*
* The program is a C wrapper for scblat1.
*
* Written by Keita Teranishi. 2/11/1998
*
*/
#include "cblas_test.h"
#include "cblas.h"
float F77_sasum(const CBLAS_INT *N, float *X, const CBLAS_INT *incX)
{
return cblas_sasum(*N, X, *incX);
}
void F77_saxpy(const CBLAS_INT *N, const float *alpha, const float *X,
const CBLAS_INT *incX, float *Y, const CBLAS_INT *incY)
{
cblas_saxpy(*N, *alpha, X, *incX, Y, *incY);
return;
}
float F77_scasum(const CBLAS_INT *N, void *X, const CBLAS_INT *incX)
{
return cblas_scasum(*N, X, *incX);
}
float F77_scnrm2(const CBLAS_INT *N, const void *X, const CBLAS_INT *incX)
{
return cblas_scnrm2(*N, X, *incX);
}
void F77_scopy(const CBLAS_INT *N, const float *X, const CBLAS_INT *incX,
float *Y, const CBLAS_INT *incY)
{
cblas_scopy(*N, X, *incX, Y, *incY);
return;
}
float F77_sdot(const CBLAS_INT *N, const float *X, const CBLAS_INT *incX,
const float *Y, const CBLAS_INT *incY)
{
return cblas_sdot(*N, X, *incX, Y, *incY);
}
float F77_snrm2(const CBLAS_INT *N, const float *X, const CBLAS_INT *incX)
{
return cblas_snrm2(*N, X, *incX);
}
void F77_srotg( float *a, float *b, float *c, float *s)
{
cblas_srotg(a,b,c,s);
return;
}
void F77_srot( const CBLAS_INT *N, float *X, const CBLAS_INT *incX, float *Y,
const CBLAS_INT *incY, const float *c, const float *s)
{
cblas_srot(*N,X,*incX,Y,*incY,*c,*s);
return;
}
void F77_sscal(const CBLAS_INT *N, const float *alpha, float *X,
const CBLAS_INT *incX)
{
cblas_sscal(*N, *alpha, X, *incX);
return;
}
void F77_sswap( const CBLAS_INT *N, float *X, const CBLAS_INT *incX,
float *Y, const CBLAS_INT *incY)
{
cblas_sswap(*N,X,*incX,Y,*incY);
return;
}
CBLAS_INT F77_isamax(const CBLAS_INT *N, const float *X, const CBLAS_INT *incX)
{
if (*N < 1 || *incX < 1) return(0);
return (cblas_isamax(*N, X, *incX)+1);
}
|
