函数风格的Scala双epsilon计算

Question

使用Java计算机器epsilon的建议方法如下，

private static float calculateMachineEpsilonFloat() {
    float machEps = 1.0f;
    do
       machEps /= 2.0f;
    while ((float) (1.0 + (machEps / 2.0)) != 1.0);
    return machEps;
}

如何计算epsilon为scala.Double使用Scala在习惯用法/功能风格从这个建议的代码?

Answer 1

可以吗?

scala> val s: Stream[Float] = 1.0f #:: s.map(f => f / 2.0f)
s: Stream[Float] = Stream(1.0, ?)
scala> val eps = s.takeWhile(e => e + 1.0f != 1.0f).last
eps: Float = 1.1920929E-7

要得到另一个epsilon(2^{-24})，可以用dropWhile(因此也可以用head)代替。

Answer 2

你可以用递归重写do-while循环:

import scala.annotation.tailrec
def calculateMachineEpsilonFloat = {
  @tailrec
  def calc(machEps: Float): Float = {
    if ((1.0 + (machEps / 2.0)).toFloat != 1.0)
       calc(machEps / 2f)
    else
      machEps
  }
  calc(1f)
}