3 #include "gromacs/utility/real.h"
4 #include "../gmx_blas.h"
7 F77_FUNC(sger
,SGER
)(int *m__
,
26 float alpha
= *alpha__
;
28 if(m
<=0 || n
<=0 || fabs(alpha
)<GMX_FLOAT_MIN
)
37 for(j
=0;j
<n
;j
++,jy
+=incy
)
38 if(fabs(y
[jy
])>GMX_FLOAT_MIN
) {
41 a
[j
*(lda
)+i
] += temp
*x
[i
];
50 for(j
=0;j
<n
;j
++,jy
+=incy
) {
51 if(fabs(y
[jy
])>GMX_FLOAT_MIN
) {
54 for(i
=0;i
<m
;i
++,ix
+=incx
)
55 a
[j
*(lda
)+i
] += temp
*x
[ix
];