| /* |
| * This file is part of FFmpeg. |
| * |
| * FFmpeg is free software; you can redistribute it and/or |
| * modify it under the terms of the GNU Lesser General Public |
| * License as published by the Free Software Foundation; either |
| * version 2.1 of the License, or (at your option) any later version. |
| * |
| * FFmpeg is distributed in the hope that it will be useful, |
| * but WITHOUT ANY WARRANTY; without even the implied warranty of |
| * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| * Lesser General Public License for more details. |
| * |
| * You should have received a copy of the GNU Lesser General Public |
| * License along with FFmpeg; if not, write to the Free Software |
| * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
| */ |
| |
| #include "libavutil/libm.h" |
| #include "libavcodec/celp_math.c" |
| |
| static inline void IsAlmostEqual(float A, float B, float epsilon) |
| { |
| float diff = fabsf(A - B); |
| float absa = fabsf(A); |
| float absb = fabsf(B); |
| float largest = (absb > absa) ? absb : absa; |
| av_assert0(diff <= largest * epsilon); |
| } |
| |
| int main(void) |
| { |
| int i; |
| const float f1[3] = {0.0, 1.1, 2.2}; |
| const float f2[3] = {3.3, 4.4, 5.5}; |
| const int16_t i1[3] = {6, 7, 8}; |
| const int16_t i2[3] = {9, 10, 11}; |
| |
| float r = ff_dot_productf(f1, f2, FF_ARRAY_ELEMS(f1)); |
| int64_t d = ff_dot_product(i1, i2, FF_ARRAY_ELEMS(i1)); |
| |
| IsAlmostEqual(16.94f, r, 0.000001f); |
| av_assert0(212 == d); |
| |
| for (i = 1024; i >= 1; i/=2) |
| av_assert0(ff_log2_q15(i) == (1<<15)*((int)log2(i))+(1<<2)); |
| |
| return 0; |
| } |