1234567891011121314151617181920212223242526272829303132 |
- import sinh from './sinh';
- import cosh from './cosh';
- export default function(pp, arg_r, arg_i) {
- var sin_arg_r = Math.sin(arg_r);
- var cos_arg_r = Math.cos(arg_r);
- var sinh_arg_i = sinh(arg_i);
- var cosh_arg_i = cosh(arg_i);
- var r = 2 * cos_arg_r * cosh_arg_i;
- var i = -2 * sin_arg_r * sinh_arg_i;
- var j = pp.length - 1;
- var hr = pp[j];
- var hi1 = 0;
- var hr1 = 0;
- var hi = 0;
- var hr2;
- var hi2;
- while (--j >= 0) {
- hr2 = hr1;
- hi2 = hi1;
- hr1 = hr;
- hi1 = hi;
- hr = -hr2 + r * hr1 - i * hi1 + pp[j];
- hi = -hi2 + i * hr1 + r * hi1;
- }
- r = sin_arg_r * cosh_arg_i;
- i = cos_arg_r * sinh_arg_i;
- return [r * hr - i * hi, r * hi + i * hr];
- }
|