PLMD::blas Namespace Reference

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

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

Function Documentation

PLMD::blas::for

(

i__

= mp1; i__ <= i__2; i__ += 6

)

PLMD::blas::for

(

)

Initial value:

=  t1 + dx[i] * dy[i]   
                +    dx[i+1] * dy[i+1] 
                +    dx[i+2] * dy[i+2] 
                +    dx[i+3] * dy[i+3]   
                +    dx[i+4] * dy[i+4]

PLMD::blas::for

(

ix

= 1;ix<=max_ix;ix+=incx

)

PLMD::blas::if

(

n<=0||incx<=

0

)

PLMD::blas::if

(

incx!

= 1

)

PLMD::blas::if

(

m!

= 0

)

PLMD::blas::if

(

n<=

0

)

PLMD::blas::if

(

tra

= ='N'

)

PLMD::blas::if

(

trb

= ='N'

)

= k

PLMD::blas::if

(

m

= =0||n==0||((fabs(alpha)< PLUMED_GMX_DOUBLE_MIN||k==0)&&fabs(beta-1.0)< PLUMED_GMX_DOUBLE_EPS)

)

PLMD::blas::if

(

)

= incx * (1 - n)

PLMD::blas::if

(

n<=0||m<=0||

fabs(alpha)< PLUMED_GMX_DOUBLE_MIN &&fabs(beta-1.0)< PLUMED_GMX_DOUBLE_EPS

)

PLMD::blas::if

(

ch

= ='N'

)

PLMD::blas::if	(	incx	,
		0
	)		= 1

PLMD::blas::if	(	incy	,
		0
	)		= 1

PLMD::blas::if	(	fabs(beta-1.0)	,
		PLUMED_GMX_DOUBLE_EPS
	)

PLMD::blas::if

(

incx

= =1

)

PLMD::blas::if

(

n

= =1

)

PLMD::blas::if

(

n<=0||

incx = =0||incy==0

)

PLMD::blas::if

(

ch2

= = 'N'

)

= n

PLMD::blas::if

(

lside

)

PLMD::blas::if

(

incx<=

0

)

PLMD::blas::if

(

xside

= ='L'

)

PLMD::blas::if

(

n<=0||m<=0||

fabs(alpha)< PLUMED_GMX_FLOAT_MIN &&fabs(beta-1.0)< PLUMED_GMX_FLOAT_EPS

)

PLMD::blas::if	(	fabs(beta-1.0)	,
		PLUMED_GMX_FLOAT_EPS
	)

double PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dasum	,
		DASUM
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	daxpy	,
		DAXPY
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dcopy	,
		DCOPY
	)

double PLMD::blas::PLUMED_BLAS_F77_FUNC	(	ddot	,
		DDOT
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dgemm	,
		DGEMM
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dgemv	,
		DGEMV
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dger	,
		DGER
	)

double PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dnrm2	,
		DNRM2
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	drot	,
		DROT
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dscal	,
		DSCAL
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dswap	,
		DSWAP
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dsymv	,
		DSYMV
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dsyr2	,
		DSYR2
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dsyr2k	,
		DSYR2K
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dtrmm	,
		DTRMM
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dtrmv	,
		DTRMV
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	dtrsm	,
		DTRSM
	)		const

int PLMD::blas::PLUMED_BLAS_F77_FUNC	(	idamax	,
		IDAMAX
	)

int PLMD::blas::PLUMED_BLAS_F77_FUNC	(	isamax	,
		ISAMAX
	)

float PLMD::blas::PLUMED_BLAS_F77_FUNC	(	sasum	,
		SASUM
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	saxpy	,
		SAXPY
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	scopy	,
		SCOPY
	)

float PLMD::blas::PLUMED_BLAS_F77_FUNC	(	sdot	,
		SDOT
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	sgemm	,
		SGEMM
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	sgemv	,
		SGEMV
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	sger	,
		SGER
	)

float PLMD::blas::PLUMED_BLAS_F77_FUNC	(	snrm2	,
		SNRM2
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	srot	,
		SROT
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	sscal	,
		SSCAL
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	sswap	,
		SSWAP
	)

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	ssymv	,
		SSYMV
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	ssyr2	,
		SSYR2
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	ssyr2k	,
		SSYR2K
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	strmm	,
		STRMM
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	strmv	,
		STRMV
	)		const

void PLMD::blas::PLUMED_BLAS_F77_FUNC	(	strsm	,
		STRSM
	)		const

return scale * PLMD::blas::sqrt

(

ssq

)

Variable Documentation

void const char const char const char int int float float * PLMD::blas::a = a_offset

PLMD::blas::a_dim1 = lda

PLMD::blas::a_offset = 1 + a_dim1

float PLMD::blas::absxi

void const char const char const char int int float * PLMD::blas::alpha = *alpha__

void const char const char const char int int float * PLMD::blas::alpha__

void const char const char const char int int float float int float * PLMD::blas::b = b_offset

