GotoBlas2之xASUM

最新推荐文章于 2024-09-20 21:00:57 发布

原创最新推荐文章于 2024-09-20 21:00:57 发布 · 1k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#vector #float #integer #numbers #function #table

C/C++ 专栏收录该内容

36 篇文章

订阅专栏

本文展示了如何使用xASUM函数计算不同类型的数组元素的绝对值之和，包括浮点数、双精度数、复数和双复数类型。通过设置不同的参数如N（元素数量）、INCX（步长）等，可以获取特定条件下数组元素的求和结果。

xASUM 绝对值求和。

#include <stdio.h> #include "gotoblas2.h" #pragma comment(lib,"libgoto2.lib") //FUNCTION xASUM (N , X, INCX) /* SASUM, DASUM, SCASUM, and DZASUM (Sum of the Magnitudes of the Elements in a Vector) On Entry n is the number of elements in vector x. Specified as: an integer; n greater than or equal to 0. x is the vector x of length n. Specified as: a one-dimensional array of (at least) length 1+(n-1)|incx|, containing numbers of the data type indicated in Table 1. incx is the stride for vector x. Specified as: an integer. It can have any value. On Return Function value is the result of the summation. SASUM and DASUM compute the sum of the absolute values of the elements of x. SCASUM and DZASUM compute the sum of the absolute values of the real and imaginary parts of the elements of x */ int main(void) { float X1[9]={1.0, 2.0, 7.0, -8.0, -5.0, -10.0, -9.0, 10.0, 6.0}; double X2[9]={1.0, 11.0, 7.0, -8.0, -5.0, -10.0, -9.0, 10.0, 6.0}; myccomplex_t X3[5]={{9.0 , 2.0}, {7.0,-8.0},{-5.0 , -10.0} , {-4.0 , 10.0}, {6.0 , 3.0}}; myzcomplex_t X4[5]={{9.0 , 2.0}, {7.0,-8.0},{-5.0 , -10.0} , {-4.0 , 10.0}, {6.0 , 3.0}}; float I1,I3; double I2,I4; blasint N; blasint INCX; N=9; INCX=1; I1=SASUM( &N , X1 , &INCX ); printf(" The IASUM is %.3f/n",I1); N=5; INCX=2; I2 = dasum( &N , X2 , &INCX ); printf(" The IASUM is %.3lf/n",I2); N=5; INCX=1; I3=scasum( &N , (float *)X3 , &INCX ); printf(" The IASUM is %.3f/n",I3); N=3; INCX=2; I4 = dzasum( &N , (double *)X4 , &INCX ); printf(" The IASUM is %.3lf/n",I4); N=5; INCX=0; I1 = sasum( &N , X1 , &INCX ); I2 = dasum( &N , X2 , &INCX ); N=5; I3 = scasum( &N , (float *)X3 , &INCX ); I4 = dzasum( &N , (double *)X4 , &INCX ); if(I1<1e-5&&I2<1e-5&&I3<1e-5&&I4<1e-5) printf("/n xASUM if incx <= 0 , Return 0/n"); N=5; INCX=0; I1 = sasum( &N , X1 , &INCX ); I2 = dasum( &N , X2 , &INCX ); N=0; INCX=1; I3 = scasum( &N , (float *)X3 , &INCX ); I4 = dzasum( &N , (double *)X4 , &INCX ); if(I1<1e-5&&I2<1e-5&&I3<1e-5&&I4<1e-5) printf(" xASUM If N <= 0 , Return 0/n"); return 0; } /*result The IASUM is 58.000 The IASUM is 28.000 The IASUM is 64.000 The IASUM is 35.000 xASUM if incx <= 0 , Return 0 xASUM If N <= 0 , Return 0 */