1 /*
2 * copyright (c) 2006 Michael Niedermayer <michaelni@gmx.at>
3 *
4 * This file is part of FFmpeg.
5 *
6 * FFmpeg is free software; you can redistribute it and/or
7 * modify it under the terms of the GNU Lesser General Public
8 * License as published by the Free Software Foundation; either
9 * version 2.1 of the License, or (at your option) any later version.
10 *
11 * FFmpeg is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 * Lesser General Public License for more details.
15 *
16 * You should have received a copy of the GNU Lesser General Public
17 * License along with FFmpeg; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19 */
20
22
23 #include <inttypes.h>
24
28
29 #include <stdio.h>
30
34
35
39 double d1, d2, d3;
40 int i, j;
42
43 d1= 1;
44 for(i= 0; i<10; i++){
45 d1= 1/(d1+1);
46 }
47 printf("test1 double=%d\n", (int)(d1 * (1<<24)));
48
49 sf1= one;
50 for(i= 0; i<10; i++){
52 }
53 printf(
"test1 sf =%d\n",
av_sf2int(sf1, 24));
54
55
56 for(i= 0; i<100; i++){
58 d1= i;
59 d2= i/100.0;
60 for(j= 0; j<1000; j++){
61 d1= (d1+1)*d2;
62 }
64 }
65 printf("test2 double=%d\n", (int)(d1 * (1<<24)));
66
67 for(i= 0; i<100; i++){
71 for(j= 0; j<1000; j++){
73 }
75 }
77
78 d1 = 0.0177764893;
79 d2 = 1374.40625;
80 d3 = 0.1249694824;
81 d2 += d1;
82 d3 += d2;
83 printf("test3 double: %.10lf\n", d3);
84
90 printf(
"test3 softfloat: %.10lf (0x%08x %d)\n", (
double)
av_sf2double(sf3), sf3.
mant, sf3.
exp);
91
93 printf(
"test4 softfloat: %.10lf (0x%08x %d)\n", (
double)
av_sf2double(sf1), sf1.
mant, sf1.
exp);
95 printf(
"test4 softfloat: %.10lf (0x%08x %d)\n", (
double)
av_sf2double(sf1), sf1.
mant, sf1.
exp);
96
98 printf(
"test4 softfloat: %.10lf (0x%08x %d)\n", (
double)
av_sf2double(sf1), sf1.
mant, sf1.
exp);
100 printf(
"test4 softfloat: %.10lf (0x%08x %d)\n", (
double)
av_sf2double(sf1), sf1.
mant, sf1.
exp);
101
102
105 printf(
"test5 softfloat: %.10lf (0x%08x %d)\n", (
double)
av_sf2double(sf1), sf1.
mant, sf1.
exp);
106
112 printf(
"test6 softfloat: %.10lf (0x%08x %d) %d %d\n", (
double)
av_sf2double(sf1), sf1.
mant, sf1.
exp, i, j);
113
114 for(i= -50; i<50; i++) {
116 for(j= -50; j<50; j++) {
121 printf("av_cmp_sf failed at %d %d as %X\n", i, j, c);
122 }
124 if ((i>j) != c) {
125 printf("av_gt_sf failed at %d %d as %X\n", i, j, c);
126 }
127 }
129 for(j = -50; j < 50; j++) {
134 printf("av_cmp_sf failed2 at %d %d as %X\n", i, j, c);
135 }
137 if ((i<j) != c) {
138 printf("av_gt_sf failed2 at %d %d as %X\n", i, j, c);
139 }
140 }
141 }
142
143
144 for(i= 0; i<4*36; i++){
146 double errs, errc;
147
149 errs = (double)s/ (1<<30) - sin(i*
M_PI/36);
150 errc = (double)c/ (1<<30) - cos(i*
M_PI/36);
151 if (fabs(errs) > 0.00000002 || fabs(errc) >0.001) {
152 printf(
"sincos FAIL %d %f %f %f %f\n", i, (
float)s/ (1<<30), (
float)c/ (1<<30), sin(i*
M_PI/36), cos(i*
M_PI/36));
153 }
154
155 }
156 return 0;
157
158 }
static const SoftFloat FLOAT_1374_40625
void av_log_set_level(int level)
Set the log level.
static av_const SoftFloat av_div_sf(SoftFloat a, SoftFloat b)
b has to be normalized and not zero.
static av_const double av_sf2double(SoftFloat v)
Convert a SoftFloat to a double precision float.
high precision timer, useful to profile code
static av_const SoftFloat av_normalize_sf(SoftFloat a)
static const SoftFloat FLOAT_0_017776489257
#define AV_LOG_DEBUG
Stuff which is only useful for libav* developers.
static av_unused void av_sincos_sf(int a, int *s, int *c)
Rounding-to-nearest used.
#define FFDIFFSIGN(x, y)
Comparator.
static av_const int av_cmp_sf(SoftFloat a, SoftFloat b)
Compares two SoftFloats.
static av_const int av_gt_sf(SoftFloat a, SoftFloat b)
Compares two SoftFloats.
static av_const int av_sf2int(SoftFloat v, int frac_bits)
Converts a SoftFloat to an integer.
static av_const SoftFloat av_add_sf(SoftFloat a, SoftFloat b)
common internal and external API header
static av_const SoftFloat av_mul_sf(SoftFloat a, SoftFloat b)
static av_const SoftFloat av_int2sf(int v, int frac_bits)
Converts a mantisse and exponent to a SoftFloat.
static const SoftFloat FLOAT_0_1249694824218