1 subroutine passb (nac,ido,ip,l1,idl1,cc,c1,c2,ch,ch2,wa)
3 1 c1(ido,l1,ip) ,wa(1) ,c2(idl1,ip),
11 if (ido .lt. l1) go
to 106
16 ch(i,k,j) = cc(i,j,k)+cc(i,jc,k)
17 ch(i,k,jc) = cc(i,j,k)-cc(i,jc,k)
31 ch(i,k,j) = cc(i,j,k)+cc(i,jc,k)
32 ch(i,k,jc) = cc(i,j,k)-cc(i,jc,k)
47 c2(ik,l) = ch2(ik,1)+wa(idl-1)*ch2(ik,2)
48 c2(ik,lc) = wa(idl)*ch2(ik,ip)
55 if (idlj .gt. idp) idlj = idlj-idp
59 c2(ik,l) = c2(ik,l)+war*ch2(ik,j)
60 c2(ik,lc) = c2(ik,lc)+wai*ch2(ik,jc)
66 ch2(ik,1) = ch2(ik,1)+ch2(ik,j)
72 ch2(ik-1,j) = c2(ik-1,j)-c2(ik,jc)
73 ch2(ik-1,jc) = c2(ik-1,j)+c2(ik,jc)
74 ch2(ik,j) = c2(ik,j)+c2(ik-1,jc)
75 ch2(ik,jc) = c2(ik,j)-c2(ik-1,jc)
79 if (ido .eq. 2)
return
90 if (idot .gt. l1) go
to 127
97 c1(i-1,k,j) = wa(idij-1)*ch(i-1,k,j)-wa(idij)*ch(i,k,j)
98 c1(i,k,j) = wa(idij-1)*ch(i,k,j)+wa(idij)*ch(i-1,k,j)
110 c1(i-1,k,j) = wa(idij-1)*ch(i-1,k,j)-wa(idij)*ch(i,k,j)
111 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 passb(nac, ido, ip, l1, idl1, cc, c1, c2, ch, ch2, wa)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)