PLMD::blas::b_dim1 = ldb

PLMD::blas::b_offset = 1 + b_dim1

void const char int int float float int float int float * PLMD::blas::beta = *beta__

void const char int int float float int float int float * PLMD::blas::beta__

void const char int int float float int float int float float * PLMD::blas::c = *c__

void float int float int float * PLMD::blas::c__

const char PLMD::blas::ch =toupper(*uplo)

PLMD::blas::ch1 = toupper(*uplo)

PLMD::blas::ch2 = toupper(*trans)

float PLMD::blas::d1

float PLMD::blas::d2

float PLMD::blas::d3

void float * PLMD::blas::da =*da_arg

void float * PLMD::blas::da_arg

void const char const char const char * PLMD::blas::diag

float PLMD::blas::dmax = fabs(dx[0])

float PLMD::blas::dtemp = 0.

int float * PLMD::blas::dx

void float int float * PLMD::blas::dy

PLMD::blas::else

Initial value:

{

    
    
    for(i=0;i<(n-4);i+=4) {
      dy[i]   += da*dx[i];
      dy[i+1] += da*dx[i+1];
      dy[i+2] += da*dx[i+2];
      dy[i+3] += da*dx[i+3];
    }
    
    for(;i<n;i++)
      dy[i]   += da*dx[i]

void float * PLMD::blas::fact = *fact__

void float * PLMD::blas::fact__

int PLMD::blas::i

int PLMD::blas::i__

PLMD::blas::i__2 = n

return PLMD::blas::idxmax = 1

int float int * PLMD::blas::incx = *incx__

void const char const char int float int float int * PLMD::blas::incx__

Initial value:

{
    int i__1, i__2

float float int * PLMD::blas::incx_arg

void int float float int float int * PLMD::blas::incy = *incy_arg

void int float float int float int * PLMD::blas::incy__

Initial value:

{
    int i,ix,iy

float float int float int * PLMD::blas::incy_arg

Initial value:

{
  int i,ix,iy

int PLMD::blas::info = 0

int PLMD::blas::ix

int PLMD::blas::iy = 0

int PLMD::blas::j

int PLMD::blas::jx = jy = kx = ky = 0

int PLMD::blas::jy = incy * (1 - n)

void const char int int * PLMD::blas::k = *k__

void const char int int * PLMD::blas::k__

int PLMD::blas::kx = 1 - (lenx -1)*(incx)

else PLMD::blas::ky = 1 - (leny -1)*(incy)

int PLMD::blas::l

void const char const char const char int int float float int * PLMD::blas::lda = *lda__

void const char const char const char int int float float int * PLMD::blas::lda__

Initial value:

{
    int ix,kx,jy

void const char const char const char int int float float int float int * PLMD::blas::ldb = *ldb__

void const char const char const char int int float float int float int * PLMD::blas::ldb__

Initial value:

{
    int a_dim1, a_offset, b_dim1, b_offset, i__1, i__2, i__3

void const char int int float float int float int float float int * PLMD::blas::ldc = *ldc__

void const char int int float float int float int float float int * PLMD::blas::ldc__

Initial value:

{
  const char tra=toupper(*transa)

int PLMD::blas::lenx

int PLMD::blas::leny = n

int PLMD::blas::lside = (*side=='L' || *side=='l')

void const char const char const char int * PLMD::blas::m = n % 6

void const char const char const char int * PLMD::blas::m__

PLMD::blas::max_ix = 1+(n-1)*(incx)

int PLMD::blas::mp1 = m + 1

void const char const char const char int int * PLMD::blas::n = *n__

void const char const char const char int int * PLMD::blas::n__

int PLMD::blas::ncola = m

int PLMD::blas::nincx

int PLMD::blas::nounit = (*diag=='N' || *diag=='n')

int PLMD::blas::nrowa = k

else PLMD::blas::nrowb = n

PLMD::blas::return

void float int float int float float * PLMD::blas::s = *s__

void float int float int float float * PLMD::blas::s__

Initial value:

{
  int i,ix,iy

float PLMD::blas::scale = 0.0

float PLMD::blas::ssq = 1.0

float PLMD::blas::t

return PLMD::blas::t1 = 0.0

float PLMD::blas::temp

float PLMD::blas::temp1

float PLMD::blas::temp2

float PLMD::blas::tmp

void const char * PLMD::blas::trans

void const char const char * PLMD::blas::transa

void const char * PLMD::blas::transb

const char PLMD::blas::trb =toupper(*transb)

void const char * PLMD::blas::uplo

int PLMD::blas::upper = (*uplo=='U' || *uplo=='u')

void const char const char int float int float * PLMD::blas::x

const char PLMD::blas::xdiag = toupper(*diag)

const char PLMD::blas::xtrans = toupper(*transa)

const char PLMD::blas::xuplo = toupper(*uplo)

void int float float int float * PLMD::blas::y