Return to tl2fit.f CVS log

Up to [Research Unix] / researchv10no / cmd / view2d / num

Annotation of researchv10no/cmd/view2d/num/tl2fit.f, revision 1.1.1.1

1.1       root        1:       subroutine tl2fit(nx,ny,nn,d,w,c)
                      2: c     planar l2 fit  - tensor grid version
                      3: c     input:
                      4: c            nx = number of points in x
                      5: c            ny = number of points in y
                      6: c            nn = nx*ny
                      7: c             d = data points
                      8: c                 on output,  d(*,*) contains residuals
                      9: c             w = workspace, of size 6*n*n
                     10: c                 internal blocking:
                     11: c                  1..3  x
                     12: c                  4     y
                     13: c                  5     qraux
                     14: c                  6     qty
                     15: c             c = coefficients    c(1) + c(2)*x + c(3)*y
                     16:       integer i, j, k, n
                     17:       real c(3), w(nn,6), d(nx,ny), z
                     18:       integer i,info
                     19:       k=1
                     20:       do 100 i=1,nx
                     21:       do 100 j=1,ny
                     22:         w(k,1)=1
                     23:         w(k,2)=(i-1)/(nx-1.)
                     24:         w(k,3)=(j-1)/(ny-1.)
                     25:         w(k,4)=d(i,j)
                     26:         k=k+1
                     27:  100    continue
                     28:       call sqrdc(w,nn,nn,3,w(1,5),z,z,0)
                     29:       call sqrsl(w,nn,nn,3,w(1,5),w(1,4),z,w(1,6),c,z,z,100,info)
                     30:       if(info.ne.0)then
                     31:         write(6,1001) info
                     32:  1001   format(' sqrsl abort. info=',i5)
                     33:         stop
                     34:         end if
                     35:       end

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.