updatecont.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	updatecont (koncont, ncont, co, vold, cg, straight, mi)

Function/Subroutine Documentation

updatecont()

subroutine updatecont	(	integer, dimension(4,*)	koncont,
		integer	ncont,
		real*8, dimension(3,*)	co,
		real*8, dimension(0:mi(2),*)	vold,
		real*8, dimension(3,*)	cg,
		real*8, dimension(16,*)	straight,
		integer, dimension(*)	mi
	)

!
!     update geometric date of the contact master surface triangulation
!
      implicit none
!
      integer koncont(4,*),ncont,i,j,k,node,mi(*)
!
      real*8 co(3,*),vold(0:mi(2),*),cg(3,*),straight(16,*),col(3,3)
!
      do i=1,ncont
         do j=1,3
            node=koncont(j,i)
            do k=1,3
               col(k,j)=co(k,node)+vold(k,node)
            enddo
         enddo
!
!        center of gravity of the triangles
!
         do k=1,3
            cg(k,i)=col(k,1)
         enddo
         do j=2,3
            do k=1,3
               cg(k,i)=cg(k,i)+col(k,j)
            enddo
         enddo
         do k=1,3
            cg(k,i)=cg(k,i)/3.d0
         enddo
!
!        calculating the equation of the triangle plane and the planes
!        perpendicular on it and through the triangle edges
!
         call straighteq3d(col,straight(1,i))
!
      enddo
!
      return