-
-
Notifications
You must be signed in to change notification settings - Fork 805
/
Copy pathzcopy_cblas.c
56 lines (52 loc) · 2.22 KB
/
zcopy_cblas.c
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
/**
* @license Apache-2.0
*
* Copyright (c) 2024 The Stdlib Authors.
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
#include "stdlib/blas/base/zcopy.h"
#include "stdlib/blas/base/zcopy_cblas.h"
#include "stdlib/blas/base/shared.h"
#include "stdlib/complex/float64/ctor.h"
#include "stdlib/strided/base/min_view_buffer_index.h"
/**
* Copies values from one complex double-precision floating-point vector to another complex double-precision floating-point vector.
*
* @param N number of indexed elements
* @param X input array
* @param strideX X stride length
* @param Y output array
* @param strideY Y stride length
*/
void API_SUFFIX(c_zcopy)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, void *Y, const CBLAS_INT strideY ) {
cblas_zcopy( N, X, strideX, Y, strideY );
}
/**
* Copies values from one complex double-precision floating-point vector to another complex double-precision floating-point vector using alternative indexing semantics.
*
* @param N number of indexed elements
* @param X input array
* @param strideX X stride length
* @param offsetX starting index for X
* @param Y output array
* @param strideY Y stride length
* @param offsetY starting index for Y
*/
void API_SUFFIX(c_zcopy_ndarray)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
stdlib_complex128_t *zx = (stdlib_complex128_t *)X;
stdlib_complex128_t *zy = (stdlib_complex128_t *)Y;
zx += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer
zy += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); // adjust array pointer
cblas_zcopy( N, (void *)zx, strideX, (void *)zy, strideY );
}