fixpoint.h
/* -*- c-basic-offset: 4; indent-tabs-mode: nil -*- */
/* ====================================================================
* Copyright (c) 2005 Carnegie Mellon University. All rights
* reserved.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions
* are met:
*
* 1. Redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer.
*
* 2. Redistributions 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.
*
* This work was supported in part by funding from the Defense Advanced
* Research Projects Agency and the National Science Foundation of the
* United States of America, and the CMU Sphinx Speech Consortium.
*
* THIS SOFTWARE IS PROVIDED BY CARNEGIE MELLON UNIVERSITY ``AS IS'' AND
* ANY EXPRESSED 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 CARNEGIE MELLON UNIVERSITY
* NOR ITS EMPLOYEES 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.
*
* ==================================================================== */
/* Fixed-point arithmetic macros.
*
* Author: David Huggins-Daines <dhuggins@cs.cmu.edu>
*/
#ifndef _FIXPOINT_H_
#define _FIXPOINT_H_
#include <limits.h>
/* Win32/WinCE DLL gunk */
#include <sphinxbase/sphinxbase_export.h>
#include <sphinxbase/prim_type.h>
#ifdef __cplusplus
extern "C" {
#endif
#if 0
/* Fool Emacs. */
}
#endif
#ifndef DEFAULT_RADIX
#define DEFAULT_RADIX 12
#endif
typedef int32 fixed32;
#define FLOAT2FIX_ANY(x,radix) \
(((x)<0.0) ? \
((fixed32)((x)*(float32)(1<<(radix)) - 0.5)) \
: ((fixed32)((x)*(float32)(1<<(radix)) + 0.5)))
#define FLOAT2FIX(x) FLOAT2FIX_ANY(x,DEFAULT_RADIX)
#define FIX2FLOAT_ANY(x,radix) ((float32)(x)/(1<<(radix)))
#define FIX2FLOAT(x) FIX2FLOAT_ANY(x,DEFAULT_RADIX)
/*
* This works on most modern ARMs but *only* in ARM mode (for obvious
* reasons), so don't use it in Thumb mode (but why are you building
* signal processing code in Thumb mode?!)
*/
#if defined(__arm__) && !defined(__thumb__)
#define FIXMUL(a,b) FIXMUL_ANY(a,b,DEFAULT_RADIX)
#define FIXMUL_ANY(a,b,r) ({ \
int cl, ch, _a = a, _b = b; \
__asm__ ("smull %0, %1, %2, %3\n" \
"mov %0, %0, lsr %4\n" \
"orr %0, %0, %1, lsl %5\n" \
: "=&r" (cl), "=&r" (ch) \
: "r" (_a), "r" (_b), "i" (r), "i" (32-(r)));\
cl; })
#elif defined(BFIN) && DEFAULT_RADIX == 16
/* Blackfin magic */
#undef FIXMUL
/* Use the accumulators for the 16.16 case (probably not as efficient as it could be). */
#define FIXMUL(a,b) ({ \
int c, _a = a, _b = b; \
__asm__("%0.L = %1.l * %2.l (FU);\n\t" \
"%0.H = %1.h * %2.h (IS);\n\t" \
"A1 = %0;\n\t" \
"A1 += %1.h * %2.l (IS, M);\n\t" \
"%0 = (A1 += %2.h * %1.l) (IS, M);\n\t" \
: "=&W" (c) \
: "d" (_a), "d" (_b) \
: "A1", "cc"); \
c; })
#define FIXMUL_ANY(a,b,radix) ((fixed32)(((int64)(a)*(b))>>(radix)))
#elif defined(_MSC_VER) || (defined(HAVE_LONG_LONG) && SIZEOF_LONG_LONG == 8)
#define FIXMUL(a,b) FIXMUL_ANY(a,b,DEFAULT_RADIX)
#define FIXMUL_ANY(a,b,radix) ((fixed32)(((int64)(a)*(b))>>(radix)))
#else /* Most general case where 'long long' doesn't exist or is slow. */
#define FIXMUL(a,b) FIXMUL_ANY(a,b,DEFAULT_RADIX)
#define FIXMUL_ANY(a,b,radix) \
(fixed32)(((((uint32)(a))&((1<<(radix))-1)) \
* (((uint32)(b))&((1<<(radix))-1)) >> (radix)) \
+ (((((int32)(a))>>(radix)) * (((int32)(b))>>(radix))) << (radix)) \
+ ((((uint32)(a))&((1<<(radix))-1)) * (((int32)(b))>>(radix))) \
+ ((((uint32)(b))&((1<<(radix))-1)) * (((int32)(a))>>(radix))))
#endif
/* Various fixed-point logarithmic functions that we need. */
#define MIN_FIXLOG -2829416 /* log(1e-300) * (1<<DEFAULT_RADIX) */
#define MIN_FIXLOG2 -4081985 /* log2(1e-300) * (1<<DEFAULT_RADIX) */
#define FIXLN_2 ((fixed32)(0.693147180559945 * (1<<DEFAULT_RADIX)))
#define FIXLN(x) (fixlog(x) - (FIXLN_2 * DEFAULT_RADIX))
int32 fixlog(uint32 x);
int32 fixlog2(uint32 x);
#ifdef __cplusplus
}
#endif
#endif /* _FIXPOINT_H_ */
prim_type.h
Basic type definitions used in Sphinx.