weight.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	weight0 (xi, yi, zi, xj, yj, zj, h, w)
subroutine	weight1 (xi, yi, zi, xj, yj, zj, h, w, wgrad)

Function/Subroutine Documentation

weight0()

subroutine weight0	(		xi,
			yi,
			zi,
			xj,
			yj,
			zj,
			h,
			w )

Definition at line 33 of file weight.F.

C-------------------------------------------------------------------------
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C----------------------------------------------------------
C   D u m m y   A r g u m e n t s   a n d   F u n c t i o n
C----------------------------------------------------------
      my_real
     .   xi,yi,zi,xj,yj,zj,h,w
C----------------------------------------------------------
      my_real
     .   r,rh,rh2,rh3,rhm,rhm3
C=======================================================================
      IF (h==zero) stop
C-----
      r =sqrt((xi-xj)*(xi-xj)+(yi-yj)*(yi-yj)+(zi-zj)*(zi-zj))
C-----
      IF (r<=h) THEN
        rh  =r/h
        rh2 =rh*rh
        rh3 =rh2*rh
        w   =(one - three_half*rh2+three_over_4*rh3)*invpi/(h*h*h)
      ELSEIF (r<=two*h) THEN
        rhm =two-r/h
        rhm =rhm/h
        rhm3=rhm*rhm*rhm
        w   =fourth*rhm3*invpi
      ELSE
        w   =zero
      ENDIF
      RETURN

weight1()

subroutine weight1	(		xi,
			yi,
			zi,
			xj,
			yj,
			zj,
			h,
			w,
			wgrad )

Definition at line 78 of file weight.F.

C-------------------------------------------------------------------------
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C----------------------------------------------------------
C   D u m m y   A r g u m e n t s   a n d   F u n c t i o n
C----------------------------------------------------------
      my_real
     .   xi,yi,zi,xj,yj,zj,h,w,wgrad(3)
C----------------------------------------------------------
      my_real
     .   r,rh,rh2,rh3,ih3,rhm,rhm2,rhm3,wprimr
C=======================================================================
      IF (h==zero) stop
C-----
      r =sqrt((xi-xj)*(xi-xj)+(yi-yj)*(yi-yj)+(zi-zj)*(zi-zj))
C-----
      IF (r<=h) THEN
        rh  =r/h
        rh2 =rh*rh
        rh3 =rh2*rh
        ih3 =invpi/(h*h*h)
        w      =(one - three_half*rh2+three_over_4*rh3)*ih3
        wprimr =(-three+twop25*rh)*ih3/(h*h)
      ELSEIF (r<=two*h) THEN
        rhm =two - r/h
        rhm =rhm/h
        rhm2=rhm*rhm*invpi
        rhm3=rhm2*rhm
        w      = fourth*rhm3
        wprimr =-three_over_4*rhm2/(h*h*r)
      ELSE
        w      =zero
        wprimr =zero
      ENDIF
      wgrad(1)=wprimr*(xi-xj)
      wgrad(2)=wprimr*(yi-yj)
      wgrad(3)=wprimr*(zi-zj)
      RETURN