Merge commit '41578f70cf8aec8e7565fba1ca7e07f3dc46c3d2'
[ffmpeg.git] / libavutil / lls.h
1 /*
2  * linear least squares model
3  *
4  * Copyright (c) 2006 Michael Niedermayer <michaelni@gmx.at>
5  *
6  * This file is part of FFmpeg.
7  *
8  * FFmpeg is free software; you can redistribute it and/or
9  * modify it under the terms of the GNU Lesser General Public
10  * License as published by the Free Software Foundation; either
11  * version 2.1 of the License, or (at your option) any later version.
12  *
13  * FFmpeg is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
16  * Lesser General Public License for more details.
17  *
18  * You should have received a copy of the GNU Lesser General Public
19  * License along with FFmpeg; if not, write to the Free Software
20  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21  */
22
23 #ifndef AVUTIL_LLS_H
24 #define AVUTIL_LLS_H
25
26 #include "version.h"
27
28 #define MAX_VARS 32
29
30 //FIXME avoid direct access to LLSModel from outside
31
32 /**
33  * Linear least squares model.
34  */
35 typedef struct LLSModel {
36     double covariance[MAX_VARS + 1][MAX_VARS + 1];
37     double coeff[MAX_VARS][MAX_VARS];
38     double variance[MAX_VARS];
39     int indep_count;
40     /**
41      * Take the outer-product of var[] with itself, and add to the covariance matrix.
42      * @param m this context
43      * @param var training samples, starting with the value to be predicted
44      */
45     void (*update_lls)(struct LLSModel *m, double *var);
46     /**
47      * Inner product of var[] and the LPC coefs.
48      * @param m this context
49      * @param var training samples, excluding the value to be predicted
50      * @param order lpc order
51      */
52     double (*evaluate_lls)(struct LLSModel *m, double *var, int order);
53 } LLSModel;
54
55 void avpriv_init_lls(LLSModel *m, int indep_count);
56 void avpriv_solve_lls(LLSModel *m, double threshold, unsigned short min_order);
57
58 #if FF_API_LLS_PRIVATE
59 void av_init_lls(LLSModel *m, int indep_count);
60 void av_update_lls(LLSModel *m, double *param, double decay);
61 void av_solve_lls(LLSModel *m, double threshold, int min_order);
62 double av_evaluate_lls(LLSModel *m, double *param, int order);
63 #endif /* FF_API_LLS_PRIVATE */
64
65 #endif /* AVUTIL_LLS_H */