|
8 | 8 |
|
9 | 9 | /*----------------------Single Precision Functions----------------------------*/ |
10 | 10 |
|
11 | | -float Fresnel_Cosine_Taylor_to_Asymptotic_Float(float x){ |
12 | | - |
| 11 | +float Fresnel_Cosine_Taylor_to_Asymptotic_Float(float x) |
| 12 | +{ |
13 | 13 | /* Variables for S(x) and powers of x, respectively. */ |
14 | 14 | float cx, arg; |
15 | 15 | arg = x*x; |
@@ -72,7 +72,8 @@ float Fresnel_Cosine_Taylor_to_Asymptotic_Float(float x){ |
72 | 72 | } |
73 | 73 | } |
74 | 74 |
|
75 | | -float Fresnel_Cosine_While_to_Asymptotic_Float(float x){ |
| 75 | +float Fresnel_Cosine_While_to_Asymptotic_Float(float x) |
| 76 | +{ |
76 | 77 | float FRESNEL_COSINE_TAYLOR_COEFFICIENTS[30] = { |
77 | 78 | FRESNEL_COSINE_TAYLOR_00, FRESNEL_COSINE_TAYLOR_01, |
78 | 79 | FRESNEL_COSINE_TAYLOR_02, FRESNEL_COSINE_TAYLOR_03, |
@@ -150,7 +151,8 @@ float Fresnel_Cosine_While_to_Asymptotic_Float(float x){ |
150 | 151 | } |
151 | 152 | } |
152 | 153 |
|
153 | | -float Fresnel_Cosine_Heald_Rational_EPS_Minus_Three_Float(float x){ |
| 154 | +float Fresnel_Cosine_Heald_Rational_EPS_Minus_Three_Float(float x) |
| 155 | +{ |
154 | 156 | float A, R, a, b, c, d, sgn_x; |
155 | 157 | sgn_x = (x>0)-(x<0); |
156 | 158 | x *= SQRT_2_BY_PI*sgn_x; |
@@ -178,7 +180,8 @@ float Fresnel_Cosine_Heald_Rational_EPS_Minus_Three_Float(float x){ |
178 | 180 | return sgn_x*(SQRT_PI_BY_8 - R*sinf(A)); |
179 | 181 | } |
180 | 182 |
|
181 | | -float Fresnel_Cosine_Heald_Rational_EPS_Minus_Four_Float(float x){ |
| 183 | +float Fresnel_Cosine_Heald_Rational_EPS_Minus_Four_Float(float x) |
| 184 | +{ |
182 | 185 | float A, R, a, b, c, d, sgn_x; |
183 | 186 | sgn_x = (x>0)-(x<0); |
184 | 187 | x *= SQRT_2_BY_PI*sgn_x; |
@@ -208,7 +211,8 @@ float Fresnel_Cosine_Heald_Rational_EPS_Minus_Four_Float(float x){ |
208 | 211 | return sgn_x*(SQRT_PI_BY_8 - R*sinf(A)); |
209 | 212 | } |
210 | 213 |
|
211 | | -float Fresnel_Cosine_Heald_Rational_EPS_Minus_Six_Float(float x){ |
| 214 | +float Fresnel_Cosine_Heald_Rational_EPS_Minus_Six_Float(float x) |
| 215 | +{ |
212 | 216 | float A, R, a, b, c, d, sgn_x; |
213 | 217 | sgn_x = (x>0)-(x<0); |
214 | 218 | x *= SQRT_2_BY_PI*sgn_x; |
@@ -242,7 +246,8 @@ float Fresnel_Cosine_Heald_Rational_EPS_Minus_Six_Float(float x){ |
242 | 246 | return sgn_x*(SQRT_PI_BY_8 - R*sinf(A)); |
243 | 247 | } |
244 | 248 |
|
245 | | -float Fresnel_Cosine_Heald_Rational_EPS_Minus_Eight_Float(float x){ |
| 249 | +float Fresnel_Cosine_Heald_Rational_EPS_Minus_Eight_Float(float x) |
| 250 | +{ |
246 | 251 | float A, R, a, b, c, d, sgn_x; |
247 | 252 | sgn_x = (x>0)-(x<0); |
248 | 253 | x *= SQRT_2_BY_PI*sgn_x; |
@@ -286,8 +291,8 @@ float Fresnel_Cosine_Heald_Rational_EPS_Minus_Eight_Float(float x){ |
286 | 291 |
|
287 | 292 | /*----------------------Double Precision Functions----------------------------*/ |
288 | 293 |
|
289 | | -double Fresnel_Cosine_Taylor_to_Asymptotic_Double(double x){ |
290 | | - |
| 294 | +double Fresnel_Cosine_Taylor_to_Asymptotic_Double(double x) |
| 295 | +{ |
291 | 296 | /* Variables for S(x) and powers of x, respectively. */ |
292 | 297 | double cx, arg; |
293 | 298 | arg = x*x; |
@@ -356,7 +361,8 @@ double Fresnel_Cosine_Taylor_to_Asymptotic_Double(double x){ |
356 | 361 | } |
357 | 362 | } |
358 | 363 |
|
359 | | -double Fresnel_Cosine_While_to_Asymptotic_Double(double x){ |
| 364 | +double Fresnel_Cosine_While_to_Asymptotic_Double(double x) |
| 365 | +{ |
360 | 366 | double FRESNEL_COSINE_TAYLOR_COEFFICIENTS[27] = { |
361 | 367 | FRESNEL_COSINE_TAYLOR_00, FRESNEL_COSINE_TAYLOR_01, |
362 | 368 | FRESNEL_COSINE_TAYLOR_02, FRESNEL_COSINE_TAYLOR_03, |
@@ -436,7 +442,8 @@ double Fresnel_Cosine_While_to_Asymptotic_Double(double x){ |
436 | 442 | } |
437 | 443 | } |
438 | 444 |
|
439 | | -double Fresnel_Cosine_Heald_Rational_EPS_Minus_Three_Double(double x){ |
| 445 | +double Fresnel_Cosine_Heald_Rational_EPS_Minus_Three_Double(double x) |
| 446 | +{ |
440 | 447 | double A, R, a, b, c, d, sgn_x; |
441 | 448 | sgn_x = (x>0)-(x<0); |
442 | 449 | x *= SQRT_2_BY_PI*sgn_x; |
@@ -464,7 +471,8 @@ double Fresnel_Cosine_Heald_Rational_EPS_Minus_Three_Double(double x){ |
464 | 471 | return sgn_x*(SQRT_PI_BY_8 - R*sin(A)); |
465 | 472 | } |
466 | 473 |
|
467 | | -double Fresnel_Cosine_Heald_Rational_EPS_Minus_Four_Double(double x){ |
| 474 | +double Fresnel_Cosine_Heald_Rational_EPS_Minus_Four_Double(double x) |
| 475 | +{ |
468 | 476 | double A, R, a, b, c, d, sgn_x; |
469 | 477 | sgn_x = (x>0)-(x<0); |
470 | 478 | x *= SQRT_2_BY_PI*sgn_x; |
@@ -494,7 +502,8 @@ double Fresnel_Cosine_Heald_Rational_EPS_Minus_Four_Double(double x){ |
494 | 502 | return sgn_x*(SQRT_PI_BY_8 - R*sin(A)); |
495 | 503 | } |
496 | 504 |
|
497 | | -double Fresnel_Cosine_Heald_Rational_EPS_Minus_Six_Double(double x){ |
| 505 | +double Fresnel_Cosine_Heald_Rational_EPS_Minus_Six_Double(double x) |
| 506 | +{ |
498 | 507 | double A, R, a, b, c, d, sgn_x; |
499 | 508 | sgn_x = (x>0)-(x<0); |
500 | 509 | x *= SQRT_2_BY_PI*sgn_x; |
@@ -528,7 +537,8 @@ double Fresnel_Cosine_Heald_Rational_EPS_Minus_Six_Double(double x){ |
528 | 537 | return sgn_x*(SQRT_PI_BY_8 - R*sin(A)); |
529 | 538 | } |
530 | 539 |
|
531 | | -double Fresnel_Cosine_Heald_Rational_EPS_Minus_Eight_Double(double x){ |
| 540 | +double Fresnel_Cosine_Heald_Rational_EPS_Minus_Eight_Double(double x) |
| 541 | +{ |
532 | 542 | double A, R, a, b, c, d, sgn_x; |
533 | 543 | sgn_x = (x>0)-(x<0); |
534 | 544 | x *= SQRT_2_BY_PI*sgn_x; |
@@ -572,8 +582,8 @@ double Fresnel_Cosine_Heald_Rational_EPS_Minus_Eight_Double(double x){ |
572 | 582 |
|
573 | 583 | /*--------------------Long Double Precision Functions-------------------------*/ |
574 | 584 |
|
575 | | -long double Fresnel_Cosine_Taylor_to_Asymptotic_Long_Double(long double x){ |
576 | | - |
| 585 | +long double Fresnel_Cosine_Taylor_to_Asymptotic_Long_Double(long double x) |
| 586 | +{ |
577 | 587 | /* Variables for S(x) and powers of x, respectively. */ |
578 | 588 | long double cx, arg; |
579 | 589 | arg = x*x; |
@@ -646,7 +656,8 @@ long double Fresnel_Cosine_Taylor_to_Asymptotic_Long_Double(long double x){ |
646 | 656 | } |
647 | 657 | } |
648 | 658 |
|
649 | | -long double Fresnel_Cosine_While_to_Asymptotic_Long_Long_Double(long double x){ |
| 659 | +long double Fresnel_Cosine_While_to_Asymptotic_Long_Long_Double(long double x) |
| 660 | +{ |
650 | 661 | double FRESNEL_COSINE_TAYLOR_COEFFICIENTS[27] = { |
651 | 662 | FRESNEL_COSINE_TAYLOR_00, FRESNEL_COSINE_TAYLOR_01, |
652 | 663 | FRESNEL_COSINE_TAYLOR_02, FRESNEL_COSINE_TAYLOR_03, |
|
0 commit comments