biod.pnpi.spb.ru
/
alexxy
/
gromacs.git
/ blob
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
history
|
raw
|
HEAD
Merge release-5-0 into master
[alexxy/gromacs.git]
/
src
/
gromacs
/
linearalgebra
/
gmx_lapack
/
dlapy2.c
1
#include <math.h>
2
#include "../gmx_lapack.h"
3
4
#include "gromacs/utility/real.h"
5
6
double
7
F77_FUNC(dlapy2,DLAPY2)(double * x, double * y)
8
{
9
double xabs,yabs;
10
double 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_DOUBLE_MIN)
24
return w;
25
else {
26
z = z/w;
27
return w*sqrt(1.0+z*z);
28
}
29
}
30