1
subroutine wwpowe(vr,vi,pr,pi,rr,ri,ierr)
3
c computes v^p with v complex and p complex
5
c subroutine wdpowe(vr,vi,p,rr,ri,ierr)
7
c double precision vr,vi,pr,pi,rr,ri
9
c vi : imaginary part of v
11
c pi : imaginary part of p
12
c rr : result's real part
13
c ri : result's imaginary part
17
c ierr=2 if 0**k with k<0
19
c Serge Steer INRIA 1996
23
double precision vr,vi,pr,pi,sr,si,rr,ri
28
call wdpowe(vr,vi,pr,rr,ri,ierr)
30
if(abs(vr)+abs(vi).ne.0.0d+0) then
31
call wlog(vr,vi,sr,si)
32
call wmul(sr,si,pr,pi,sr,si)