181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray/* 281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * cblas_cgemv.c 381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * The program is a C interface of cgemv 481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * Keita Teranishi 5/20/98 681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * 781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray */ 881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include <stdio.h> 981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include <stdlib.h> 1081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas.h" 1181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas_f77.h" 1281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murrayvoid cblas_cgemv(const enum CBLAS_ORDER order, 1381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const enum CBLAS_TRANSPOSE TransA, const int M, const int N, 1481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const void *alpha, const void *A, const int lda, 1581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const void *X, const int incX, const void *beta, 1681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray void *Y, const int incY) 1781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray{ 1881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray char TA; 1981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_CHAR 2081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_CHAR F77_TA; 2181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_TA &TA 2381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 2481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_INT 2581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_INT F77_M=M, F77_N=N, F77_lda=lda, F77_incX=incX, F77_incY=incY; 2681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else 2781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_M M 2881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_N N 2981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_lda lda 3081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_incX incx 3181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #define F77_incY incY 3281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif 3381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 3481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray int n=0, i=0, incx=incX; 3581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const float *xx= (const float *)X; 3681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray float ALPHA[2],BETA[2]; 3781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray int tincY, tincx; 3881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray float *x=(float *)X, *y=(float *)Y, *st=0, *tx=0; 3981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray const float *stx = x; 4081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int CBLAS_CallFromC; 4181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray extern int RowMajorStrg; 4281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 4381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 4481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 1; 4581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 4681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (order == CblasColMajor) 4781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 4881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (TransA == CblasNoTrans) TA = 'N'; 4981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (TransA == CblasTrans) TA = 'T'; 5081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (TransA == CblasConjTrans) TA = 'C'; 5181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 5281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 5381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_cgemv","Illegal TransA setting, %d\n", TransA); 5481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 5581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 5681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 5781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 5881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 5981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TA = C2F_CHAR(&TA); 6081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 6181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_cgemv(F77_TA, &F77_M, &F77_N, alpha, A, &F77_lda, X, &F77_incX, 6281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray beta, Y, &F77_incY); 6381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 6481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (order == CblasRowMajor) 6581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 6681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 1; 6781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 6881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (TransA == CblasNoTrans) TA = 'T'; 6981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (TransA == CblasTrans) TA = 'N'; 7081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else if (TransA == CblasConjTrans) 7181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 7281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray ALPHA[0]= *( (const float *) alpha ); 7381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray ALPHA[1]= -( *( (const float *) alpha+1) ); 7481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray BETA[0]= *( (const float *) beta ); 7581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray BETA[1]= -( *( (const float *) beta+1 ) ); 7681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray TA = 'N'; 7781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (M > 0) 7881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 7981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray n = M << 1; 8081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x = malloc(n*sizeof(float)); 8181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray tx = x; 8281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if( incX > 0 ) { 8381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray i = incX << 1 ; 8481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray tincx = 2; 8581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray st= x+n; 8681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } else { 8781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray i = incX *(-2); 8881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray tincx = -2; 8981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray st = x-2; 9081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x +=(n-2); 9181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 9281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 9381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray do 9481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 9581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *x = *xx; 9681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x[1] = -xx[1]; 9781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x += tincx ; 9881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray xx += i; 9981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 10081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray while (x != st); 10181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray x=tx; 10281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 10381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_incX = 1; 10481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 10581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if(incY > 0) 10681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray tincY = incY; 10781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 10881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray tincY = -incY; 10981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 11081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y++; 11181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 11281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (N > 0) 11381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 11481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray i = tincY << 1; 11581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray n = i * N ; 11681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray st = y + n; 11781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray do { 11881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *y = -(*y); 11981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y += i; 12081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } while(y != st); 12181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y -= n; 12281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 12381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray stx = x; 12481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 12581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else stx = (const float *)X; 12681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 12781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 12881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 12981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray cblas_xerbla(2, "cblas_cgemv","Illegal TransA setting, %d\n", TransA); 13081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 13181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 13281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 13381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 13481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #ifdef F77_CHAR 13581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_TA = C2F_CHAR(&TA); 13681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray #endif 13781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (TransA == CblasConjTrans) 13881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_cgemv(F77_TA, &F77_N, &F77_M, ALPHA, A, &F77_lda, stx, 13981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray &F77_incX, BETA, Y, &F77_incY); 14081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else 14181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray F77_cgemv(F77_TA, &F77_N, &F77_M, alpha, A, &F77_lda, x, 14281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray &F77_incX, beta, Y, &F77_incY); 14381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray 14481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (TransA == CblasConjTrans) 14581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 14681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (x != (const float *)X) free(x); 14781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray if (N > 0) 14881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 14981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray do 15081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray { 15181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *y = -(*y); 15281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray y += i; 15381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 15481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray while (y != st); 15581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 15681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 15781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray } 15881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray else cblas_xerbla(1, "cblas_cgemv", "Illegal Order setting, %d\n", order); 15981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray CBLAS_CallFromC = 0; 16081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray RowMajorStrg = 0; 16181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray return; 16281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray} 163