1
subroutine wdpowe(vr,vi,p,rr,ri,ierr)
3
c computes v^p with v complex and p real
5
c subroutine wdpowe(vr,vi,p,rr,ri,ierr)
7
c double precision vr,vi,p,rr,ri
9
c vi : imaginary part of v
10
c rr : result's real part
11
c ri : result's imaginary part
15
c ierr=2 if 0**k with k<0
17
c Serge Steer INRIA 1996
22
double precision vr,vi,p,sr,si,rr,ri
26
if(dble(int(p)).eq.p) then
27
call wipowe(vr,vi,int(p),rr,ri,ierr)
29
if(abs(vr)+abs(vi).ne.0.0d+0) then
30
call wlog(vr,vi,sr,si)
39
elseif(p.lt.0.0d+0) then