evalshapefunc.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	evalshapefunc (xil, etl, xl2, nopes, p)

Function/Subroutine Documentation

evalshapefunc()

subroutine evalshapefunc	(	real*8	xil,
		real*8	etl,
		real*8, dimension(3,8)	xl2,
		integer	nopes,
		real*8, dimension(3)	p
	)

!     
      implicit none
!     
      integer i,j,nopes,iflag
!     
      real*8  xl2(3,8),xs2(3,2),shp2(7,8),p(3),xsj2(3),
     &  xil,etl
!     
      iflag=1
      do j=1,3
         p(j)=0.0
      enddo 
      if(nopes.eq.8)then
         call shape8q(xil,etl,xl2,xsj2,xs2,shp2,iflag)
      elseif(nopes.eq.4)then
         call shape4q(xil,etl,xl2,xsj2,xs2,shp2,iflag)
      elseif(nopes.eq.6)then
         call shape6tri(xil,etl,xl2,xsj2,xs2,shp2,iflag)
      else
         call shape3tri(xil,etl,xl2,xsj2,xs2,shp2,iflag)
      endif 
      do i=1,nopes
         do j=1,3
            p(j)=p(j)+xl2(j,i)*shp2(4,i)
         enddo
      enddo 
!     
      return