3 * This source code is part of
7 * GROningen MAchine for Chemical Simulations
9 * Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.
10 * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
11 * Copyright (c) 2001-2008, The GROMACS development team,
12 * check out http://www.gromacs.org for more information.
14 * This program is free software; you can redistribute it and/or
15 * modify it under the terms of the GNU General Public License
16 * as published by the Free Software Foundation; either version 2
17 * of the License, or (at your option) any later version.
19 * If you want to redistribute modifications, please consider that
20 * scientific software is very special. Version control is crucial -
21 * bugs must be traceable. We will be happy to consider code for
22 * inclusion in the official distribution, but derived work must not
23 * be called official GROMACS. Details are found in the README & COPYING
24 * files - if they are missing, get the official version at www.gromacs.org.
26 * To help us fund GROMACS development, we humbly ask that you cite
27 * the papers on the package - you can find them in the top README file.
29 * For more info, check our website at http://www.gromacs.org
32 * Gallium Rubidium Oxygen Manganese Argon Carbon Silicon
38 #include "types/simple.h"
40 /* Suppress Cygwin compiler warnings from using newlib version of
49 * @brief Header definitions for the standard BLAS library.
51 * This is the subset of BLAS routines used for the
52 * linear algebra operations in Gromacs.
53 * Do NOT use this for other purposes - we only provide this as a
54 * simple fallback/reference implementation when no optimized BLAS
55 * is present. If you need an implementation for your own code
56 * there are several much faster versions out there.
58 * All routines are compatible with the BLAS reference implementation,
59 * meaning they assume fortran-style matrix row/column organization.
61 * There is plenty of documentation for these routines available
62 * at http://www.netlib.org/blas , so there is no point in repeating
74 #define F77_FUNC(name,NAME) name ## _
78 /* Double precision versions */
80 F77_FUNC(dasum
,DASUM
)(int *n
, double *dx
, int *incx
);
83 F77_FUNC(daxpy
,DAXPY
)(int *n
, double *da
, double *dx
, int *incx
, double *dy
, int *incy
);
86 F77_FUNC(dcopy
,DCOPY
)(int *n
, double *dx
, int *incx
, double *dy
, int *incy
);
89 F77_FUNC(ddot
,DDOT
)(int *n
, double *dx
, int *incx
, double *dy
, int *incy
);
92 F77_FUNC(dgemm
,DGEMM
)(const char *transa
, const char *transb
, int *m
, int *n
, int *k
,
93 double *alpha
, double *a
, int *lda
, double *b
, int *ldb
,
94 double *beta
, double *c
, int *ldc
);
97 F77_FUNC(dgemv
,DGEMV
)(const char *trans
, int *m
, int *n
, double *alpha
, double *a
, int *lda
,
98 double *x
, int *incx
, double *beta
, double *y
, int *incy
);
101 F77_FUNC(dger
,DGER
)(int *m
, int *n
, double *alpha
, double *x
, int *incx
,
102 double *y
, int *incy
, double *a
, int *lda
);
105 F77_FUNC(dnrm2
,DNRM2
)(int *n
, double *x
, int *incx
);
108 F77_FUNC(drot
,DROT
)(int *n
, double *dx
, int *incx
,
109 double *dy
, int *incy
, double *c
, double *s
);
112 F77_FUNC(dscal
,DSCAL
)(int *n
, double *fact
, double *dx
, int *incx
);
115 F77_FUNC(dswap
,DSWAP
)(int *n
, double *dx
, int *incx
, double *dy
, int *incy
);
118 F77_FUNC(dsymv
,DSYMV
)(const char *uplo
, int *n
, double *alpha
, double *a
, int *lda
,
119 double *x
, int *incx
, double *beta
, double *y
, int *incy
);
122 F77_FUNC(dsyr2
,DSYR2
)(const char *uplo
, int *n
, double *alpha
, double *x
, int *incx
,
123 double *y
, int *incy
, double *a
, int *lda
);
126 F77_FUNC(dsyr2k
,DSYR2K
)(const char *uplo
, const char *trans
, int *n
, int *k
, double *alpha
, double *a
,
127 int *lda
, double *b
, int *ldb
, double *beta
, double *c
, int *ldc
);
130 F77_FUNC(dtrmm
,DTRMM
)(const char *side
, const char *uplo
, const char *transa
, const char *diag
, int *m
, int *n
,
131 double *alpha
, double *a
, int *lda
, double *b
, int *ldb
);
134 F77_FUNC(dtrmv
,DTRMV
)(const char *uplo
, const char *trans
, const char *diag
, int *n
,
135 double *a
, int *lda
, double *x
, int *incx
);
138 F77_FUNC(dtrsm
,DTRSM
)(const char *side
, const char *uplo
, const char *transa
, const char *diag
, int *m
, int *n
,
139 double *alpha
, double *a
,int *lda
, double *b
, int *ldb
);
142 F77_FUNC(idamax
,IDAMAX
)(int *n
, double *dx
, int *incx
);
146 /* Single precision versions */
148 F77_FUNC(sasum
,SASUM
)(int *n
, float *dx
, int *incx
);
151 F77_FUNC(saxpy
,SAXPY
)(int *n
, float *da
, float *dx
, int *incx
, float *dy
, int *incy
);
154 F77_FUNC(scopy
,SCOPY
)(int *n
, float *dx
, int *incx
, float *dy
, int *incy
);
157 F77_FUNC(sdot
,SDOT
)(int *n
, float *dx
, int *incx
, float *dy
, int *incy
);
160 F77_FUNC(sgemm
,SGEMM
)(const char *transa
, const char *transb
, int *m
, int *n
, int *k
,
161 float *alpha
, float *a
, int *lda
, float *b
, int *ldb
,
162 float *beta
, float *c
, int *ldc
);
165 F77_FUNC(sgemv
,SGEMV
)(const char *trans
, int *m
, int *n
, float *alpha
, float *a
, int *lda
,
166 float *x
, int *incx
, float *beta
, float *y
, int *incy
);
169 F77_FUNC(sger
,SGER
)(int *m
, int *n
, float *alpha
, float *x
, int *incx
,
170 float *y
, int *incy
, float *a
, int *lda
);
173 F77_FUNC(snrm2
,SNRM2
)(int *n
, float *x
, int *incx
);
176 F77_FUNC(srot
,SROT
)(int *n
, float *dx
, int *incx
,
177 float *dy
, int *incy
, float *c
, float *s
);
180 F77_FUNC(sscal
,SSCAL
)(int *n
, float *fact
, float *dx
, int *incx
);
183 F77_FUNC(sswap
,SSWAP
)(int *n
, float *dx
, int *incx
, float *dy
, int *incy
);
186 F77_FUNC(ssymv
,SSYMV
)(const char *uplo
, int *n
, float *alpha
, float *a
, int *lda
,
187 float *x
, int *incx
, float *beta
, float *y
, int *incy
);
190 F77_FUNC(ssyr2
,SSYR2
)(const char *uplo
, int *n
, float *alpha
, float *x
, int *incx
,
191 float *y
, int *incy
, float *a
, int *lda
);
194 F77_FUNC(ssyr2k
,SSYR2K
)(const char *uplo
, const char *trans
, int *n
, int *k
, float *alpha
, float *a
,
195 int *lda
, float *b
, int *ldb
, float *beta
, float *c
, int *ldc
);
198 F77_FUNC(strmm
,STRMM
)(const char *side
, const char *uplo
, const char *transa
, const char *diag
, int *m
, int *n
,
199 float *alpha
, float *a
, int *lda
, float *b
, int *ldb
);
202 F77_FUNC(strmv
,STRMV
)(const char *uplo
, const char *trans
, const char *diag
, int *n
,
203 float *a
, int *lda
, float *x
, int *incx
);
206 F77_FUNC(strsm
,STRSM
)(const char *side
, const char *uplo
, const char *transa
, const char *diag
, int *m
, int *n
,
207 float *alpha
, float *a
,int *lda
, float *b
, int *ldb
);
210 F77_FUNC(isamax
,ISAMAX
)(int *n
, float *dx
, int *incx
);
219 #endif /* _BLAS_H_ */