Namespaces | Functions | Variables
blas.cpp File Reference
#include <math.h>
#include "blas.h"
#include <ctype.h>
#include "simple.h"

Namespaces

 PLMD
 
 PLMD::blas
 

Functions

double PLMD::blas::PLUMED_BLAS_F77_FUNC (dasum, DASUM)(int *n__
 
 PLMD::blas::if (n<=0||incx<=0)
 
 PLMD::blas::if (incx!=1)
 
 PLMD::blas::if (m!=0)
 
 PLMD::blas::for (i__=mp1;i__<=i__2;i__+=6)
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (daxpy, DAXPY)(int *n_arg
 
 PLMD::blas::if (n<=0) return
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dcopy, DCOPY)(int *n__
 
double PLMD::blas::PLUMED_BLAS_F77_FUNC (ddot, DDOT)(int *n_arg
 
 PLMD::blas::for (i=0;i< m;i++) t1+
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dgemm, DGEMM)(const char *transa
 
 PLMD::blas::if (tra=='N')
 
 PLMD::blas::if (trb=='N') nrowb = k
 
 PLMD::blas::if (m==0||n==0||((fabs(alpha)< PLUMED_GMX_DOUBLE_MIN||k==0)&&fabs(beta-1.0)< PLUMED_GMX_DOUBLE_EPS)) return
 
 PLMD::blas::if (fabs(alpha)< PLUMED_GMX_DOUBLE_MIN) = incx * (1 - n)
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dgemv, DGEMV)(const char *trans
 
 PLMD::blas::if (n<=0||m<=0||(fabs(alpha)< PLUMED_GMX_DOUBLE_MIN &&fabs(beta-1.0)< PLUMED_GMX_DOUBLE_EPS)) return
 
 PLMD::blas::if (ch=='N')
 
 PLMD::blas::if (incx >0) kx = 1
 
 PLMD::blas::if (incy >0) ky = 1
 
 PLMD::blas::if (fabs(beta-1.0)>PLUMED_GMX_DOUBLE_EPS)
 
 PLMD::blas::if (incx==1)
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dger, DGER)(int *m__
 
double PLMD::blas::PLUMED_BLAS_F77_FUNC (dnrm2, DNRM2)(int *n__
 
else PLMD::blas::if (n==1)
 
 PLMD::blas::for (ix=1;ix<=max_ix;ix+=incx)
 
return scale * PLMD::blas::sqrt (ssq)
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (drot, DROT)(int *n__
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dscal, DSCAL)(int *n__
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dswap, DSWAP)(int *n__
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dsymv, DSYMV)(const char *uplo
 
 PLMD::blas::if (n<=0||incx==0||incy==0) return
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dsyr2, DSYR2)(const char *uplo
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dsyr2k, DSYR2K)(const char *uplo
 
 PLMD::blas::if (ch2== 'N') nrowa = n
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dtrmm, DTRMM)(const char *side
 
 PLMD::blas::if (lside)
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dtrmv, DTRMV)(const char *uplo
 
 PLMD::blas::if (incx<=0)
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (dtrsm, DTRSM)(const char *side
 
 PLMD::blas::if (xside=='L')
 
int PLMD::blas::PLUMED_BLAS_F77_FUNC (idamax, IDAMAX)(int *n__
 
int PLMD::blas::PLUMED_BLAS_F77_FUNC (isamax, ISAMAX)(int *n__
 
float PLMD::blas::PLUMED_BLAS_F77_FUNC (sasum, SASUM)(int *n__
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (saxpy, SAXPY)(int *n_arg
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (scopy, SCOPY)(int *n__
 
float PLMD::blas::PLUMED_BLAS_F77_FUNC (sdot, SDOT)(int *n_arg
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (sgemm, SGEMM)(const char *transa
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (sgemv, SGEMV)(const char *trans
 
 PLMD::blas::if (n<=0||m<=0||(fabs(alpha)< PLUMED_GMX_FLOAT_MIN &&fabs(beta-1.0)< PLUMED_GMX_FLOAT_EPS)) return
 
 PLMD::blas::if (fabs(beta-1.0)>PLUMED_GMX_FLOAT_EPS)
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (sger, SGER)(int *m__
 
float PLMD::blas::PLUMED_BLAS_F77_FUNC (snrm2, SNRM2)(int *n__
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (srot, SROT)(int *n__
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (sscal, SSCAL)(int *n__
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (sswap, SSWAP)(int *n__
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (ssymv, SSYMV)(const char *uplo
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (ssyr2, SSYR2)(const char *uplo
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (ssyr2k, SSYR2K)(const char *uplo
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (strmm, STRMM)(const char *side
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (strmv, STRMV)(const char *uplo
 
void PLMD::blas::PLUMED_BLAS_F77_FUNC (strsm, STRSM)(const char *side
 

Variables

double double * PLMD::blas::dx
 
double double int * PLMD::blas::incx__
 
int PLMD::blas::i__
 
int PLMD::blas::m = n % 6
 
int PLMD::blas::mp1 = m + 1
 
double PLMD::blas::dtemp = 0.
 
int PLMD::blas::nincx
 
int PLMD::blas::n = *n__
 
int PLMD::blas::incx = *incx__
 
 PLMD::blas::i__2 = n
 
void double * PLMD::blas::da_arg
 
void double double int * PLMD::blas::incx_arg
 
void double double int double * PLMD::blas::dy
 
void double double int double int * PLMD::blas::incy_arg
 
double PLMD::blas::da =*da_arg
 
int PLMD::blas::incy = *incy_arg
 
 PLMD::blas::else
 
void double int double int * PLMD::blas::incy__
 
double PLMD::blas::t1 = 0.0
 
void const char * PLMD::blas::transb
 
void const char int * PLMD::blas::m__
 
void const char int int * PLMD::blas::n__
 
void const char int int int * PLMD::blas::k__
 
void const char int int int double * PLMD::blas::alpha__
 
void const char int int int double double * PLMD::blas::a = a_offset
 
void const char int int int double double int * PLMD::blas::lda__
 
void const char int int int double double int double * PLMD::blas::b = b_offset
 
void const char int int int double double int double int * PLMD::blas::ldb__
 
void const char int int int double double int double int double * PLMD::blas::beta__
 
void const char int int int double double int double int double double * PLMD::blas::c = *c__
 
void const char int int int double double int double int double double int * PLMD::blas::ldc__
 
const char PLMD::blas::trb =toupper(*transb)
 
double PLMD::blas::temp
 
int PLMD::blas::i
 
int PLMD::blas::j
 
int PLMD::blas::l
 
int PLMD::blas::nrowa = k
 
int PLMD::blas::ncola = m
 
int PLMD::blas::nrowb = n
 
int PLMD::blas::k = *k__
 
int PLMD::blas::lda = *lda__
 
int PLMD::blas::ldb = *ldb__
 
int PLMD::blas::ldc = *ldc__
 
double PLMD::blas::alpha = *alpha__
 
double PLMD::blas::beta = *beta__
 
void int int double double int double * PLMD::blas::x
 
void int int double double int double int double double * PLMD::blas::y
 
int PLMD::blas::lenx
 
int PLMD::blas::leny = n
 
int PLMD::blas::kx = 1 - (lenx -1)*(incx)
 
int PLMD::blas::ky = 1 - (leny -1)*(incy)
 
int PLMD::blas::jx = jy = kx = ky = 0
 
int PLMD::blas::jy = incy * (1 - n)
 
int PLMD::blas::ix
 
int PLMD::blas::iy = 0
 
 PLMD::blas::return
 
double PLMD::blas::ssq = 1.0
 
double PLMD::blas::scale = 0.0
 
double PLMD::blas::absxi
 
double PLMD::blas::t
 
 PLMD::blas::max_ix = 1+(n-1)*(incx)
 
void double int double int double * PLMD::blas::c__
 
void double int double int double double * PLMD::blas::s__
 
double PLMD::blas::s = *s__
 
void double * PLMD::blas::fact__
 
double PLMD::blas::fact = *fact__
 
double PLMD::blas::d1
 
double PLMD::blas::d2
 
double PLMD::blas::d3
 
double PLMD::blas::temp1
 
double PLMD::blas::temp2
 
const char PLMD::blas::ch =toupper(*uplo)
 
void const char * PLMD::blas::trans
 
 PLMD::blas::ch1 = toupper(*uplo)
 
 PLMD::blas::ch2 = toupper(*trans)
 
void const char * PLMD::blas::uplo
 
void const char const char * PLMD::blas::transa
 
void const char const char const char * PLMD::blas::diag
 
int PLMD::blas::info = 0
 
int PLMD::blas::lside = (*side=='L' || *side=='l')
 
int PLMD::blas::upper = (*uplo=='U' || *uplo=='u')
 
int PLMD::blas::nounit = (*diag=='N' || *diag=='n')
 
 PLMD::blas::a_dim1 = lda
 
 PLMD::blas::a_offset = 1 + a_dim1
 
 PLMD::blas::b_dim1 = ldb
 
 PLMD::blas::b_offset = 1 + b_dim1
 
const char PLMD::blas::xuplo = toupper(*uplo)
 
const char PLMD::blas::xtrans = toupper(*transa)
 
const char PLMD::blas::xdiag = toupper(*diag)
 
double PLMD::blas::dmax = fabs(dx[0])
 
double PLMD::blas::tmp
 
 PLMD::blas::idxmax = 1