Числовая устойчивость аргумента комплексного числа / ответвлений - PullRequest
0 голосов
/ 23 октября 2018

Я выполняю числовую оценку некоторых аналитических выражений, которые включают такие факторы, как exp(1i*arg(z) / 2), где z в принципе комплексное число, которое иногда бывает почти действительным (т. Е. С плавающей запятой).точность, например 4.440892098500626e-16j).

Я реализовал свои вычисления в Python и C ++ и обнаружил, что иногда результаты не совпадают, так как небольшая мнимая часть «почти вещественных» чисел немного отличается по знаку, а затем сокращается ветвьПоведение arg(z) (то есть arg(-1+0j) = pi, но arg(-1-0j) = -pi) значительно меняет результат ... Мне было интересно, существует ли какой-либо обычно используемый протокол для смягчения этих проблем?

Большое спасибо заранее.

...