biod.pnpi.spb.ru
/
alexxy
/
gromacs.git
/ blob
commit
grep
author
committer
pickaxe
?
search:
re
252d2c9c01bacf98d1b28ecb396ac5d21d71589d
[alexxy/gromacs.git]
/
src
/
gmxlib
/
gmx_lapack
/
slapy2.c
1
#include <math.h>
2
#include "gmx_lapack.h"
3
4
#include <types/simple.h>
5
6
float
7
F77_FUNC(slapy2,SLAPY2)(float * x, float * y)
8
{
9
float xabs,yabs;
10
float w,z;
11
12
xabs = fabs(*x);
13
yabs = fabs(*y);
14
15
if(xabs>yabs) {
16
w = xabs;
17
z = yabs;
18
} else {
19
w = yabs;
20
z = xabs;
21
}
22
23
if( fabs(z)<GMX_FLOAT_MIN)
24
return w;
25
else {
26
z = z/w;
27
return w*sqrt(1.0+z*z);
28
}
29
}
30