plane_eq.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	plane_eq (x1, y1, z1, x2, y2, z2, x3, y3, z3, x0, y0, output)

Function/Subroutine Documentation

plane_eq()

subroutine plane_eq	(	real*8	x1,
		real*8	y1,
		real*8	z1,
		real*8	x2,
		real*8	y2,
		real*8	z2,
		real*8	x3,
		real*8	y3,
		real*8	z3,
		real*8	x0,
		real*8	y0,
		real*8	output
	)

!
      implicit none
!
      real*8 x1,y1,z1,x2,y2,z2,x3,y3,z3,x0,y0,output,
     &  a,b,c,d
!
      d=x1*y2*z3+y1*z2*x3+z1*x2*y3-x1*z2*y3-y1*x2*z3-z1*y2*x3
      if(d.ne.0.d0) then
         a=1.d0/d*(y2*z3+y1*z2+z1*y3-z2*y3-y1*z3-z1*y2)
      endif  
      if(d.ne.0.d0) then
         b=1.d0/d*(x1*z3+z2*x3+z1*x2-x1*z2-x2*z3-z1*x3)
      endif  
      if(d.ne.0.d0) then
         c=1.d0/d*(x1*y2+y1*x3+x2*y3-x1*y3-y1*x2-y2*x3)
      endif  
      if(d.ne.0.d0) then
         output=1.d0/c*(1.d0-a*x0-b*y0)
      else
         output=0.d0
      endif
      return