/freebsd/contrib/arm-optimized-routines/math/tools/ |
H A D | v_log.sollya | 6 deg = 6; // poly degree 12 deg = deg-1; // because of /x 20 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 25 for i from 1 to deg do { 34 for i from 0 to deg do coeff(poly,i);
|
H A D | v_sin.sollya | 6 deg = 15; // polynomial degree 13 deg = deg-1; 22 return remez(f(x)-poly(x), deg-d, [a;b], x^d, 1e-10); 27 for i from 1 to deg/2 do { 36 for i from 0 to deg do coeff(poly,i);
|
H A D | sin.sollya | 6 deg = 7; // polynomial degree 13 deg = deg-1; 22 return remez(f(x)-poly(x), deg-d, [a;b], x^d, 1e-10); 27 for i from 1 to deg/2 do { 37 for i from 0 to deg do coeff(poly,i);
|
H A D | log.sollya | 6 deg = 12; // poly degree 13 deg = deg-1; // because of /x 21 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 26 for i from 1 to deg do { 35 for i from 0 to deg do coeff(poly,i);
|
H A D | v_log2f.sollya | 6 deg = 9; // poly degree 15 deg = deg-1; // because of /x 24 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 29 for i from 1 to deg do { 38 for i from 0 to deg do coeff(poly,i);
|
H A D | v_log10.sollya | 6 deg = 6; // poly degree 12 deg = deg-1; // because of /x 20 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 25 for i from 1 to deg do { 38 for i from 0 to deg do double(coeff(poly,i));
|
H A D | log10f.sollya | 8 deg = 4; // poly degree 15 deg = deg-1; // because of /x 23 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 28 for i from 1 to deg do { 37 for i from 0 to deg do double(coeff(poly,i));
|
H A D | log2.sollya | 6 deg = 11; // poly degree 17 deg = deg-1; // because of /x 26 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 31 for i from 1 to deg do { 42 for i from 0 to deg do coeff(poly,i);
|
H A D | v_log10f.sollya | 6 deg = 9; // poly degree 19 deg = deg-1; // because of /x 28 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 33 for i from 1 to deg do { 42 for i from 0 to deg do single(coeff(poly,i));
|
H A D | tanpi.sollya | 11 if (method == 0) then { deg = 5; } 12 else if (method == 1) then { deg = 3; }; 14 if (method == 0) then { deg = 13; } 15 else if (method == 1) then { deg = 8; }; 31 poly = fpminimax(f, deg, [|dtype ...|], [a*a;b*b]); 40 for i from 0 to deg do coeff(poly, i);
|
H A D | log10.sollya | 6 deg = 6; // poly degree 17 deg = deg-1; // because of /x 26 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 31 for i from 1 to deg do { 41 for i from 0 to deg do coeff(poly,i);
|
H A D | exp2.sollya | 7 deg = 3; // poly degree 13 //deg = 5; // poly degree 24 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 28 return remez(f(x) - poly(x), deg-d, [a;b], x^d, 1e-10); 33 for i from 1 to deg do { 48 for i from 0 to deg do coeff(poly,i);
|
H A D | exp10.sollya | 7 deg = 5; // poly degree 14 //deg = 4; // poly degree - bump to 5 for ~1 ULP 27 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 31 return remez(f(x) - poly(x), deg-d, [a;b], x^d, 1e-10); 36 for i from 1 to deg do { 47 for i from 0 to deg do coeff(poly,i);
|
H A D | tanf.sollya | 9 deg = 5; // poly degree 13 // deg = 3; // poly degree 43 poly = fpminimax(f, deg, [|dtype ...|], [a*a;b*b]); 56 return fpminimax(func - poly / x^-(deg-d), 0, [|dtype|], [a;b], absolute, floating); 59 poly = fpminimax(f, [|0,...,deg|], [|dtype ...|], [a;b], absolute, floating); 70 for i from 0 to deg do coeff(poly, i);
|
H A D | asinhf.sollya | 6 deg = 9; 16 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 20 for i from 2 to deg do { 29 for i from 2 to deg do coeff(poly,i);
|
H A D | atan.sollya | 12 deg = 20; 13 mons = [|1,...,deg|]; 14 for i from 0 to deg-1 do mons[i] = mons[i] * 2 + 1; 23 for i from 0 to deg-1 do coeff(poly,mons[i]);
|
H A D | log1p.sollya | 6 deg = 20; 16 return remez(1 - poly(x)/f(x), deg-d, [a;b], x^d/f(x), 1e-10); 20 for i from 2 to deg do { 28 for i from 2 to deg do coeff(poly,i);
|
H A D | asinf.sollya | 12 deg = 4; 23 return remez(1 - poly(x) / forward(f, x), deg - d, [a;b], x^d/forward(f, x), 1e-16); 27 for i from 0 to deg do { 36 for i from 0 to deg do print(coeff(poly, i));
|
H A D | cos.sollya | 6 deg = 8; // polynomial degree 16 return remez(f(x)-poly(x), deg-d, [a;b], x^d, 1e-10); 21 for i from 1 to deg/2 do { 31 for i from 0 to deg do coeff(poly,i);
|
H A D | v_exp.sollya | 6 deg = 4; // poly degree 15 return remez(exp(x)-poly(x), deg-d, [a;b], x^d, 1e-10); 20 for i from 2 to deg do { 30 for i from 0 to deg do coeff(poly,i);
|
H A D | sincosf.sollya | 8 deg = 8; // polynomial degree 18 return remez(f(x)-poly(x), deg-d, [a;b], x^d, 1e-10); 23 for i from 1 to deg/2 do { 33 for i from 0 to deg do coeff(poly,i);
|
H A D | sinpi.sollya | 6 deg = 19; // polynomial degree 19 return remez(f(x)-poly(x), deg-d, [a;b], x^d, 1e-10); 24 for i from 0 to (deg-1)/2 do { 33 for i from 0 to deg do coeff(poly,i);
|
H A D | sincos.sollya | 8 deg = 14; // polynomial degree 18 return remez(f(x)-poly(x), deg-d, [a;b], x^d, 1e-10); 23 for i from 1 to deg/2 do { 33 for i from 0 to deg do coeff(poly,i);
|
/freebsd/contrib/arm-optimized-routines/math/ |
H A D | erf_data.c | 43 /* Generated using Sollya::remez(f(c*x+d), deg, [(a-d)/c;(b-d)/c], 1, 1e-16), [|D ...|] with deg=15… 54 /* Generated using Sollya::remez(f(c*x+d), deg, [(a-d)/c;(b-d)/c], 1, 1e-16), [|D ...|] with deg=17… 65 /* Generated using Sollya::remez(f(c*x+d), deg, [(a-d)/c;(b-d)/c], 1, 1e-16), [|D ...|] with deg=13… 75 /* Generated using Sollya::remez(f(c*x+d), deg, [(a-d)/c;(b-d)/c], 1, 1e-16), [|D ...|] with deg=16…
|
/freebsd/usr.bin/pom/ |
H A D | pom.c | 211 dtor(double deg) in dtor() argument 214 return(deg * PI / 180); in dtor() 219 * adjust value so 0 <= deg <= 360 222 adj360(double *deg) in adj360() argument 226 if (*deg < 0) in adj360() 227 *deg += 360; in adj360() 228 else if (*deg > 360) in adj360() 229 *deg -= 360; in adj360()
|