---
title: kml_sparse_?dotci_sub
description: "x · y = xˆH * y，即x向量中的元素取其共轭复数后与y向量中对应元素相乘然后累加（复数a + bi的共轭为a - bi），计算公式表示为：conjg(x[0]) * y[indx[0]] + ... + conjg(x[nz-1]) * y[indx[nz-1]]。"
url: https://www.hikunpeng.com/document/detail/zh/kunpengaccel/math-lib/devg-kml/kunpengaccel_kml_0022.html
sourcePath: /source/zh/kunpengaccel/math-lib/devg-kml/kunpengaccel_kml_0022.html
indexId: 8f0f540eec78b1264a76e4129ca440c7f8f0271dc7146b1a33d357248e1f45ed68
---
# kml_sparse_?dotci_sub

x · y = xˆH * y，即x向量中的元素取其共轭复数后与y向量中对应元素相乘然后累加（复数a + bi的共轭为a - bi），计算公式表示为：conjg(x[0]) * y[indx[0]] + ... + conjg(x[nz-1]) * y[indx[nz-1]]。

#### 接口定义

C interface：

kml_sparse_status_t kml_sparse_cdotci_sub(const KML_INT nz, const KML_Complex8 *x, const KML_INT *indx, const KML_Complex8 *y, KML_Complex8 *dotci);

kml_sparse_status_t kml_sparse_zdotci_sub(const KML_INT nz, const KML_Complex16 *x, const KML_INT *indx, const KML_Complex16 *y, KML_Complex16 *dotci);

Fortran interface：

RES = KML_SPARSE_CDOTCI_SUB(NZ, X, INDX, Y, DOTCI);

RES = KML_SPARSE_ZDOTCI_SUB(NZ, X, INDX, Y, DOTCI);


#### 参数

| 参数名 | 类型 | 描述 | 输入/输出 |
| --- | --- | --- | --- |
| nz | 整型数 | x及indx数组中元素的个数。 | 输入 |
| x | 在cdotci中，x是单精度浮点类型数组。 在zdotci中，x是双精度浮点类型数组。 | 存储非零元素的数组x，大于或等于nz。 | 输入 |
| indx | 整型数组 | indx[i]表示x数组中第i个元素在稠密向量中的序列号，数组大于或等于nz。 | 输入 |
| y | 在cdotci中，y是单精度浮点类型数组。 在zdotci中，y是双精度浮点类型数组。 | 数组y，大于或等于max(indx[i])。 | 输入 |
| dotci | 在cdotci中，dotci是单精度浮点类型。 在zdotci中，dotci是双精度浮点类型。 | 返回共轭x与y的点积。 | 输入/输出 |


#### 返回值

函数执行状态，枚举类型kml_sparse_status_t。


#### 依赖

C: "kspblas.h"

Fortran: "kspblas.f03"


#### 示例

C interface：

```
KML_INT nz = 2; 
    KML_Complex8 dotci = {0, 0}; 
    KML_INT indx[2] = {1, 2}; 
    KML_Complex8 x[2] = {{1, 2}, {3, 4}}; 
    KML_Complex8 y[4] = {{-1, 1}, {5, 3}, {-2, -3}, {4, 1}}; 
    kml_sparse_status_t status = kml_sparse_cdotci_sub(nz, x, indx, y, &dotci); 
 
    /* 
     *  Output dotci: 
     *     {-7.000000, -8.000000} 
     */
```

Fortran interface：

```
INTEGER(C_INT) :: NZ = 2
INTEGER(C_INT) :: STATUS
INTEGER(C_INT) :: I
TYPE(KML_COMPLEX8), TARGET :: DOTCI = KML_COMPLEX8(0, 0)
INTEGER(C_INT) :: INDX(2)
TYPE(KML_COMPLEX8) :: CX(2), Y(4)
REAL(C_FLOAT) :: X(4), Y(8)
TYPE(C_PTR) :: PDOTCI
PDOTCI = C_LOC(DOTCI)
DATA INDX/1, 2/
DATA X/1, 2, 3, 4/
DATA Y /-1, 1, 5, 3, -2, -3, 4, 1/
I = 1
DO WHILE(I <= 2)
CX(I)%REAL = X(2 * I - 1)
CX(I)%IMAG = X(2 * I)
END DO
I = 1
DO WHILE(I <= 4)
CY(I)%REAL = Y(2 * I - 1)
CY(I)%IMAG = Y(2 * I)
END DO
STATUS = KML_SPARSE_CDOTCI_SUB(NZ, CX, INDX, CY, PDOTCI)
!
!  OUTPUT DOTCI:
!     {-7.000000, -8.000000}
!
```