Lines Matching refs:BDiagl
29 dimension BDiagl(numel,nshape,nflow,nflow),
53 call local (BDiag, BDiagl(iel:inum,:,:,:), abs(mien(iblk)%p),
70 & - BDiagl(n+iv,inode,2,1) * EGmass(n+iv,i+1,j))
73 & - BDiagl(n+iv,inode,3,1) * EGmass(n+iv,i+1,j)
74 & - BDiagl(n+iv,inode,3,2) * EGmass(n+iv,i+2,j))
77 & - BDiagl(n+iv,inode,4,1) * EGmass(n+iv,i+1,j)
78 & - BDiagl(n+iv,inode,4,2) * EGmass(n+iv,i+2,j)
79 & - BDiagl(n+iv,inode,4,3) * EGmass(n+iv,i+3,j))
82 & - BDiagl(n+iv,inode,5,1) * EGmass(n+iv,i+1,j)
83 & - BDiagl(n+iv,inode,5,2) * EGmass(n+iv,i+2,j)
84 & - BDiagl(n+iv,inode,5,3) * EGmass(n+iv,i+3,j)
85 & - BDiagl(n+iv,inode,5,4) * EGmass(n+iv,i+4,j))
100 EGmass(n+iv,j,i+1) = BDiagl(n+iv,inode,1,1) *
103 EGmass(n+iv,j,i+2) = BDiagl(n+iv,inode,2,2) * (
105 & - BDiagl(n+iv,inode,1,2) * EGmass(n+iv,j,i+1))
107 EGmass(n+iv,j,i+3) = BDiagl(n+iv,inode,3,3) * (
109 & - BDiagl(n+iv,inode,1,3) * EGmass(n+iv,j,i+1)
110 & - BDiagl(n+iv,inode,2,3) * EGmass(n+iv,j,i+2))
112 EGmass(n+iv,j,i+4) = BDiagl(n+iv,inode,4,4) * (
114 & - BDiagl(n+iv,inode,1,4) * EGmass(n+iv,j,i+1)
115 & - BDiagl(n+iv,inode,2,4) * EGmass(n+iv,j,i+2)
116 & - BDiagl(n+iv,inode,3,4) * EGmass(n+iv,j,i+3))
118 EGmass(n+iv,j,i+5) = BDiagl(n+iv,inode,5,5) * (
120 & - BDiagl(n+iv,inode,1,5) * EGmass(n+iv,j,i+1)
121 & - BDiagl(n+iv,inode,2,5) * EGmass(n+iv,j,i+2)
122 & - BDiagl(n+iv,inode,3,5) * EGmass(n+iv,j,i+3)
123 & - BDiagl(n+iv,inode,4,5) * EGmass(n+iv,j,i+4))