181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray/*
281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * cblas_saxpy.c
381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *
481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * The program is a C interface to saxpy.
581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * It calls the fortran wrapper before calling saxpy.
681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *
781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray * Written by Keita Teranishi.  2/11/1998
881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray *
981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray */
1081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas.h"
1181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#include "cblas_f77.h"
1281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murrayvoid cblas_saxpy( const int N, const float alpha, const float *X,
1381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray                       const int incX, float *Y, const int incY)
1481253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray{
1581253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#ifdef F77_INT
1681253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray   F77_INT F77_N=N, F77_incX=incX, F77_incY=incY;
1781253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#else
1881253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray   #define F77_N N
1981253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray   #define F77_incX incX
2081253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray   #define F77_incY incY
2181253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray#endif
2281253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray   F77_saxpy( &F77_N, &alpha, X, &F77_incX, Y, &F77_incY);
2381253e9afbc34b99e9adc22ddf33a1bfac56c697Tim Murray}
24