Quelle  spl_custom_float.c   Sprache: C

// SPDX-License-Identifier: MIT
//
// Copyright 2024 Advanced Micro Devices, Inc.

#include "spl_debug.h"
#include "spl_custom_float.h"

static bool spl_build_custom_float(struct spl_fixed31_32 value,
          const struct spl_custom_float_format *format,
          bool *negative,
          uint32_t *mantissa,
          uint32_t *exponenta)
{
 uint32_t exp_offset = (1 << (format->exponenta_bits - 1)) - 1;

 const struct spl_fixed31_32 mantissa_constant_plus_max_fraction =
  spl_fixpt_from_fraction((1LL << (format->mantissa_bits + 1)) - 1,
           1LL << format->mantissa_bits);

 struct spl_fixed31_32 mantiss;

 if (spl_fixpt_eq(value, spl_fixpt_zero)) {
  *negative = false;
  *mantissa = 0;
  *exponenta = 0;
  return true;
 }

 if (spl_fixpt_lt(value, spl_fixpt_zero)) {
  *negative = format->sign;
  value = spl_fixpt_neg(value);
 } else {
  *negative = false;
 }

 if (spl_fixpt_lt(value, spl_fixpt_one)) {
  uint32_t i = 1;

  do {
   value = spl_fixpt_shl(value, 1);
   ++i;
  } while (spl_fixpt_lt(value, spl_fixpt_one));

  --i;

  if (exp_offset <= i) {
   *mantissa = 0;
   *exponenta = 0;
   return true;
  }

  *exponenta = exp_offset - i;
 } else if (spl_fixpt_le(mantissa_constant_plus_max_fraction, value)) {
  uint32_t i = 1;

  do {
   value = spl_fixpt_shr(value, 1);
   ++i;
  } while (spl_fixpt_lt(mantissa_constant_plus_max_fraction, value));

  *exponenta = exp_offset + i - 1;
 } else {
  *exponenta = exp_offset;
 }

 mantiss = spl_fixpt_sub(value, spl_fixpt_one);

 if (spl_fixpt_lt(mantiss, spl_fixpt_zero) ||
     spl_fixpt_lt(spl_fixpt_one, mantiss))
  mantiss = spl_fixpt_zero;
 else
  mantiss = spl_fixpt_shl(mantiss, format->mantissa_bits);

 *mantissa = spl_fixpt_floor(mantiss);

 return true;
}

static bool spl_setup_custom_float(const struct spl_custom_float_format *format,
          bool negative,
          uint32_t mantissa,
          uint32_t exponenta,
          uint32_t *result)
{
 uint32_t i = 0;
 uint32_t j = 0;
 uint32_t value = 0;

 /* verification code:
 * once calculation is ok we can remove it
 */

 const uint32_t mantissa_mask =
  (1 << (format->mantissa_bits + 1)) - 1;

 const uint32_t exponenta_mask =
  (1 << (format->exponenta_bits + 1)) - 1;

 if (mantissa & ~mantissa_mask) {
  SPL_BREAK_TO_DEBUGGER();
  mantissa = mantissa_mask;
 }

 if (exponenta & ~exponenta_mask) {
  SPL_BREAK_TO_DEBUGGER();
  exponenta = exponenta_mask;
 }

 /* end of verification code */

 while (i < format->mantissa_bits) {
  uint32_t mask = 1 << i;

  if (mantissa & mask)
   value |= mask;

  ++i;
 }

 while (j < format->exponenta_bits) {
  uint32_t mask = 1 << j;

  if (exponenta & mask)
   value |= mask << i;

  ++j;
 }

 if (negative && format->sign)
  value |= 1 << (i + j);

 *result = value;

 return true;
}

bool spl_convert_to_custom_float_format(struct spl_fixed31_32 value,
        const struct spl_custom_float_format *format,
        uint32_t *result)
{
 uint32_t mantissa;
 uint32_t exponenta;
 bool negative;

 return spl_build_custom_float(value, format, &negative, &mantissa, &exponenta) &&
      spl_setup_custom_float(format,
           negative,
           mantissa,
           exponenta,
           result);
}