minifloat.c 1.9 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162
  1. /*
  2. * Copyright (C) 2014 The Android Open Source Project
  3. *
  4. * Licensed under the Apache License, Version 2.0 (the "License");
  5. * you may not use this file except in compliance with the License.
  6. * You may obtain a copy of the License at
  7. *
  8. * http://www.apache.org/licenses/LICENSE-2.0
  9. *
  10. * Unless required by applicable law or agreed to in writing, software
  11. * distributed under the License is distributed on an "AS IS" BASIS,
  12. * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  13. * See the License for the specific language governing permissions and
  14. * limitations under the License.
  15. */
  16. #include <math.h>
  17. #include <audio_utils/minifloat.h>
  18. #define EXPONENT_BITS 3
  19. #define EXPONENT_MAX ((1 << EXPONENT_BITS) - 1)
  20. #define EXCESS ((1 << EXPONENT_BITS) - 2)
  21. #define MANTISSA_BITS 13
  22. #define MANTISSA_MAX ((1 << MANTISSA_BITS) - 1)
  23. #define HIDDEN_BIT (1 << MANTISSA_BITS)
  24. #define ONE_FLOAT ((float) (1 << (MANTISSA_BITS + 1)))
  25. #define MINIFLOAT_MAX ((EXPONENT_MAX << MANTISSA_BITS) | MANTISSA_MAX)
  26. #if EXPONENT_BITS + MANTISSA_BITS != 16
  27. #error EXPONENT_BITS and MANTISSA_BITS must sum to 16
  28. #endif
  29. gain_minifloat_t gain_from_float(float v)
  30. {
  31. if (isnan(v) || v <= 0.0f) {
  32. return 0;
  33. }
  34. if (v >= 2.0f) {
  35. return MINIFLOAT_MAX;
  36. }
  37. int exp;
  38. float r = frexpf(v, &exp);
  39. if ((exp += EXCESS) > EXPONENT_MAX) {
  40. return MINIFLOAT_MAX;
  41. }
  42. if (-exp >= MANTISSA_BITS) {
  43. return 0;
  44. }
  45. int mantissa = (int) (r * ONE_FLOAT);
  46. return exp > 0 ? (exp << MANTISSA_BITS) | (mantissa & ~HIDDEN_BIT) :
  47. (mantissa >> (1 - exp)) & MANTISSA_MAX;
  48. }
  49. float float_from_gain(gain_minifloat_t a)
  50. {
  51. int mantissa = a & MANTISSA_MAX;
  52. int exponent = (a >> MANTISSA_BITS) & EXPONENT_MAX;
  53. return ldexpf((exponent > 0 ? HIDDEN_BIT | mantissa : mantissa << 1) / ONE_FLOAT,
  54. exponent - EXCESS);
  55. }