MAKE_MASS_MATRIX.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_mass_matrix (nn, xq, wq, dd, rho, aa11, aa12, aa13, aa21, aa22, aa23, aa31, aa32, aa33, bb11, bb12, bb13, bb21, bb22, bb23, bb31, bb32, bb33, gg1, gg2, gg3, dd1, dd2, dd3, mass)
	Makes mass matrix.

Function/Subroutine Documentation

make_mass_matrix()

subroutine make_mass_matrix	(	integer*4	nn,
		real*8, dimension(nn)	xq,
		real*8, dimension(nn)	wq,
		real*8, dimension(nn,nn)	dd,
		real*8, dimension(nn,nn,nn)	rho,
		real*8	aa11,
		real*8	aa12,
		real*8	aa13,
		real*8	aa21,
		real*8	aa22,
		real*8	aa23,
		real*8	aa31,
		real*8	aa32,
		real*8	aa33,
		real*8	bb11,
		real*8	bb12,
		real*8	bb13,
		real*8	bb21,
		real*8	bb22,
		real*8	bb23,
		real*8	bb31,
		real*8	bb32,
		real*8	bb33,
		real*8	gg1,
		real*8	gg2,
		real*8	gg3,
		real*8	dd1,
		real*8	dd2,
		real*8	dd3,
		real*8, dimension(nn,nn,nn)	mass
	)

Makes mass matrix.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nn	number of 1-D GLL nodes
[in]	xq	GLL nodes
[in]	wq	GLL weights
[in]	dd	spectral derivatives matrix
[in]	rho	mass density
[in]	A11	costant values for the bilinear map
[in]	A12	costant values for the bilinear map
[in]	A13	costant values for the bilinear map
[in]	A21	costant values for the bilinear map
[in]	A22	costant values for the bilinear map
[in]	A23	costant values for the bilinear map
[in]	A31	costant values for the bilinear map
[in]	A32	costant values for the bilinear map
[in]	A33	costant values for the bilinear map
[in]	B11	costant values for the bilinear map
[in]	B12	costant values for the bilinear map
[in]	B13	costant values for the bilinear map
[in]	B21	costant values for the bilinear map
[in]	B22	costant values for the bilinear map
[in]	B23	costant values for the bilinear map
[in]	B31	costant values for the bilinear map
[in]	B32	costant values for the bilinear map
[in]	B33	costant values for the bilinear map
[in]	GG1	costant values for the bilinear map
[in]	GG2	costant values for the bilinear map
[in]	GG3	costant values for the bilinear map
[in]	DD1	costant values for the bilinear map
[in]	DD2	costant values for the bilinear map
[in]	DD3	costant values for the bilinear map
[out]	mass	diagonal mass matrix

Definition at line 54 of file MAKE_MASS_MATRIX.f90.

      
      
      implicit none
      
      integer*4 :: nn
      integer*4 :: i,j,k
 
      real*8 :: aa11,aa12,aa13,aa21,aa22,aa23,aa31,aa32,aa33
      real*8 :: bb11,bb12,bb13,bb21,bb22,bb23,bb31,bb32,bb33
      real*8 :: gg1,gg2,gg3,dd1,dd2,dd3
      real*8 :: dxdx,dxdy,dxdz,dydx,dydy,dydz,dzdx,dzdy,dzdz,det_j
 
      real*8, dimension(nn) :: xq, wq
 
      real*8, dimension(nn,nn) :: dd
 
      real*8, dimension(nn,nn,nn) :: rho 
      real*8, dimension(nn,nn,nn) :: mass
      
      
!     ELEMENT NODAL MASS CALCULATION
      
      do k = 1,nn
         do j = 1,nn
            do i = 1,nn
               dxdx = aa11 + bb12*xq(k) + bb13*xq(j) + gg1*xq(j)*xq(k)
               dydx = aa21 + bb22*xq(k) + bb23*xq(j) + gg2*xq(j)*xq(k)
               dzdx = aa31 + bb32*xq(k) + bb33*xq(j) + gg3*xq(j)*xq(k)
               
               dxdy = aa12 + bb11*xq(k) + bb13*xq(i) + gg1*xq(k)*xq(i)
               dydy = aa22 + bb21*xq(k) + bb23*xq(i) + gg2*xq(k)*xq(i)
               dzdy = aa32 + bb31*xq(k) + bb33*xq(i) + gg3*xq(k)*xq(i)
               
               dxdz = aa13 + bb11*xq(j) + bb12*xq(i) + gg1*xq(i)*xq(j)
               dydz = aa23 + bb21*xq(j) + bb22*xq(i) + gg2*xq(i)*xq(j)
               dzdz = aa33 + bb31*xq(j) + bb32*xq(i) + gg3*xq(i)*xq(j)
               
               det_j = dxdz * (dydx*dzdy - dzdx*dydy) &
                     - dydz * (dxdx*dzdy - dzdx*dxdy) &
                     + dzdz * (dxdx*dydy - dydx*dxdy)
               
               mass(i,j,k) = rho(i,j,k) * det_j * wq(i)*wq(j)*wq(k)
            enddo
         enddo
      enddo
      
      return