Lines Matching refs:shp
2 subroutine trishp(el,xl,ndm,iord, xsj,shp) argument
37 real*8 el(3),xl(ndm,*), shp(3,*) local
121 shp(1,1) = (y2-y3)*xsjr
122 shp(2,1) = (x3-x2)*xsjr
123 shp(3,1) = el(1)
125 shp(1,2) = (y3-y1)*xsjr
126 shp(2,2) = (x1-x3)*xsjr
127 shp(3,2) = el(2)
129 shp(1,3) = (y1-y2)*xsjr
130 shp(2,3) = (x2-x1)*xsjr
131 shp(3,3) = el(3)
137 shp(1,4) = shp(1,1)*fel2 + shp(1,2)*fel1
138 shp(2,4) = shp(2,1)*fel2 + shp(2,2)*fel1
139 shp(3,4) = el(1)*fel2
141 shp(1,5) = shp(1,2)*fel3 + shp(1,3)*fel2
142 shp(2,5) = shp(2,2)*fel3 + shp(2,3)*fel2
143 shp(3,5) = el(2)*fel3
145 shp(1,6) = shp(1,3)*fel1 + shp(1,1)*fel3
146 shp(2,6) = shp(2,3)*fel1 + shp(2,1)*fel3
147 shp(3,6) = el(3)*fel1
152 shp(1,7) = shp(1,1)*fel23 + shp(1,2)*fel31 + shp(1,3)*fel12
153 shp(2,7) = shp(2,1)*fel23 + shp(2,2)*fel31 + shp(2,3)*fel12
154 shp(3,7) = fel12*el(3)
156 shp(1,7) = (shp(1,1)*fel23*2.d0
157 & + shp(1,2)*fel31
158 & + shp(1,3)*fel12)*el(1)
159 shp(2,7) = (shp(2,1)*fel23*2.d0
160 & + shp(2,2)*fel31
161 & + shp(2,3)*fel12)*el(1)
162 shp(3,7) = fel12*fel31
164 shp(1,8) = (shp(1,1)*fel23
165 & + shp(1,2)*fel31*2.d0
166 & + shp(1,3)*fel12)*el(2)
167 shp(2,8) = (shp(2,1)*fel23
168 & + shp(2,2)*fel31*2.d0
169 & + shp(2,3)*fel12)*el(2)
170 shp(3,8) = fel12*fel23
172 shp(1,9) = (shp(1,1)*fel23
173 & + shp(1,2)*fel31
174 & + shp(1,3)*fel12*2.d0)*el(3)
175 shp(2,9) = (shp(2,1)*fel23
176 & + shp(2,2)*fel31
177 & + shp(2,3)*fel12*2.d0)*el(3)
178 shp(3,9) = fel31*fel23
189 shp(i,1) = shp(i,1) - one3*shp(i,7)
190 shp(i,2) = shp(i,2) - one3*shp(i,7)
191 shp(i,3) = shp(i,3) - one3*shp(i,7)
192 shp(i,4) = shp(i,4) - four9*shp(i,7)
193 shp(i,5) = shp(i,5) - four9*shp(i,7)
194 shp(i,6) = shp(i,6) - four9*shp(i,7)
201 shp(i,1) = shp(i,1) - 0.5d0*(shp(i,4) + shp(i,6))
202 shp(i,2) = shp(i,2) - 0.5d0*(shp(i,5) + shp(i,4))
203 shp(i,3) = shp(i,3) - 0.5d0*(shp(i,6) + shp(i,5))