pow_ci.c 412 B

1234567891011121314151617181920212223242526
  1. #include "f2c.h"
  2. #ifdef __cplusplus
  3. extern "C" {
  4. #endif
  5. #ifdef KR_headers
  6. VOID pow_ci(p, a, b) /* p = a**b */
  7. complex *p, *a; integer *b;
  8. #else
  9. extern void pow_zi(doublecomplex*, doublecomplex*, integer*);
  10. void pow_ci(complex *p, complex *a, integer *b) /* p = a**b */
  11. #endif
  12. {
  13. doublecomplex p1, a1;
  14. a1.r = a->r;
  15. a1.i = a->i;
  16. pow_zi(&p1, &a1, b);
  17. p->r = p1.r;
  18. p->i = p1.i;
  19. }
  20. #ifdef __cplusplus
  21. }
  22. #endif