+++ /dev/null
-/*M///////////////////////////////////////////////////////////////////////////////////////
-//
-// IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
-//
-// By downloading, copying, installing or using the software you agree to this license.
-// If you do not agree to this license, do not download, install,
-// copy or use the software.
-//
-//
-// Intel License Agreement
-// For Open Source Computer Vision Library
-//
-// Copyright (C) 2000, Intel Corporation, all rights reserved.
-// Third party copyrights are property of their respective owners.
-//
-// Redistribution and use in source and binary forms, with or without modification,
-// are permitted provided that the following conditions are met:
-//
-// * Redistribution's of source code must retain the above copyright notice,
-// this list of conditions and the following disclaimer.
-//
-// * Redistribution's in binary form must reproduce the above copyright notice,
-// this list of conditions and the following disclaimer in the documentation
-// and/or other materials provided with the distribution.
-//
-// * The name of Intel Corporation may not be used to endorse or promote products
-// derived from this software without specific prior written permission.
-//
-// This software is provided by the copyright holders and contributors "as is" and
-// any express or implied warranties, including, but not limited to, the implied
-// warranties of merchantability and fitness for a particular purpose are disclaimed.
-// In no event shall the Intel Corporation or contributors be liable for any direct,
-// indirect, incidental, special, exemplary, or consequential damages
-// (including, but not limited to, procurement of substitute goods or services;
-// loss of use, data, or profits; or business interruption) however caused
-// and on any theory of liability, whether in contract, strict liability,
-// or tort (including negligence or otherwise) arising in any way out of
-// the use of this software, even if advised of the possibility of such damage.
-//
-//M*/
-
-#include "_cv.h"
-
-#define ICV_DEF_ACC_FUNC( name, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, \
-name,( const srctype *src, int srcstep, dsttype *dst, \
- int dststep, CvSize size ), (src, srcstep, dst, dststep, size )) \
- \
-{ \
- srcstep /= sizeof(src[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src += srcstep, dst += dststep ) \
- { \
- int x; \
- for( x = 0; x <= size.width - 4; x += 4 ) \
- { \
- dsttype t0 = dst[x] + cvtmacro(src[x]); \
- dsttype t1 = dst[x + 1] + cvtmacro(src[x + 1]); \
- dst[x] = t0; dst[x + 1] = t1; \
- \
- t0 = dst[x + 2] + cvtmacro(src[x + 2]); \
- t1 = dst[x + 3] + cvtmacro(src[x + 3]); \
- dst[x + 2] = t0; dst[x + 3] = t1; \
- } \
- \
- for( ; x < size.width; x++ ) \
- dst[x] += cvtmacro(src[x]); \
- } \
- \
- return CV_OK; \
-}
-
-
-ICV_DEF_ACC_FUNC( icvAdd_8u32f_C1IR, uchar, float, CV_8TO32F )
-ICV_DEF_ACC_FUNC( icvAdd_32f_C1IR, float, float, CV_NOP )
-ICV_DEF_ACC_FUNC( icvAddSquare_8u32f_C1IR, uchar, float, CV_8TO32F_SQR )
-ICV_DEF_ACC_FUNC( icvAddSquare_32f_C1IR, float, float, CV_SQR )
-
-
-#define ICV_DEF_ACCPROD_FUNC( flavor, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, icvAddProduct_##flavor##_C1IR, \
-( const srctype *src1, int step1, const srctype *src2, int step2, \
- dsttype *dst, int dststep, CvSize size ), \
- (src1, step1, src2, step2, dst, dststep, size) ) \
-{ \
- step1 /= sizeof(src1[0]); \
- step2 /= sizeof(src2[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src1 += step1, src2 += step2, dst += dststep ) \
- { \
- int x; \
- for( x = 0; x <= size.width - 4; x += 4 ) \
- { \
- dsttype t0 = dst[x] + cvtmacro(src1[x])*cvtmacro(src2[x]); \
- dsttype t1 = dst[x+1] + cvtmacro(src1[x+1])*cvtmacro(src2[x+1]);\
- dst[x] = t0; dst[x + 1] = t1; \
- \
- t0 = dst[x + 2] + cvtmacro(src1[x + 2])*cvtmacro(src2[x + 2]); \
- t1 = dst[x + 3] + cvtmacro(src1[x + 3])*cvtmacro(src2[x + 3]); \
- dst[x + 2] = t0; dst[x + 3] = t1; \
- } \
- \
- for( ; x < size.width; x++ ) \
- dst[x] += cvtmacro(src1[x])*cvtmacro(src2[x]); \
- } \
- \
- return CV_OK; \
-}
-
-
-ICV_DEF_ACCPROD_FUNC( 8u32f, uchar, float, CV_8TO32F )
-ICV_DEF_ACCPROD_FUNC( 32f, float, float, CV_NOP )
-
-
-#define ICV_DEF_ACCWEIGHT_FUNC( flavor, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, icvAddWeighted_##flavor##_C1IR, \
-( const srctype *src, int srcstep, dsttype *dst, int dststep, \
- CvSize size, dsttype alpha ), (src, srcstep, dst, dststep, size, alpha) )\
-{ \
- dsttype beta = (dsttype)(1 - alpha); \
- srcstep /= sizeof(src[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src += srcstep, dst += dststep ) \
- { \
- int x; \
- for( x = 0; x <= size.width - 4; x += 4 ) \
- { \
- dsttype t0 = dst[x]*beta + cvtmacro(src[x])*alpha; \
- dsttype t1 = dst[x+1]*beta + cvtmacro(src[x+1])*alpha; \
- dst[x] = t0; dst[x + 1] = t1; \
- \
- t0 = dst[x + 2]*beta + cvtmacro(src[x + 2])*alpha; \
- t1 = dst[x + 3]*beta + cvtmacro(src[x + 3])*alpha; \
- dst[x + 2] = t0; dst[x + 3] = t1; \
- } \
- \
- for( ; x < size.width; x++ ) \
- dst[x] = dst[x]*beta + cvtmacro(src[x])*alpha; \
- } \
- \
- return CV_OK; \
-}
-
-
-ICV_DEF_ACCWEIGHT_FUNC( 8u32f, uchar, float, CV_8TO32F )
-ICV_DEF_ACCWEIGHT_FUNC( 32f, float, float, CV_NOP )
-
-
-#define ICV_DEF_ACCMASK_FUNC_C1( name, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, \
-name,( const srctype *src, int srcstep, const uchar* mask, int maskstep,\
- dsttype *dst, int dststep, CvSize size ), \
- (src, srcstep, mask, maskstep, dst, dststep, size )) \
-{ \
- srcstep /= sizeof(src[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src += srcstep, \
- dst += dststep, mask += maskstep ) \
- { \
- int x; \
- for( x = 0; x <= size.width - 2; x += 2 ) \
- { \
- if( mask[x] ) \
- dst[x] += cvtmacro(src[x]); \
- if( mask[x+1] ) \
- dst[x+1] += cvtmacro(src[x+1]); \
- } \
- \
- for( ; x < size.width; x++ ) \
- if( mask[x] ) \
- dst[x] += cvtmacro(src[x]); \
- } \
- \
- return CV_OK; \
-}
-
-
-ICV_DEF_ACCMASK_FUNC_C1( icvAdd_8u32f_C1IMR, uchar, float, CV_8TO32F )
-ICV_DEF_ACCMASK_FUNC_C1( icvAdd_32f_C1IMR, float, float, CV_NOP )
-ICV_DEF_ACCMASK_FUNC_C1( icvAddSquare_8u32f_C1IMR, uchar, float, CV_8TO32F_SQR )
-ICV_DEF_ACCMASK_FUNC_C1( icvAddSquare_32f_C1IMR, float, float, CV_SQR )
-
-
-#define ICV_DEF_ACCPRODUCTMASK_FUNC_C1( flavor, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, icvAddProduct_##flavor##_C1IMR, \
-( const srctype *src1, int step1, const srctype* src2, int step2, \
- const uchar* mask, int maskstep, dsttype *dst, int dststep, CvSize size ),\
- (src1, step1, src2, step2, mask, maskstep, dst, dststep, size )) \
-{ \
- step1 /= sizeof(src1[0]); \
- step2 /= sizeof(src2[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src1 += step1, src2 += step2, \
- dst += dststep, mask += maskstep ) \
- { \
- int x; \
- for( x = 0; x <= size.width - 2; x += 2 ) \
- { \
- if( mask[x] ) \
- dst[x] += cvtmacro(src1[x])*cvtmacro(src2[x]); \
- if( mask[x+1] ) \
- dst[x+1] += cvtmacro(src1[x+1])*cvtmacro(src2[x+1]); \
- } \
- \
- for( ; x < size.width; x++ ) \
- if( mask[x] ) \
- dst[x] += cvtmacro(src1[x])*cvtmacro(src2[x]); \
- } \
- \
- return CV_OK; \
-}
-
-
-ICV_DEF_ACCPRODUCTMASK_FUNC_C1( 8u32f, uchar, float, CV_8TO32F )
-ICV_DEF_ACCPRODUCTMASK_FUNC_C1( 32f, float, float, CV_NOP )
-
-#define ICV_DEF_ACCWEIGHTMASK_FUNC_C1( flavor, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, icvAddWeighted_##flavor##_C1IMR, \
-( const srctype *src, int srcstep, const uchar* mask, int maskstep, \
- dsttype *dst, int dststep, CvSize size, dsttype alpha ), \
- (src, srcstep, mask, maskstep, dst, dststep, size, alpha )) \
-{ \
- dsttype beta = (dsttype)(1 - alpha); \
- srcstep /= sizeof(src[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src += srcstep, \
- dst += dststep, mask += maskstep ) \
- { \
- int x; \
- for( x = 0; x <= size.width - 2; x += 2 ) \
- { \
- if( mask[x] ) \
- dst[x] = dst[x]*beta + cvtmacro(src[x])*alpha; \
- if( mask[x+1] ) \
- dst[x+1] = dst[x+1]*beta + cvtmacro(src[x+1])*alpha; \
- } \
- \
- for( ; x < size.width; x++ ) \
- if( mask[x] ) \
- dst[x] = dst[x]*beta + cvtmacro(src[x])*alpha; \
- } \
- \
- return CV_OK; \
-}
-
-ICV_DEF_ACCWEIGHTMASK_FUNC_C1( 8u32f, uchar, float, CV_8TO32F )
-ICV_DEF_ACCWEIGHTMASK_FUNC_C1( 32f, float, float, CV_NOP )
-
-
-#define ICV_DEF_ACCMASK_FUNC_C3( name, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, \
-name,( const srctype *src, int srcstep, const uchar* mask, int maskstep,\
- dsttype *dst, int dststep, CvSize size ), \
- (src, srcstep, mask, maskstep, dst, dststep, size )) \
-{ \
- srcstep /= sizeof(src[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src += srcstep, \
- dst += dststep, mask += maskstep ) \
- { \
- int x; \
- for( x = 0; x < size.width; x++ ) \
- if( mask[x] ) \
- { \
- dsttype t0, t1, t2; \
- t0 = dst[x*3] + cvtmacro(src[x*3]); \
- t1 = dst[x*3+1] + cvtmacro(src[x*3+1]); \
- t2 = dst[x*3+2] + cvtmacro(src[x*3+2]); \
- dst[x*3] = t0; \
- dst[x*3+1] = t1; \
- dst[x*3+2] = t2; \
- } \
- } \
- \
- return CV_OK; \
-}
-
-
-ICV_DEF_ACCMASK_FUNC_C3( icvAdd_8u32f_C3IMR, uchar, float, CV_8TO32F )
-ICV_DEF_ACCMASK_FUNC_C3( icvAdd_32f_C3IMR, float, float, CV_NOP )
-ICV_DEF_ACCMASK_FUNC_C3( icvAddSquare_8u32f_C3IMR, uchar, float, CV_8TO32F_SQR )
-ICV_DEF_ACCMASK_FUNC_C3( icvAddSquare_32f_C3IMR, float, float, CV_SQR )
-
-
-#define ICV_DEF_ACCPRODUCTMASK_FUNC_C3( flavor, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, icvAddProduct_##flavor##_C3IMR, \
-( const srctype *src1, int step1, const srctype* src2, int step2, \
- const uchar* mask, int maskstep, dsttype *dst, int dststep, CvSize size ),\
- (src1, step1, src2, step2, mask, maskstep, dst, dststep, size )) \
-{ \
- step1 /= sizeof(src1[0]); \
- step2 /= sizeof(src2[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src1 += step1, src2 += step2, \
- dst += dststep, mask += maskstep ) \
- { \
- int x; \
- for( x = 0; x < size.width; x++ ) \
- if( mask[x] ) \
- { \
- dsttype t0, t1, t2; \
- t0 = dst[x*3]+cvtmacro(src1[x*3])*cvtmacro(src2[x*3]); \
- t1 = dst[x*3+1]+cvtmacro(src1[x*3+1])*cvtmacro(src2[x*3+1]);\
- t2 = dst[x*3+2]+cvtmacro(src1[x*3+2])*cvtmacro(src2[x*3+2]);\
- dst[x*3] = t0; \
- dst[x*3+1] = t1; \
- dst[x*3+2] = t2; \
- } \
- } \
- \
- return CV_OK; \
-}
-
-
-ICV_DEF_ACCPRODUCTMASK_FUNC_C3( 8u32f, uchar, float, CV_8TO32F )
-ICV_DEF_ACCPRODUCTMASK_FUNC_C3( 32f, float, float, CV_NOP )
-
-
-#define ICV_DEF_ACCWEIGHTMASK_FUNC_C3( flavor, srctype, dsttype, cvtmacro ) \
-IPCVAPI_IMPL( CvStatus, icvAddWeighted_##flavor##_C3IMR, \
-( const srctype *src, int srcstep, const uchar* mask, int maskstep, \
- dsttype *dst, int dststep, CvSize size, dsttype alpha ), \
- (src, srcstep, mask, maskstep, dst, dststep, size, alpha )) \
-{ \
- dsttype beta = (dsttype)(1 - alpha); \
- srcstep /= sizeof(src[0]); \
- dststep /= sizeof(dst[0]); \
- \
- for( ; size.height--; src += srcstep, \
- dst += dststep, mask += maskstep ) \
- { \
- int x; \
- for( x = 0; x < size.width; x++ ) \
- if( mask[x] ) \
- { \
- dsttype t0, t1, t2; \
- t0 = dst[x*3]*beta + cvtmacro(src[x*3])*alpha; \
- t1 = dst[x*3+1]*beta + cvtmacro(src[x*3+1])*alpha; \
- t2 = dst[x*3+2]*beta + cvtmacro(src[x*3+2])*alpha; \
- dst[x*3] = t0; \
- dst[x*3+1] = t1; \
- dst[x*3+2] = t2; \
- } \
- } \
- \
- return CV_OK; \
-}
-
-ICV_DEF_ACCWEIGHTMASK_FUNC_C3( 8u32f, uchar, float, CV_8TO32F )
-ICV_DEF_ACCWEIGHTMASK_FUNC_C3( 32f, float, float, CV_NOP )
-
-
-#define ICV_DEF_INIT_ACC_TAB( FUNCNAME ) \
-static void icvInit##FUNCNAME##Table( CvFuncTable* tab, CvBigFuncTable* masktab ) \
-{ \
- tab->fn_2d[CV_8U] = (void*)icv##FUNCNAME##_8u32f_C1IR; \
- tab->fn_2d[CV_32F] = (void*)icv##FUNCNAME##_32f_C1IR; \
- \
- masktab->fn_2d[CV_8UC1] = (void*)icv##FUNCNAME##_8u32f_C1IMR; \
- masktab->fn_2d[CV_32FC1] = (void*)icv##FUNCNAME##_32f_C1IMR; \
- \
- masktab->fn_2d[CV_8UC3] = (void*)icv##FUNCNAME##_8u32f_C3IMR; \
- masktab->fn_2d[CV_32FC3] = (void*)icv##FUNCNAME##_32f_C3IMR; \
-}
-
-
-ICV_DEF_INIT_ACC_TAB( Add )
-ICV_DEF_INIT_ACC_TAB( AddSquare )
-ICV_DEF_INIT_ACC_TAB( AddProduct )
-ICV_DEF_INIT_ACC_TAB( AddWeighted )
-
-
-CV_IMPL void
-cvAcc( const void* arr, void* sumarr, const void* maskarr )
-{
- static CvFuncTable acc_tab;
- static CvBigFuncTable accmask_tab;
- static int inittab = 0;
-
- CV_FUNCNAME( "cvAcc" );
-
- __BEGIN__;
-
- int type, sumdepth;
- int mat_step, sum_step, mask_step = 0;
- CvSize size;
- CvMat stub, *mat = (CvMat*)arr;
- CvMat sumstub, *sum = (CvMat*)sumarr;
- CvMat maskstub, *mask = (CvMat*)maskarr;
-
- if( !inittab )
- {
- icvInitAddTable( &acc_tab, &accmask_tab );
- inittab = 1;
- }
-
- if( !CV_IS_MAT( mat ) || !CV_IS_MAT( sum ))
- {
- int coi1 = 0, coi2 = 0;
- CV_CALL( mat = cvGetMat( mat, &stub, &coi1 ));
- CV_CALL( sum = cvGetMat( sum, &sumstub, &coi2 ));
- if( coi1 + coi2 != 0 )
- CV_ERROR( CV_BadCOI, "" );
- }
-
- if( CV_MAT_DEPTH( sum->type ) != CV_32F )
- CV_ERROR( CV_BadDepth, "" );
-
- if( !CV_ARE_CNS_EQ( mat, sum ))
- CV_ERROR( CV_StsUnmatchedFormats, "" );
-
- sumdepth = CV_MAT_DEPTH( sum->type );
- if( sumdepth != CV_32F && (maskarr != 0 || sumdepth != CV_64F))
- CV_ERROR( CV_BadDepth, "Bad accumulator type" );
-
- if( !CV_ARE_SIZES_EQ( mat, sum ))
- CV_ERROR( CV_StsUnmatchedSizes, "" );
-
- size = cvGetMatSize( mat );
- type = CV_MAT_TYPE( mat->type );
-
- mat_step = mat->step;
- sum_step = sum->step;
-
- if( !mask )
- {
- CvFunc2D_2A func=(CvFunc2D_2A)acc_tab.fn_2d[CV_MAT_DEPTH(type)];
-
- if( !func )
- CV_ERROR( CV_StsUnsupportedFormat, "Unsupported type combination" );
-
- size.width *= CV_MAT_CN(type);
- if( CV_IS_MAT_CONT( mat->type & sum->type ))
- {
- size.width *= size.height;
- mat_step = sum_step = CV_STUB_STEP;
- size.height = 1;
- }
-
- IPPI_CALL( func( mat->data.ptr, mat_step, sum->data.ptr, sum_step, size ));
- }
- else
- {
- CvFunc2D_3A func = (CvFunc2D_3A)accmask_tab.fn_2d[type];
-
- if( !func )
- CV_ERROR( CV_StsUnsupportedFormat, "" );
-
- CV_CALL( mask = cvGetMat( mask, &maskstub ));
-
- if( !CV_IS_MASK_ARR( mask ))
- CV_ERROR( CV_StsBadMask, "" );
-
- if( !CV_ARE_SIZES_EQ( mat, mask ))
- CV_ERROR( CV_StsUnmatchedSizes, "" );
-
- mask_step = mask->step;
-
- if( CV_IS_MAT_CONT( mat->type & sum->type & mask->type ))
- {
- size.width *= size.height;
- mat_step = sum_step = mask_step = CV_STUB_STEP;
- size.height = 1;
- }
-
- IPPI_CALL( func( mat->data.ptr, mat_step, mask->data.ptr, mask_step,
- sum->data.ptr, sum_step, size ));
- }
-
- __END__;
-}
-
-
-CV_IMPL void
-cvSquareAcc( const void* arr, void* sq_sum, const void* maskarr )
-{
- static CvFuncTable acc_tab;
- static CvBigFuncTable accmask_tab;
- static int inittab = 0;
-
- CV_FUNCNAME( "cvSquareAcc" );
-
- __BEGIN__;
-
- int coi1, coi2;
- int type;
- int mat_step, sum_step, mask_step = 0;
- CvSize size;
- CvMat stub, *mat = (CvMat*)arr;
- CvMat sumstub, *sum = (CvMat*)sq_sum;
- CvMat maskstub, *mask = (CvMat*)maskarr;
-
- if( !inittab )
- {
- icvInitAddSquareTable( &acc_tab, &accmask_tab );
- inittab = 1;
- }
-
- CV_CALL( mat = cvGetMat( mat, &stub, &coi1 ));
- CV_CALL( sum = cvGetMat( sum, &sumstub, &coi2 ));
-
- if( coi1 != 0 || coi2 != 0 )
- CV_ERROR( CV_BadCOI, "" );
-
- if( !CV_ARE_CNS_EQ( mat, sum ))
- CV_ERROR( CV_StsUnmatchedFormats, "" );
-
- if( CV_MAT_DEPTH( sum->type ) != CV_32F )
- CV_ERROR( CV_BadDepth, "" );
-
- if( !CV_ARE_SIZES_EQ( mat, sum ))
- CV_ERROR( CV_StsUnmatchedSizes, "" );
-
- size = cvGetMatSize( mat );
- type = CV_MAT_TYPE( mat->type );
-
- mat_step = mat->step;
- sum_step = sum->step;
-
- if( !mask )
- {
- CvFunc2D_2A func = (CvFunc2D_2A)acc_tab.fn_2d[CV_MAT_DEPTH(type)];
-
- if( !func )
- CV_ERROR( CV_StsUnsupportedFormat, "" );
-
- size.width *= CV_MAT_CN(type);
-
- if( CV_IS_MAT_CONT( mat->type & sum->type ))
- {
- size.width *= size.height;
- mat_step = sum_step = CV_STUB_STEP;;
- size.height = 1;
- }
-
- IPPI_CALL( func( mat->data.ptr, mat_step, sum->data.ptr, sum_step, size ));
- }
- else
- {
- CvFunc2D_3A func = (CvFunc2D_3A)accmask_tab.fn_2d[type];
-
- if( !func )
- CV_ERROR( CV_StsUnsupportedFormat, "" );
-
- CV_CALL( mask = cvGetMat( mask, &maskstub ));
-
- if( !CV_IS_MASK_ARR( mask ))
- CV_ERROR( CV_StsBadMask, "" );
-
- if( !CV_ARE_SIZES_EQ( mat, mask ))
- CV_ERROR( CV_StsUnmatchedSizes, "" );
-
- mask_step = mask->step;
-
- if( CV_IS_MAT_CONT( mat->type & sum->type & mask->type ))
- {
- size.width *= size.height;
- mat_step = sum_step = mask_step = CV_STUB_STEP;
- size.height = 1;
- }
-
- IPPI_CALL( func( mat->data.ptr, mat_step, mask->data.ptr, mask_step,
- sum->data.ptr, sum_step, size ));
- }
-
- __END__;
-}
-
-
-CV_IMPL void
-cvMultiplyAcc( const void* arrA, const void* arrB,
- void* acc, const void* maskarr )
-{
- static CvFuncTable acc_tab;
- static CvBigFuncTable accmask_tab;
- static int inittab = 0;
-
- CV_FUNCNAME( "cvMultiplyAcc" );
-
- __BEGIN__;
-
- int coi1, coi2, coi3;
- int type;
- int mat1_step, mat2_step, sum_step, mask_step = 0;
- CvSize size;
- CvMat stub1, *mat1 = (CvMat*)arrA;
- CvMat stub2, *mat2 = (CvMat*)arrB;
- CvMat sumstub, *sum = (CvMat*)acc;
- CvMat maskstub, *mask = (CvMat*)maskarr;
-
- if( !inittab )
- {
- icvInitAddProductTable( &acc_tab, &accmask_tab );
- inittab = 1;
- }
-
- CV_CALL( mat1 = cvGetMat( mat1, &stub1, &coi1 ));
- CV_CALL( mat2 = cvGetMat( mat2, &stub2, &coi2 ));
- CV_CALL( sum = cvGetMat( sum, &sumstub, &coi3 ));
-
- if( coi1 != 0 || coi2 != 0 || coi3 != 0 )
- CV_ERROR( CV_BadCOI, "" );
-
- if( !CV_ARE_CNS_EQ( mat1, mat2 ) || !CV_ARE_CNS_EQ( mat1, sum ))
- CV_ERROR( CV_StsUnmatchedFormats, "" );
-
- if( CV_MAT_DEPTH( sum->type ) != CV_32F )
- CV_ERROR( CV_BadDepth, "" );
-
- if( !CV_ARE_SIZES_EQ( mat1, sum ) || !CV_ARE_SIZES_EQ( mat2, sum ))
- CV_ERROR( CV_StsUnmatchedSizes, "" );
-
- size = cvGetMatSize( mat1 );
- type = CV_MAT_TYPE( mat1->type );
-
- mat1_step = mat1->step;
- mat2_step = mat2->step;
- sum_step = sum->step;
-
- if( !mask )
- {
- CvFunc2D_3A func = (CvFunc2D_3A)acc_tab.fn_2d[CV_MAT_DEPTH(type)];
-
- if( !func )
- CV_ERROR( CV_StsUnsupportedFormat, "" );
-
- size.width *= CV_MAT_CN(type);
-
- if( CV_IS_MAT_CONT( mat1->type & mat2->type & sum->type ))
- {
- size.width *= size.height;
- mat1_step = mat2_step = sum_step = CV_STUB_STEP;
- size.height = 1;
- }
-
- IPPI_CALL( func( mat1->data.ptr, mat1_step, mat2->data.ptr, mat2_step,
- sum->data.ptr, sum_step, size ));
- }
- else
- {
- CvFunc2D_4A func = (CvFunc2D_4A)accmask_tab.fn_2d[type];
-
- if( !func )
- CV_ERROR( CV_StsUnsupportedFormat, "" );
-
- CV_CALL( mask = cvGetMat( mask, &maskstub ));
-
- if( !CV_IS_MASK_ARR( mask ))
- CV_ERROR( CV_StsBadMask, "" );
-
- if( !CV_ARE_SIZES_EQ( mat1, mask ))
- CV_ERROR( CV_StsUnmatchedSizes, "" );
-
- mask_step = mask->step;
-
- if( CV_IS_MAT_CONT( mat1->type & mat2->type & sum->type & mask->type ))
- {
- size.width *= size.height;
- mat1_step = mat2_step = sum_step = mask_step = CV_STUB_STEP;
- size.height = 1;
- }
-
- IPPI_CALL( func( mat1->data.ptr, mat1_step, mat2->data.ptr, mat2_step,
- mask->data.ptr, mask_step,
- sum->data.ptr, sum_step, size ));
- }
-
- __END__;
-}
-
-
-typedef CvStatus (CV_STDCALL *CvAddWeightedFunc)( const void* src, int srcstep,
- void* dst, int dststep,
- CvSize size, float alpha );
-
-typedef CvStatus (CV_STDCALL *CvAddWeightedMaskFunc)( const void* src, int srcstep,
- void* dst, int dststep,
- const void* mask, int maskstep,
- CvSize size, float alpha );
-
-CV_IMPL void
-cvRunningAvg( const void* arrY, void* arrU,
- double alpha, const void* maskarr )
-{
- static CvFuncTable acc_tab;
- static CvBigFuncTable accmask_tab;
- static int inittab = 0;
-
- CV_FUNCNAME( "cvRunningAvg" );
-
- __BEGIN__;
-
- int coi1, coi2;
- int type;
- int mat_step, sum_step, mask_step = 0;
- CvSize size;
- CvMat stub, *mat = (CvMat*)arrY;
- CvMat sumstub, *sum = (CvMat*)arrU;
- CvMat maskstub, *mask = (CvMat*)maskarr;
-
- if( !inittab )
- {
- icvInitAddWeightedTable( &acc_tab, &accmask_tab );
- inittab = 1;
- }
-
- CV_CALL( mat = cvGetMat( mat, &stub, &coi1 ));
- CV_CALL( sum = cvGetMat( sum, &sumstub, &coi2 ));
-
- if( coi1 != 0 || coi2 != 0 )
- CV_ERROR( CV_BadCOI, "" );
-
- if( !CV_ARE_CNS_EQ( mat, sum ))
- CV_ERROR( CV_StsUnmatchedFormats, "" );
-
- if( CV_MAT_DEPTH( sum->type ) != CV_32F )
- CV_ERROR( CV_BadDepth, "" );
-
- if( !CV_ARE_SIZES_EQ( mat, sum ))
- CV_ERROR( CV_StsUnmatchedSizes, "" );
-
- size = cvGetMatSize( mat );
- type = CV_MAT_TYPE( mat->type );
-
- mat_step = mat->step;
- sum_step = sum->step;
-
- if( !mask )
- {
- CvAddWeightedFunc func = (CvAddWeightedFunc)acc_tab.fn_2d[CV_MAT_DEPTH(type)];
-
- if( !func )
- CV_ERROR( CV_StsUnsupportedFormat, "" );
-
- size.width *= CV_MAT_CN(type);
- if( CV_IS_MAT_CONT( mat->type & sum->type ))
- {
- size.width *= size.height;
- mat_step = sum_step = CV_STUB_STEP;
- size.height = 1;
- }
-
- IPPI_CALL( func( mat->data.ptr, mat_step,
- sum->data.ptr, sum_step, size, (float)alpha ));
- }
- else
- {
- CvAddWeightedMaskFunc func = (CvAddWeightedMaskFunc)accmask_tab.fn_2d[type];
-
- if( !func )
- CV_ERROR( CV_StsUnsupportedFormat, "" );
-
- CV_CALL( mask = cvGetMat( mask, &maskstub ));
-
- if( !CV_IS_MASK_ARR( mask ))
- CV_ERROR( CV_StsBadMask, "" );
-
- if( !CV_ARE_SIZES_EQ( mat, mask ))
- CV_ERROR( CV_StsUnmatchedSizes, "" );
-
- mask_step = mask->step;
-
- if( CV_IS_MAT_CONT( mat->type & sum->type & mask->type ))
- {
- size.width *= size.height;
- mat_step = sum_step = mask_step = CV_STUB_STEP;
- size.height = 1;
- }
-
- IPPI_CALL( func( mat->data.ptr, mat_step, mask->data.ptr, mask_step,
- sum->data.ptr, sum_step, size, (float)alpha ));
- }
-
- __END__;
-}
-
-
-/* End of file. */