filter.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	filter (dgdxglob, nobject, nk, nodedesi, ndesi, objectset, xo, yo, zo, x, y, z, nx, ny, nz, neighbor, r, ndesia, ndesib)

Function/Subroutine Documentation

filter()

subroutine filter	(	real*8, dimension(2,nk,nobject)	dgdxglob,
		integer	nobject,
		integer	nk,
		integer, dimension(*)	nodedesi,
		integer	ndesi,
		character*81, dimension(4,*)	objectset,
		real*8, dimension(ndesi)	xo,
		real*8, dimension(ndesi)	yo,
		real*8, dimension(ndesi)	zo,
		real*8, dimension(ndesi)	x,
		real*8, dimension(ndesi)	y,
		real*8, dimension(ndesi)	z,
		integer, dimension(ndesi)	nx,
		integer, dimension(ndesi)	ny,
		integer, dimension(ndesi)	nz,
		integer, dimension(ndesi+6)	neighbor,
		real*8, dimension(ndesi+6)	r,
		integer	ndesia,
		integer	ndesib
	)

!
!     Filtering of sensitivities      
!
      implicit none
!
      character*81 objectset(4,*)

      integer nobject,nk,nodedesi(*),nnodesinside,i,
     &        ndesi,j,m,neighbor(ndesi+6),nx(ndesi),
     &        ny(ndesi),nz(ndesi),istat,ndesia,ndesib
!
      real*8 dgdxglob(2,nk,nobject),xo(ndesi),yo(ndesi),zo(ndesi),
     &       x(ndesi),y(ndesi),z(ndesi),filterrad,r(ndesi+6),
     &       filterval(ndesi),nominator,denominator
!
!     Calculate filtered sensitivities
!   
!     Assign filter radius (taken from first defined object function)
!
      read(objectset(2,1)(21:40),'(f20.0)',iostat=istat) filterrad     
!        
      do j=ndesia,ndesib
!     
         call near3d_se(xo,yo,zo,x,y,z,nx,ny,nz,xo(j),yo(j),zo(j),
     &        ndesi,neighbor,r,nnodesinside,filterrad)
!  
!        Calculate function value of the filterfunction (linear) 
!
         if(objectset(2,1)(1:6).eq.'LINEAR') then
          do i=1,nnodesinside
               filterval(i)=filterrad-dsqrt(r(i))
            enddo
!
!        Calculate function value of the filterfunction (gauss) 
!
         elseif(objectset(2,1)(1:5).eq.'GAUSS') then
            do i=1,nnodesinside
               filterval(i)=dexp(-(2*r(i)/filterrad**2))
            enddo 
!
!        Calculate function value of the filterfunction (cubic) 
!
         elseif(objectset(2,1)(1:5).eq.'CUBIC') then
            do i=1,nnodesinside
               filterval(i)=filterrad-(1/filterrad**2)*dsqrt(r(i))**3
            enddo
         endif
!  
!        Calculate filtered sensitivity
 
         do m=1,nobject
            if(objectset(1,m)(1:9).eq.'THICKNESS') cycle             
            nominator=0.d0
            denominator=0.d0
            do i=1,nnodesinside
               nominator=nominator+filterval(i)*
     &             dgdxglob(1,nodedesi(neighbor(i)),m)
               denominator=denominator+filterval(i)
            enddo 
            dgdxglob(2,nodedesi(j),m)=nominator/denominator
         enddo
      enddo
!
      return