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
|
/*
* c_zblas1.c
*
* The program is a C wrapper for zcblat1.
*
* Written by Keita Teranishi. 2/11/1998
*
*/
#include "cblas_test.h"
#include "cblas.h"
void F77_zaxpy(const CBLAS_INT *N, const void *alpha, void *X,
const CBLAS_INT *incX, void *Y, const CBLAS_INT *incY)
{
cblas_zaxpy(*N, alpha, X, *incX, Y, *incY);
return;
}
void F77_zcopy(const CBLAS_INT *N, void *X, const CBLAS_INT *incX,
void *Y, const CBLAS_INT *incY)
{
cblas_zcopy(*N, X, *incX, Y, *incY);
return;
}
void F77_zdotc(const CBLAS_INT *N, const void *X, const CBLAS_INT *incX,
const void *Y, const CBLAS_INT *incY,void *dotc)
{
cblas_zdotc_sub(*N, X, *incX, Y, *incY, dotc);
return;
}
void F77_zdotu(const CBLAS_INT *N, void *X, const CBLAS_INT *incX,
void *Y, const CBLAS_INT *incY,void *dotu)
{
cblas_zdotu_sub(*N, X, *incX, Y, *incY, dotu);
return;
}
void F77_zdscal(const CBLAS_INT *N, const double *alpha, void *X,
const CBLAS_INT *incX)
{
cblas_zdscal(*N, *alpha, X, *incX);
return;
}
void F77_zscal(const CBLAS_INT *N, const void * *alpha, void *X,
const CBLAS_INT *incX)
{
cblas_zscal(*N, alpha, X, *incX);
return;
}
void F77_zswap( const CBLAS_INT *N, void *X, const CBLAS_INT *incX,
void *Y, const CBLAS_INT *incY)
{
cblas_zswap(*N,X,*incX,Y,*incY);
return;
}
CBLAS_INT F77_izamax(const CBLAS_INT *N, const void *X, const CBLAS_INT *incX)
{
if (*N < 1 || *incX < 1) return(0);
return(cblas_izamax(*N, X, *incX)+1);
}
double F77_dznrm2(const CBLAS_INT *N, const void *X, const CBLAS_INT *incX)
{
return cblas_dznrm2(*N, X, *incX);
}
double F77_dzasum(const CBLAS_INT *N, void *X, const CBLAS_INT *incX)
{
return cblas_dzasum(*N, X, *incX);
}
|
