void Mjoin( PATL, f77dotc_sub ) ( const int N, const TYPE * X, const int INCX, const TYPE * Y, const int INCY, TYPE * DOT ) { #ifdef ATL_FunkyInts const F77_INTEGER F77N = N, F77incx = INCX, F77incy = INCY; #else #define F77N N #define F77incx INCX #define F77incy INCY #endif if( INCX < 0 ) X -= ( ( 1 - N ) * INCX ) SHIFT; if( INCY < 0 ) Y -= ( ( 1 - N ) * INCY ) SHIFT; #ifdef RETURN_BY_STACK F77dotc(DOT, &F77N, X, &F77incx, Y, &F77incy); #else *((TYPE _Complex *)(DOT)) = F77dotc( &F77N, X, &F77incx, Y, &F77incy); #endif }
void Mjoin( PATL, f77dotc_sub ) ( const int N, const TYPE * X, const int INCX, const TYPE * Y, const int INCY, TYPE * DOT ) { #ifdef ATL_FunkyInts const F77_INTEGER F77N = N, F77incx = INCX, F77incy = INCY; #else #define F77N N #define F77incx INCX #define F77incy INCY #endif if( INCX < 0 ) X -= ( ( 1 - N ) * INCX ) SHIFT; if( INCY < 0 ) Y -= ( ( 1 - N ) * INCY ) SHIFT; F77dotc( &F77N, X, &F77incx, Y, &F77incy, DOT ); }