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
|
// Copyright (c) 2017-2023, University of Tennessee. All rights reserved.
// SPDX-License-Identifier: BSD-3-Clause
// This program is free software: you can redistribute it and/or modify it under
// the terms of the BSD 3-Clause license. See the accompanying LICENSE file.
#include <stdio.h>
#if defined(BLAS_HAVE_MKL)
#if defined(BLAS_ILP64) && ! defined(MKL_ILP64)
#define MKL_ILP64
#endif
#include <mkl_cblas.h>
#elif defined(BLAS_HAVE_ESSL)
#if defined(BLAS_ILP64) && ! defined(_ESV6464)
#define _ESV6464
#endif
#include <essl.h>
#elif defined(BLAS_HAVE_ACCELERATE)
// On macOS, the official way to include cblas is via Accelerate.h.
// Unfortunately with Xcode 10.3 and GNU g++ 9.3, that doesn't compile.
// If we can find cblas.h, use it, otherwise use Accelerate.h.
#ifdef BLAS_HAVE_ACCELERATE_CBLAS_H
#include <cblas.h>
#else
#include <Accelerate/Accelerate.h>
#endif
#else
#ifdef __cplusplus
// Some ancient cblas.h don't include extern C. It's okay to nest.
extern "C" {
#include <cblas.h>
}
#else
#include <cblas.h>
#endif
#endif
//------------------------------------------------------------------------------
int main()
{
int n = 5;
double x[] = { 1, 2, 3, 4, 5 };
double y[] = { 5, 4, 3, 2, 1 };
for (int i = 0; i < n; ++i) {
printf( "x[ %d ] = %.1f; y[ %d ] = %.1f\n",
i, x[ i ],
i, y[ i ] );
}
double result = cblas_ddot( n, x, 1, y, 1 );
printf( "result = %.1f; should be 35.0\n", result );
bool okay = (result == 35);
printf( "%s\n", okay ? "ok" : "failed" );
return ! okay;
}
|
