math_utils.cpp
上传用户:gelin96
上传日期:2017-01-08
资源大小:20993k
文件大小:5k
- /*****************************************************************************
- * Copyright Statement:
- * --------------------
- * This software is protected by Copyright and the information contained
- * herein is confidential. The software may not be copied and the information
- * contained herein may not be used or disclosed except with the written
- * permission of MediaTek Inc. (C) 2005
- *
- * BY OPENING THIS FILE, BUYER HEREBY UNEQUIVOCALLY ACKNOWLEDGES AND AGREES
- * THAT THE SOFTWARE/FIRMWARE AND ITS DOCUMENTATIONS ("MEDIATEK SOFTWARE")
- * RECEIVED FROM MEDIATEK AND/OR ITS REPRESENTATIVES ARE PROVIDED TO BUYER ON
- * AN "AS-IS" BASIS ONLY. MEDIATEK EXPRESSLY DISCLAIMS ANY AND ALL WARRANTIES,
- * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE IMPLIED WARRANTIES OF
- * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT.
- * NEITHER DOES MEDIATEK PROVIDE ANY WARRANTY WHATSOEVER WITH RESPECT TO THE
- * SOFTWARE OF ANY THIRD PARTY WHICH MAY BE USED BY, INCORPORATED IN, OR
- * SUPPLIED WITH THE MEDIATEK SOFTWARE, AND BUYER AGREES TO LOOK ONLY TO SUCH
- * THIRD PARTY FOR ANY WARRANTY CLAIM RELATING THERETO. MEDIATEK SHALL ALSO
- * NOT BE RESPONSIBLE FOR ANY MEDIATEK SOFTWARE RELEASES MADE TO BUYER'S
- * SPECIFICATION OR TO CONFORM TO A PARTICULAR STANDARD OR OPEN FORUM.
- *
- * BUYER'S SOLE AND EXCLUSIVE REMEDY AND MEDIATEK'S ENTIRE AND CUMULATIVE
- * LIABILITY WITH RESPECT TO THE MEDIATEK SOFTWARE RELEASED HEREUNDER WILL BE,
- * AT MEDIATEK'S OPTION, TO REVISE OR REPLACE THE MEDIATEK SOFTWARE AT ISSUE,
- * OR REFUND ANY SOFTWARE LICENSE FEES OR SERVICE CHARGE PAID BY BUYER TO
- * MEDIATEK FOR SUCH MEDIATEK SOFTWARE AT ISSUE.
- *
- * THE TRANSACTION CONTEMPLATED HEREUNDER SHALL BE CONSTRUED IN ACCORDANCE
- * WITH THE LAWS OF THE STATE OF CALIFORNIA, USA, EXCLUDING ITS CONFLICT OF
- * LAWS PRINCIPLES. ANY DISPUTES, CONTROVERSIES OR CLAIMS ARISING THEREOF AND
- * RELATED THERETO SHALL BE SETTLED BY ARBITRATION IN SAN FRANCISCO, CA, UNDER
- * THE RULES OF THE INTERNATIONAL CHAMBER OF COMMERCE (ICC).
- *
- *****************************************************************************/
- /*****************************************************************************
- *
- * Filename:
- * ---------
- * math_utils.cpp
- *
- * Project:
- * --------
- * Maui META APP
- *
- * Description:
- * ------------
- * Math utility source
- *
- * Author:
- * -------
- * Andy Ueng (mtk00490)
- *
- *============================================================================
- * HISTORY
- * Below this line, this part is controlled by PVCS VM. DO NOT MODIFY!!
- *------------------------------------------------------------------------------
- * $Revision$
- * $Modtime$
- * $Log$
- *
- *------------------------------------------------------------------------------
- * Upper this line, this part is controlled by PVCS VM. DO NOT MODIFY!!
- *============================================================================
- ****************************************************************************/
- #include <math.h>
- #include <math.hpp>
- #include "visatype.h"
- #pragma hdrstop
- #ifndef _MATH_UTILS_H_
- #include "math_utils.h"
- #endif
- //------------------------------------------------------------------------------
- double Round_double( double d_in )
- {
- double d_diff1, d_diff2, d_out;
- d_diff1 = d_in - (double)Ceil(d_in);
- d_diff2 = d_in - (double)Floor(d_in);
- if( Abs_double(d_diff1) < Abs_double(d_diff2) )
- {
- d_out = (double) Ceil(d_in);
- }
- else
- {
- d_out = (double) Floor(d_in);
- }
- return d_out;
- }
- //------------------------------------------------------------------------------
- double Abs_double(double d_in)
- {
- double d_out;
- if( d_in>0 )
- {
- d_out = d_in;
- }
- else
- {
- d_out = -d_in;
- }
- return d_out;
- }
- //------------------------------------------------------------------------------
- bool AbsOfComplexArray( double *d_A, double *d_B, int size, double *d_ABS )
- {
- for( int i=0; i<size; i++ )
- {
- *(d_ABS+i) = sqrt( pow(*(d_A+i), 2.0) + pow(*(d_B+i), 2.0) );
- }
- return true;
- }
- //------------------------------------------------------------------------------
- bool MaxOfDoubleArray( double *d_A, int size, double *d_max )
- {
- double d_temp_max = -100000;
- for( int i=0; i<size; i++ )
- {
- if( *(d_A+i) > d_temp_max )
- {
- d_temp_max = *(d_A+i);
- }
- }
- *d_max = d_temp_max;
- return true;
- }
- //------------------------------------------------------------------------------
- double SumOfDoubleArray( double *d_A, int size )
- {
- double d_sum = 0;
- for(int i=0; i<size; i++)
- {
- d_sum += *(d_A + i);
- }
- return d_sum;
- }
- //------------------------------------------------------------------------------
- int SortByIncreasing( const void *v1, const void *v2 )
- {
- // need to use pointer to ViReal64 to access values
- const ViReal64 *ViReal1 = (ViReal64 *) v1;
- const ViReal64 *ViReal2 = (ViReal64 *) v2;
- if( *ViReal1 < *ViReal2 )
- return -1;
- else if ( *ViReal1 == *ViReal2 )
- return 0;
- else
- return 1;
- }