For input of a complex number, Z = X + iY, the complex arccosine is given by,
acos(Z) = acos(B) - i alog(A + sqrt(A2 - 1)) if Y >= 0
acos(Z) = acos(B) + i alog(A + sqrt(A2 - 1)) if Y < 0
A = 0.5 sqrt((X + 1)2 + Y2) + 0.5 sqrt((X - 1)2 + Y2)
B = 0.5 sqrt((X + 1)2 + Y2) - 0.5 sqrt((X - 1)2 + Y2)
The separation of the two formulas at Y = 0 takes into account the branch-cut discontinuity along the real axis from -∞ to -1 and +1 to +∞, and ensures that cos(acos(Z)) is equal to Z. For reference, see formulas 4.4.37-39 in Abramowitz, M. and Stegun, I.A., 1964: Handbook of Mathematical Functions (Washington: National Bureau of Standards).
Find the angle whose cosine is 0.707 and print the result in degrees by entering:
Find the complex arccosine of 2 + i and print the result by entering:
( 0.507356, -1.46935)
See the ATAN function for an example of visualizing the complex arccosine.
Result = ACOS(X)
Returns the angle, expressed in radians, whose cosine is X (i.e., the arc-cosine).
The cosine of the desired angle. For real input, X should be in the range -1 to +1. If X is double-precision floating or complex, the result is of the same type. All other types are converted to single-precision floating-point and yield floating-point results. If X is an array, the result has the same structure, with each element containing the arc-cosine of the corresponding element of X.
Thread Pool Keywords
This routine is written to make use of IDL’s thread pool, which can increase execution speed on systems with multiple CPUs. The values stored in the !CPU system variable control whether IDL uses the thread pool for a given computation. In addition, you can use the thread pool keywords TPOOL_MAX_ELTS, TPOOL_MIN_ELTS, and TPOOL_NOTHREAD to override the defaults established by !CPU for a single invocation of this routine. See Thread Pool Keywords for details.
This page has no comments yet. Be the first one!