Забавная арифметика.
Есть фунция на Julia
function pconcat(x,y,arc)
xe = exp(x)
xm = 1.0/xe
ye = exp(y)
ym = 1.0/ye
c = cos(arc)
xc = xe+xm
xs = xe-xm
yc = ye+ym
ys = ye-ym
l = (xc*yc + xs*ys*c)/4
larc = (xc*l - yc)/(xs*(sqrt((l-1)*(l+1))))
(acosh(l), acos(larc)*sign(arc))
end
Вольфрам альфа говорит, что pconcat(x,y,0) должен быть равен (x+y,0). Чаще всего так оно и есть, но если вдруг x и y окажутся близки к N*1.09054966350708617, иногда наступает военное время и косинус оказывается больше единицы. Там надо явную проверку вставить, или есть какой-то более правильный способ это починить?