i8_nearest_seg.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i8_nearest_seg (x, is, m1, m2, m3, bmin)

Function/Subroutine Documentation

i8_nearest_seg()

subroutine i8_nearest_seg	(	dimension(3,*), intent(in)	x,
		integer, intent(in)	is,
		integer, intent(in)	m1,
		integer, intent(in)	m2,
		integer, intent(in)	m3,
		intent(out)	bmin )

Definition at line 28 of file i8_nearest_seg.F.

C  
C  COMPUTE BARYCENTRIC COORDINATES OF PROJECTION OF NODE IS ON TRIANGLE M1 M2 M3
C  WOLFGANG HEIDRICH, 2005
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
C-----------------------------------------------
      INTEGER, INTENT(IN)  :: IS, M1, M2, M3
      my_real, INTENT(IN)  :: x(3,*)
      my_real, INTENT(OUT) :: bmin
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real n(3), u(3), v(3), w(3), uw(3), wv(3)
      my_real b0, b1, b2, fac
C---------------------------------------------------
      u(1)=x(1,m2)-x(1,m1)
      u(2)=x(2,m2)-x(2,m1)
      u(3)=x(3,m2)-x(3,m1)
C
      v(1)=x(1,m3)-x(1,m1)
      v(2)=x(2,m3)-x(2,m1)
      v(3)=x(3,m3)-x(3,m1)
C
      n(1)=u(2)*v(3)-u(3)*v(2)
      n(2)=u(3)*v(1)-u(1)*v(3)
      n(3)=u(1)*v(2)-u(2)*v(1)
C
      w(1)=x(1,is)-x(1,m1)
      w(2)=x(2,is)-x(2,m1)
      w(3)=x(3,is)-x(3,m1)
C
      uw(1)=u(2)*w(3)-u(3)*w(2)
      uw(2)=u(3)*w(1)-u(1)*w(3)
      uw(3)=u(1)*w(2)-u(2)*w(1)
C
      wv(1)=w(2)*v(3)-w(3)*v(2)
      wv(2)=w(3)*v(1)-w(1)*v(3)
      wv(3)=w(1)*v(2)-w(2)*v(1)
C
      fac=one/(n(1)*n(1)+n(2)*n(2)+n(3)*n(3))
      b2 =(uw(1)*n(1)+uw(2)*n(2)+uw(3)*n(3))*fac
      b1 =(wv(1)*n(1)+wv(2)*n(2)+wv(3)*n(3))*fac
      b0 = one - b1 - b2
      bmin=min(b0,b1,b2)    
C BMIN > 0 Projection est a l'interieur du triangle
C BMIN = 0 Projection est sur un cote du triangle
C BMIN < 0 Projection est a l'exterieur du triangle
      RETURN