Return to l2fit.f CVS log

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

researchv10 Norman

      subroutine l2fit(n,nn,ld,d,w,c)
c     planar l2 fit
c     input:
c             n = number of points in each direction
c            nn = n**2
c            ld = declared dimension of d
c             d = data points
c                 on output,  d(3,*,*) contains residuals
c             w = workspace, of size 6*n*n
c                 internal blocking:
c                  1..3  x
c                  4     y
c                  5     qraux
c                  6     qty
c             c = coefficients    c(1) + c(2)*x + c(3)*y
      integer i, j, k, n
      real c(3), w(nn,6), d(3,ld,ld), z
      integer i,info
      k=1
      do 100 i=1,n
      do 100 j=1,n
        w(k,1)=1
        w(k,2)=d(1,i,j)
        w(k,3)=d(2,i,j)
        w(k,4)=d(3,i,j)
        k=k+1
 100    continue
      call sqrdc(w,nn,nn,3,w(1,5),z,z,0)
      call sqrsl(w,nn,nn,3,w(1,5),w(1,4),z,w(1,6),c,z,z,100,info)
      if(info.ne.0)then
        write(6,1001) info
 1001   format(' sqrsl abort. info=',i5)
        stop
        end if
      end

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