math_utils.cpp
上传用户:gelin96
上传日期:2017-01-08
资源大小:20993k
文件大小:5k
源码类别:

MTK

开发平台:

C++ Builder

  1. /*****************************************************************************
  2. *  Copyright Statement:
  3. *  --------------------
  4. *  This software is protected by Copyright and the information contained
  5. *  herein is confidential. The software may not be copied and the information
  6. *  contained herein may not be used or disclosed except with the written
  7. *  permission of MediaTek Inc. (C) 2005
  8. *
  9. *  BY OPENING THIS FILE, BUYER HEREBY UNEQUIVOCALLY ACKNOWLEDGES AND AGREES
  10. *  THAT THE SOFTWARE/FIRMWARE AND ITS DOCUMENTATIONS ("MEDIATEK SOFTWARE")
  11. *  RECEIVED FROM MEDIATEK AND/OR ITS REPRESENTATIVES ARE PROVIDED TO BUYER ON
  12. *  AN "AS-IS" BASIS ONLY. MEDIATEK EXPRESSLY DISCLAIMS ANY AND ALL WARRANTIES,
  13. *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE IMPLIED WARRANTIES OF
  14. *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT.
  15. *  NEITHER DOES MEDIATEK PROVIDE ANY WARRANTY WHATSOEVER WITH RESPECT TO THE
  16. *  SOFTWARE OF ANY THIRD PARTY WHICH MAY BE USED BY, INCORPORATED IN, OR
  17. *  SUPPLIED WITH THE MEDIATEK SOFTWARE, AND BUYER AGREES TO LOOK ONLY TO SUCH
  18. *  THIRD PARTY FOR ANY WARRANTY CLAIM RELATING THERETO. MEDIATEK SHALL ALSO
  19. *  NOT BE RESPONSIBLE FOR ANY MEDIATEK SOFTWARE RELEASES MADE TO BUYER'S
  20. *  SPECIFICATION OR TO CONFORM TO A PARTICULAR STANDARD OR OPEN FORUM.
  21. *
  22. *  BUYER'S SOLE AND EXCLUSIVE REMEDY AND MEDIATEK'S ENTIRE AND CUMULATIVE
  23. *  LIABILITY WITH RESPECT TO THE MEDIATEK SOFTWARE RELEASED HEREUNDER WILL BE,
  24. *  AT MEDIATEK'S OPTION, TO REVISE OR REPLACE THE MEDIATEK SOFTWARE AT ISSUE,
  25. *  OR REFUND ANY SOFTWARE LICENSE FEES OR SERVICE CHARGE PAID BY BUYER TO
  26. *  MEDIATEK FOR SUCH MEDIATEK SOFTWARE AT ISSUE.
  27. *
  28. *  THE TRANSACTION CONTEMPLATED HEREUNDER SHALL BE CONSTRUED IN ACCORDANCE
  29. *  WITH THE LAWS OF THE STATE OF CALIFORNIA, USA, EXCLUDING ITS CONFLICT OF
  30. *  LAWS PRINCIPLES.  ANY DISPUTES, CONTROVERSIES OR CLAIMS ARISING THEREOF AND
  31. *  RELATED THERETO SHALL BE SETTLED BY ARBITRATION IN SAN FRANCISCO, CA, UNDER
  32. *  THE RULES OF THE INTERNATIONAL CHAMBER OF COMMERCE (ICC).
  33. *
  34. *****************************************************************************/
  35. /*****************************************************************************
  36.  *
  37.  * Filename:
  38.  * ---------
  39.  *   math_utils.cpp
  40.  *
  41.  * Project:
  42.  * --------
  43.  *   Maui META APP
  44.  *
  45.  * Description:
  46.  * ------------
  47.  *  Math utility source
  48.  *
  49.  * Author:
  50.  * -------
  51.  *  Andy Ueng (mtk00490)
  52.  *
  53.  *============================================================================
  54.  *             HISTORY
  55.  * Below this line, this part is controlled by PVCS VM. DO NOT MODIFY!!
  56.  *------------------------------------------------------------------------------
  57.  * $Revision$
  58.  * $Modtime$
  59.  * $Log$
  60.  *
  61.  *------------------------------------------------------------------------------
  62.  * Upper this line, this part is controlled by PVCS VM. DO NOT MODIFY!!
  63.  *============================================================================
  64.  ****************************************************************************/
  65. #include <math.h>
  66. #include <math.hpp>
  67. #include "visatype.h"
  68. #pragma hdrstop
  69. #ifndef  _MATH_UTILS_H_
  70. #include "math_utils.h"
  71. #endif
  72. //------------------------------------------------------------------------------
  73. double  Round_double( double d_in )
  74. {
  75.     double d_diff1, d_diff2, d_out;
  76.     d_diff1 = d_in - (double)Ceil(d_in);
  77.     d_diff2 = d_in - (double)Floor(d_in);
  78.     if( Abs_double(d_diff1) < Abs_double(d_diff2) )
  79.     {
  80.         d_out = (double) Ceil(d_in);
  81.     }
  82.     else
  83.     {
  84.         d_out = (double) Floor(d_in);
  85.     }
  86.     return  d_out;
  87. }
  88. //------------------------------------------------------------------------------
  89. double Abs_double(double d_in)
  90. {
  91.     double  d_out;
  92.     if( d_in>0 )
  93.     {
  94.        d_out = d_in;
  95.     }
  96.     else
  97.     {
  98.        d_out = -d_in;
  99.     }
  100.     return  d_out;
  101. }
  102. //------------------------------------------------------------------------------
  103. bool AbsOfComplexArray( double *d_A, double *d_B, int size, double *d_ABS )
  104. {
  105.     for( int i=0; i<size; i++ )
  106.     {
  107.         *(d_ABS+i) = sqrt( pow(*(d_A+i), 2.0) + pow(*(d_B+i), 2.0) );
  108.     }
  109.     return  true;
  110. }
  111. //------------------------------------------------------------------------------
  112. bool MaxOfDoubleArray( double *d_A, int size, double *d_max )
  113. {
  114.     double d_temp_max = -100000;
  115.     for( int i=0; i<size; i++ )
  116.     {
  117.          if( *(d_A+i) > d_temp_max )
  118.          {
  119.               d_temp_max =  *(d_A+i);
  120.          }
  121.     }
  122.     *d_max = d_temp_max;
  123.     return  true;
  124. }
  125. //------------------------------------------------------------------------------
  126. double SumOfDoubleArray( double *d_A, int size )
  127. {
  128.     double d_sum = 0;
  129.     for(int i=0; i<size; i++)
  130.     {
  131.         d_sum += *(d_A + i);
  132.     }
  133.     return d_sum;
  134. }
  135. //------------------------------------------------------------------------------
  136. int SortByIncreasing( const void *v1, const void *v2 )
  137. {
  138.     // need to use pointer to ViReal64 to access values
  139.     const ViReal64 *ViReal1 = (ViReal64 *) v1;
  140.     const ViReal64 *ViReal2 = (ViReal64 *) v2;
  141.     if( *ViReal1 < *ViReal2 )
  142.         return -1;
  143.     else if ( *ViReal1 == *ViReal2 )
  144.         return 0;
  145.     else
  146.         return 1;
  147. }