Lines Matching refs:ax
53 double x, y, xh, yh, w, ax, ay;
65 /* force ax = |x| ~>~ ay = |y| */
67 ax = fabs(y);
76 ax = fabs(x);
88 return ((ax == ay)? ay : ax);
90 return ((ay == ax)? ax : ay);
92 return (ax * ay);
97 ax += ay;
98 return (ax);
100 ax *= twom768;
110 return (ax);
112 return (ax + ay);
114 ax *= two768;
117 if (ax == zero) /* guard subnormal flush to zero */
118 return (ax);
119 ix = ((int *)&ax)[HIWORD];
125 return (ax * twom768);
135 return ((ax + ay) * twom768);
141 return (ax + ay);
145 * Medium range ax and ay with max{|ax/ay|,|ay/ax|} bounded by 2^32.
151 ax += ay / (ax + sqrt(ax * ax + ay));
155 w = sqrt(ax * ax + ay * ay);
156 ax += ay / ((ax + w) / ay);
159 w = ax - ay;
163 ay = ay * ay + (ax - xh) * (ax + xh);
164 ax = sqrt(xh * xh + ay);
166 ax = ax + ax;
171 ay = w * w + ((ax - xh) * yh + (ay - yh) * ax);
172 ax = sqrt(xh * yh + ay);
177 ax *= twom768;
179 ax *= two768; /* must generate side effect here */
181 return (ax);