File: dipowe.f

package info (click to toggle)
scilab 2.6-4
  • links: PTS
  • area: non-free
  • in suites: woody
  • size: 54,632 kB
  • ctags: 40,267
  • sloc: ansic: 267,851; fortran: 166,549; sh: 10,005; makefile: 4,119; tcl: 1,070; cpp: 233; csh: 143; asm: 135; perl: 130; java: 39
file content (45 lines) | stat: -rw-r--r-- 870 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
 
      subroutine dipowe(v,p,r,ierr)
c!purpose
c     computes r=v^p where v double precision and b integer
c!calling sequence
c     subroutine dipowe(v,p,r,ierr)
c     integer p ,ierr
c     double precision v,r
c
c     ierr : error flag
c            ierr=0 if ok
c            ierr=1 if 0**0
c            ierr=2 if  0**k with k<0
c!origin
c Serge Steer INRIA 1996
c!
c     Copyright INRIA
      integer p ,ierr
      double precision v,r
c
      ierr=0
c
      if(p.eq.1) then
         r=v
      elseif(p.eq.0) then
c     .  v^0
c         if(v.ne.0.0d+0) then
            r=1.0d+0
c         else

c            ierr=1
c         endif
      else if(p.lt.0) then
         if(v.ne.0.0d+0) then
            r=v**p
         else
            r=0.0d0
            r=infinity(r)
            ierr=2
         endif
      else
         r=v**p
      endif
c     
      return
      end