!************************************************************** !* AceGen 6.405 Windows (22 Jun 15) * !* Co. J. Korelc 2013 24 Jun 15 07:30:51 * !************************************************************** ! User : USER ! Notebook : Remote Session ! Evaluation time : 30 s Mode : Optimal ! Number of formulae : 488 Method: Automatic ! Subroutine : SKR size: 11019 ! Total size of Mathematica code : 11019 subexpressions ! Total size of Fortran code : 39103 bytes ! Element: Standard 3D serendipity 20 node element, isoparametric, displacement based (H2S). ! Constitutive model: Strain driven hyperelastic, ! Unified element code = SED3H2SDFHYH2SNeoHooke ! Bibliography = ! Variational formulation = Standard displacement formulation Hexahedron ! Integration = Standard numerical integration ! Memory requirement/element: ! DOF = 60 ! History = 0 ! Code structure = {"SE", "D3", "H2S", "DF", "HY", "H2S", "NeoHooke"} ! Problem type = Stationar Nonlinear ! Group data: ! 1 "E -elastic modulus" ! 2 "\[Nu] -Poisson ratio" ! Parameter sensitivity: None ! Material point postprocessing quantities: ! Nodal point postprocessing quantities: ! ! User element declaration ! *User element, type=U1, nodes=20, 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="SED3H2SDFHYH2SNeoHooke") integer ngdof,ncdof,hlen integer maxdof,dofglobal(20),nnodes,nalldof,ndim double precision ul(3,20),ul0(3,20) integer i,j,jj,ll,ii,k,kk,i1,i2,i3,icode,iow double precision vv(17608),gpost(64,0),npost(20,0) integer ipordl(54) data (ipordl(i),i=1,54)/1,12,4,11,3,10,2,9,1,1,9,2,14,6,17, & 5,13,1,2,10,3,15,7,18,6,14,2,4,16,8, & 19,7,15,3,11,4,1,13,5,20,8,16,4,12,1, & 5,17,6,18,7,19,8,20,5/ data (dofglobal(i),i=1,20)/3,3,3,3,3,3,3,3,3,3,3,3,3,3,3, & 3,3,3,3,3/ 1235 format(i3,">",3g17.9) 1236 format(" >",5g17.9) 1237 format(i3,"> ",60g12.5) 1238 format(i3,"> ",60g12.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=8 ncdof=0 ngdof=60 nalldof=60 symmetric=.true. maxdof=3 nnodes=20 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)="E -elastic modulus" datades(2)="[Nu] -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 i3,i4,i469,i488 LOGICAL b1,b2,b484,b486 DOUBLE PRECISION v(17608),d(2),ul(3,20),ul0(3,20),xl(3 &,20),s(60,60),p(60),ht(*),hp(*) v(643)=xl(1,1) v(642)=xl(2,1) v(641)=xl(3,1) v(640)=xl(1,2) v(639)=xl(2,2) v(638)=xl(3,2) v(637)=xl(1,3) v(636)=xl(2,3) v(635)=xl(3,3) v(634)=xl(1,4) v(633)=xl(2,4) v(632)=xl(3,4) v(631)=xl(1,5) v(630)=xl(2,5) v(629)=xl(3,5) v(628)=xl(1,6) v(627)=xl(2,6) v(626)=xl(3,6) v(625)=xl(1,7) v(624)=xl(2,7) v(623)=xl(3,7) v(622)=xl(1,8) v(621)=xl(2,8) v(620)=xl(3,8) v(619)=xl(1,9) v(618)=xl(2,9) v(617)=xl(3,9) v(616)=xl(1,10) v(615)=xl(2,10) v(614)=xl(3,10) v(613)=xl(1,11) v(612)=xl(2,11) v(611)=xl(3,11) v(610)=xl(1,12) v(609)=xl(2,12) v(608)=xl(3,12) v(607)=xl(1,13) v(606)=xl(2,13) v(605)=xl(3,13) v(604)=xl(1,14) v(603)=xl(2,14) v(602)=xl(3,14) v(601)=xl(1,15) v(600)=xl(2,15) v(599)=xl(3,15) v(598)=xl(1,16) v(597)=xl(2,16) v(596)=xl(3,16) v(595)=xl(1,17) v(594)=xl(2,17) v(593)=xl(3,17) v(592)=xl(1,18) v(591)=xl(2,18) v(590)=xl(3,18) v(589)=xl(1,19) v(588)=xl(2,19) v(587)=xl(3,19) v(586)=xl(1,20) v(585)=xl(2,20) v(584)=xl(3,20) v(583)=2.*d(2) v(582)=1d0/(1.+d(2)) v(581)=ul(1,1) v(580)=ul(2,1) v(579)=ul(3,1) v(578)=ul(1,2) v(577)=ul(2,2) v(576)=ul(3,2) v(575)=ul(1,3) v(574)=ul(2,3) v(573)=ul(3,3) v(572)=ul(1,4) v(571)=ul(2,4) v(570)=ul(3,4) v(569)=ul(1,5) v(568)=ul(2,5) v(567)=ul(3,5) v(566)=ul(1,6) v(565)=ul(2,6) v(564)=ul(3,6) v(563)=ul(1,7) v(562)=ul(2,7) v(561)=ul(3,7) v(560)=ul(1,8) v(559)=ul(2,8) v(558)=ul(3,8) v(557)=ul(1,9) v(556)=ul(2,9) v(555)=ul(3,9) v(554)=ul(1,10) v(553)=ul(2,10) v(552)=ul(3,10) v(551)=ul(1,11) v(550)=ul(2,11) v(549)=ul(3,11) v(548)=ul(1,12) v(547)=ul(2,12) v(546)=ul(3,12) v(545)=ul(1,13) v(544)=ul(2,13) v(543)=ul(3,13) v(542)=ul(1,14) v(541)=ul(2,14) v(540)=ul(3,14) v(539)=ul(1,15) v(538)=ul(2,15) v(537)=ul(3,15) v(536)=ul(1,16) v(535)=ul(2,16) v(534)=ul(3,16) v(533)=ul(1,17) v(532)=ul(2,17) v(531)=ul(3,17) v(530)=ul(1,18) v(529)=ul(2,18) v(528)=ul(3,18) v(527)=ul(1,19) v(526)=ul(2,19) v(525)=ul(3,19) v(524)=ul(1,20) v(523)=ul(2,20) v(522)=ul(3,20) v(451)=(d(1)*v(582))/2. v(449)=(v(451)*v(583))/(1.-v(583)) b1=int(idata(ID_SkipResidual)).eq.1 b2=int(idata(ID_SkipTangent)).eq.1 i3=int(ngpo) DO i4=1,i3 v(5)=gp(1,i4) v(182)=-2.*v(5) v(44)=1.-(v(5)*v(5)) v(191)=v(44)/4. v(17)=1.-v(5) v(13)=1.+v(5) v(6)=gp(2,i4) v(183)=-2.*v(6) v(46)=1.-(v(6)*v(6)) v(189)=v(46)/4. v(245)=-(v(17)*v(189)) v(241)=-(v(13)*v(189)) v(25)=1.+v(6) v(658)=v(25)/4. v(243)=-(v(191)*v(25)) v(193)=v(25)/8. v(20)=1.-v(6) v(659)=v(20)/4. v(239)=-(v(191)*v(20)) v(194)=v(20)/8. v(7)=gp(3,i4) v(184)=-2.*v(7) v(644)=v(184)/4. v(187)=v(20)*v(644) v(249)=v(13)*v(187) v(247)=v(17)*v(187) v(186)=v(25)*v(644) v(251)=v(17)*v(186) v(250)=v(13)*v(186) v(40)=1.-(v(7)*v(7)) v(185)=v(40)/4. v(248)=-(v(13)*v(185)) v(246)=-(v(17)*v(185)) v(42)=v(185)*v(25) v(38)=v(185)*v(20) v(29)=-2.+v(6)+v(7) v(233)=v(29)-v(5) v(646)=v(17)*v(233) v(228)=v(29)+v(5) v(647)=v(13)*v(228) v(26)=1.+v(7) v(645)=v(26)/4. v(192)=v(26)/8. v(190)=v(182)*v(645) v(254)=v(190)*v(25) v(252)=v(190)*v(20) v(188)=v(183)*v(645) v(255)=v(17)*v(188) v(253)=v(13)*v(188) v(50)=v(46)*v(645) v(48)=v(44)*v(645) v(28)=v(192)*v(25) v(235)=-(v(17)*v(28)) v(237)=-v(235)+v(193)*v(646) v(236)=-v(235)+v(192)*v(646) v(234)=v(235)-v(233)*v(28) v(230)=v(13)*v(28) v(232)=v(230)+v(193)*v(647) v(231)=v(230)+v(192)*v(647) v(229)=v(230)+v(228)*v(28) v(23)=v(183)+v(29) v(223)=v(23)+v(5) v(648)=v(13)*v(223) v(218)=v(23)-v(5) v(649)=v(17)*v(218) v(22)=v(192)*v(20) v(225)=v(13)*v(22) v(227)=v(225)+v(194)*v(648) v(226)=-v(225)-v(192)*v(648) v(224)=v(22)*v(223)+v(225) v(220)=-(v(17)*v(22)) v(222)=-v(220)+v(194)*v(649) v(221)=v(220)-v(192)*v(649) v(219)=-(v(218)*v(22))+v(220) v(18)=v(184)+v(29) v(213)=v(18)-v(5) v(651)=v(17)*v(213) v(208)=v(18)+v(5) v(652)=v(13)*v(208) v(14)=1.-v(7) v(650)=v(14)/4. v(197)=v(14)/8. v(196)=v(182)*v(650) v(242)=v(196)*v(25) v(238)=v(196)*v(20) v(195)=v(183)*v(650) v(244)=v(17)*v(195) v(240)=v(13)*v(195) v(35)=v(46)*v(650) v(33)=v(44)*v(650) v(16)=v(197)*v(25) v(215)=-(v(16)*v(17)) v(217)=v(215)-v(193)*v(651) v(216)=-v(215)+v(197)*v(651) v(214)=-(v(16)*v(213))+v(215) v(210)=v(13)*v(16) v(212)=-v(210)-v(193)*v(652) v(211)=v(210)+v(197)*v(652) v(209)=v(16)*v(208)+v(210) v(11)=-2.-v(6)-v(7) v(203)=v(11)+v(5) v(653)=-(v(13)*v(203)) v(198)=v(11)-v(5) v(654)=-(v(17)*v(198)) v(10)=v(197)*v(20) v(205)=v(10)*v(13) v(207)=-v(205)+v(194)*v(653) v(206)=-v(205)+v(197)*v(653) v(204)=v(10)*v(203)+v(205) v(200)=-(v(10)*v(17)) v(202)=v(200)+v(194)*v(654) v(201)=v(200)+v(197)*v(654) v(199)=-(v(10)*v(198))+v(200) v(8)=gp(4,i4) v(256)=v(254)*v(589)+v(50)*(-v(586)+v(592))+v(252)*v(595)+v(42 & )*(-v(598)+v(601))+v(38)*(v(604)-v(607))+v(242)*v(613)+v(35)*( & -v(610)+v(616))+v(238)*v(619)+v(234)*v(622)+v(229)*v(625)+v & (224)*v(628)+v(219)*v(631)+v(214)*v(634)+v(209)*v(637)+v(204 & )*v(640)+v(199)*v(643) v(257)=v(254)*v(588)+v(50)*(-v(585)+v(591))+v(252)*v(594)+v(42 & )*(-v(597)+v(600))+v(38)*(v(603)-v(606))+v(242)*v(612)+v(35)*( & -v(609)+v(615))+v(238)*v(618)+v(234)*v(621)+v(229)*v(624)+v & (224)*v(627)+v(219)*v(630)+v(214)*v(633)+v(209)*v(636)+v(204 & )*v(639)+v(199)*v(642) v(258)=v(254)*v(587)+v(50)*(-v(584)+v(590))+v(252)*v(593)+v(42 & )*(-v(596)+v(599))+v(38)*(v(602)-v(605))+v(242)*v(611)+v(35)*( & -v(608)+v(614))+v(238)*v(617)+v(234)*v(620)+v(229)*v(623)+v & (224)*v(626)+v(219)*v(629)+v(214)*v(632)+v(209)*v(635)+v(204 & )*v(638)+v(199)*v(641) v(259)=v(255)*v(586)+v(253)*v(592)+v(48)*(v(589)-v(595))+v(248 & )*(-v(601)+v(604))+v(246)*(-v(598)+v(607))+v(244)*v(610)+v(240 & )*v(616)+v(33)*(v(613)-v(619))+v(236)*v(622)+v(231)*v(625)+v & (226)*v(628)+v(221)*v(631)+v(216)*v(634)+v(211)*v(637)+v(206 & )*v(640)+v(201)*v(643) v(260)=v(255)*v(585)+v(253)*v(591)+v(48)*(v(588)-v(594))+v(248 & )*(-v(600)+v(603))+v(246)*(-v(597)+v(606))+v(244)*v(609)+v(240 & )*v(615)+v(33)*(v(612)-v(618))+v(236)*v(621)+v(231)*v(624)+v & (226)*v(627)+v(221)*v(630)+v(216)*v(633)+v(211)*v(636)+v(206 & )*v(639)+v(201)*v(642) v(657)=-(v(257)*v(259))+v(256)*v(260) v(261)=v(255)*v(584)+v(253)*v(590)+v(48)*(v(587)-v(593))+v(248 & )*(-v(599)+v(602))+v(246)*(-v(596)+v(605))+v(244)*v(608)+v(240 & )*v(614)+v(33)*(v(611)-v(617))+v(236)*v(620)+v(231)*v(623)+v & (226)*v(626)+v(221)*v(629)+v(216)*v(632)+v(211)*v(635)+v(206 & )*v(638)+v(201)*v(641) v(656)=v(258)*v(259)-v(256)*v(261) v(655)=-(v(258)*v(260))+v(257)*v(261) v(262)=v(251)*v(598)+v(250)*v(601)+v(249)*v(604)+v(247)*v(607) & +v(245)*(-v(586)+v(610))+v(243)*(-v(589)+v(613))+v(241)*(-v & (592)+v(616))+v(239)*(-v(595)+v(619))+v(237)*v(622)+v(232)*v & (625)+v(227)*v(628)+v(222)*v(631)+v(217)*v(634)+v(212)*v(637) & +v(207)*v(640)+v(202)*v(643) v(263)=v(251)*v(597)+v(250)*v(600)+v(249)*v(603)+v(247)*v(606) & +v(245)*(-v(585)+v(609))+v(243)*(-v(588)+v(612))+v(241)*(-v & (591)+v(615))+v(239)*(-v(594)+v(618))+v(237)*v(621)+v(232)*v & (624)+v(227)*v(627)+v(222)*v(630)+v(217)*v(633)+v(212)*v(636) & +v(207)*v(639)+v(202)*v(642) v(264)=v(251)*v(596)+v(250)*v(599)+v(249)*v(602)+v(247)*v(605) & +v(245)*(-v(584)+v(608))+v(243)*(-v(587)+v(611))+v(241)*(-v & (590)+v(614))+v(239)*(-v(593)+v(617))+v(237)*v(620)+v(232)*v & (623)+v(227)*v(626)+v(222)*v(629)+v(217)*v(632)+v(212)*v(635) & +v(207)*v(638)+v(202)*v(641) v(265)=v(262)*v(655)+v(263)*v(656)+v(264)*v(657) v(266)=(-(v(261)*v(263))+v(260)*v(264))/v(265) v(663)=v(16)*v(266) v(660)=-(v(10)*v(266)) v(431)=v(266)*v(50) v(413)=v(266)*v(42) v(404)=-(v(266)*v(38)) v(395)=v(266)*v(35) v(267)=(v(261)*v(262)-v(259)*v(264))/v(265) v(664)=v(16)*v(267) v(661)=-(v(10)*v(267)) v(433)=v(267)*v(50) v(415)=v(267)*v(42) v(406)=-(v(267)*v(38)) v(397)=v(267)*v(35) v(268)=(-(v(260)*v(262))+v(259)*v(263))/v(265) v(665)=v(16)*v(268) v(662)=-(v(10)*v(268)) v(435)=v(268)*v(50) v(417)=v(268)*v(42) v(408)=-(v(268)*v(38)) v(399)=v(268)*v(35) v(269)=(v(258)*v(263)-v(257)*v(264))/v(265) v(425)=-(v(269)*v(48)) v(389)=-(v(269)*v(33)) v(338)=v(197)*v(269) v(314)=v(192)*v(269) v(293)=v(185)*v(269) v(270)=(-(v(258)*v(262))+v(256)*v(264))/v(265) v(427)=-(v(270)*v(48)) v(391)=-(v(270)*v(33)) v(341)=v(197)*v(270) v(316)=v(192)*v(270) v(295)=v(185)*v(270) v(271)=(v(257)*v(262)-v(256)*v(263))/v(265) v(429)=-(v(271)*v(48)) v(393)=-(v(271)*v(33)) v(344)=v(197)*v(271) v(318)=v(192)*v(271) v(297)=v(185)*v(271) v(272)=v(655)/v(265) v(339)=v(194)*v(272) v(332)=v(193)*v(272) v(326)=v(191)*v(272) v(320)=v(189)*v(272) v(273)=v(656)/v(265) v(342)=v(194)*v(273) v(334)=v(193)*v(273) v(328)=v(191)*v(273) v(322)=v(189)*v(273) v(274)=v(657)/v(265) v(345)=v(194)*v(274) v(336)=v(193)*v(274) v(330)=v(191)*v(274) v(324)=v(189)*v(274) v(281)=v(182)*v(266) v(282)=v(182)*v(267) v(283)=v(182)*v(268) v(284)=v(183)*v(269) v(285)=v(183)*v(270) v(286)=v(183)*v(271) v(287)=v(184)*v(272) v(288)=v(184)*v(273) v(289)=v(184)*v(274) v(290)=v(293)+v(287)*v(658) v(291)=v(295)+v(288)*v(658) v(292)=v(297)+v(289)*v(658) v(294)=-v(293)+v(287)*v(659) v(296)=-v(295)+v(288)*v(659) v(298)=-v(297)+v(289)*v(659) v(299)=v(269)+v(272) v(377)=v(266)-v(299) v(371)=-v(266)-v(299) v(300)=v(270)+v(273) v(379)=v(267)-v(300) v(373)=-v(267)-v(300) v(301)=v(271)+v(274) v(381)=v(268)-v(301) v(375)=-v(268)-v(301) v(302)=v(320)+v(284)*v(645) v(303)=v(322)+v(285)*v(645) v(304)=v(324)+v(286)*v(645) v(305)=v(326)+v(281)*v(645) v(306)=v(328)+v(282)*v(645) v(307)=v(330)+v(283)*v(645) v(308)=v(314)+v(332) v(309)=v(316)+v(334) v(310)=v(318)+v(336) v(311)=-v(269)+v(272) v(365)=-v(266)-v(311) v(359)=v(266)-v(311) v(312)=-v(270)+v(273) v(367)=-v(267)-v(312) v(361)=v(267)-v(312) v(313)=-v(271)+v(274) v(369)=-v(268)-v(313) v(363)=v(268)-v(313) v(315)=-v(314)+v(339) v(317)=-v(316)+v(342) v(319)=-v(318)+v(345) v(321)=-v(320)+v(284)*v(650) v(323)=-v(322)+v(285)*v(650) v(325)=-v(324)+v(286)*v(650) v(327)=-v(326)+v(281)*v(650) v(329)=-v(328)+v(282)*v(650) v(331)=-v(330)+v(283)*v(650) v(333)=-v(332)+v(338) v(335)=-v(334)+v(341) v(337)=-v(336)+v(344) v(340)=-v(338)-v(339) v(343)=-v(341)-v(342) v(346)=-v(344)-v(345) v(347)=-(v(200)*v(371))+v(198)*(v(17)*v(340)+v(660)) v(348)=-(v(200)*v(373))+v(198)*(v(17)*v(343)+v(661)) v(349)=-(v(200)*v(375))+v(198)*(v(17)*v(346)+v(662)) v(350)=v(205)*v(377)+v(203)*(v(13)*v(340)-v(660)) v(351)=v(205)*v(379)+v(203)*(v(13)*v(343)-v(661)) v(352)=v(205)*v(381)+v(203)*(v(13)*v(346)-v(662)) v(353)=v(210)*v(359)+v(208)*(v(13)*v(333)+v(663)) v(354)=v(210)*v(361)+v(208)*(v(13)*v(335)+v(664)) v(355)=v(210)*v(363)+v(208)*(v(13)*v(337)+v(665)) v(356)=-(v(215)*v(365))+v(213)*(v(17)*v(333)-v(663)) v(357)=-(v(215)*v(367))+v(213)*(v(17)*v(335)-v(664)) v(358)=-(v(215)*v(369))+v(213)*(v(17)*v(337)-v(665)) v(360)=-(v(22)*(v(218)*v(266)+v(17)*v(359)))+v(315)*v(649) v(362)=-(v(22)*(v(218)*v(267)+v(17)*v(361)))+v(317)*v(649) v(364)=-(v(22)*(v(218)*v(268)+v(17)*v(363)))+v(319)*v(649) v(366)=v(223)*(v(22)*v(266)+v(13)*v(315))-v(225)*v(365) v(368)=v(223)*(v(22)*v(267)+v(13)*v(317))-v(225)*v(367) v(370)=v(223)*(v(22)*v(268)+v(13)*v(319))-v(225)*v(369) v(372)=v(228)*(v(266)*v(28)+v(13)*v(308))-v(230)*v(371) v(374)=v(228)*(v(267)*v(28)+v(13)*v(309))-v(230)*v(373) v(376)=v(228)*(v(268)*v(28)+v(13)*v(310))-v(230)*v(375) v(378)=-(v(28)*(v(233)*v(266)+v(17)*v(377)))+v(308)*v(646) v(380)=-(v(28)*(v(233)*v(267)+v(17)*v(379)))+v(309)*v(646) v(382)=-(v(28)*(v(233)*v(268)+v(17)*v(381)))+v(310)*v(646) v(383)=v(20)*v(327)+v(389) v(384)=v(20)*v(329)+v(391) v(385)=v(20)*v(331)+v(393) v(386)=v(13)*v(321)+v(395) v(387)=v(13)*v(323)+v(397) v(388)=v(13)*v(325)+v(399) v(390)=v(25)*v(327)-v(389) v(392)=v(25)*v(329)-v(391) v(394)=v(25)*v(331)-v(393) v(396)=v(17)*v(321)-v(395) v(398)=v(17)*v(323)-v(397) v(400)=v(17)*v(325)-v(399) v(401)=v(17)*v(294)+v(404) v(402)=v(17)*v(296)+v(406) v(403)=v(17)*v(298)+v(408) v(405)=v(13)*v(294)-v(404) v(407)=v(13)*v(296)-v(406) v(409)=v(13)*v(298)-v(408) v(410)=v(13)*v(290)+v(413) v(411)=v(13)*v(291)+v(415) v(412)=v(13)*v(292)+v(417) v(414)=v(17)*v(290)-v(413) v(416)=v(17)*v(291)-v(415) v(418)=v(17)*v(292)-v(417) v(419)=v(20)*v(305)+v(425) v(420)=v(20)*v(306)+v(427) v(421)=v(20)*v(307)+v(429) v(422)=v(13)*v(302)+v(431) v(423)=v(13)*v(303)+v(433) v(424)=v(13)*v(304)+v(435) v(426)=v(25)*v(305)-v(425) v(428)=v(25)*v(306)-v(427) v(430)=v(25)*v(307)-v(429) v(432)=v(17)*v(302)-v(431) v(17367)=0. v(17368)=0. v(17369)=v(347) v(17370)=0. v(17371)=0. v(17372)=v(350) v(17373)=0. v(17374)=0. v(17375)=v(353) v(17376)=0. v(17377)=0. v(17378)=v(356) v(17379)=0. v(17380)=0. v(17381)=v(360) v(17382)=0. v(17383)=0. v(17384)=v(366) v(17385)=0. v(17386)=0. v(17387)=v(372) v(17388)=0. v(17389)=0. v(17390)=v(378) v(17391)=0. v(17392)=0. v(17393)=v(383) v(17394)=0. v(17395)=0. v(17396)=v(386) v(17397)=0. v(17398)=0. v(17399)=v(390) v(17400)=0. v(17401)=0. v(17402)=v(396) v(17403)=0. v(17404)=0. v(17405)=v(401) v(17406)=0. v(17407)=0. v(17408)=v(405) v(17409)=0. v(17410)=0. v(17411)=v(410) v(17412)=0. v(17413)=0. v(17414)=v(414) v(17415)=0. v(17416)=0. v(17417)=v(419) v(17418)=0. v(17419)=0. v(17420)=v(422) v(17421)=0. v(17422)=0. v(17423)=v(426) v(17424)=0. v(17425)=0. v(17426)=v(432) v(17187)=0. v(17188)=v(347) v(17189)=0. v(17190)=0. v(17191)=v(350) v(17192)=0. v(17193)=0. v(17194)=v(353) v(17195)=0. v(17196)=0. v(17197)=v(356) v(17198)=0. v(17199)=0. v(17200)=v(360) v(17201)=0. v(17202)=0. v(17203)=v(366) v(17204)=0. v(17205)=0. v(17206)=v(372) v(17207)=0. v(17208)=0. v(17209)=v(378) v(17210)=0. v(17211)=0. v(17212)=v(383) v(17213)=0. v(17214)=0. v(17215)=v(386) v(17216)=0. v(17217)=0. v(17218)=v(390) v(17219)=0. v(17220)=0. v(17221)=v(396) v(17222)=0. v(17223)=0. v(17224)=v(401) v(17225)=0. v(17226)=0. v(17227)=v(405) v(17228)=0. v(17229)=0. v(17230)=v(410) v(17231)=0. v(17232)=0. v(17233)=v(414) v(17234)=0. v(17235)=0. v(17236)=v(419) v(17237)=0. v(17238)=0. v(17239)=v(422) v(17240)=0. v(17241)=0. v(17242)=v(426) v(17243)=0. v(17244)=0. v(17245)=v(432) v(17246)=0. v(17007)=v(347) v(17008)=0. v(17009)=0. v(17010)=v(350) v(17011)=0. v(17012)=0. v(17013)=v(353) v(17014)=0. v(17015)=0. v(17016)=v(356) v(17017)=0. v(17018)=0. v(17019)=v(360) v(17020)=0. v(17021)=0. v(17022)=v(366) v(17023)=0. v(17024)=0. v(17025)=v(372) v(17026)=0. v(17027)=0. v(17028)=v(378) v(17029)=0. v(17030)=0. v(17031)=v(383) v(17032)=0. v(17033)=0. v(17034)=v(386) v(17035)=0. v(17036)=0. v(17037)=v(390) v(17038)=0. v(17039)=0. v(17040)=v(396) v(17041)=0. v(17042)=0. v(17043)=v(401) v(17044)=0. v(17045)=0. v(17046)=v(405) v(17047)=0. v(17048)=0. v(17049)=v(410) v(17050)=0. v(17051)=0. v(17052)=v(414) v(17053)=0. v(17054)=0. v(17055)=v(419) v(17056)=0. v(17057)=0. v(17058)=v(422) v(17059)=0. v(17060)=0. v(17061)=v(426) v(17062)=0. v(17063)=0. v(17064)=v(432) v(17065)=0. v(17066)=0. v(434)=v(17)*v(303)-v(433) v(17427)=0. v(17428)=0. v(17429)=v(348) v(17430)=0. v(17431)=0. v(17432)=v(351) v(17433)=0. v(17434)=0. v(17435)=v(354) v(17436)=0. v(17437)=0. v(17438)=v(357) v(17439)=0. v(17440)=0. v(17441)=v(362) v(17442)=0. v(17443)=0. v(17444)=v(368) v(17445)=0. v(17446)=0. v(17447)=v(374) v(17448)=0. v(17449)=0. v(17450)=v(380) v(17451)=0. v(17452)=0. v(17453)=v(384) v(17454)=0. v(17455)=0. v(17456)=v(387) v(17457)=0. v(17458)=0. v(17459)=v(392) v(17460)=0. v(17461)=0. v(17462)=v(398) v(17463)=0. v(17464)=0. v(17465)=v(402) v(17466)=0. v(17467)=0. v(17468)=v(407) v(17469)=0. v(17470)=0. v(17471)=v(411) v(17472)=0. v(17473)=0. v(17474)=v(416) v(17475)=0. v(17476)=0. v(17477)=v(420) v(17478)=0. v(17479)=0. v(17480)=v(423) v(17481)=0. v(17482)=0. v(17483)=v(428) v(17484)=0. v(17485)=0. v(17486)=v(434) v(17247)=0. v(17248)=v(348) v(17249)=0. v(17250)=0. v(17251)=v(351) v(17252)=0. v(17253)=0. v(17254)=v(354) v(17255)=0. v(17256)=0. v(17257)=v(357) v(17258)=0. v(17259)=0. v(17260)=v(362) v(17261)=0. v(17262)=0. v(17263)=v(368) v(17264)=0. v(17265)=0. v(17266)=v(374) v(17267)=0. v(17268)=0. v(17269)=v(380) v(17270)=0. v(17271)=0. v(17272)=v(384) v(17273)=0. v(17274)=0. v(17275)=v(387) v(17276)=0. v(17277)=0. v(17278)=v(392) v(17279)=0. v(17280)=0. v(17281)=v(398) v(17282)=0. v(17283)=0. v(17284)=v(402) v(17285)=0. v(17286)=0. v(17287)=v(407) v(17288)=0. v(17289)=0. v(17290)=v(411) v(17291)=0. v(17292)=0. v(17293)=v(416) v(17294)=0. v(17295)=0. v(17296)=v(420) v(17297)=0. v(17298)=0. v(17299)=v(423) v(17300)=0. v(17301)=0. v(17302)=v(428) v(17303)=0. v(17304)=0. v(17305)=v(434) v(17306)=0. v(17067)=v(348) v(17068)=0. v(17069)=0. v(17070)=v(351) v(17071)=0. v(17072)=0. v(17073)=v(354) v(17074)=0. v(17075)=0. v(17076)=v(357) v(17077)=0. v(17078)=0. v(17079)=v(362) v(17080)=0. v(17081)=0. v(17082)=v(368) v(17083)=0. v(17084)=0. v(17085)=v(374) v(17086)=0. v(17087)=0. v(17088)=v(380) v(17089)=0. v(17090)=0. v(17091)=v(384) v(17092)=0. v(17093)=0. v(17094)=v(387) v(17095)=0. v(17096)=0. v(17097)=v(392) v(17098)=0. v(17099)=0. v(17100)=v(398) v(17101)=0. v(17102)=0. v(17103)=v(402) v(17104)=0. v(17105)=0. v(17106)=v(407) v(17107)=0. v(17108)=0. v(17109)=v(411) v(17110)=0. v(17111)=0. v(17112)=v(416) v(17113)=0. v(17114)=0. v(17115)=v(420) v(17116)=0. v(17117)=0. v(17118)=v(423) v(17119)=0. v(17120)=0. v(17121)=v(428) v(17122)=0. v(17123)=0. v(17124)=v(434) v(17125)=0. v(17126)=0. v(436)=v(17)*v(304)-v(435) v(17487)=0. v(17488)=0. v(17489)=v(349) v(17490)=0. v(17491)=0. v(17492)=v(352) v(17493)=0. v(17494)=0. v(17495)=v(355) v(17496)=0. v(17497)=0. v(17498)=v(358) v(17499)=0. v(17500)=0. v(17501)=v(364) v(17502)=0. v(17503)=0. v(17504)=v(370) v(17505)=0. v(17506)=0. v(17507)=v(376) v(17508)=0. v(17509)=0. v(17510)=v(382) v(17511)=0. v(17512)=0. v(17513)=v(385) v(17514)=0. v(17515)=0. v(17516)=v(388) v(17517)=0. v(17518)=0. v(17519)=v(394) v(17520)=0. v(17521)=0. v(17522)=v(400) v(17523)=0. v(17524)=0. v(17525)=v(403) v(17526)=0. v(17527)=0. v(17528)=v(409) v(17529)=0. v(17530)=0. v(17531)=v(412) v(17532)=0. v(17533)=0. v(17534)=v(418) v(17535)=0. v(17536)=0. v(17537)=v(421) v(17538)=0. v(17539)=0. v(17540)=v(424) v(17541)=0. v(17542)=0. v(17543)=v(430) v(17544)=0. v(17545)=0. v(17546)=v(436) v(17307)=0. v(17308)=v(349) v(17309)=0. v(17310)=0. v(17311)=v(352) v(17312)=0. v(17313)=0. v(17314)=v(355) v(17315)=0. v(17316)=0. v(17317)=v(358) v(17318)=0. v(17319)=0. v(17320)=v(364) v(17321)=0. v(17322)=0. v(17323)=v(370) v(17324)=0. v(17325)=0. v(17326)=v(376) v(17327)=0. v(17328)=0. v(17329)=v(382) v(17330)=0. v(17331)=0. v(17332)=v(385) v(17333)=0. v(17334)=0. v(17335)=v(388) v(17336)=0. v(17337)=0. v(17338)=v(394) v(17339)=0. v(17340)=0. v(17341)=v(400) v(17342)=0. v(17343)=0. v(17344)=v(403) v(17345)=0. v(17346)=0. v(17347)=v(409) v(17348)=0. v(17349)=0. v(17350)=v(412) v(17351)=0. v(17352)=0. v(17353)=v(418) v(17354)=0. v(17355)=0. v(17356)=v(421) v(17357)=0. v(17358)=0. v(17359)=v(424) v(17360)=0. v(17361)=0. v(17362)=v(430) v(17363)=0. v(17364)=0. v(17365)=v(436) v(17366)=0. v(17127)=v(349) v(17128)=0. v(17129)=0. v(17130)=v(352) v(17131)=0. v(17132)=0. v(17133)=v(355) v(17134)=0. v(17135)=0. v(17136)=v(358) v(17137)=0. v(17138)=0. v(17139)=v(364) v(17140)=0. v(17141)=0. v(17142)=v(370) v(17143)=0. v(17144)=0. v(17145)=v(376) v(17146)=0. v(17147)=0. v(17148)=v(382) v(17149)=0. v(17150)=0. v(17151)=v(385) v(17152)=0. v(17153)=0. v(17154)=v(388) v(17155)=0. v(17156)=0. v(17157)=v(394) v(17158)=0. v(17159)=0. v(17160)=v(400) v(17161)=0. v(17162)=0. v(17163)=v(403) v(17164)=0. v(17165)=0. v(17166)=v(409) v(17167)=0. v(17168)=0. v(17169)=v(412) v(17170)=0. v(17171)=0. v(17172)=v(418) v(17173)=0. v(17174)=0. v(17175)=v(421) v(17176)=0. v(17177)=0. v(17178)=v(424) v(17179)=0. v(17180)=0. v(17181)=v(430) v(17182)=0. v(17183)=0. v(17184)=v(436) v(17185)=0. v(17186)=0. v(452)=1.+v(432)*v(524)+v(426)*v(527)+v(422)*v(530)+v(419)*v & (533)+v(414)*v(536)+v(410)*v(539)+v(405)*v(542)+v(401)*v(545) & +v(396)*v(548)+v(390)*v(551)+v(386)*v(554)+v(383)*v(557)+v(378 & )*v(560)+v(372)*v(563)+v(366)*v(566)+v(360)*v(569)+v(356)*v & (572)+v(353)*v(575)+v(350)*v(578)+v(347)*v(581) v(453)=v(434)*v(524)+v(428)*v(527)+v(423)*v(530)+v(420)*v(533) & +v(416)*v(536)+v(411)*v(539)+v(407)*v(542)+v(402)*v(545)+v(398 & )*v(548)+v(392)*v(551)+v(387)*v(554)+v(384)*v(557)+v(380)*v & (560)+v(374)*v(563)+v(368)*v(566)+v(362)*v(569)+v(357)*v(572) & +v(354)*v(575)+v(351)*v(578)+v(348)*v(581) v(454)=v(436)*v(524)+v(430)*v(527)+v(424)*v(530)+v(421)*v(533) & +v(418)*v(536)+v(412)*v(539)+v(409)*v(542)+v(403)*v(545)+v(400 & )*v(548)+v(394)*v(551)+v(388)*v(554)+v(385)*v(557)+v(382)*v & (560)+v(376)*v(563)+v(370)*v(566)+v(364)*v(569)+v(358)*v(572) & +v(355)*v(575)+v(352)*v(578)+v(349)*v(581) v(455)=v(432)*v(523)+v(426)*v(526)+v(422)*v(529)+v(419)*v(532) & +v(414)*v(535)+v(410)*v(538)+v(405)*v(541)+v(401)*v(544)+v(396 & )*v(547)+v(390)*v(550)+v(386)*v(553)+v(383)*v(556)+v(378)*v & (559)+v(372)*v(562)+v(366)*v(565)+v(360)*v(568)+v(356)*v(571) & +v(353)*v(574)+v(350)*v(577)+v(347)*v(580) v(456)=1.+v(434)*v(523)+v(428)*v(526)+v(423)*v(529)+v(420)*v & (532)+v(416)*v(535)+v(411)*v(538)+v(407)*v(541)+v(402)*v(544) & +v(398)*v(547)+v(392)*v(550)+v(387)*v(553)+v(384)*v(556)+v(380 & )*v(559)+v(374)*v(562)+v(368)*v(565)+v(362)*v(568)+v(357)*v & (571)+v(354)*v(574)+v(351)*v(577)+v(348)*v(580) v(508)=-(v(453)*v(455))+v(452)*v(456) v(457)=v(436)*v(523)+v(430)*v(526)+v(424)*v(529)+v(421)*v(532) & +v(418)*v(535)+v(412)*v(538)+v(409)*v(541)+v(403)*v(544)+v(400 & )*v(547)+v(394)*v(550)+v(388)*v(553)+v(385)*v(556)+v(382)*v & (559)+v(376)*v(562)+v(370)*v(565)+v(364)*v(568)+v(358)*v(571) & +v(355)*v(574)+v(352)*v(577)+v(349)*v(580) v(507)=v(454)*v(455)-v(452)*v(457) v(506)=-(v(454)*v(456))+v(453)*v(457) v(458)=v(432)*v(522)+v(426)*v(525)+v(422)*v(528)+v(419)*v(531) & +v(414)*v(534)+v(410)*v(537)+v(405)*v(540)+v(401)*v(543)+v(396 & )*v(546)+v(390)*v(549)+v(386)*v(552)+v(383)*v(555)+v(378)*v & (558)+v(372)*v(561)+v(366)*v(564)+v(360)*v(567)+v(356)*v(570) & +v(353)*v(573)+v(350)*v(576)+v(347)*v(579) v(459)=v(434)*v(522)+v(428)*v(525)+v(423)*v(528)+v(420)*v(531) & +v(416)*v(534)+v(411)*v(537)+v(407)*v(540)+v(402)*v(543)+v(398 & )*v(546)+v(392)*v(549)+v(387)*v(552)+v(384)*v(555)+v(380)*v & (558)+v(374)*v(561)+v(368)*v(564)+v(362)*v(567)+v(357)*v(570) & +v(354)*v(573)+v(351)*v(576)+v(348)*v(579) v(505)=v(453)*v(458)-v(452)*v(459) v(502)=-(v(456)*v(458))+v(455)*v(459) v(460)=1.+v(436)*v(522)+v(430)*v(525)+v(424)*v(528)+v(421)*v & (531)+v(418)*v(534)+v(412)*v(537)+v(409)*v(540)+v(403)*v(543) & +v(400)*v(546)+v(394)*v(549)+v(388)*v(552)+v(385)*v(555)+v(382 & )*v(558)+v(376)*v(561)+v(370)*v(564)+v(364)*v(567)+v(358)*v & (570)+v(355)*v(573)+v(352)*v(576)+v(349)*v(579) v(504)=-(v(454)*v(458))+v(452)*v(460) v(503)=v(454)*v(459)-v(453)*v(460) v(501)=v(457)*v(458)-v(455)*v(460) v(500)=-(v(457)*v(459))+v(456)*v(460) v(467)=v(458)*v(506)+v(459)*v(507)+v(460)*v(508) v(675)=v(451)/v(467)**2 v(472)=v(449)*(-1.+v(467))-v(451)/v(467) v(473)=v(451)*v(460)+v(472)*v(508) v(474)=v(451)*v(459)+v(472)*v(507) v(475)=v(451)*v(458)+v(472)*v(506) v(476)=v(451)*v(457)+v(472)*v(505) v(477)=v(451)*v(456)+v(472)*v(504) v(478)=v(451)*v(455)+v(472)*v(503) v(479)=v(451)*v(454)+v(472)*v(502) v(480)=v(451)*v(453)+v(472)*v(501) v(481)=v(451)*v(452)+v(472)*v(500) v(16943)=v(349)*v(479)+v(348)*v(480)+v(347)*v(481) v(16944)=v(349)*v(476)+v(348)*v(477)+v(347)*v(478) v(16945)=v(349)*v(473)+v(348)*v(474)+v(347)*v(475) v(16946)=v(352)*v(479)+v(351)*v(480)+v(350)*v(481) v(16947)=v(352)*v(476)+v(351)*v(477)+v(350)*v(478) v(16948)=v(352)*v(473)+v(351)*v(474)+v(350)*v(475) v(16949)=v(355)*v(479)+v(354)*v(480)+v(353)*v(481) v(16950)=v(355)*v(476)+v(354)*v(477)+v(353)*v(478) v(16951)=v(355)*v(473)+v(354)*v(474)+v(353)*v(475) v(16952)=v(358)*v(479)+v(357)*v(480)+v(356)*v(481) v(16953)=v(358)*v(476)+v(357)*v(477)+v(356)*v(478) v(16954)=v(358)*v(473)+v(357)*v(474)+v(356)*v(475) v(16955)=v(364)*v(479)+v(362)*v(480)+v(360)*v(481) v(16956)=v(364)*v(476)+v(362)*v(477)+v(360)*v(478) v(16957)=v(364)*v(473)+v(362)*v(474)+v(360)*v(475) v(16958)=v(370)*v(479)+v(368)*v(480)+v(366)*v(481) v(16959)=v(370)*v(476)+v(368)*v(477)+v(366)*v(478) v(16960)=v(370)*v(473)+v(368)*v(474)+v(366)*v(475) v(16961)=v(376)*v(479)+v(374)*v(480)+v(372)*v(481) v(16962)=v(376)*v(476)+v(374)*v(477)+v(372)*v(478) v(16963)=v(376)*v(473)+v(374)*v(474)+v(372)*v(475) v(16964)=v(382)*v(479)+v(380)*v(480)+v(378)*v(481) v(16965)=v(382)*v(476)+v(380)*v(477)+v(378)*v(478) v(16966)=v(382)*v(473)+v(380)*v(474)+v(378)*v(475) v(16967)=v(385)*v(479)+v(384)*v(480)+v(383)*v(481) v(16968)=v(385)*v(476)+v(384)*v(477)+v(383)*v(478) v(16969)=v(385)*v(473)+v(384)*v(474)+v(383)*v(475) v(16970)=v(388)*v(479)+v(387)*v(480)+v(386)*v(481) v(16971)=v(388)*v(476)+v(387)*v(477)+v(386)*v(478) v(16972)=v(388)*v(473)+v(387)*v(474)+v(386)*v(475) v(16973)=v(394)*v(479)+v(392)*v(480)+v(390)*v(481) v(16974)=v(394)*v(476)+v(392)*v(477)+v(390)*v(478) v(16975)=v(394)*v(473)+v(392)*v(474)+v(390)*v(475) v(16976)=v(400)*v(479)+v(398)*v(480)+v(396)*v(481) v(16977)=v(400)*v(476)+v(398)*v(477)+v(396)*v(478) v(16978)=v(400)*v(473)+v(398)*v(474)+v(396)*v(475) v(16979)=v(403)*v(479)+v(402)*v(480)+v(401)*v(481) v(16980)=v(403)*v(476)+v(402)*v(477)+v(401)*v(478) v(16981)=v(403)*v(473)+v(402)*v(474)+v(401)*v(475) v(16982)=v(409)*v(479)+v(407)*v(480)+v(405)*v(481) v(16983)=v(409)*v(476)+v(407)*v(477)+v(405)*v(478) v(16984)=v(409)*v(473)+v(407)*v(474)+v(405)*v(475) v(16985)=v(412)*v(479)+v(411)*v(480)+v(410)*v(481) v(16986)=v(412)*v(476)+v(411)*v(477)+v(410)*v(478) v(16987)=v(412)*v(473)+v(411)*v(474)+v(410)*v(475) v(16988)=v(418)*v(479)+v(416)*v(480)+v(414)*v(481) v(16989)=v(418)*v(476)+v(416)*v(477)+v(414)*v(478) v(16990)=v(418)*v(473)+v(416)*v(474)+v(414)*v(475) v(16991)=v(421)*v(479)+v(420)*v(480)+v(419)*v(481) v(16992)=v(421)*v(476)+v(420)*v(477)+v(419)*v(478) v(16993)=v(421)*v(473)+v(420)*v(474)+v(419)*v(475) v(16994)=v(424)*v(479)+v(423)*v(480)+v(422)*v(481) v(16995)=v(424)*v(476)+v(423)*v(477)+v(422)*v(478) v(16996)=v(424)*v(473)+v(423)*v(474)+v(422)*v(475) v(16997)=v(430)*v(479)+v(428)*v(480)+v(426)*v(481) v(16998)=v(430)*v(476)+v(428)*v(477)+v(426)*v(478) v(16999)=v(430)*v(473)+v(428)*v(474)+v(426)*v(475) v(17000)=v(436)*v(479)+v(434)*v(480)+v(432)*v(481) v(17001)=v(436)*v(476)+v(434)*v(477)+v(432)*v(478) v(17002)=v(436)*v(473)+v(434)*v(474)+v(432)*v(475) DO i469=1,60 v(491)=v(265)*v(17006+i469) v(492)=v(265)*v(17066+i469) v(493)=v(265)*v(17126+i469) v(494)=v(265)*v(17186+i469) v(495)=v(265)*v(17246+i469) v(496)=v(265)*v(17306+i469) v(497)=v(265)*v(17366+i469) v(498)=v(265)*v(17426+i469) v(499)=v(265)*v(17486+i469) v(509)=(v(491)*v(500)+v(492)*v(501)+v(493)*v(502)+v(494)*v & (503)+v(495)*v(504)+v(496)*v(505)+v(497)*v(506)+v(498)*v(507) & +v(499)*v(508))*(v(449)+v(675)) v(510)=v(472)*(v(456)*v(491)-v(455)*v(492)-v(453)*v(494)+v & (452)*v(495))+v(451)*v(499)+v(508)*v(509) v(511)=v(472)*(-(v(457)*v(491))+v(455)*v(493)+v(454)*v(494)-v & (452)*v(496))+v(451)*v(498)+v(507)*v(509) v(512)=v(472)*(v(457)*v(492)-v(456)*v(493)-v(454)*v(495)+v & (453)*v(496))+v(451)*v(497)+v(506)*v(509) v(513)=v(451)*v(496)+v(472)*(-(v(459)*v(491))+v(458)*v(492)+v & (453)*v(497)-v(452)*v(498))+v(505)*v(509) v(514)=v(451)*v(495)+v(472)*(v(460)*v(491)-v(458)*v(493)-v & (454)*v(497)+v(452)*v(499))+v(504)*v(509) v(515)=v(451)*v(494)+v(472)*(-(v(460)*v(492))+v(459)*v(493)+v & (454)*v(498)-v(453)*v(499))+v(503)*v(509) v(516)=v(451)*v(493)+v(472)*(v(459)*v(494)-v(458)*v(495)-v & (456)*v(497)+v(455)*v(498))+v(502)*v(509) v(517)=v(451)*v(492)+v(472)*(-(v(460)*v(494))+v(458)*v(496)+v & (457)*v(497)-v(455)*v(499))+v(501)*v(509) v(518)=v(451)*v(491)+v(472)*(v(460)*v(495)-v(459)*v(496)-v & (457)*v(498)+v(456)*v(499))+v(500)*v(509) v(17547)=v(349)*v(516)+v(348)*v(517)+v(347)*v(518) v(17548)=v(349)*v(513)+v(348)*v(514)+v(347)*v(515) v(17549)=v(349)*v(510)+v(348)*v(511)+v(347)*v(512) v(17550)=v(352)*v(516)+v(351)*v(517)+v(350)*v(518) v(17551)=v(352)*v(513)+v(351)*v(514)+v(350)*v(515) v(17552)=v(352)*v(510)+v(351)*v(511)+v(350)*v(512) v(17553)=v(355)*v(516)+v(354)*v(517)+v(353)*v(518) v(17554)=v(355)*v(513)+v(354)*v(514)+v(353)*v(515) v(17555)=v(355)*v(510)+v(354)*v(511)+v(353)*v(512) v(17556)=v(358)*v(516)+v(357)*v(517)+v(356)*v(518) v(17557)=v(358)*v(513)+v(357)*v(514)+v(356)*v(515) v(17558)=v(358)*v(510)+v(357)*v(511)+v(356)*v(512) v(17559)=v(364)*v(516)+v(362)*v(517)+v(360)*v(518) v(17560)=v(364)*v(513)+v(362)*v(514)+v(360)*v(515) v(17561)=v(364)*v(510)+v(362)*v(511)+v(360)*v(512) v(17562)=v(370)*v(516)+v(368)*v(517)+v(366)*v(518) v(17563)=v(370)*v(513)+v(368)*v(514)+v(366)*v(515) v(17564)=v(370)*v(510)+v(368)*v(511)+v(366)*v(512) v(17565)=v(376)*v(516)+v(374)*v(517)+v(372)*v(518) v(17566)=v(376)*v(513)+v(374)*v(514)+v(372)*v(515) v(17567)=v(376)*v(510)+v(374)*v(511)+v(372)*v(512) v(17568)=v(382)*v(516)+v(380)*v(517)+v(378)*v(518) v(17569)=v(382)*v(513)+v(380)*v(514)+v(378)*v(515) v(17570)=v(382)*v(510)+v(380)*v(511)+v(378)*v(512) v(17571)=v(385)*v(516)+v(384)*v(517)+v(383)*v(518) v(17572)=v(385)*v(513)+v(384)*v(514)+v(383)*v(515) v(17573)=v(385)*v(510)+v(384)*v(511)+v(383)*v(512) v(17574)=v(388)*v(516)+v(387)*v(517)+v(386)*v(518) v(17575)=v(388)*v(513)+v(387)*v(514)+v(386)*v(515) v(17576)=v(388)*v(510)+v(387)*v(511)+v(386)*v(512) v(17577)=v(394)*v(516)+v(392)*v(517)+v(390)*v(518) v(17578)=v(394)*v(513)+v(392)*v(514)+v(390)*v(515) v(17579)=v(394)*v(510)+v(392)*v(511)+v(390)*v(512) v(17580)=v(400)*v(516)+v(398)*v(517)+v(396)*v(518) v(17581)=v(400)*v(513)+v(398)*v(514)+v(396)*v(515) v(17582)=v(400)*v(510)+v(398)*v(511)+v(396)*v(512) v(17583)=v(403)*v(516)+v(402)*v(517)+v(401)*v(518) v(17584)=v(403)*v(513)+v(402)*v(514)+v(401)*v(515) v(17585)=v(403)*v(510)+v(402)*v(511)+v(401)*v(512) v(17586)=v(409)*v(516)+v(407)*v(517)+v(405)*v(518) v(17587)=v(409)*v(513)+v(407)*v(514)+v(405)*v(515) v(17588)=v(409)*v(510)+v(407)*v(511)+v(405)*v(512) v(17589)=v(412)*v(516)+v(411)*v(517)+v(410)*v(518) v(17590)=v(412)*v(513)+v(411)*v(514)+v(410)*v(515) v(17591)=v(412)*v(510)+v(411)*v(511)+v(410)*v(512) v(17592)=v(418)*v(516)+v(416)*v(517)+v(414)*v(518) v(17593)=v(418)*v(513)+v(416)*v(514)+v(414)*v(515) v(17594)=v(418)*v(510)+v(416)*v(511)+v(414)*v(512) v(17595)=v(421)*v(516)+v(420)*v(517)+v(419)*v(518) v(17596)=v(421)*v(513)+v(420)*v(514)+v(419)*v(515) v(17597)=v(421)*v(510)+v(420)*v(511)+v(419)*v(512) v(17598)=v(424)*v(516)+v(423)*v(517)+v(422)*v(518) v(17599)=v(424)*v(513)+v(423)*v(514)+v(422)*v(515) v(17600)=v(424)*v(510)+v(423)*v(511)+v(422)*v(512) v(17601)=v(430)*v(516)+v(428)*v(517)+v(426)*v(518) v(17602)=v(430)*v(513)+v(428)*v(514)+v(426)*v(515) v(17603)=v(430)*v(510)+v(428)*v(511)+v(426)*v(512) v(17604)=v(436)*v(516)+v(434)*v(517)+v(432)*v(518) v(17605)=v(436)*v(513)+v(434)*v(514)+v(432)*v(515) v(17606)=v(436)*v(510)+v(434)*v(511)+v(432)*v(512) IF(.not.(b1)) THEN p(i469)=p(i469)+v(265)*v(8)*v(16942+i469) ELSE ENDIF IF(b2) THEN cycle ELSE ENDIF DO i488=i469,60 s(i469,i488)=s(i469,i488)+v(8)*v(17546+i488) ENDDO ENDDO ENDDO END