!************************************************************** !* AceGen 1.1 Windows (13 Aug 07) * !* Co. J. Korelc 2007 13 Aug 07 15:10:16* !************************************************************** ! User : USER ! Evaluation time : 11 s Mode : Optimal ! Number of formulae : 163 Method: Automatic ! Subroutine : SKR size :4249 ! Total size of Mathematica code : 4249 subexpressions ! Total size of Fortran code : 13805 bytes ! ! User element declaration ! *User element, type=U1, nodes=8, coordinates=3, properties=2, variables=1 ! 1,2,3 ! include 'SMSUtility.for' SUBROUTINE UEL(p,s,SVARS,ENERGY,NDOFEL,NRHS,NSVARS, 1 PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME, 2 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF, 3 NPREDF,LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP, 4 PERIOD) INCLUDE 'ABA_PARAM.INC' include 'sms.h' DIMENSION p(MLVARX,*),s(NDOFEL,NDOFEL), 1 SVARS(NSVARS),ENERGY(8),PROPS(*),COORDS(MCRD,NNODE), 2 U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2), 3 PARAMS(3),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*), 4 DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS(*), 5 JPROPS(*) logical DEBUG,symmetric,DEBUGELEMENT,cedebug,firststep character*50 SELEM,datades(2),postdes(0) parameter (DEBUG=.false.,DEBUGELEMENT=-1, # SELEM="ExamplesHypersolidABAQUS") integer ngdof,ncdof,hlen integer maxdof,dofglobal(8),nnodes,nalldof,ndim double precision ul(3,8),ul0(3,8) integer i,j,jj,ll,ii,k,kk,i1,i2,i3,icode,iow double precision vv(5301),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,"> ",24g12.5) 1238 format(i3,"> ",24g12.5) 1239 format(i3,"> ",3g12.5) if(TIME(2).eq. 0.0 .and. uecheck.ne.JTYPE) then do i=1,IData_Last idata(i)=0 enddo do i=1,RData_Last rdata(i)=0.d0 enddo endif !.... integration code (if changed, the number of state variables can change also) icode=7 ncdof=0 ngdof=24 nalldof=24 symmetric=.true. maxdof=3 nnodes=8 ndim=3 !.... write input data to .dat file iow=6 idata(ID_OutputFile)=iow idata(ID_DebugElement)=DEBUGELEMENT cedebug=DEBUG.and.(DEBUGELEMENT.eq.-1.or.JELEM.eq.DEBUGELEMENT) rdata(RD_SubIterationTolerance)=1d-9 rdata(RD_Multiplier)=TIME(1) rdata(RD_MultiplierIncrement)=DTIME rdata(RD_Time)=TIME(1) rdata(RD_TimeIncrement)=DTIME idata(ID_CurrentElement)=JELEM idata(ID_MaxMessages)=100 if(KINC.ne.ckinc.or.KSTEP.ne.ckstep) then idata(ID_NoMessages)=0 idata(ID_TotalIteration)=idata(ID_TotalIteration)+1 ckinc=KINC ckstep=KSTEP endif if(idata(ID_LastIntCode).ne.icode) then call SMSIntPoints(icode,ngpo,gp) idata(ID_LastIntCode)=icode endif hlen=0 if(TIME(2).eq. 0.0 .and. uecheck.ne.JTYPE) then write(iow,*)"AceGen - initialization" write(iow,*)"Element:",SELEM if(NDOFEL.ne.ngdof) then write(iow,*)"USER ELEMENT ERROR!!!:", # "wrong number of d.o.f per element" write(iow,*)"Number of d.o.f defined in input file: ",NDOFEL write(iow,*)"Number of d.o.f required by the element: ",ngdof endif if(NSVARS.gt.1.and.NSVARS.ne.hlen) then write(iow,*)"USER ELEMENT ERROR!!!:", # "wrong number of state variables per element" write(iow,*)"No. of state variables", # " defined in input file: ",NSVARS write(iow,*)"No. of state variables", # " required by the element: ",hlen endif !....Description of the input data datades(1)="Elastic modulus" datades(2)="Poisson ratio" write(iow,"(10x,f15.5,A3,A50)") # (PROPS(i)," = ",datades(i),i=1,2) if(NPROPS.ne.2) then write(iow,*)"ERROR in AceGen element", # " - number uf user properties" write(iow,*)"Supplied: ",NPROPS, " Required:",2 endif uecheck=JTYPE endif !... write debug data to .msg file iow=7 idata(ID_OutputFile)=iow firststep=.true. k=0 do i=1,nnodes do j=1,dofglobal(i) k=k+1 ul(j,i)=U(k) ul0(j,i)=ul(j,i)-DU(k,1) firststep=firststep.and. abs(DU(k,1)) .lt. 1d-10 enddo enddo if(firststep) then idata(ID_Iteration)=1 else idata(ID_Iteration)=2 endif if(cedebug) then write(iow,*)"#################################################" write(iow,*)SELEM write(iow,"(A,i3,A,f15.5,A,f15.5)")"ELEMENT=",JELEM, # " TIME=",TIME(2)," INCREMENT=",DTIME write(iow,"(A,i3,A,i3,A,i3)")"ITERATION=",idata(ID_Iteration), # " TOTAL_ITERATION=",idata(ID_TotalIteration)," TASK",LFLAGS(3) write(iow,"(A,5i4)")"LFLAGS",(LFLAGS(i),i=1,5) write(iow,*)"NODE X" write(iow,1239)(i,(COORDS(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),j=1,maxdof),i=1,nnodes) write(iow,*)"ELEMENT hp" write(iow,1236)(SVARS(i),i=1,NSVARS) endif SELECT CASE(LFLAGS(3)) CASE DEFAULT write(iow,*)"User switch ",LFLAGS(3)," not implemented" return CASE (4) !... initial mass matrix do i=1,ngdof do j=1,ngdof s(i,j)=0 enddo s(i,i)=1. enddo if(cedebug) then write(iow,*)"INITIAL MASS MATRIX" write(iow,*)"MatrixLocal=" write(iow,1237)( i,(s(i,j),j=1,ngdof),i=1,ngdof) endif return CASE (1) !... tangent and residuum call SMSZero(s,ngdof*ngdof) call SMSZero(p,ngdof) call SKR(vv,PROPS,ul,ul0,COORDS,s,p,SVARS,SVARS) if(symmetric) then do i=2,nalldof do j=1,i-1 s(i,j)=s(j,i) enddo enddo endif if(cedebug) then write(iow,*)"TANGENT AND RESIDUAL SUBROUTINE" write(iow,*)"ELEMENT ht" write(iow,1236)(SVARS(i),i=1,NSVARS) 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,1),i=1,ngdof) endif return CASE (2) !....only tangent if(cedebug) then write(iow,*)"TANGENT for ",SELEM endif idata(ID_SkipResidual)=1 call SMSZero(s,ngdof*ngdof) call SMSZero(p,ngdof) call SKR(vv,PROPS,ul,ul0,COORDS,s,p,SVARS,SVARS) if(symmetric) then do i=2,nalldof do j=1,i-1 s(i,j)=s(j,i) enddo enddo endif idata(ID_SkipResidual)=0 CASE (5) !....only residual if(cedebug) then write(iow,*)"RESIDUAL for ",SELEM endif idata(ID_SkipTangent)=1 call SMSZero(s,ngdof*ngdof) call SMSZero(p,ngdof) call SKR(vv,PROPS,ul,ul0,COORDS,s,p,SVARS,SVARS) if(symmetric) then do i=2,nalldof do j=1,i-1 s(i,j)=s(j,i) enddo enddo endif idata(ID_SkipTangent)=0 END SELECT End !******************* S U B R O U T I N E ********************** SUBROUTINE SKR(v,d,ul,ul0,xl,s,p,ht,hp) IMPLICIT NONE include 'sms.h' INTEGER i1,i206,i221 DOUBLE PRECISION v(5301),d(2),ul(3,8),ul0(3,8),xl(3,8),s(24,24) &,p(24),ht(*),hp(*) v(255)=2d0*d(2) v(203)=1d0/(1d0+d(2)) v(204)=(d(1)*v(203))/2d0 v(202)=(v(204)*v(255))/(1d0-v(255)) 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(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(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(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(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(258)=-(v(90)*v(92))+v(89)*v(93) 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(257)=v(91)*v(92)-v(89)*v(94) v(256)=-(v(91)*v(93))+v(90)*v(94) 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(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(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(98)=v(256)*v(95)+v(257)*v(96)+v(258)*v(97) v(260)=v(55)*v(98) 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(256)/v(98) v(127)=v(108)*v(79) v(120)=v(108)*v(78) v(109)=v(257)/v(98) v(130)=v(109)*v(79) v(122)=v(109)*v(78) v(110)=v(258)/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(5204)=v(135) v(5205)=0d0 v(5206)=0d0 v(5207)=v(139) v(5208)=0d0 v(5209)=0d0 v(5210)=v(144) v(5211)=0d0 v(5212)=0d0 v(5213)=v(148) v(5214)=0d0 v(5215)=0d0 v(5216)=v(153) v(5217)=0d0 v(5218)=0d0 v(5219)=v(157) v(5220)=0d0 v(5221)=0d0 v(5222)=v(162) v(5223)=0d0 v(5224)=0d0 v(5225)=v(166) v(5226)=0d0 v(5227)=0d0 v(5084)=0d0 v(5085)=0d0 v(5086)=v(135) v(5087)=0d0 v(5088)=0d0 v(5089)=v(139) v(5090)=0d0 v(5091)=0d0 v(5092)=v(144) v(5093)=0d0 v(5094)=0d0 v(5095)=v(148) v(5096)=0d0 v(5097)=0d0 v(5098)=v(153) v(5099)=0d0 v(5100)=0d0 v(5101)=v(157) v(5102)=0d0 v(5103)=0d0 v(5104)=v(162) v(5105)=0d0 v(5106)=0d0 v(5107)=v(166) v(5060)=0d0 v(5061)=v(135) v(5062)=0d0 v(5063)=0d0 v(5064)=v(139) v(5065)=0d0 v(5066)=0d0 v(5067)=v(144) v(5068)=0d0 v(5069)=0d0 v(5070)=v(148) v(5071)=0d0 v(5072)=0d0 v(5073)=v(153) v(5074)=0d0 v(5075)=0d0 v(5076)=v(157) v(5077)=0d0 v(5078)=0d0 v(5079)=v(162) v(5080)=0d0 v(5081)=0d0 v(5082)=v(166) v(5083)=0d0 v(168)=-v(167)+v(112)*v(63) v(5228)=0d0 v(5229)=v(136) v(5230)=0d0 v(5231)=0d0 v(5232)=v(141) v(5233)=0d0 v(5234)=0d0 v(5235)=v(145) v(5236)=0d0 v(5237)=0d0 v(5238)=v(150) v(5239)=0d0 v(5240)=0d0 v(5241)=v(154) v(5242)=0d0 v(5243)=0d0 v(5244)=v(159) v(5245)=0d0 v(5246)=0d0 v(5247)=v(163) v(5248)=0d0 v(5249)=0d0 v(5250)=v(168) v(5251)=0d0 v(5132)=0d0 v(5133)=0d0 v(5134)=v(136) v(5135)=0d0 v(5136)=0d0 v(5137)=v(141) v(5138)=0d0 v(5139)=0d0 v(5140)=v(145) v(5141)=0d0 v(5142)=0d0 v(5143)=v(150) v(5144)=0d0 v(5145)=0d0 v(5146)=v(154) v(5147)=0d0 v(5148)=0d0 v(5149)=v(159) v(5150)=0d0 v(5151)=0d0 v(5152)=v(163) v(5153)=0d0 v(5154)=0d0 v(5155)=v(168) v(5108)=v(136) v(5109)=0d0 v(5110)=0d0 v(5111)=v(141) v(5112)=0d0 v(5113)=0d0 v(5114)=v(145) v(5115)=0d0 v(5116)=0d0 v(5117)=v(150) v(5118)=0d0 v(5119)=0d0 v(5120)=v(154) v(5121)=0d0 v(5122)=0d0 v(5123)=v(159) v(5124)=0d0 v(5125)=0d0 v(5126)=v(163) v(5127)=0d0 v(5128)=0d0 v(5129)=v(168) v(5130)=0d0 v(5131)=0d0 v(170)=-v(169)+v(113)*v(63) v(5252)=0d0 v(5253)=0d0 v(5254)=v(137) v(5255)=0d0 v(5256)=0d0 v(5257)=v(143) v(5258)=0d0 v(5259)=0d0 v(5260)=v(146) v(5261)=0d0 v(5262)=0d0 v(5263)=v(152) v(5264)=0d0 v(5265)=0d0 v(5266)=v(155) v(5267)=0d0 v(5268)=0d0 v(5269)=v(161) v(5270)=0d0 v(5271)=0d0 v(5272)=v(164) v(5273)=0d0 v(5274)=0d0 v(5275)=v(170) v(5180)=0d0 v(5181)=v(137) v(5182)=0d0 v(5183)=0d0 v(5184)=v(143) v(5185)=0d0 v(5186)=0d0 v(5187)=v(146) v(5188)=0d0 v(5189)=0d0 v(5190)=v(152) v(5191)=0d0 v(5192)=0d0 v(5193)=v(155) v(5194)=0d0 v(5195)=0d0 v(5196)=v(161) v(5197)=0d0 v(5198)=0d0 v(5199)=v(164) v(5200)=0d0 v(5201)=0d0 v(5202)=v(170) v(5203)=0d0 v(5156)=v(137) v(5157)=0d0 v(5158)=0d0 v(5159)=v(143) v(5160)=0d0 v(5161)=0d0 v(5162)=v(146) v(5163)=0d0 v(5164)=0d0 v(5165)=v(152) v(5166)=0d0 v(5167)=0d0 v(5168)=v(155) v(5169)=0d0 v(5170)=0d0 v(5171)=v(161) v(5172)=0d0 v(5173)=0d0 v(5174)=v(164) v(5175)=0d0 v(5176)=0d0 v(5177)=v(170) v(5178)=0d0 v(5179)=0d0 v(172)=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(173)=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(174)=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(176)=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(177)=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(178)=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(192)=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(238)=v(173)*v(174)-v(176)*v(192) v(236)=v(172)*v(177)-v(178)*v(192) v(193)=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(241)=-(v(172)*v(174))+v(192)*v(193) v(237)=v(172)*v(176)-v(173)*v(193) v(234)=v(174)*v(178)-v(177)*v(193) v(194)=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(240)=-(v(173)*v(177))+v(192)*v(194) v(239)=-(v(176)*v(178))+v(193)*v(194) v(235)=v(173)*v(178)-v(172)*v(194) v(233)=v(176)*v(177)-v(174)*v(194) v(201)=v(177)*v(237)+v(178)*v(238)+v(194)*v(241) v(209)=((-1d0)+v(201))*v(202)-v(204)/v(201) DO i206=1,24 v(224)=v(260)*v(5059+i206) v(225)=v(260)*v(5083+i206) v(226)=v(260)*v(5107+i206) v(227)=v(260)*v(5131+i206) v(228)=v(260)*v(5155+i206) v(229)=v(260)*v(5179+i206) v(230)=v(260)*v(5203+i206) v(231)=v(260)*v(5227+i206) v(232)=v(260)*v(5251+i206) v(242)=(v(202)+v(204)/v(201)**2)*(v(224)*v(233)+v(225)*v(234) & +v(226)*v(235)+v(227)*v(236)+v(228)*v(237)+v(229)*v(238)+v & (230)*v(239)+v(231)*v(240)+v(232)*v(241)) v(243)=v(209)*(-(v(174)*v(224))-v(172)*v(226)+v(193)*v(230)+v & (192)*v(231))+v(204)*v(232)+v(241)*v(242) v(244)=v(204)*v(231)+v(209)*(-(v(177)*v(225))-v(173)*v(228)+v & (194)*v(230)+v(192)*v(232))+v(240)*v(242) v(245)=v(204)*v(230)+v(209)*(-(v(178)*v(227))-v(176)*v(229)+v & (194)*v(231)+v(193)*v(232))+v(239)*v(242) v(246)=v(204)*v(229)+v(209)*(v(174)*v(225)+v(173)*v(226)-v & (192)*v(227)-v(176)*v(230))+v(238)*v(242) v(247)=v(204)*v(228)+v(209)*(v(176)*v(224)-v(193)*v(225)+v & (172)*v(227)-v(173)*v(231))+v(237)*v(242) v(248)=v(204)*v(227)+v(209)*(v(177)*v(224)+v(172)*v(228)-v & (192)*v(229)-v(178)*v(230))+v(236)*v(242) v(249)=v(204)*v(226)+v(209)*(-(v(194)*v(224))+v(178)*v(225)+v & (173)*v(229)-v(172)*v(232))+v(235)*v(242) v(250)=v(204)*v(225)+v(209)*(v(178)*v(226)-v(193)*v(228)+v & (174)*v(229)-v(177)*v(231))+v(234)*v(242) v(251)=v(204)*v(224)+v(209)*(-(v(194)*v(226))+v(177)*v(227)+v & (176)*v(228)-v(174)*v(232))+v(233)*v(242) v(5276)=v(135)*v(245)+v(137)*v(247)+v(136)*v(249) v(5277)=v(136)*v(244)+v(137)*v(246)+v(135)*v(251) v(5278)=v(137)*v(243)+v(136)*v(248)+v(135)*v(250) v(5279)=v(139)*v(245)+v(143)*v(247)+v(141)*v(249) v(5280)=v(141)*v(244)+v(143)*v(246)+v(139)*v(251) v(5281)=v(143)*v(243)+v(141)*v(248)+v(139)*v(250) v(5282)=v(144)*v(245)+v(146)*v(247)+v(145)*v(249) v(5283)=v(145)*v(244)+v(146)*v(246)+v(144)*v(251) v(5284)=v(146)*v(243)+v(145)*v(248)+v(144)*v(250) v(5285)=v(148)*v(245)+v(152)*v(247)+v(150)*v(249) v(5286)=v(150)*v(244)+v(152)*v(246)+v(148)*v(251) v(5287)=v(152)*v(243)+v(150)*v(248)+v(148)*v(250) v(5288)=v(153)*v(245)+v(155)*v(247)+v(154)*v(249) v(5289)=v(154)*v(244)+v(155)*v(246)+v(153)*v(251) v(5290)=v(155)*v(243)+v(154)*v(248)+v(153)*v(250) v(5291)=v(157)*v(245)+v(161)*v(247)+v(159)*v(249) v(5292)=v(159)*v(244)+v(161)*v(246)+v(157)*v(251) v(5293)=v(161)*v(243)+v(159)*v(248)+v(157)*v(250) v(5294)=v(162)*v(245)+v(164)*v(247)+v(163)*v(249) v(5295)=v(163)*v(244)+v(164)*v(246)+v(162)*v(251) v(5296)=v(164)*v(243)+v(163)*v(248)+v(162)*v(250) v(5297)=v(166)*v(245)+v(170)*v(247)+v(168)*v(249) v(5298)=v(168)*v(244)+v(170)*v(246)+v(166)*v(251) v(5299)=v(170)*v(243)+v(168)*v(248)+v(166)*v(250) DO i221=i206,24 s(i206,i221)=s(i206,i221)+v(5275+i221) ENDDO ENDDO ENDDO END