evaluate complex polynomial and derivative */
| 21 | } |
| 22 | /* evaluate complex polynomial and derivative */ |
| 23 | COMPLEX |
| 24 | pj_zpolyd1(COMPLEX z, const COMPLEX *C, int n, COMPLEX *der) { |
| 25 | COMPLEX a, b; |
| 26 | double t; |
| 27 | int first = 1; |
| 28 | |
| 29 | a = *(C += n); |
| 30 | b = a; |
| 31 | while (n-- > 0) { |
| 32 | if (first) { |
| 33 | first = 0; |
| 34 | } else { |
| 35 | b.r = a.r + z.r * (t = b.r) - z.i * b.i; |
| 36 | b.i = a.i + z.r * b.i + z.i * t; |
| 37 | } |
| 38 | a.r = (--C)->r + z.r * (t = a.r) - z.i * a.i; |
| 39 | a.i = C->i + z.r * a.i + z.i * t; |
| 40 | } |
| 41 | b.r = a.r + z.r * (t = b.r) - z.i * b.i; |
| 42 | b.i = a.i + z.r * b.i + z.i * t; |
| 43 | a.r = z.r * (t = a.r) - z.i * a.i; |
| 44 | a.i = z.r * a.i + z.i * t; |
| 45 | *der = b; |
| 46 | return a; |
| 47 | } |
no outgoing calls
no test coverage detected