moment.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	moment (intrect, badpts, npt, usage)
	Compute enclosing rectangle around flagged points. More...

Function/Subroutine Documentation

moment()

subroutine moment	(	dimension(nsize)	intrect,
		dimension(2,npt)	badpts,
			npt,
		save	usage
	)

Compute enclosing rectangle around flagged points.

save some info., even tho. usage might be low and rect. scrapped.

Parameters

[in]	badpts	x,y coords of cells around which a rectangle is to be made and compute efficiency on
[in]	npt	num. of badpts. in the cluster.
[out]	usage	ratio of flagged to unflagged badpts. in new grid measures goodness of fit and clustering
[out]	intrect	stores some info. for grid created herein. sometimes rect = rnode, sometimes = temp. array. sometimes intrect = node. depending on calling prog. (grdfit or expand)

Definition at line 17 of file moment.f.

References grdfit(), amr_module::ndihi, amr_module::ndilo, amr_module::ndjhi, amr_module::ndjlo, amr_module::node, amr_module::nsize, and amr_module::rnode.

Referenced by grdfit(), and smartbis().

c
      use amr_module
      implicit double precision (a-h,o-z)


      dimension     intrect(nsize),badpts(2,npt)
c
c :::::::::::::::::::::::: MOMENT ::::::::::::::::::::::::::::::::::
c  moment = compute enclosing rectangle around flagged points.
c  save some info., even tho. usage might be low and rect. scrapped.

c input parameters:
c     badpts      = x,y coords of flagged badpts grouped into clusters
c                   are in the first two rows
c     npt         = num. of badpts. in the cluster.
c
c output parameters:
c     usage       = ratio of flagged to unflagged badpts. in new grid
c                   measures goodness of fit and clustering
c    intrect( )    = stores some info. for grid created herein.
c                   sometimes rect = rnode, sometimes = temp. array.
c                   sometimes intrect = node.
c                   depending on calling prog. (grdfit or expand)
c
c
c :::::::::::::::::::::::: MOMENT ::::::::::::::::::::::::::::::::::
c
      rn = dble(npt)
c
c compute length of enclosing rectangles to include all flagged badpts.
c
      emx1 = badpts(1,1)
      emn1 = emx1
      emx2 = badpts(2,1)
      emn2 = emx2
      do 80 ipt = 1, npt
          if (badpts(1,ipt) .gt. emx1) emx1 = badpts(1,ipt)
          if (badpts(1,ipt) .lt. emn1) emn1 = badpts(1,ipt)
          if (badpts(2,ipt) .gt. emx2) emx2 = badpts(2,ipt)
          if (badpts(2,ipt) .lt. emn2) emn2 = badpts(2,ipt)
 80   continue
c
c from length of the sides determine rect. corners.
c transform to cell numbers (subtract .5)
c
      intrect(ndilo) = nint(emn1 - .5)
      intrect(ndjlo) = nint(emn2 - .5)
      intrect(ndihi) = nint(emx1 - .5)
      intrect(ndjhi) = nint(emx2 - .5)
c
c compute usage
c
      iside1 = intrect(ndihi) - intrect(ndilo) + 1
      iside2 = intrect(ndjhi) - intrect(ndjlo) + 1
      gpall  = iside1 * iside2
      usage  = rn / gpall
c
      return

Here is the call graph for this function:

Here is the caller graph for this function: