s8slke3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s8slke3 (kl, hh, vol, dn_y, dn_z, ksi, eta, zeta, dn_r, dn_s, dn_t, invj, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, a11, a12, a13, a21, a22, a23, a31, a32, a33, nel)

Function/Subroutine Documentation

s8slke3()

subroutine s8slke3	(	double precision, dimension(24,24,mvsiz)	kl,
		double precision, dimension(2,mvsiz)	hh,
			vol,
		double precision, dimension(mvsiz,8)	dn_y,
		double precision, dimension(mvsiz,8)	dn_z,
			ksi,
			eta,
			zeta,
		double precision, dimension(8)	dn_r,
		double precision, dimension(8)	dn_s,
		double precision, dimension(8)	dn_t,
		double precision, dimension(9,mvsiz)	invj,
			x1,
			x2,
			x3,
			x4,
			x5,
			x6,
			x7,
			x8,
			y1,
			y2,
			y3,
			y4,
			y5,
			y6,
			y7,
			y8,
			z1,
			z2,
			z3,
			z4,
			z5,
			z6,
			z7,
			z8,
			a11,
			a12,
			a13,
			a21,
			a22,
			a23,
			a31,
			a32,
			a33,
		integer, intent(in)	nel )

Definition at line 31 of file s8slke3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   vol(mvsiz),ksi,eta,zeta
      double precision
     .   kl(24,24,mvsiz), hh(2,mvsiz), 
     .   dn_x(mvsiz,8),dn_y(mvsiz,8),dn_z(mvsiz,8),ktemp(6,24),b(6,24),
     .   hv(6,6),dn_r(8),dn_s(8),dn_t(8),invj(9,mvsiz)!,KLL(24,24)
      my_real
     .   x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
     .   x5(mvsiz), x6(mvsiz), x7(mvsiz), x8(mvsiz),
     .   y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
     .   y5(mvsiz), y6(mvsiz), y7(mvsiz), y8(mvsiz),
     .   z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
     .   z5(mvsiz), z6(mvsiz), z7(mvsiz), z8(mvsiz) 
       my_real
     .   a11(mvsiz), a12(mvsiz), a13(mvsiz), 
     .   a21(mvsiz), a22(mvsiz), a23(mvsiz), 
     .   a31(mvsiz), a32(mvsiz), a33(mvsiz)      
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
#if defined(MYREAL8) && !defined(WITHOUT_LINALG)
      INTEGER I, J,JJ,K,JK, IS
      b(1:6,1:24) = zero
      hv(1:6,1:6) = zero
 
      DO i=1,nel
! Full integration
       !DO J=1,8
       ! JJ=3*(J-1)
       ! B(1,JJ+1) = DN_X(I,J)
       ! B(2,JJ+2) = DN_Y(I,J)
       ! B(3,JJ+3) = DN_Z(I,J)
       ! B(4,JJ+1) = DN_Y(I,J)
       ! B(4,JJ+2) = DN_X(I,J)
       ! B(5,JJ+1) = DN_Z(I,J)
       ! B(5,JJ+3) = DN_X(I,J)
       ! B(6,JJ+2) = DN_Z(I,J)
       ! B(6,JJ+3) = DN_Y(I,J)
       !ENDDO
! ANS formulation
       CALL s8sansb(
     .    x1(i),x2(i),x3(i),x4(i),x5(i),x6(i),x7(i),x8(i),
     .    y1(i),y2(i),y3(i),y4(i),y5(i),y6(i),y7(i),y8(i),                                                                
     .    z1(i),z2(i),z3(i),z4(i),z5(i),z6(i),z7(i),z8(i), 
     .    invj(1,i),ksi,eta,zeta,
     .    a11(i),a12(i),a13(i),
     .    a21(i),a22(i),a23(i),
     .    a31(i),a32(i),a33(i),
     .    dn_r, dn_s, dn_t,b)
       
       hv(1,1) = (hh(1,i)+two*hh(2,i))*vol(i)
       hv(2,1) = hh(1,i)*vol(i)
       hv(3,1) = hh(1,i)*vol(i)
       hv(1,2) = hh(1,i)*vol(i)
       hv(2,2) = (hh(1,i)+two*hh(2,i))*vol(i)
       hv(3,2) = hh(1,i)*vol(i)
       hv(1,3) = hh(1,i)*vol(i)
       hv(2,3) = hh(1,i)*vol(i)
       hv(3,3) = (hh(1,i)+two*hh(2,i))*vol(i)
       hv(4,4) = hh(2,i)*vol(i)
       hv(5,5) = hh(2,i)*vol(i)
       hv(6,6) = hh(2,i)*vol(i)
       
       ktemp(1:6,1:24)=zero
       !DO J=1,6
       ! DO K=1,24
       !  DO JK=1,6
       !   KTEMP(J,K) = KTEMP(J,K)+HV(J,JK)*B(JK,K)
       !  ENDDO
       ! ENDDO
       !ENDDO
 
       ! CALL DGEMM('N','N',6,24,6,ONE,HV,6,B,6,ZERO,KTEMP,6)
        CALL dsymm('L','U',6,24,one,hv,6,b,6,zero,ktemp,6)
       !DO J=1,24
       ! DO K=1,24
       !  DO JK=1,6
       !   KL(J,K,I) = KL(J,K,I)+B(JK,J)*KTEMP(JK,K)
       !  ENDDO
       ! ENDDO
       !ENDDO
       CALL dgemm('T','N',24,24,6,one,b,6,ktemp,6,one,kl(1,1,i),24)
      ENDDO
#else
      CALL arret(5)
#endif
      RETURN