X-Git-Url: https://vcs.maemo.org/git/?a=blobdiff_plain;f=cv%2Fsrc%2Fcvaccum.cpp;fp=cv%2Fsrc%2Fcvaccum.cpp;h=0000000000000000000000000000000000000000;hb=e4c14cdbdf2fe805e79cd96ded236f57e7b89060;hp=21b715f472b49a4743a2bea382b255259611fa5f;hpb=454138ff8a20f6edb9b65a910101403d8b520643;p=opencv diff --git a/cv/src/cvaccum.cpp b/cv/src/cvaccum.cpp deleted file mode 100644 index 21b715f..0000000 --- a/cv/src/cvaccum.cpp +++ /dev/null @@ -1,786 +0,0 @@ -/*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. */