File: cuproi.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 (15 lines) | stat: -rw-r--r-- 326 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
	subroutine cuproi(n,wr,wi)
c      Utility fct: cumulated product, complex argument
c	Copyright INRIA
	double precision wr(*),wi(*),tr,ti,wwr
	tr=1.0d0
        ti=0.0d0
	do 1 k=1,n
c	w(k)=t*w(k)
           wwr=wr(k)
           wr(k)=tr*wwr-ti*wi(k)
           wi(k)=tr*wi(k)+ti*wwr
	tr=wr(k)
        ti=wi(k)
 1	continue
	end