Lines Matching refs:shp

33       subroutine shape8qtilde(xi,et,xl,xsj,xs,shp,iflag)  argument
56       real*8 shp(7,8),xs(3,7),xsi(2,3),xl(3,8),sh(3),xsj(3),xi,  local
77       shp(4,1)=1.0*shpold(4,1)+alpha*shpold(4,5)+alpha*shpold(4,8)
78       shp(4,2)=1.0*shpold(4,2)+alpha*shpold(4,5)+alpha*shpold(4,6)
79       shp(4,3)=1.0*shpold(4,3)+alpha*shpold(4,6)+alpha*shpold(4,7)
80       shp(4,4)=1.0*shpold(4,4)+alpha*shpold(4,7)+alpha*shpold(4,8)
81       shp(4,5)=(1.0-2.0*alpha)*shpold(4,5)
82       shp(4,6)=(1.0-2.0*alpha)*shpold(4,6)
83       shp(4,7)=(1.0-2.0*alpha)*shpold(4,7)
84       shp(4,8)=(1.0-2.0*alpha)*shpold(4,8)
94       shp(1,1)=(1.d0-et)*(2.d0*xi+et)/4.d0
95       shp(1,2)=(1.d0-et)*(2.d0*xi-et)/4.d0
96       shp(1,3)=(1.d0+et)*(2.d0*xi+et)/4.d0
97       shp(1,4)=(1.d0+et)*(2.d0*xi-et)/4.d0
98       shp(1,5)=-xi*(1.d0-et)
99       shp(1,6)=(1.d0-et*et)/2.d0
100       shp(1,7)=-xi*(1.d0+et)
101       shp(1,8)=-(1.d0-et*et)/2.d0
105       shp(2,1)=(1.d0-xi)*(2.d0*et+xi)/4.d0
106       shp(2,2)=(1.d0+xi)*(2.d0*et-xi)/4.d0
107       shp(2,3)=(1.d0+xi)*(2.d0*et+xi)/4.d0
108       shp(2,4)=(1.d0-xi)*(2.d0*et-xi)/4.d0
109       shp(2,5)=-(1.d0-xi*xi)/2.d0
110       shp(2,6)=-et*(1.d0+xi)
111       shp(2,7)=(1.d0-xi*xi)/2.d0
112       shp(2,8)=-et*(1.d0-xi)
121                xs(i,j)=xs(i,j)+xl(i,k)*shp(j,k)
177                sh(j)=shp(1,k)*xsi(1,j)+shp(2,k)*xsi(2,j)
180                shp(j,k)=sh(j)
188          shp(5,1)=(1.d0-et)/2.d0
189          shp(5,2)=(1.d0-et)/2.d0
190          shp(5,3)=(1.d0+et)/2.d0
191          shp(5,4)=(1.d0+et)/2.d0
192          shp(5,5)=-(1.d0-et)
193          shp(5,6)=0.d0
194          shp(5,7)=-(1.d0+et)
195          shp(5,8)=0.d0
197          shp(5,1)=1.0*shpold(5,1)+alpha*shpold(5,5)+alpha*shpold(5,8)
198          shp(5,2)=1.0*shpold(5,2)+alpha*shpold(5,5)+alpha*shpold(5,6)
199          shp(5,3)=1.0*shpold(5,3)+alpha*shpold(5,6)+alpha*shpold(5,7)
200          shp(5,4)=1.0*shpold(5,4)+alpha*shpold(5,7)+alpha*shpold(5,8)
201          shp(5,5)=(1.0-2.0*alpha)*shpold(5,5)
202          shp(5,6)=(1.0-2.0*alpha)*shpold(5,6)
203          shp(5,7)=(1.0-2.0*alpha)*shpold(5,7)
204          shp(5,8)=(1.0-2.0*alpha)*shpold(5,8)
208          shp(6,1)=(1.d0-2.d0*(xi+et))/4.d0
209          shp(6,2)=(-1.d0-2.d0*(xi-et))/4.d0
210          shp(6,3)=(1.d0+2.d0*(xi+et))/4.d0
211          shp(6,4)=(-1.d0-2.d0*(xi+et))/4.d0
212          shp(6,5)=xi
213          shp(6,6)=-et
214          shp(6,7)=-xi
215          shp(6,8)=et
217          shp(6,1)=1.0*shpold(6,1)+alpha*shpold(6,5)+alpha*shpold(6,8)
218          shp(6,2)=1.0*shpold(6,2)+alpha*shpold(6,5)+alpha*shpold(6,6)
219          shp(6,3)=1.0*shpold(6,3)+alpha*shpold(6,6)+alpha*shpold(6,7)
220          shp(6,4)=1.0*shpold(6,4)+alpha*shpold(6,7)+alpha*shpold(6,8)
221          shp(6,5)=(1.0-2.0*alpha)*shpold(6,5)
222          shp(6,6)=(1.0-2.0*alpha)*shpold(6,6)
223          shp(6,7)=(1.0-2.0*alpha)*shpold(6,7)
224          shp(6,8)=(1.0-2.0*alpha)*shpold(6,8)
228          shp(7,1)=(1.d0-xi)/2.d0
229          shp(7,2)=(1.d0+xi)/2.d0
230          shp(7,3)=(1.d0+xi)/2.d0
231          shp(7,4)=(1.d0-xi)/2.d0
232          shp(7,5)=0.d0
233          shp(7,6)=-(1.d0+xi)
234          shp(7,7)=0.d0
235          shp(7,8)=-(1.d0-xi)
237          shp(7,1)=1.0*shpold(7,1)+alpha*shpold(7,5)+alpha*shpold(7,8)
238          shp(7,2)=1.0*shpold(7,2)+alpha*shpold(7,5)+alpha*shpold(7,6)
239          shp(7,3)=1.0*shpold(7,3)+alpha*shpold(7,6)+alpha*shpold(7,7)
240          shp(7,4)=1.0*shpold(7,4)+alpha*shpold(7,7)+alpha*shpold(7,8)
241          shp(7,5)=(1.0-2.0*alpha)*shpold(7,5)
242          shp(7,6)=(1.0-2.0*alpha)*shpold(7,6)
243          shp(7,7)=(1.0-2.0*alpha)*shpold(7,7)
244          shp(7,8)=(1.0-2.0*alpha)*shpold(7,8)
253                   xs(i,j)=xs(i,j)+xl(i,k)*shp(j,k)