v?sqrt
计算输入的平方根,输入为向量。
接口定义
C interface:
void vssqrt(const int len, const float* src, float* dst);
void vdsqrt(const int len, const double* src, double* dst);
void vcsqrt(const int len, const float complex *src, float complex *dst);
void vzsqrt(const int len, const double complex *src, double complex *dst);
Fortran interface:
CALL VSSQRT(LEN, SRC, DST);
CALL VDSQRT(LEN, SRC, DST);
参数
参数名 |
类型 |
描述 |
输入/输出 |
---|---|---|---|
len |
整型数 |
表示输入向量的元素个数。 len≤0时会提示len无效并返回。 |
输入 |
src |
|
输入向量src,向量长度为len。 若为空指针,会提示空指针错误并返回。 |
输入 |
dst |
|
输出dst,向量长度为len。 若为空指针,会提示空指针错误并返回。 |
输出 |
输出结果
- 每一个运算值返回数的平方根,取值范围为[0, +INF)。
- 对于vssqrt和vdsqrt,有以下输出情况。
输入值(src)
输出值(dst)
±0
±0
+inf
+inf
< 0
nan
nan
nan
- 对于vcsqrt和vzsqrt,使用欧拉公式和三角函数二倍角公式计算:
特殊取值如下表所示。
RE(z)
i * IM(z)
-∞
-X
-0
+0
+X
+∞
NAN
+i * ∞
+∞+i*∞
+∞+i*∞
+∞+i*∞
+∞+i*∞
+∞+i*∞
+∞+i*∞
+∞+i*∞
+i * Y
+0+i*∞
-
-
-
-
+∞+i*0
QNAN+i*QNAN
+i * 0
+0+i*∞
-
+0+i*0
+0+i*0
-
+∞+i*0
QNAN+i*QNAN
-i * 0
+0-i*∞
-
+0+i*0
+0+i*0
-
+∞-i*0
QNAN+i*QNAN
-i * Y
+0-i*∞
-
-
-
-
+∞-i*0
QNAN+i*QNAN
-i * ∞
+∞-i*∞
+∞-i*∞
+∞-i*∞
+∞-i*∞
+∞-i*∞
+∞-i*∞
+∞-i*∞
+i * NAN
QNAN+i*∞
QNAN+i*QNAN
QNAN+i*QNAN
QNAN+i*QNAN
QNAN+i*QNAN
+∞+i*QNAN
QNAN+i*QNAN
依赖
C: "kvml.h"
示例
C interface:
int i, len = 4; float src[len] = {0.0f, 0.25f, 16.0f, 1.0f}; float* dst = (float*)malloc(sizeof(float) * len); if (dst == NULL) { printf("Malloc Failed!\n"); return 0; } vssqrt(len, src, dst); /** * Output dst: * 0.0 0.5 4.0 1.0 * */
Fortran interface:
INTEGER :: LEN = 4 REAL(4) SRC(4) REAL(4) DST(4) DATA SRC /1, 2, 3, 4/ CALL VSSQRT(LEN, SRC1, SRC2, DST) ! ! OUTPUT DST: ! 1.000000000, 1.414213538, 1.732050776, 2.000000000 !