!************************************************************** !* AceGen 1.1 Windows (13 Aug 07) * !* Co. J. Korelc 2007 13 Aug 07 15:04:56* !************************************************************** ! User : USER ! Evaluation time : 18 s Mode : Optimal ! Number of formulae : 236 Method: Automatic ! Subroutine : SKR10 size :6605 ! Total size of Mathematica code : 6605 subexpressions ! Total size of Fortran code : 21026 bytes C subroutine elmt10(d,ul,xl,ix,tl,s,p,ndfe,ndme,nste,isw) implicit none include 'sms.h' logical DEBUG,symmetric character*50 SELEM,datades(2),postdes(0) parameter (DEBUG=.false., # SELEM="ExamplesHypersolidFEAP") integer ix(nen),ndme,ndfe,nste,isw,ngdof,ncdof integer maxdof,dofglobal(8),nnodes,nalldof,ndim double precision xl(ndme,nen),d(*),ul(ndfe,nen,*) double precision s(nste,nste),p(nste),tl(nen),sxd(24) double precision sall(33,33),pall(33),du(24),dh(9) double precision ul0(ndfe,nen),sg(24),sg0(24) integer i,j,jj,ll,ii,k,kk,i1,i2,i3,hlen,icode integer cdata1,cdata2,cdata3,cdata4 double precision w,v(5409),gpost(64,0),npost(8,0) integer ipordl(30) data (ipordl(i),i=1,30)/1,4,3,2,1,2,6,5,1,2,1,5,8,4,1, & 4,8,7,3,4,3,7,6,2,3,7,8,5,6,7/ data (dofglobal(i),i=1,8)/3,3,3,3,3,3,3,3/ 1235 format(i3,">",3g17.9) 1236 format(" >",5g17.9) 1237 format(i3,"> ",24g11.5) 1238 format(i3,"> ",33g11.5) 1239 format(i3,"> ",3g11.5) icode=7 ncdof=9 ngdof=24 nalldof=33 symmetric=.true. maxdof=3 nnodes=8 ndim=3 rdata(RD_SubIterationTolerance)=1d-9 rdata(RD_Multiplier)=ttim rdata(RD_MultiplierIncrement)=dt rdata(RD_Time)=ttim rdata(RD_TimeIncrement)=dt idata(ID_MaxMessages)=20 if(isw.ne.1) then if(idata(ID_LastIntCode).ne.icode) then call SMSIntPoints(icode,ngpo,gp) idata(ID_LastIntCode)=icode endif endif If(DEBUG) then hlen=234+9*idata(ID_NoSensParameters) endif cdata1=1 cdata2=10 cdata3=19 cdata4=235 do i=1,ndfe do j=1,nen ul0(i,j)=ul(i,j,1)-ul(i,j,2) enddo enddo idata(ID_Iteration)=niter+1 if(ckiter.ne.idata(ID_Iteration)) then ckiter=idata(ID_Iteration) idata(ID_NoMessages)=0 endif if(nstep.eq.1) then idata(ID_TotalIteration)=niter+1 else idata(ID_TotalIteration)=1000 endif go to(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17, # 18,19,20,21), isw c c.... input record 1 of material properties 1 call dinput(d(1),2) call SMSIntPoints(icode,ngpo,gp) write(*,*) SELEM write(iow,*)SELEM c.....Description of the input data datades(1)="Elastic modulus" datades(2)="Poisson ratio" write(iow,"(10x,f15.5,A3,A50)") # (d(i)," = ",datades(i),i=1,2) c.... number of history variables idata(ID_NoSensParameters)=0 nsenpa=idata(ID_NoSensParameters) mct=234+9*idata(ID_NoSensParameters) c.....number of data for TECPLOT ntecdata=6 c.... define node numbering for plot mesh routine, see pltord inord(10) = 30 do ii = 1,30 ipord(ii,10) = ipordl(ii) end do c.... number of projected postprocess quantities istv=0 c.... description of the postprocessing data idata(ID_OutputFile)=iow idata(ID_PostIteration)=0 idata(ID_SkipTangent)=0 idata(ID_SkipResidual)=0 idata(ID_SkipSubIteration)=0; return 2 write(*,*)"User switch 2 not implemented" return 3 continue c.... tangent and residuum if(DEBUG) then write(iow,*)"#################################################" write(iow,*)"TANGENT AND RESIDUAL for ",SELEM write(iow,*)"NODE X" write(iow,1239)(i,(xl(j,i),j=1,ndim),i=1,nnodes) write(iow,*)"NODE ap" write(iow,1235)(i,(ul0(j,i),j=1,maxdof),i=1,nnodes) write(iow,*)"NODE at" write(iow,1235)(i,(ul(j,i,1),j=1,maxdof),i=1,nnodes) write(iow,*)"NODE da" write(iow,1235)(i,(ul(j,i,3),j=1,maxdof),i=1,nnodes) write(iow,*)"ELEMENT hp" write(iow,1236)(hr(nh1+i-1),i=1,hlen) endif if(ncdof.gt.0) then k=0 do i=1,nnodes do j=1,dofglobal(i) k=k+1 du(k)=ul(j,i,3) enddo enddo c..... recover internal modes if(rnmax.gt.0.0d0) then call SMSSRecover(hr(nh2-1+cdata2),hr(nh2-1+cdata3), # du,dh,ngdof,ncdof) do i=1,ncdof hr(nh2-1+cdata1+i-1)=hr(nh2-1+cdata1+i-1)+dh(i) enddo endif call SMSZero(sall,nalldof*nalldof) call SMSZero(pall,nalldof) call SKR10(v,d,ul,ul0,xl,sall,pall, # hr(nh2),hr(nh1)) if(symmetric) then do i=2,nalldof do j=1,i-1 sall(i,j)=sall(j,i) enddo enddo endif call SMSSCondense(sall,s,pall,p,hr(nh2-1+cdata2), # hr(nh2-1+cdata3),ngdof,ncdof) if(DEBUG) then write(iow,*)"STATIC CONDENSATION" write(iow,*)"MODE INCREMENT dh" write(iow,1236)(dh(i),i=1,ncdof) write(iow,*)"FULL MatrixLocal=" write(iow,1238)( i,(s(i,j),j=1,nalldof),i=1,nalldof) write(iow,*)"FULL VectorLocal=" write(iow,1238)1,(p(i),i=1,nalldof) endif else call SMSZero(s,ngdof*ngdof) call SMSZero(p,ngdof) call SKR10(v,d,ul,ul0,xl,s,p, # hr(nh2),hr(nh1)) if(symmetric) then do i=2,nalldof do j=1,i-1 s(i,j)=s(j,i) enddo enddo endif endif if(DEBUG) then write(iow,*)"ELEMENT ht" write(iow,1236)(hr(nh2+i-1),i=1,hlen) write(iow,*)"MatrixLocal=" write(iow,1237)( i,(s(i,j),j=1,ngdof),i=1,ngdof) write(iow,*)"VectorLocal=" write(iow,1237)1,(p(i),i=1,ngdof) endif return 4 continue goto 8 5 write(*,*)"User switch 5 not implemented" return 6 goto 3 7 write(*,*)"User switch 7 not implemented" return c.... postprocessing 8 continue c.....Description of the post-processing data return 9 continue 10 continue 11 continue 12 continue 13 continue c...... initialize history 14 return 15 continue 16 continue 17 continue 18 continue 19 continue write(*,*)"User switch ",isw," not implemented" return c.....sensitivity analysis - external 20 continue return c..... internal sensitivity 21 continue return End !******************* S U B R O U T I N E ********************** SUBROUTINE SKR10(v,d,ul,ul0,xl,s,p,ht,hp) IMPLICIT NONE include 'sms.h' INTEGER i1,i244,i269 DOUBLE PRECISION v(5409),d(2),ul(3,8),ul0(3,8),xl(3,8),s(33,33) &,p(33),ht(*),hp(*) v(321)=-xl(3,1)+xl(3,7) v(320)=-xl(3,4)+xl(3,6) v(319)=-xl(3,2)+xl(3,8) v(318)=-xl(2,1)+xl(2,7) v(317)=-xl(2,4)+xl(2,6) v(316)=-xl(2,2)+xl(2,8) v(315)=-xl(1,1)+xl(1,7) v(314)=-xl(1,4)+xl(1,6) v(313)=-xl(1,2)+xl(1,8) v(312)=2d0*d(2) v(241)=1d0/(1d0+d(2)) v(242)=(d(1)*v(241))/2d0 v(240)=(v(242)*v(312))/(1d0-v(312)) v(187)=v(321)+xl(3,3)-xl(3,5) v(185)=v(318)+xl(2,3)-xl(2,5) v(183)=v(315)+xl(1,3)-xl(1,5) DO i1=1,int(ngpo) v(52)=gp(1,i1) v(63)=1d0-v(52) v(59)=1d0+v(52) v(53)=gp(2,i1) v(69)=1d0+v(53) v(78)=v(69)/8d0 v(86)=-(v(63)*v(78)) v(85)=-(v(59)*v(78)) v(65)=1d0-v(53) v(79)=v(65)/8d0 v(84)=-(v(59)*v(79)) v(82)=-(v(63)*v(79)) v(54)=gp(3,i1) v(70)=1d0+v(54) v(77)=v(70)/8d0 v(88)=-(v(59)*v(77)) v(87)=-(v(63)*v(77)) v(72)=v(69)*v(77) v(67)=v(65)*v(77) v(60)=1d0-v(54) v(80)=v(60)/8d0 v(83)=-(v(59)*v(80)) v(81)=-(v(63)*v(80)) v(62)=v(69)*v(80) v(57)=v(65)*v(80) v(55)=gp(4,i1) v(89)=v(57)*(-xl(1,1)+xl(1,2))+v(62)*(xl(1,3)-xl(1,4))+v(67)*( & -xl(1,5)+xl(1,6))+v(72)*(xl(1,7)-xl(1,8)) v(180)=(v(183)+v(314)+xl(1,2)-xl(1,8))/8d0 v(90)=v(57)*(-xl(2,1)+xl(2,2))+v(62)*(xl(2,3)-xl(2,4))+v(67)*( & -xl(2,5)+xl(2,6))+v(72)*(xl(2,7)-xl(2,8)) v(181)=(v(185)+v(317)+xl(2,2)-xl(2,8))/8d0 v(91)=v(57)*(-xl(3,1)+xl(3,2))+v(62)*(xl(3,3)-xl(3,4))+v(67)*( & -xl(3,5)+xl(3,6))+v(72)*(xl(3,7)-xl(3,8)) v(182)=(v(187)+v(320)+xl(3,2)-xl(3,8))/8d0 v(92)=v(83)*(xl(1,2)-xl(1,3))+v(81)*(xl(1,1)-xl(1,4))+v(88)* & (xl(1,6)-xl(1,7))+v(87)*(xl(1,5)-xl(1,8)) v(184)=(v(183)+v(313)+xl(1,4)-xl(1,6))/8d0 v(93)=v(83)*(xl(2,2)-xl(2,3))+v(81)*(xl(2,1)-xl(2,4))+v(88)* & (xl(2,6)-xl(2,7))+v(87)*(xl(2,5)-xl(2,8)) v(325)=-(v(90)*v(92))+v(89)*v(93) v(186)=(v(185)+v(316)+xl(2,4)-xl(2,6))/8d0 v(222)=-(v(181)*v(184))+v(180)*v(186) v(94)=v(83)*(xl(3,2)-xl(3,3))+v(81)*(xl(3,1)-xl(3,4))+v(88)* & (xl(3,6)-xl(3,7))+v(87)*(xl(3,5)-xl(3,8)) v(324)=v(91)*v(92)-v(89)*v(94) v(323)=-(v(91)*v(93))+v(90)*v(94) v(188)=(v(187)+v(319)+xl(3,4)-xl(3,6))/8d0 v(218)=v(182)*v(184)-v(180)*v(188) v(211)=-(v(182)*v(186))+v(181)*v(188) v(95)=v(82)*(xl(1,1)-xl(1,5))+v(84)*(xl(1,2)-xl(1,6))+v(85)* & (xl(1,3)-xl(1,7))+v(86)*(xl(1,4)-xl(1,8)) v(189)=(v(313)+v(314)+v(315)-xl(1,3)+xl(1,5))/8d0 v(96)=v(82)*(xl(2,1)-xl(2,5))+v(84)*(xl(2,2)-xl(2,6))+v(85)* & (xl(2,3)-xl(2,7))+v(86)*(xl(2,4)-xl(2,8)) v(190)=(v(316)+v(317)+v(318)-xl(2,3)+xl(2,5))/8d0 v(221)=v(181)*v(189)-v(180)*v(190) v(220)=-(v(186)*v(189))+v(184)*v(190) v(97)=v(82)*(xl(3,1)-xl(3,5))+v(84)*(xl(3,2)-xl(3,6))+v(85)* & (xl(3,3)-xl(3,7))+v(86)*(xl(3,4)-xl(3,8)) v(191)=(v(319)+v(320)+v(321)-xl(3,3)+xl(3,5))/8d0 v(216)=-(v(182)*v(189))+v(180)*v(191) v(214)=v(188)*v(189)-v(184)*v(191) v(210)=v(182)*v(190)-v(181)*v(191) v(209)=-(v(188)*v(190))+v(186)*v(191) v(202)=v(189)*v(211)+v(190)*v(218)+v(191)*v(222) v(328)=v(52)/v(202) v(327)=v(53)/v(202) v(326)=v(54)/v(202) v(98)=v(323)*v(95)+v(324)*v(96)+v(325)*v(97) v(330)=v(55)*v(98) v(322)=v(202)/v(98) v(265)=v(218)*v(322) v(264)=v(222)*v(322) v(263)=v(211)*v(322) v(262)=v(216)*v(322) v(261)=v(221)*v(322) v(260)=v(210)*v(322) v(259)=v(214)*v(322) v(258)=v(220)*v(322) v(257)=v(209)*v(322) v(102)=(-(v(94)*v(96))+v(93)*v(97))/v(98) v(165)=v(102)*v(72) v(156)=-(v(102)*v(67)) v(147)=v(102)*v(62) v(138)=-(v(102)*v(57)) v(103)=(v(94)*v(95)-v(92)*v(97))/v(98) v(167)=v(103)*v(72) v(158)=-(v(103)*v(67)) v(149)=v(103)*v(62) v(140)=-(v(103)*v(57)) v(104)=(-(v(93)*v(95))+v(92)*v(96))/v(98) v(169)=v(104)*v(72) v(160)=-(v(104)*v(67)) v(151)=v(104)*v(62) v(142)=-(v(104)*v(57)) v(105)=(v(91)*v(96)-v(90)*v(97))/v(98) v(126)=v(105)*v(80) v(114)=v(105)*v(77) v(106)=(-(v(91)*v(95))+v(89)*v(97))/v(98) v(129)=v(106)*v(80) v(116)=v(106)*v(77) v(107)=(v(90)*v(95)-v(89)*v(96))/v(98) v(132)=v(107)*v(80) v(118)=v(107)*v(77) v(108)=v(323)/v(98) v(127)=v(108)*v(79) v(120)=v(108)*v(78) v(109)=v(324)/v(98) v(130)=v(109)*v(79) v(122)=v(109)*v(78) v(110)=v(325)/v(98) v(133)=v(110)*v(79) v(124)=v(110)*v(78) v(111)=v(114)+v(120) v(112)=v(116)+v(122) v(113)=v(118)+v(124) v(115)=-v(114)+v(127) v(117)=-v(116)+v(130) v(119)=-v(118)+v(133) v(121)=-v(120)+v(126) v(123)=-v(122)+v(129) v(125)=-v(124)+v(132) v(128)=-v(126)-v(127) v(131)=-v(129)-v(130) v(134)=-v(132)-v(133) v(135)=v(138)+v(128)*v(63) v(136)=v(140)+v(131)*v(63) v(137)=v(142)+v(134)*v(63) v(139)=-v(138)+v(128)*v(59) v(141)=-v(140)+v(131)*v(59) v(143)=-v(142)+v(134)*v(59) v(144)=v(147)+v(121)*v(59) v(145)=v(149)+v(123)*v(59) v(146)=v(151)+v(125)*v(59) v(148)=-v(147)+v(121)*v(63) v(150)=-v(149)+v(123)*v(63) v(152)=-v(151)+v(125)*v(63) v(153)=v(156)+v(115)*v(63) v(154)=v(158)+v(117)*v(63) v(155)=v(160)+v(119)*v(63) v(157)=-v(156)+v(115)*v(59) v(159)=-v(158)+v(117)*v(59) v(161)=-v(160)+v(119)*v(59) v(162)=v(165)+v(111)*v(59) v(163)=v(167)+v(112)*v(59) v(164)=v(169)+v(113)*v(59) v(166)=-v(165)+v(111)*v(63) v(168)=-v(167)+v(112)*v(63) v(170)=-v(169)+v(113)*v(63) v(285)=v(265)*v(326) v(284)=v(262)*v(327) v(283)=v(259)*v(328) v(5309)=0d0 v(5310)=v(136) v(5311)=0d0 v(5312)=0d0 v(5313)=v(141) v(5314)=0d0 v(5315)=0d0 v(5316)=v(145) v(5317)=0d0 v(5318)=0d0 v(5319)=v(150) v(5320)=0d0 v(5321)=0d0 v(5322)=v(154) v(5323)=0d0 v(5324)=0d0 v(5325)=v(159) v(5326)=0d0 v(5327)=0d0 v(5328)=v(163) v(5329)=0d0 v(5330)=0d0 v(5331)=v(168) v(5332)=0d0 v(5333)=0d0 v(5334)=0d0 v(5335)=0d0 v(5336)=v(283) v(5337)=v(284) v(5338)=v(285) v(5339)=0d0 v(5340)=0d0 v(5341)=0d0 v(5243)=0d0 v(5244)=0d0 v(5245)=v(136) v(5246)=0d0 v(5247)=0d0 v(5248)=v(141) v(5249)=0d0 v(5250)=0d0 v(5251)=v(145) v(5252)=0d0 v(5253)=0d0 v(5254)=v(150) v(5255)=0d0 v(5256)=0d0 v(5257)=v(154) v(5258)=0d0 v(5259)=0d0 v(5260)=v(159) v(5261)=0d0 v(5262)=0d0 v(5263)=v(163) v(5264)=0d0 v(5265)=0d0 v(5266)=v(168) v(5267)=0d0 v(5268)=0d0 v(5269)=0d0 v(5270)=0d0 v(5271)=0d0 v(5272)=0d0 v(5273)=v(283) v(5274)=v(284) v(5275)=v(285) v(5078)=v(136) v(5079)=0d0 v(5080)=0d0 v(5081)=v(141) v(5082)=0d0 v(5083)=0d0 v(5084)=v(145) v(5085)=0d0 v(5086)=0d0 v(5087)=v(150) v(5088)=0d0 v(5089)=0d0 v(5090)=v(154) v(5091)=0d0 v(5092)=0d0 v(5093)=v(159) v(5094)=0d0 v(5095)=0d0 v(5096)=v(163) v(5097)=0d0 v(5098)=0d0 v(5099)=v(168) v(5100)=0d0 v(5101)=0d0 v(5102)=v(283) v(5103)=v(284) v(5104)=v(285) v(5105)=0d0 v(5106)=0d0 v(5107)=0d0 v(5108)=0d0 v(5109)=0d0 v(5110)=0d0 v(281)=v(263)*v(326) v(280)=v(260)*v(327) v(279)=v(257)*v(328) v(5276)=v(135) v(5277)=0d0 v(5278)=0d0 v(5279)=v(139) v(5280)=0d0 v(5281)=0d0 v(5282)=v(144) v(5283)=0d0 v(5284)=0d0 v(5285)=v(148) v(5286)=0d0 v(5287)=0d0 v(5288)=v(153) v(5289)=0d0 v(5290)=0d0 v(5291)=v(157) v(5292)=0d0 v(5293)=0d0 v(5294)=v(162) v(5295)=0d0 v(5296)=0d0 v(5297)=v(166) v(5298)=0d0 v(5299)=0d0 v(5300)=v(279) v(5301)=v(280) v(5302)=v(281) v(5303)=0d0 v(5304)=0d0 v(5305)=0d0 v(5306)=0d0 v(5307)=0d0 v(5308)=0d0 v(5210)=0d0 v(5211)=0d0 v(5212)=v(135) v(5213)=0d0 v(5214)=0d0 v(5215)=v(139) v(5216)=0d0 v(5217)=0d0 v(5218)=v(144) v(5219)=0d0 v(5220)=0d0 v(5221)=v(148) v(5222)=0d0 v(5223)=0d0 v(5224)=v(153) v(5225)=0d0 v(5226)=0d0 v(5227)=v(157) v(5228)=0d0 v(5229)=0d0 v(5230)=v(162) v(5231)=0d0 v(5232)=0d0 v(5233)=v(166) v(5234)=0d0 v(5235)=0d0 v(5236)=0d0 v(5237)=0d0 v(5238)=0d0 v(5239)=0d0 v(5240)=v(279) v(5241)=v(280) v(5242)=v(281) v(5144)=0d0 v(5145)=v(135) v(5146)=0d0 v(5147)=0d0 v(5148)=v(139) v(5149)=0d0 v(5150)=0d0 v(5151)=v(144) v(5152)=0d0 v(5153)=0d0 v(5154)=v(148) v(5155)=0d0 v(5156)=0d0 v(5157)=v(153) v(5158)=0d0 v(5159)=0d0 v(5160)=v(157) v(5161)=0d0 v(5162)=0d0 v(5163)=v(162) v(5164)=0d0 v(5165)=0d0 v(5166)=v(166) v(5167)=0d0 v(5168)=0d0 v(5169)=0d0 v(5170)=0d0 v(5171)=v(279) v(5172)=v(280) v(5173)=v(281) v(5174)=0d0 v(5175)=0d0 v(5176)=0d0 v(277)=v(264)*v(326) v(276)=v(261)*v(327) v(275)=v(258)*v(328) v(5342)=0d0 v(5343)=0d0 v(5344)=v(137) v(5345)=0d0 v(5346)=0d0 v(5347)=v(143) v(5348)=0d0 v(5349)=0d0 v(5350)=v(146) v(5351)=0d0 v(5352)=0d0 v(5353)=v(152) v(5354)=0d0 v(5355)=0d0 v(5356)=v(155) v(5357)=0d0 v(5358)=0d0 v(5359)=v(161) v(5360)=0d0 v(5361)=0d0 v(5362)=v(164) v(5363)=0d0 v(5364)=0d0 v(5365)=v(170) v(5366)=0d0 v(5367)=0d0 v(5368)=0d0 v(5369)=0d0 v(5370)=0d0 v(5371)=0d0 v(5372)=v(275) v(5373)=v(276) v(5374)=v(277) v(5177)=0d0 v(5178)=v(137) v(5179)=0d0 v(5180)=0d0 v(5181)=v(143) v(5182)=0d0 v(5183)=0d0 v(5184)=v(146) v(5185)=0d0 v(5186)=0d0 v(5187)=v(152) v(5188)=0d0 v(5189)=0d0 v(5190)=v(155) v(5191)=0d0 v(5192)=0d0 v(5193)=v(161) v(5194)=0d0 v(5195)=0d0 v(5196)=v(164) v(5197)=0d0 v(5198)=0d0 v(5199)=v(170) v(5200)=0d0 v(5201)=0d0 v(5202)=0d0 v(5203)=0d0 v(5204)=v(275) v(5205)=v(276) v(5206)=v(277) v(5207)=0d0 v(5208)=0d0 v(5209)=0d0 v(5111)=v(137) v(5112)=0d0 v(5113)=0d0 v(5114)=v(143) v(5115)=0d0 v(5116)=0d0 v(5117)=v(146) v(5118)=0d0 v(5119)=0d0 v(5120)=v(152) v(5121)=0d0 v(5122)=0d0 v(5123)=v(155) v(5124)=0d0 v(5125)=0d0 v(5126)=v(161) v(5127)=0d0 v(5128)=0d0 v(5129)=v(164) v(5130)=0d0 v(5131)=0d0 v(5132)=v(170) v(5133)=0d0 v(5134)=0d0 v(5135)=v(275) v(5136)=v(276) v(5137)=v(277) v(5138)=0d0 v(5139)=0d0 v(5140)=0d0 v(5141)=0d0 v(5142)=0d0 v(5143)=0d0 v(227)=ht(9)*v(326) v(226)=ht(8)*v(327) v(225)=ht(7)*v(328) v(217)=ht(6)*v(326) v(215)=ht(5)*v(327) v(213)=ht(4)*v(328) v(206)=ht(3)*v(326) v(205)=ht(2)*v(327) v(204)=ht(1)*v(328) v(207)=ul(1,1)*v(136)+ul(1,2)*v(141)+ul(1,3)*v(145)+ul(1,4)*v & (150)+ul(1,5)*v(154)+ul(1,6)*v(159)+ul(1,7)*v(163)+ul(1,8)*v & (168)+(v(204)*v(214)+v(205)*v(216)+v(206)*v(218))*v(322) v(208)=ul(1,1)*v(137)+ul(1,2)*v(143)+ul(1,3)*v(146)+ul(1,4)*v & (152)+ul(1,5)*v(155)+ul(1,6)*v(161)+ul(1,7)*v(164)+ul(1,8)*v & (170)+(v(204)*v(220)+v(205)*v(221)+v(206)*v(222))*v(322) v(212)=ul(2,1)*v(135)+ul(2,2)*v(139)+ul(2,3)*v(144)+ul(2,4)*v & (148)+ul(2,5)*v(153)+ul(2,6)*v(157)+ul(2,7)*v(162)+ul(2,8)*v & (166)+(v(209)*v(213)+v(210)*v(215)+v(211)*v(217))*v(322) v(223)=ul(2,1)*v(137)+ul(2,2)*v(143)+ul(2,3)*v(146)+ul(2,4)*v & (152)+ul(2,5)*v(155)+ul(2,6)*v(161)+ul(2,7)*v(164)+ul(2,8)*v & (170)+(v(213)*v(220)+v(215)*v(221)+v(217)*v(222))*v(322) v(224)=ul(3,1)*v(135)+ul(3,2)*v(139)+ul(3,3)*v(144)+ul(3,4)*v & (148)+ul(3,5)*v(153)+ul(3,6)*v(157)+ul(3,7)*v(162)+ul(3,8)*v & (166)+(v(209)*v(225)+v(210)*v(226)+v(211)*v(227))*v(322) v(228)=ul(3,1)*v(136)+ul(3,2)*v(141)+ul(3,3)*v(145)+ul(3,4)*v & (150)+ul(3,5)*v(154)+ul(3,6)*v(159)+ul(3,7)*v(163)+ul(3,8)*v & (168)+(v(214)*v(225)+v(216)*v(226)+v(218)*v(227))*v(322) v(230)=1d0+ul(1,1)*v(135)+ul(1,2)*v(139)+ul(1,3)*v(144)+ul(1,4 & )*v(148)+ul(1,5)*v(153)+ul(1,6)*v(157)+ul(1,7)*v(162)+ul(1,8 & )*v(166)+(v(204)*v(209)+v(205)*v(210)+v(206)*v(211))*v(322) v(295)=v(208)*v(212)-v(223)*v(230) v(293)=v(207)*v(224)-v(228)*v(230) v(231)=1d0+ul(2,1)*v(136)+ul(2,2)*v(141)+ul(2,3)*v(145)+ul(2,4 & )*v(150)+ul(2,5)*v(154)+ul(2,6)*v(159)+ul(2,7)*v(163)+ul(2,8 & )*v(168)+(v(213)*v(214)+v(215)*v(216)+v(217)*v(218))*v(322) v(298)=-(v(207)*v(212))+v(230)*v(231) v(294)=v(207)*v(223)-v(208)*v(231) v(291)=v(212)*v(228)-v(224)*v(231) v(232)=1d0+ul(3,1)*v(137)+ul(3,2)*v(143)+ul(3,3)*v(146)+ul(3,4 & )*v(152)+ul(3,5)*v(155)+ul(3,6)*v(161)+ul(3,7)*v(164)+ul(3,8 & )*v(170)+(v(220)*v(225)+v(221)*v(226)+v(222)*v(227))*v(322) v(297)=-(v(208)*v(224))+v(230)*v(232) v(296)=-(v(223)*v(228))+v(231)*v(232) v(292)=v(208)*v(228)-v(207)*v(232) v(290)=v(223)*v(224)-v(212)*v(232) v(239)=v(224)*v(294)+v(228)*v(295)+v(232)*v(298) v(247)=((-1d0)+v(239))*v(240)-v(242)/v(239) v(248)=v(232)*v(242)+v(247)*v(298) v(249)=v(231)*v(242)+v(247)*v(297) v(250)=v(230)*v(242)+v(247)*v(296) v(251)=v(228)*v(242)+v(247)*v(295) v(252)=v(224)*v(242)+v(247)*v(294) v(253)=v(223)*v(242)+v(247)*v(293) v(254)=v(212)*v(242)+v(247)*v(292) v(255)=v(208)*v(242)+v(247)*v(291) v(256)=v(207)*v(242)+v(247)*v(290) v(5041)=v(135)*v(250)+v(137)*v(255)+v(136)*v(256) v(5042)=v(136)*v(249)+v(137)*v(253)+v(135)*v(254) v(5043)=v(137)*v(248)+v(136)*v(251)+v(135)*v(252) v(5044)=v(139)*v(250)+v(143)*v(255)+v(141)*v(256) v(5045)=v(141)*v(249)+v(143)*v(253)+v(139)*v(254) v(5046)=v(143)*v(248)+v(141)*v(251)+v(139)*v(252) v(5047)=v(144)*v(250)+v(146)*v(255)+v(145)*v(256) v(5048)=v(145)*v(249)+v(146)*v(253)+v(144)*v(254) v(5049)=v(146)*v(248)+v(145)*v(251)+v(144)*v(252) v(5050)=v(148)*v(250)+v(152)*v(255)+v(150)*v(256) v(5051)=v(150)*v(249)+v(152)*v(253)+v(148)*v(254) v(5052)=v(152)*v(248)+v(150)*v(251)+v(148)*v(252) v(5053)=v(153)*v(250)+v(155)*v(255)+v(154)*v(256) v(5054)=v(154)*v(249)+v(155)*v(253)+v(153)*v(254) v(5055)=v(155)*v(248)+v(154)*v(251)+v(153)*v(252) v(5056)=v(157)*v(250)+v(161)*v(255)+v(159)*v(256) v(5057)=v(159)*v(249)+v(161)*v(253)+v(157)*v(254) v(5058)=v(161)*v(248)+v(159)*v(251)+v(157)*v(252) v(5059)=v(162)*v(250)+v(164)*v(255)+v(163)*v(256) v(5060)=v(163)*v(249)+v(164)*v(253)+v(162)*v(254) v(5061)=v(164)*v(248)+v(163)*v(251)+v(162)*v(252) v(5062)=v(166)*v(250)+v(170)*v(255)+v(168)*v(256) v(5063)=v(168)*v(249)+v(170)*v(253)+v(166)*v(254) v(5064)=v(170)*v(248)+v(168)*v(251)+v(166)*v(252) v(5065)=(v(250)*v(257)+v(255)*v(258)+v(256)*v(259))*v(328) v(5066)=(v(250)*v(260)+v(255)*v(261)+v(256)*v(262))*v(327) v(5067)=(v(250)*v(263)+v(255)*v(264)+v(256)*v(265))*v(326) v(5068)=(v(254)*v(257)+v(253)*v(258)+v(249)*v(259))*v(328) v(5069)=(v(254)*v(260)+v(253)*v(261)+v(249)*v(262))*v(327) v(5070)=(v(254)*v(263)+v(253)*v(264)+v(249)*v(265))*v(326) v(5071)=(v(252)*v(257)+v(248)*v(258)+v(251)*v(259))*v(328) v(5072)=(v(252)*v(260)+v(248)*v(261)+v(251)*v(262))*v(327) v(5073)=(v(252)*v(263)+v(248)*v(264)+v(251)*v(265))*v(326) DO i244=1,33 v(272)=v(330)*v(5077+i244) v(273)=v(330)*v(5110+i244) v(274)=v(330)*v(5143+i244) v(278)=v(330)*v(5176+i244) v(282)=v(330)*v(5209+i244) v(286)=v(330)*v(5242+i244) v(287)=v(330)*v(5275+i244) v(288)=v(330)*v(5308+i244) v(289)=v(330)*v(5341+i244) v(299)=(v(240)+v(242)/v(239)**2)*(v(272)*v(290)+v(273)*v(291) & +v(274)*v(292)+v(278)*v(293)+v(282)*v(294)+v(286)*v(295)+v & (287)*v(296)+v(288)*v(297)+v(289)*v(298)) v(300)=v(247)*(-(v(212)*v(272))-v(207)*v(274)+v(231)*v(287)+v & (230)*v(288))+v(242)*v(289)+v(298)*v(299) v(301)=v(242)*v(288)+v(247)*(-(v(224)*v(273))-v(208)*v(282)+v & (232)*v(287)+v(230)*v(289))+v(297)*v(299) v(302)=v(242)*v(287)+v(247)*(-(v(228)*v(278))-v(223)*v(286)+v & (232)*v(288)+v(231)*v(289))+v(296)*v(299) v(303)=v(242)*v(286)+v(247)*(v(212)*v(273)+v(208)*v(274)-v & (230)*v(278)-v(223)*v(287))+v(295)*v(299) v(304)=v(242)*v(282)+v(247)*(v(223)*v(272)-v(231)*v(273)+v & (207)*v(278)-v(208)*v(288))+v(294)*v(299) v(305)=v(242)*v(278)+v(247)*(v(224)*v(272)+v(207)*v(282)-v & (230)*v(286)-v(228)*v(287))+v(293)*v(299) v(306)=v(242)*v(274)+v(247)*(-(v(232)*v(272))+v(228)*v(273)+v & (208)*v(286)-v(207)*v(289))+v(292)*v(299) v(307)=v(242)*v(273)+v(247)*(v(228)*v(274)-v(231)*v(282)+v & (212)*v(286)-v(224)*v(288))+v(291)*v(299) v(308)=v(242)*v(272)+v(247)*(-(v(232)*v(274))+v(224)*v(278)+v & (223)*v(282)-v(212)*v(289))+v(290)*v(299) v(5375)=v(135)*v(302)+v(137)*v(307)+v(136)*v(308) v(5376)=v(136)*v(301)+v(137)*v(305)+v(135)*v(306) v(5377)=v(137)*v(300)+v(136)*v(303)+v(135)*v(304) v(5378)=v(139)*v(302)+v(143)*v(307)+v(141)*v(308) v(5379)=v(141)*v(301)+v(143)*v(305)+v(139)*v(306) v(5380)=v(143)*v(300)+v(141)*v(303)+v(139)*v(304) v(5381)=v(144)*v(302)+v(146)*v(307)+v(145)*v(308) v(5382)=v(145)*v(301)+v(146)*v(305)+v(144)*v(306) v(5383)=v(146)*v(300)+v(145)*v(303)+v(144)*v(304) v(5384)=v(148)*v(302)+v(152)*v(307)+v(150)*v(308) v(5385)=v(150)*v(301)+v(152)*v(305)+v(148)*v(306) v(5386)=v(152)*v(300)+v(150)*v(303)+v(148)*v(304) v(5387)=v(153)*v(302)+v(155)*v(307)+v(154)*v(308) v(5388)=v(154)*v(301)+v(155)*v(305)+v(153)*v(306) v(5389)=v(155)*v(300)+v(154)*v(303)+v(153)*v(304) v(5390)=v(157)*v(302)+v(161)*v(307)+v(159)*v(308) v(5391)=v(159)*v(301)+v(161)*v(305)+v(157)*v(306) v(5392)=v(161)*v(300)+v(159)*v(303)+v(157)*v(304) v(5393)=v(162)*v(302)+v(164)*v(307)+v(163)*v(308) v(5394)=v(163)*v(301)+v(164)*v(305)+v(162)*v(306) v(5395)=v(164)*v(300)+v(163)*v(303)+v(162)*v(304) v(5396)=v(166)*v(302)+v(170)*v(307)+v(168)*v(308) v(5397)=v(168)*v(301)+v(170)*v(305)+v(166)*v(306) v(5398)=v(170)*v(300)+v(168)*v(303)+v(166)*v(304) v(5399)=(v(257)*v(302)+v(258)*v(307)+v(259)*v(308))*v(328) v(5400)=(v(260)*v(302)+v(261)*v(307)+v(262)*v(308))*v(327) v(5401)=(v(263)*v(302)+v(264)*v(307)+v(265)*v(308))*v(326) v(5402)=(v(259)*v(301)+v(258)*v(305)+v(257)*v(306))*v(328) v(5403)=(v(262)*v(301)+v(261)*v(305)+v(260)*v(306))*v(327) v(5404)=(v(265)*v(301)+v(264)*v(305)+v(263)*v(306))*v(326) v(5405)=(v(258)*v(300)+v(259)*v(303)+v(257)*v(304))*v(328) v(5406)=(v(261)*v(300)+v(262)*v(303)+v(260)*v(304))*v(327) v(5407)=(v(264)*v(300)+v(265)*v(303)+v(263)*v(304))*v(326) p(i244)=p(i244)-v(330)*v(5040+i244) DO i269=i244,33 s(i244,i269)=s(i244,i269)+v(5374+i269) ENDDO ENDDO ENDDO END