Stable Cotangente

Est-il plus stable de la mise en œuvre de la cotangente de la fonction de retour 1.0/tan(x);?

  • Je suppose que vous êtes problème avec 1/tan est que ses pas défini à pi/2 lorsque la fonction de retour à zéro?
  • Faire l'inversion est douteux quand près de pi/2 oui... je me demandais si il y avait une meilleure façon de le faire.
InformationsquelleAutor Meh. | 2010-09-17
  1. 37

    cot(x) = cos(x)/sin(x) devrait être plus numériquement stable à proximité de π/2 que cot(x) = 1/tan(x). Vous pouvez mettre en œuvre de façon efficace à l'aide de sincos sur les plates-formes qui ont.

    Une autre possibilité est cot(x) = tan(M_PI_2 - x). Cela devrait être plus rapide que le dessus (même si sincos est disponible), mais il peut également être moins précis, car M_PI_2 est évidemment qu'une approximation de la transcendance du nombre π/2, de sorte que la différence M_PI_2 - x ne sera pas précis à la largeur d'une double mantisse -- en fait, si vous êtes malchanceux, il peut n'avoir que quelques bits significatifs.

    • Ok. Je suppose que c'est le mieux que je puisse faire. Aussi j'ai appris beaucoup sur sincos, qui je dois l'avouer, je n'avais jamais rencontré avant!
    • Viens de me rappeler un trig d'identité qui permet de faire beaucoup mieux & édité ma réponse.
    • +1. Certainement utiliser tan(pi/2 - x).
    InformationsquelleAutor zwol