/* * -- High Performance Computing Linpack Benchmark (HPL) * HPL - 1.0 - September 27, 2000 * Antoine P. Petitet * University of Tennessee, Knoxville * Innovative Computing Laboratories * (C) Copyright 2000 All Rights Reserved * * -- Copyright notice and Licensing terms: * * Redistribution and use in source and binary forms, with or without * modification, are permitted provided that the following conditions * are met: * * 1. Redistributions of source code must retain the above copyright * notice, this list of conditions and the following disclaimer. * * 2. Redistributions in binary form must reproduce the above copyright * notice, this list of conditions, and the following disclaimer in the * documentation and/or other materials provided with the distribution. * * 3. All advertising materials mentioning features or use of this * software must display the following acknowledgement: * This product includes software developed at the University of * Tennessee, Knoxville, Innovative Computing Laboratories. * * 4. The name of the University, the name of the Laboratory, or the * names of its contributors may not be used to endorse or promote * products derived from this software without specific written * permission. * * -- Disclaimer: * * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR * A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE UNIVERSITY * OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, * DATA OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. */ #ifndef HPL_BLAS_H #define HPL_BLAS_H /* * --------------------------------------------------------------------- * Include files * --------------------------------------------------------------------- */ #include "hpl_misc.h" /* * --------------------------------------------------------------------- * typedef definitions * --------------------------------------------------------------------- */ enum HPL_ORDER { HplRowMajor = 101, HplColumnMajor = 102 }; enum HPL_TRANS { HplNoTrans = 111, HplTrans = 112, HplConjTrans = 113 }; enum HPL_UPLO { HplUpper = 121, HplLower = 122 }; enum HPL_DIAG { HplNonUnit = 131, HplUnit = 132 }; enum HPL_SIDE { HplLeft = 141, HplRight = 142 }; #ifdef HPL_CALL_CBLAS /* * --------------------------------------------------------------------- * The C interface of the BLAS is available ... * --------------------------------------------------------------------- * #define macro constants * --------------------------------------------------------------------- */ #define CBLAS_INDEX int #define CBLAS_ORDER HPL_ORDER #define CblasRowMajor HplRowMajor #define CblasColMajor HplColMajor #define CBLAS_TRANSPOSE HPL_TRANS #define CblasNoTrans HplNoTrans #define CblasTrans HplTrans #define CblasConjTrans HplConjTrans #define CBLAS_UPLO HPL_UPLO #define CblasUpper HplUpper #define CblasLower HplLower #define CBLAS_DIAG HPL_DIAG #define CblasNonUnit HplNonUnit #define CblasUnit HplUnit #define CBLAS_SIDE HPL_SIDE #define CblasLeft HplLeft #define CblasRight HplRight /* * --------------------------------------------------------------------- * CBLAS Function prototypes * --------------------------------------------------------------------- */ CBLAS_INDEX cblas_idamax STDC_ARGS( ( const int, const double *, const int ) ); void cblas_dswap STDC_ARGS( ( const int, double *, const int, double *, const int ) ); void cblas_dcopy STDC_ARGS( ( const int, const double *, const int, double *, const int ) ); void cblas_daxpy STDC_ARGS( ( const int, const double, const double *, const int, double *, const int ) ); void cblas_dscal STDC_ARGS( ( const int, const double, double *, const int ) ); void cblas_dgemv STDC_ARGS( ( const enum CBLAS_ORDER, const enum CBLAS_TRANSPOSE, const int, const int, const double, const double *, const int, const double *, const int, const double, double *, const int ) ); void cblas_dger STDC_ARGS( ( const enum CBLAS_ORDER, const int, const int, const double, const double *, const int, const double *, const int, double *, const int ) ); void cblas_dtrsv STDC_ARGS( ( const enum CBLAS_ORDER, const enum CBLAS_UPLO, const enum CBLAS_TRANSPOSE, const enum CBLAS_DIAG, const int, const double *, const int, double *, const int ) ); void cblas_dgemm STDC_ARGS( ( const enum CBLAS_ORDER, const enum CBLAS_TRANSPOSE, const enum CBLAS_TRANSPOSE, const int, const int, const int, const double, const double *, const int, const double *, const int, const double, double *, const int ) ); void cblas_dtrsm STDC_ARGS( ( const enum CBLAS_ORDER, const enum CBLAS_SIDE, const enum CBLAS_UPLO, const enum CBLAS_TRANSPOSE, const enum CBLAS_DIAG, const int, const int, const double, const double *, const int, double *, const int ) ); /* * --------------------------------------------------------------------- * HPL C BLAS macro definition * --------------------------------------------------------------------- */ #define HPL_dswap cblas_dswap #define HPL_dcopy cblas_dcopy #define HPL_daxpy cblas_daxpy #define HPL_dscal cblas_dscal #define HPL_idamax cblas_idamax #define HPL_dgemv cblas_dgemv #define HPL_dtrsv cblas_dtrsv #define HPL_dger cblas_dger #define HPL_dgemm cblas_dgemm #define HPL_dtrsm cblas_dtrsm #endif #ifdef HPL_CALL_FBLAS /* * --------------------------------------------------------------------- * Use the Fortran 77 interface of the BLAS ... * --------------------------------------------------------------------- * Defaults: Add_, F77_INTEGER=int, StringSunStyle * --------------------------------------------------------------------- */ #ifndef NoChange /* #ifndef Upcase */ #ifndef UpCase #ifndef Add__ #ifndef Add_ #define Add_ #endif #endif #endif #endif #ifndef F77_INTEGER #define F77_INTEGER int #else #define HPL_USE_F77_INTEGER_DEF #endif #ifndef StringCrayStyle #ifndef StringStructVal #ifndef StringStructPtr #ifndef StringSunStyle #define StringSunStyle #endif #endif #endif #endif /* * --------------------------------------------------------------------- * Fortran 77 <-> C interface * --------------------------------------------------------------------- * * These macros identifies how Fortran routines will be called. * * Add_ : the Fortran compiler expects the name of C functions to be * in all lower case and to have an underscore postfixed it (Suns, Intel * compilers expect this). * * NoChange : the Fortran compiler expects the name of C functions to be * in all lower case (IBM RS6K compilers do this). * * UpCase : the Fortran compiler expects the name of C functions to be * in all upcase. (Cray compilers expect this). * * Add__ : the Fortran compiler in use is f2c, a Fortran to C conver- * ter. */ #ifdef NoChange /* * These defines set up the naming scheme required to have a FORTRAN * routine called by a C routine with the following FORTRAN to C inter- * face: * * FORTRAN DECLARATION C CALL * SUBROUTINE DGEMM(...) dgemm(...) */ #define F77dswap dswap #define F77dscal dscal #define F77dcopy dcopy #define F77daxpy daxpy #define F77idamax idamax #define F77dgemv dgemv #define F77dtrsv dtrsv #define F77dger dger #define F77dgemm dgemm #define F77dtrsm dtrsm #endif #ifdef UpCase /* * These defines set up the naming scheme required to have a FORTRAN * routine called by a C routine with the following FORTRAN to C inter- * face: * * FORTRAN DECLARATION C CALL * SUBROUTINE DGEMM(...) DGEMM(...) */ /* ORIGINAL: #define F77dswap DSWAP #define F77dscal DSCAL #define F77dcopy DCOPY #define F77daxpy DAXPY #define F77idamax IDAMAX #define F77dgemv DGEMV #define F77dtrsv DTRSV #define F77dger DGER #define F77dgemm DGEMM #define F77dtrsm DTRSM */ //Changed for Piotr's suggesion (email Jan 26, 04) #define F77dswap SSWAP #define F77dscal SSCAL #define F77dcopy SCOPY #define F77daxpy SAXPY #define F77idamax ISAMAX #define F77dgemv SGEMV #define F77dtrsv STRSV #define F77dger SGER #define F77dgemm SGEMM #define F77dtrsm STRSM #endif #ifdef Add_ /* * These defines set up the naming scheme required to have a FORTRAN * routine called by a C routine with the following FORTRAN to C inter- * face: * * FORTRAN DECLARATION C CALL * SUBROUTINE DGEMM(...) dgemm_(...) */ #define F77dswap dswap_ #define F77dscal dscal_ #define F77dcopy dcopy_ #define F77daxpy daxpy_ #define F77idamax idamax_ #define F77dgemv dgemv_ #define F77dtrsv dtrsv_ #define F77dger dger_ #define F77dgemm dgemm_ #define F77dtrsm dtrsm_ #endif #ifdef Add__ /* * These defines set up the naming scheme required to have a FORTRAN * routine called by a C routine with the following FORTRAN to C inter- * face: * * FORTRAN DECLARATION C CALL * SUBROUTINE DGEMM(...) dgemm_(...) */ #define F77dswap dswap_ #define F77dscal dscal_ #define F77dcopy dcopy_ #define F77daxpy daxpy_ #define F77idamax idamax_ #define F77dgemv dgemv_ #define F77dtrsv dtrsv_ #define F77dger dger_ #define F77dgemm dgemm_ #define F77dtrsm dtrsm_ #endif /* * --------------------------------------------------------------------- * Typedef definitions and conversion utilities * --------------------------------------------------------------------- */ #ifdef StringCrayStyle #include /* Type of character argument in a FORTRAN call */ #define F77_CHAR _fcd /* Character conversion utilities */ #define HPL_F2C_CHAR(c) (*(_fcdtocp(c) )) #define HPL_C2F_CHAR(c) (_cptofcd(&(c), 1)) #define F77_CHAR_DECL F77_CHAR /* input CHARACTER*1 */ #endif /* ------------------------------------------------------------------ */ #ifdef StringStructVal /* Type of character argument in a FORTRAN call */ typedef struct { char *cp; F77_INTEGER len; } F77_CHAR; /* Character conversion utilities */ #define HPL_F2C_CHAR(c) (*(c.cp)) #define F77_CHAR_DECL F77_CHAR /* input CHARACTER*1 */ #endif /* ------------------------------------------------------------------ */ #ifdef StringStructPtr /* Type of character argument in a FORTRAN call */ typedef struct { char *cp; F77_INTEGER len; } F77_CHAR; /* Character conversion utilities */ #define HPL_F2C_CHAR(c) (*(c->cp)) #define F77_CHAR_DECL F77_CHAR * /* input CHARACTER*1 */ #endif /* ------------------------------------------------------------------ */ #ifdef StringSunStyle /* Type of character argument in a FORTRAN call */ #define F77_CHAR char * /* Character conversion utilities */ #define HPL_F2C_CHAR(c) (*(c)) #define HPL_C2F_CHAR(c) (&(c)) #define F77_CHAR_DECL F77_CHAR /* input CHARACTER*1 */ #define F77_1_CHAR , F77_INTEGER #define F77_2_CHAR F77_1_CHAR F77_1_CHAR #define F77_3_CHAR F77_2_CHAR F77_1_CHAR #define F77_4_CHAR F77_3_CHAR F77_1_CHAR #endif /* ------------------------------------------------------------------ */ #ifndef F77_1_CHAR #define F77_1_CHAR #define F77_2_CHAR #define F77_3_CHAR #define F77_4_CHAR #endif #define F77_INT_DECL const F77_INTEGER * /* input integer */ #define F77_SIN_DECL const double * /* input scalar */ #define F77_VIN_DECL const double * /* input vector */ #define F77_VINOUT_DECL double * /* input/output matrix */ #define F77_MIN_DECL const double * /* input matrix */ #define F77_MINOUT_DECL double * /* input/output matrix */ //#ifdef CRAY /* Type of FORTRAN functions */ #if defined(CRAY) && !defined(__crayx1) && !defined(_CRAYSV1) && !defined(_CRAYSV2) #define F77_VOID_FUN extern fortran void /* subroutine */ #define F77_INT_FUN extern fortran int /* integer function */ #else #define F77_VOID_FUN extern void /* subroutine */ #define F77_INT_FUN extern int /* integer function */ #endif /* * --------------------------------------------------------------------- * Fortran 77 BLAS function prototypes * --------------------------------------------------------------------- */ F77_VOID_FUN F77dswap STDC_ARGS( ( F77_INT_DECL, F77_VINOUT_DECL, F77_INT_DECL, F77_VINOUT_DECL, F77_INT_DECL ) ); F77_VOID_FUN F77dscal STDC_ARGS( ( F77_INT_DECL, F77_SIN_DECL, F77_VINOUT_DECL, F77_INT_DECL ) ); F77_VOID_FUN F77dcopy STDC_ARGS( ( F77_INT_DECL, F77_VIN_DECL, F77_INT_DECL, F77_VINOUT_DECL, F77_INT_DECL ) ); F77_VOID_FUN F77daxpy STDC_ARGS( ( F77_INT_DECL, F77_SIN_DECL, F77_VIN_DECL, F77_INT_DECL, F77_VINOUT_DECL, F77_INT_DECL ) ); F77_INT_FUN F77idamax STDC_ARGS( ( F77_INT_DECL, F77_VIN_DECL, F77_INT_DECL ) ); F77_VOID_FUN F77dgemv STDC_ARGS( ( F77_CHAR_DECL, F77_INT_DECL, F77_INT_DECL, F77_SIN_DECL, F77_MIN_DECL, F77_INT_DECL, F77_VIN_DECL, F77_INT_DECL, F77_SIN_DECL, F77_VINOUT_DECL, F77_INT_DECL F77_1_CHAR ) ); F77_VOID_FUN F77dger STDC_ARGS( ( F77_INT_DECL, F77_INT_DECL, F77_SIN_DECL, F77_VIN_DECL, F77_INT_DECL, F77_VIN_DECL, F77_INT_DECL, F77_MINOUT_DECL, F77_INT_DECL ) ); F77_VOID_FUN F77dtrsv STDC_ARGS( ( F77_CHAR_DECL, F77_CHAR_DECL, F77_CHAR_DECL, F77_INT_DECL, F77_MIN_DECL, F77_INT_DECL, F77_VINOUT_DECL, F77_INT_DECL F77_3_CHAR ) ); F77_VOID_FUN F77dgemm STDC_ARGS( ( F77_CHAR_DECL, F77_CHAR_DECL, F77_INT_DECL, F77_INT_DECL, F77_INT_DECL, F77_SIN_DECL, F77_MIN_DECL, F77_INT_DECL, F77_MIN_DECL, F77_INT_DECL, F77_SIN_DECL, F77_MINOUT_DECL, F77_INT_DECL F77_2_CHAR ) ); F77_VOID_FUN F77dtrsm STDC_ARGS( ( F77_CHAR_DECL, F77_CHAR_DECL, F77_CHAR_DECL, F77_CHAR_DECL, F77_INT_DECL, F77_INT_DECL, F77_SIN_DECL, F77_MIN_DECL, F77_INT_DECL, F77_MINOUT_DECL, F77_INT_DECL F77_4_CHAR ) ); #endif /* * --------------------------------------------------------------------- * HPL BLAS Function prototypes * --------------------------------------------------------------------- */ #ifndef HPL_CALL_CBLAS int HPL_idamax STDC_ARGS( ( const int, const double *, const int ) ); void HPL_daxpy STDC_ARGS( ( const int, const double, const double *, const int, double *, const int ) ); void HPL_dcopy STDC_ARGS( ( const int, const double *, const int, double *, const int ) ); void HPL_dscal STDC_ARGS( ( const int, const double, double *, const int ) ); void HPL_dswap STDC_ARGS( ( const int, double *, const int, double *, const int ) ); void HPL_dgemv STDC_ARGS( ( const enum HPL_ORDER, const enum HPL_TRANS, const int, const int, const double, const double *, const int, const double *, const int, const double, double *, const int ) ); void HPL_dger STDC_ARGS( ( const enum HPL_ORDER, const int, const int, const double, const double *, const int, double *, const int, double *, const int ) ); void HPL_dtrsv STDC_ARGS( ( const enum HPL_ORDER, const enum HPL_UPLO, const enum HPL_TRANS, const enum HPL_DIAG, const int, const double *, const int, double *, const int ) ); void HPL_dgemm STDC_ARGS( ( const enum HPL_ORDER, const enum HPL_TRANS, const enum HPL_TRANS, const int, const int, const int, const double, const double *, const int, const double *, const int, const double, double *, const int ) ); void HPL_dtrsm STDC_ARGS( ( const enum HPL_ORDER, const enum HPL_SIDE, const enum HPL_UPLO, const enum HPL_TRANS, const enum HPL_DIAG, const int, const int, const double, const double *, const int, double *, const int ) ); #endif #endif /* * hpl_blas.h */