Есть ли Java-эквивалент frexp? - PullRequest
5 голосов
/ 12 октября 2009

Существует ли Java-эквивалент функции C / C ++ с именем frexp ? Если вы не знакомы, frexp определен в Википедии , чтобы «разбить число с плавающей точкой на мантиссу и экспоненту».

Я ищу реализацию как по скорости, так и по точности, но я предпочел бы иметь точность, если бы мог выбрать только одну.

Это пример кода из первой ссылки. Это должно сделать контракт frexp немного более понятным:

/* frexp example */
#include <stdio.h>
#include <math.h>

int main ()
{
  double param, result;
  int n;

  param = 8.0;
  result = frexp (param , &n);
  printf ("%lf * 2^%d = %f\n", result, n, param);
  return 0;
}

/* Will produce: 0.500000 * 2^4 = 8.000000 */

Ответы [ 6 ]

3 голосов
/ 16 октября 2010

Как это?

public static class FRexpResult
{
   public int exponent = 0;
   public double mantissa = 0.;
}

public static FRexpResult frexp(double value)
{
   final FRexpResult result = new FRexpResult();
   long bits = Double.doubleToLongBits(value);
   double realMant = 1.;

   // Test for NaN, infinity, and zero.
   if (Double.isNaN(value) || 
       value + value == value || 
       Double.isInfinite(value))
   {
      result.exponent = 0;
      result.mantissa = value;
   }
   else
   {

      boolean neg = (bits < 0);
      int exponent = (int)((bits >> 52) & 0x7ffL);
      long mantissa = bits & 0xfffffffffffffL;

      if(exponent == 0)
      {
         exponent++;
      }
      else
      {
         mantissa = mantissa | (1L<<52);
      }

      // bias the exponent - actually biased by 1023.
      // we are treating the mantissa as m.0 instead of 0.m
      //  so subtract another 52.
      exponent -= 1075;
      realMant = mantissa;

      // normalize
      while(realMant > 1.0) 
      {
         mantissa >>= 1;
         realMant /= 2.;
         exponent++;
      }

      if(neg)
      {
         realMant = realMant * -1;
      }

      result.exponent = exponent;
      result.mantissa = realMant;
   }
   return result;
}

Это "вдохновлено" или фактически почти идентично скопировано из ответа на аналогичный вопрос C #. Он работает с битами, а затем делает число мантиссы между 1,0 и 0,0.

1 голос
/ 12 октября 2009

См. Float.floatToIntBits и Double.doubleToLongBits. Вам все еще нужна дополнительная логика для декодирования плавающих точек IEEE 754.

0 голосов
/ 14 октября 2009

Если я правильно читаю ...

public class Frexp {
  public static void main (String[] args)
  {
    double param, result;
    int n;

    param = 8.0;
    n = Math.getExponent(param);
    //result = ??

    System.out.printf ("%f * 2^%d = %f\n", result, n, param);
  }
}

К сожалению, не существует встроенного метода для получения мантиссы без предварительного преобразования ее в BigDecimal (или просто деления: result = param / Math.pow(2,n).

Как ни странно, scalb делает прямо противоположное: возьмите мантиссу и экспоненту и сгенерируйте из нее новое число.

0 голосов
/ 14 октября 2009

Нет, в ядре Java или в Commons Lang нет текущей реализации (скорее всего, в другом месте, где ее можно найти), которая обладает точно такой же функциональностью и простотой frexp ; что я знаю Если он существует, то, вероятно, он не используется широко.

0 голосов
/ 14 октября 2009

Это делает то, что вы хотите.

public class Test {
  public class FRex {

    public FRexPHolder frexp (double value) {
      FRexPHolder ret = new FRexPHolder();

      ret.exponent = 0;
      ret.mantissa = 0;

      if (value == 0.0 || value == -0.0) {
        return ret;
      }

      if (Double.isNaN(value)) {
        ret.mantissa = Double.NaN;
        ret.exponent = -1;
        return ret;
      }

      if (Double.isInfinite(value)) {
        ret.mantissa = value;
        ret.exponent = -1;
        return ret;
      }

      ret.mantissa = value;
      ret.exponent = 0;
      int sign = 1;

      if (ret.mantissa < 0f) {
        sign--;
        ret.mantissa = -(ret.mantissa);
      }
      while (ret.mantissa < 0.5f) {
        ret.mantissa *= 2.0f;
        ret.exponent -= 1;
      }
      while (ret.mantissa >= 1.0f) {
        ret.mantissa *= 0.5f;
        ret.exponent++;
      }
      ret.mantissa *= sign;
      return ret;
    }
  }

  public class FRexPHolder {
    int exponent;
    double mantissa;
  }

  public static void main(String args[]) {
    new Test();
  }

  public Test() {
    double value = 8.0;
    //double value = 0.0;
    //double value = -0.0;
    //double value = Double.NaN;
    //double value = Double.NEGATIVE_INFINITY;
    //double value = Double.POSITIVE_INFINITY;

    FRex test = new FRex();
    FRexPHolder frexp = test.frexp(value);
    System.out.println("Mantissa: " + frexp.mantissa);
    System.out.println("Exponent: " + frexp.exponent);
    System.out.println("Original value was: " + value);
    System.out.println(frexp.mantissa+" * 2^" + frexp.exponent + " = ");
    System.out.println(frexp.mantissa*(1<<frexp.exponent));
  }
}
0 голосов
/ 12 октября 2009

Я не знаком с функцией frexp, но думаю, вам нужно посмотреть на BigDecimal 'масштабированные и немасштабированные значения. «немасштабированный» - точная мантисса, масштаб - показатель степени. В psuedocode: значение = unscaledValue 10 ^ (- масштаб)

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...