fortran acos函数参数健壮性

Question

使用fortran时，我有一些例程，可以多次调用acos()和asin()以及atan()等。在特定的例程中，当使用单精度变量时，参数的计算有时会变得不精确(即acos(1.0000001))。使用双精度，这种情况会消失，但速度会减慢约50 -70%。我能想到的唯一另一种方法是限制参数：

arg3= arg1/arg2
if (angle > 1) then
  arg3= 1
else if (angle < -1) then
  arg3= -1
end if
angle= acos(arg3)

尽管这看起来是一种缓慢的方法。有没有更快的方法？例如，预先限制angle，以便它只接受值，在本例中，介于-1和1之间？

编辑:参数是实数，结果需要是实数(非复数)。

Answer 1

我可能会试试

angle = acos(sign(min(abs(arg1/arg2),1.0),arg1/arg2))

并且希望我的编译器不会两次执行arg1/arg2。如果是这样，那么我会尝试预先计算它，并使用它两次。