performance - 有办法避免这些嵌套循环吗？

我需要使用每个变量（wndspeed，xh，yh，A和B）的值进行函数调用
有谁知道如何在不使用这些嵌套循环的情况下做到这一点？

  do i=1,windsteps
     wndspeed =windinc*float(i-1)+windstart

     write(*,*)i,wndspeed

     do j=1,ixsteps
        xh=xinc*float(j-1)
        xh=xh/height

        do k=1,iysteps
           yh=10.*float(k-1)*1e3
           yh=yinc*float(k-1)
           yh=yh/height

           do l=1,iasteps
              A=(ainc*float(l-1)+50)*1e3


              do m=1,ibsteps
                 B=binc*float(m-1)+1



                 Pdelt=wndspeed*wndspeed*rho*exp(1.)/B

                 call trackslopes(elevang,wndspeed,coxmkcor,skewpeak,
 1                    numslops,slopary,tanumax,xh,yh,A,B,Pdelt,height,numchips,
 2                    bkgrdwnd,antgain)
                 skewpeak=0

                 call trackconvolve(height,tanumax,numslops,slopary,numchips,
 1                    stppchip,deadband,elevang,convary)


enddo
enddo
enddo
enddo

最佳答案

除非可以将tracklopes / trackconvolve过程简化为简单的方法，否则无法（非平凡地）删除嵌套循环。

最好的观察方法是每个子循环都取决于其至少一个超级循环的值，因此，除非最内层的计算能够使整个过程在数学上得以简化，否则它已经处于最简单的形式。

如果您不希望看到这么多的嵌套层次，则可以尝试使用索引数组而不是5个不同的索引来展平5-D空间，然后必须编写一个难看的1-D循环它检查5个索引中每个索引的边界，然后适当地扩展索引数组。