1 subroutine zpassf (nac,ido,ip,l1,idl1,cc,c1,c2,ch,ch2,wa)
2 implicit double precision (a-h,o-z)
4 1 c1(ido,l1,ip) ,wa(1) ,c2(idl1,ip),
12 if (ido .lt. l1) go
to 106
17 ch(i,k,j) = cc(i,j,k)+cc(i,jc,k)
18 ch(i,k,jc) = cc(i,j,k)-cc(i,jc,k)
32 ch(i,k,j) = cc(i,j,k)+cc(i,jc,k)
33 ch(i,k,jc) = cc(i,j,k)-cc(i,jc,k)
48 c2(ik,l) = ch2(ik,1)+wa(idl-1)*ch2(ik,2)
49 c2(ik,lc) = -wa(idl)*ch2(ik,ip)
56 if (idlj .gt. idp) idlj = idlj-idp
60 c2(ik,l) = c2(ik,l)+war*ch2(ik,j)
61 c2(ik,lc) = c2(ik,lc)-wai*ch2(ik,jc)
67 ch2(ik,1) = ch2(ik,1)+ch2(ik,j)
73 ch2(ik-1,j) = c2(ik-1,j)-c2(ik,jc)
74 ch2(ik-1,jc) = c2(ik-1,j)+c2(ik,jc)
75 ch2(ik,j) = c2(ik,j)+c2(ik-1,jc)
76 ch2(ik,jc) = c2(ik,j)-c2(ik-1,jc)
80 if (ido .eq. 2)
return
91 if (idot .gt. l1) go
to 127
98 c1(i-1,k,j) = wa(idij-1)*ch(i-1,k,j)+wa(idij)*ch(i,k,j)
99 c1(i,k,j) = wa(idij-1)*ch(i,k,j)-wa(idij)*ch(i-1,k,j)
111 c1(i-1,k,j) = wa(idij-1)*ch(i-1,k,j)+wa(idij)*ch(i,k,j)
112 c1(i,k,j) = wa(idij-1)*ch(i,k,j)-wa(idij)*ch(i-1,k,j)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
subroutine zpassf(nac, ido, ip, l1, idl1, cc, c1, c2, ch, ch2, wa)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)