FFmpeg
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
float_dsp.c
Go to the documentation of this file.
1 /*
2  * Copyright 2005 Balatoni Denes
3  * Copyright 2006 Loren Merritt
4  *
5  * This file is part of FFmpeg.
6  *
7  * FFmpeg is free software; you can redistribute it and/or
8  * modify it under the terms of the GNU Lesser General Public
9  * License as published by the Free Software Foundation; either
10  * version 2.1 of the License, or (at your option) any later version.
11  *
12  * FFmpeg is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15  * Lesser General Public License for more details.
16  *
17  * You should have received a copy of the GNU Lesser General Public
18  * License along with FFmpeg; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21 
22 #include "config.h"
23 #include "attributes.h"
24 #include "float_dsp.h"
25 #include "mem.h"
26 
27 static void vector_fmul_c(float *dst, const float *src0, const float *src1,
28  int len)
29 {
30  int i;
31  for (i = 0; i < len; i++)
32  dst[i] = src0[i] * src1[i];
33 }
34 
35 static void vector_fmac_scalar_c(float *dst, const float *src, float mul,
36  int len)
37 {
38  int i;
39  for (i = 0; i < len; i++)
40  dst[i] += src[i] * mul;
41 }
42 
43 static void vector_fmul_scalar_c(float *dst, const float *src, float mul,
44  int len)
45 {
46  int i;
47  for (i = 0; i < len; i++)
48  dst[i] = src[i] * mul;
49 }
50 
51 static void vector_dmul_scalar_c(double *dst, const double *src, double mul,
52  int len)
53 {
54  int i;
55  for (i = 0; i < len; i++)
56  dst[i] = src[i] * mul;
57 }
58 
59 static void vector_fmul_window_c(float *dst, const float *src0,
60  const float *src1, const float *win, int len)
61 {
62  int i, j;
63 
64  dst += len;
65  win += len;
66  src0 += len;
67 
68  for (i = -len, j = len - 1; i < 0; i++, j--) {
69  float s0 = src0[i];
70  float s1 = src1[j];
71  float wi = win[i];
72  float wj = win[j];
73  dst[i] = s0 * wj - s1 * wi;
74  dst[j] = s0 * wi + s1 * wj;
75  }
76 }
77 
78 static void vector_fmul_add_c(float *dst, const float *src0, const float *src1,
79  const float *src2, int len){
80  int i;
81 
82  for (i = 0; i < len; i++)
83  dst[i] = src0[i] * src1[i] + src2[i];
84 }
85 
86 static void vector_fmul_reverse_c(float *dst, const float *src0,
87  const float *src1, int len)
88 {
89  int i;
90 
91  src1 += len-1;
92  for (i = 0; i < len; i++)
93  dst[i] = src0[i] * src1[-i];
94 }
95 
96 static void butterflies_float_c(float *av_restrict v1, float *av_restrict v2,
97  int len)
98 {
99  int i;
100 
101  for (i = 0; i < len; i++) {
102  float t = v1[i] - v2[i];
103  v1[i] += v2[i];
104  v2[i] = t;
105  }
106 }
107 
108 float avpriv_scalarproduct_float_c(const float *v1, const float *v2, int len)
109 {
110  float p = 0.0;
111  int i;
112 
113  for (i = 0; i < len; i++)
114  p += v1[i] * v2[i];
115 
116  return p;
117 }
118 
120 {
122  if (!fdsp)
123  return NULL;
124 
125  fdsp->vector_fmul = vector_fmul_c;
134 
135  if (ARCH_AARCH64)
137  if (ARCH_ARM)
138  ff_float_dsp_init_arm(fdsp);
139  if (ARCH_PPC)
140  ff_float_dsp_init_ppc(fdsp, bit_exact);
141  if (ARCH_X86)
142  ff_float_dsp_init_x86(fdsp);
143  if (ARCH_MIPS)
145  return fdsp;
146 }
#define NULL
Definition: coverity.c:32
void ff_float_dsp_init_x86(AVFloatDSPContext *fdsp)
static void vector_fmul_scalar_c(float *dst, const float *src, float mul, int len)
Definition: float_dsp.c:43
Memory handling functions.
float(* scalarproduct_float)(const float *v1, const float *v2, int len)
Calculate the scalar product of two vectors of floats.
Definition: float_dsp.h:159
void(* vector_fmul_reverse)(float *dst, const float *src0, const float *src1, int len)
Calculate the entry wise product of two vectors of floats, and store the result in a vector of floats...
Definition: float_dsp.h:138
void * av_mallocz(size_t size)
Allocate a memory block with alignment suitable for all memory accesses (including vectors if availab...
Definition: mem.c:222
#define src
Definition: vp8dsp.c:254
void ff_float_dsp_init_mips(AVFloatDSPContext *fdsp)
Macro definitions for various function/variable attributes.
void(* vector_fmac_scalar)(float *dst, const float *src, float mul, int len)
Multiply a vector of floats by a scalar float and add to destination vector.
Definition: float_dsp.h:54
float avpriv_scalarproduct_float_c(const float *v1, const float *v2, int len)
Return the scalar product of two vectors.
Definition: float_dsp.c:108
void(* vector_fmul_window)(float *dst, const float *src0, const float *src1, const float *win, int len)
Overlap/add with window function.
Definition: float_dsp.h:103
av_cold void ff_float_dsp_init_aarch64(AVFloatDSPContext *fdsp)
av_cold void ff_float_dsp_init_arm(AVFloatDSPContext *fdsp)
#define av_cold
Definition: attributes.h:82
static void butterflies_float_c(float *av_restrict v1, float *av_restrict v2, int len)
Definition: float_dsp.c:96
void(* vector_fmul)(float *dst, const float *src0, const float *src1, int len)
Calculate the entry wise product of two vectors of floats and store the result in a vector of floats...
Definition: float_dsp.h:38
void(* butterflies_float)(float *av_restrict v1, float *av_restrict v2, int len)
Calculate the sum and difference of two vectors of floats.
Definition: float_dsp.h:148
#define s0
Definition: regdef.h:37
static void vector_fmul_add_c(float *dst, const float *src0, const float *src1, const float *src2, int len)
Definition: float_dsp.c:78
void(* vector_dmul_scalar)(double *dst, const double *src, double mul, int len)
Multiply a vector of double by a scalar double.
Definition: float_dsp.h:84
static void vector_dmul_scalar_c(double *dst, const double *src, double mul, int len)
Definition: float_dsp.c:51
void(* vector_fmul_scalar)(float *dst, const float *src, float mul, int len)
Multiply a vector of floats by a scalar float.
Definition: float_dsp.h:69
#define src1
Definition: h264pred.c:139
void ff_float_dsp_init_ppc(AVFloatDSPContext *fdsp, int strict)
void(* vector_fmul_add)(float *dst, const float *src0, const float *src1, const float *src2, int len)
Calculate the entry wise product of two vectors of floats, add a third vector of floats and store the...
Definition: float_dsp.h:121
#define src0
Definition: h264pred.c:138
#define s1
Definition: regdef.h:38
av_cold AVFloatDSPContext * avpriv_float_dsp_alloc(int bit_exact)
Allocate a float DSP context.
Definition: float_dsp.c:119
static void vector_fmul_window_c(float *dst, const float *src0, const float *src1, const float *win, int len)
Definition: float_dsp.c:59
static void vector_fmul_c(float *dst, const float *src0, const float *src1, int len)
Definition: float_dsp.c:27
static void vector_fmac_scalar_c(float *dst, const float *src, float mul, int len)
Definition: float_dsp.c:35
int len
static void vector_fmul_reverse_c(float *dst, const float *src0, const float *src1, int len)
Definition: float_dsp.c:86