public
/
Ack-5.5.git
/ blob
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
history
|
raw
|
HEAD
Pristine Ack-5.5
[Ack-5.5.git]
/
lang
/
fortran
/
lib
/
libF77
/
pow_zz.c
1
#include "f2c.h"
2
3
VOID pow_zz(r,a,b)
4
doublecomplex *r, *a, *b;
5
{
6
double logr, logi, x, y;
7
double log(), exp(), cos(), sin(), atan2(), cabs();
8
9
logr = log( cabs(a->r, a->i) );
10
logi = atan2(a->i, a->r);
11
12
x = exp( logr * b->r - logi * b->i );
13
y = logr * b->i + logi * b->r;
14
15
r->r = x * cos(y);
16
r->i = x * sin(y);
17
}