!************************************************************** !* AceGen 6.405 Windows (22 Jun 15) * !* Co. J. Korelc 2013 24 Jun 15 07:27:26 * !************************************************************** ! User : USER ! Notebook : Remote Session ! Evaluation time : 56 s Mode : Optimal ! Number of formulae : 942 Method: Automatic ! Subroutine : SKR10 size: 11019 ! Subroutine : SPP10 size: 7489 ! Total size of Mathematica code : 18508 subexpressions ! Total size of Fortran code : 64137 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: ! 1 "Sxx" ! 2 "Sxy" ! 3 "Sxz" ! 4 "Syx" ! 5 "Syy" ! 6 "Syz" ! 7 "Szx" ! 8 "Szy" ! 9 "Szz" ! 10 "Exx" ! 11 "Exy" ! 12 "Exz" ! 13 "Eyx" ! 14 "Eyy" ! 15 "Eyz" ! 16 "Ezx" ! 17 "Ezy" ! 18 "Ezz" ! 19 "Mises stress" ! Nodal point postprocessing quantities: ! 1 "DeformedMeshX" ! 2 "DeformedMeshY" ! 3 "DeformedMeshZ" ! 4 "u" ! 5 "v" ! 6 "w" 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(19) parameter (DEBUG=.false., # SELEM="SED3H2SDFHYH2SNeoHooke") integer ix(nen),ndme,ndfe,nste,isw,ngdof,ncdof integer maxdof,dofglobal(20),nnodes,nalldof,ndim double precision xl(ndme,nen),d(*),ul(ndfe,nen,*) double precision s(nste,nste),p(nste),tl(nen),sxd(60) double precision sall(60,60),pall(60),du(60),dh(1) double precision ul0(ndfe,nen),sg(60),sg0(60) integer i,j,jj,ll,ii,k,kk,i1,i2,i3,hlen,icode integer cdata1,cdata2,cdata3,cdata4 double precision w,v(17608),gpost(64,19),npost(20,6) 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,"> ",60g11.5) 1238 format(i3,"> ",60g11.5) 1239 format(i3,"> ",3g11.5) icode=8 ncdof=0 ngdof=60 nalldof=60 symmetric=.true. maxdof=3 nnodes=20 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=0 endif 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 if(isw.eq.1)then c c.... input record 1 of material properties call dinput(d(1),2) call SMSIntPoints(icode,ngpo,gp) write(*,*) SELEM write(iow,*)SELEM c.....Description of the input data datades(1)="E -elastic modulus" datades(2)="[Nu] -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 j=3 call dinput(d(j),1) idata(ID_NoSensParameters)=int(d(j)) if(idata(ID_NoSensParameters).gt.0) then k=0 do i=10,idata(ID_NoSensParameters),10 call dinput(d(j),10) j=j+10 k=k+1 enddo call dinput(d(j),idata(ID_NoSensParameters)-10*k) write(iow,*)'Sensitivity parameters:' write(iow,"(5(i5,2x,f15.5))") # (i,d(2+i),i=1,idata(ID_NoSensParameters)) endif nsenpa=idata(ID_NoSensParameters) mct=0 c.... number of projected postprocess quantities istv=19 c.... description of the postprocessing data postdes(1)="Sxx" postdes(2)="Sxy" postdes(3)="Sxz" postdes(4)="Syx" postdes(5)="Syy" postdes(6)="Syz" postdes(7)="Szx" postdes(8)="Szy" postdes(9)="Szz" postdes(10)="Exx" postdes(11)="Exy" postdes(12)="Exz" postdes(13)="Eyx" postdes(14)="Eyy" postdes(15)="Eyz" postdes(16)="Ezx" postdes(17)="Ezy" postdes(18)="Ezz" postdes(19)="Mises stress" write(iow,"(A)")"Numbering of the postprocessing data:" write(iow,"(10x,i15,A3,A50)") # (i," = ",postdes(i),i=1,19) write(iow,2004)(postdes(i)(1:5),i=1,19) 2004 format(" ",6(A5,6X) # ,/,12X,6(A5,6X),/,12X,6(A5,6X),/,12X,6(A5,6X)) idata(ID_OutputFile)=iow idata(ID_PostIteration)=0 idata(ID_SkipTangent)=0 idata(ID_SkipResidual)=0 idata(ID_SkipSubIteration)=0; return elseif(isw.eq.2) then write(*,*)"User switch 2 not implemented" return elseif(isw.eq.3.or.isw.eq.6) then 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 elseif(isw.eq.5) then write(*,*)"User switch 5 not implemented" return elseif(isw.eq.7) then write(*,*)"User switch 7 not implemented" return c.... postprocessing elseif(isw.eq.4.or.isw.eq.8) then c.....Description of the post-processing data call SPP10(v,d,ul,ul0,xl,s,p,hr(nh2),hr(nh1), # sg,sg0,sxd,gpost,npost) Do j=1,ngpo if(isw.eq.4) then c..... textual output mct = mct - 2 if(mct.le.0) then if(ior.lt.0) then write(*,2001)o,head,(postdes(i)(1:5),i=1,19) endif mct = 60 endif write(iow,2003) n,j,(gpost(j,i),i=1,19) if(ior.lt.0) then write(*,2002) n,j,(gpost(j,i),i=1,19) endif 2001 format(a1,20a4//," Elm. Ptn. ",6(A5,6X) # ,/,12X,6(A5,6X),/,12X,6(A5,6X),/,12X,6(A5,6X)) 2002 format(2i4,6e11.4/8x,6e11.4/8x,6e11.4) 2003 format("@E/ip",i5,i5,6e12.4/15x,6e12.4/15x,6e12.4/15x,6e12.4) else c..... graphics postprocessing do jj=1,nnodes ll = ix(jj) if(ll.gt.0) then c.. w - nodal weight for post-processing w=1.d0 hr(nph+ll-1)=hr(nph+ll-1)+w hr(ner+ll-1)=hr(ner+ll-1)+erav*w do ii=1,19 k=nph+numnp+numnp*(ii-1)+ll-1 hr(k)=hr(k)+w*gpost(j,ii) enddo endif enddo endif enddo if(isw.eq.4) write(iow,"(a)")"@" return c.....sensitivity analysis - external elseif(isw.eq.20) then return elseif(isw.eq.21) then c..... internal sensitivity return endif 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 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(1135)=xl(1,1) v(1134)=xl(2,1) v(1133)=xl(3,1) v(1132)=xl(1,2) v(1131)=xl(2,2) v(1130)=xl(3,2) v(1129)=xl(1,3) v(1128)=xl(2,3) v(1127)=xl(3,3) v(1126)=xl(1,4) v(1125)=xl(2,4) v(1124)=xl(3,4) v(1123)=xl(1,5) v(1122)=xl(2,5) v(1121)=xl(3,5) v(1120)=xl(1,6) v(1119)=xl(2,6) v(1118)=xl(3,6) v(1117)=xl(1,7) v(1116)=xl(2,7) v(1115)=xl(3,7) v(1114)=xl(1,8) v(1113)=xl(2,8) v(1112)=xl(3,8) v(1111)=xl(1,9) v(1110)=xl(2,9) v(1109)=xl(3,9) v(1108)=xl(1,10) v(1107)=xl(2,10) v(1106)=xl(3,10) v(1105)=xl(1,11) v(1104)=xl(2,11) v(1103)=xl(3,11) v(1102)=xl(1,12) v(1101)=xl(2,12) v(1100)=xl(3,12) v(1099)=xl(1,13) v(1098)=xl(2,13) v(1097)=xl(3,13) v(1096)=xl(1,14) v(1095)=xl(2,14) v(1094)=xl(3,14) v(1093)=xl(1,15) v(1092)=xl(2,15) v(1091)=xl(3,15) v(1090)=xl(1,16) v(1089)=xl(2,16) v(1088)=xl(3,16) v(1087)=xl(1,17) v(1086)=xl(2,17) v(1085)=xl(3,17) v(1084)=xl(1,18) v(1083)=xl(2,18) v(1082)=xl(3,18) v(1081)=xl(1,19) v(1080)=xl(2,19) v(1079)=xl(3,19) v(1078)=xl(1,20) v(1077)=xl(2,20) v(1076)=xl(3,20) v(1075)=2.*d(2) v(1074)=1d0/(1.+d(2)) v(1073)=ul(1,1) v(1072)=ul(2,1) v(1071)=ul(3,1) v(1070)=ul(1,2) v(1069)=ul(2,2) v(1068)=ul(3,2) v(1067)=ul(1,3) v(1066)=ul(2,3) v(1065)=ul(3,3) v(1064)=ul(1,4) v(1063)=ul(2,4) v(1062)=ul(3,4) v(1061)=ul(1,5) v(1060)=ul(2,5) v(1059)=ul(3,5) v(1058)=ul(1,6) v(1057)=ul(2,6) v(1056)=ul(3,6) v(1055)=ul(1,7) v(1054)=ul(2,7) v(1053)=ul(3,7) v(1052)=ul(1,8) v(1051)=ul(2,8) v(1050)=ul(3,8) v(1049)=ul(1,9) v(1048)=ul(2,9) v(1047)=ul(3,9) v(1046)=ul(1,10) v(1045)=ul(2,10) v(1044)=ul(3,10) v(1043)=ul(1,11) v(1042)=ul(2,11) v(1041)=ul(3,11) v(1040)=ul(1,12) v(1039)=ul(2,12) v(1038)=ul(3,12) v(1037)=ul(1,13) v(1036)=ul(2,13) v(1035)=ul(3,13) v(1034)=ul(1,14) v(1033)=ul(2,14) v(1032)=ul(3,14) v(1031)=ul(1,15) v(1030)=ul(2,15) v(1029)=ul(3,15) v(1028)=ul(1,16) v(1027)=ul(2,16) v(1026)=ul(3,16) v(1025)=ul(1,17) v(1024)=ul(2,17) v(1023)=ul(3,17) v(1022)=ul(1,18) v(1021)=ul(2,18) v(1020)=ul(3,18) v(1019)=ul(1,19) v(1018)=ul(2,19) v(1017)=ul(3,19) v(1016)=ul(1,20) v(1015)=ul(2,20) v(1014)=ul(3,20) v(451)=(d(1)*v(1074))/2. v(449)=(v(1075)*v(451))/(1.-v(1075)) 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(1150)=v(25)/4. v(243)=-(v(191)*v(25)) v(193)=v(25)/8. v(20)=1.-v(6) v(1151)=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(1136)=v(184)/4. v(187)=v(1136)*v(20) v(249)=v(13)*v(187) v(247)=v(17)*v(187) v(186)=v(1136)*v(25) 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(1138)=v(17)*v(233) v(228)=v(29)+v(5) v(1139)=v(13)*v(228) v(26)=1.+v(7) v(1137)=v(26)/4. v(192)=v(26)/8. v(190)=v(1137)*v(182) v(254)=v(190)*v(25) v(252)=v(190)*v(20) v(188)=v(1137)*v(183) v(255)=v(17)*v(188) v(253)=v(13)*v(188) v(50)=v(1137)*v(46) v(48)=v(1137)*v(44) v(28)=v(192)*v(25) v(235)=-(v(17)*v(28)) v(237)=v(1138)*v(193)-v(235) v(236)=v(1138)*v(192)-v(235) v(234)=v(235)-v(233)*v(28) v(230)=v(13)*v(28) v(232)=v(1139)*v(193)+v(230) v(231)=v(1139)*v(192)+v(230) v(229)=v(230)+v(228)*v(28) v(23)=v(183)+v(29) v(223)=v(23)+v(5) v(1140)=v(13)*v(223) v(218)=v(23)-v(5) v(1141)=v(17)*v(218) v(22)=v(192)*v(20) v(225)=v(13)*v(22) v(227)=v(1140)*v(194)+v(225) v(226)=-(v(1140)*v(192))-v(225) v(224)=v(22)*v(223)+v(225) v(220)=-(v(17)*v(22)) v(222)=v(1141)*v(194)-v(220) v(221)=-(v(1141)*v(192))+v(220) v(219)=-(v(218)*v(22))+v(220) v(18)=v(184)+v(29) v(213)=v(18)-v(5) v(1143)=v(17)*v(213) v(208)=v(18)+v(5) v(1144)=v(13)*v(208) v(14)=1.-v(7) v(1142)=v(14)/4. v(197)=v(14)/8. v(196)=v(1142)*v(182) v(242)=v(196)*v(25) v(238)=v(196)*v(20) v(195)=v(1142)*v(183) v(244)=v(17)*v(195) v(240)=v(13)*v(195) v(35)=v(1142)*v(46) v(33)=v(1142)*v(44) v(16)=v(197)*v(25) v(215)=-(v(16)*v(17)) v(217)=-(v(1143)*v(193))+v(215) v(216)=v(1143)*v(197)-v(215) v(214)=-(v(16)*v(213))+v(215) v(210)=v(13)*v(16) v(212)=-(v(1144)*v(193))-v(210) v(211)=v(1144)*v(197)+v(210) v(209)=v(16)*v(208)+v(210) v(11)=-2.-v(6)-v(7) v(203)=v(11)+v(5) v(1145)=-(v(13)*v(203)) v(198)=v(11)-v(5) v(1146)=-(v(17)*v(198)) v(10)=v(197)*v(20) v(205)=v(10)*v(13) v(207)=v(1145)*v(194)-v(205) v(206)=v(1145)*v(197)-v(205) v(204)=v(10)*v(203)+v(205) v(200)=-(v(10)*v(17)) v(202)=v(1146)*v(194)+v(200) v(201)=v(1146)*v(197)+v(200) v(199)=-(v(10)*v(198))+v(200) v(8)=gp(4,i4) v(256)=v(1135)*v(199)+v(1132)*v(204)+v(1129)*v(209)+v(1126)*v & (214)+v(1123)*v(219)+v(1120)*v(224)+v(1117)*v(229)+v(1114)*v & (234)+v(1111)*v(238)+v(1105)*v(242)+v(1087)*v(252)+v(1081)*v & (254)+(-v(1102)+v(1108))*v(35)+(v(1096)-v(1099))*v(38)+(-v & (1090)+v(1093))*v(42)+(-v(1078)+v(1084))*v(50) v(257)=v(1134)*v(199)+v(1131)*v(204)+v(1128)*v(209)+v(1125)*v & (214)+v(1122)*v(219)+v(1119)*v(224)+v(1116)*v(229)+v(1113)*v & (234)+v(1110)*v(238)+v(1104)*v(242)+v(1086)*v(252)+v(1080)*v & (254)+(-v(1101)+v(1107))*v(35)+(v(1095)-v(1098))*v(38)+(-v & (1089)+v(1092))*v(42)+(-v(1077)+v(1083))*v(50) v(258)=v(1133)*v(199)+v(1130)*v(204)+v(1127)*v(209)+v(1124)*v & (214)+v(1121)*v(219)+v(1118)*v(224)+v(1115)*v(229)+v(1112)*v & (234)+v(1109)*v(238)+v(1103)*v(242)+v(1085)*v(252)+v(1079)*v & (254)+(-v(1100)+v(1106))*v(35)+(v(1094)-v(1097))*v(38)+(-v & (1088)+v(1091))*v(42)+(-v(1076)+v(1082))*v(50) v(259)=v(1135)*v(201)+v(1132)*v(206)+v(1129)*v(211)+v(1126)*v & (216)+v(1123)*v(221)+v(1120)*v(226)+v(1117)*v(231)+v(1114)*v & (236)+v(1108)*v(240)+v(1102)*v(244)+(-v(1090)+v(1099))*v(246)+ & (-v(1093)+v(1096))*v(248)+v(1084)*v(253)+v(1078)*v(255)+(v & (1105)-v(1111))*v(33)+(v(1081)-v(1087))*v(48) v(260)=v(1134)*v(201)+v(1131)*v(206)+v(1128)*v(211)+v(1125)*v & (216)+v(1122)*v(221)+v(1119)*v(226)+v(1116)*v(231)+v(1113)*v & (236)+v(1107)*v(240)+v(1101)*v(244)+(-v(1089)+v(1098))*v(246)+ & (-v(1092)+v(1095))*v(248)+v(1083)*v(253)+v(1077)*v(255)+(v & (1104)-v(1110))*v(33)+(v(1080)-v(1086))*v(48) v(1149)=-(v(257)*v(259))+v(256)*v(260) v(261)=v(1133)*v(201)+v(1130)*v(206)+v(1127)*v(211)+v(1124)*v & (216)+v(1121)*v(221)+v(1118)*v(226)+v(1115)*v(231)+v(1112)*v & (236)+v(1106)*v(240)+v(1100)*v(244)+(-v(1088)+v(1097))*v(246)+ & (-v(1091)+v(1094))*v(248)+v(1082)*v(253)+v(1076)*v(255)+(v & (1103)-v(1109))*v(33)+(v(1079)-v(1085))*v(48) v(1148)=v(258)*v(259)-v(256)*v(261) v(1147)=-(v(258)*v(260))+v(257)*v(261) v(262)=v(1135)*v(202)+v(1132)*v(207)+v(1129)*v(212)+v(1126)*v & (217)+v(1123)*v(222)+v(1120)*v(227)+v(1117)*v(232)+v(1114)*v & (237)+(-v(1087)+v(1111))*v(239)+(-v(1084)+v(1108))*v(241)+(-v & (1081)+v(1105))*v(243)+(-v(1078)+v(1102))*v(245)+v(1099)*v(247 & )+v(1096)*v(249)+v(1093)*v(250)+v(1090)*v(251) v(263)=v(1134)*v(202)+v(1131)*v(207)+v(1128)*v(212)+v(1125)*v & (217)+v(1122)*v(222)+v(1119)*v(227)+v(1116)*v(232)+v(1113)*v & (237)+(-v(1086)+v(1110))*v(239)+(-v(1083)+v(1107))*v(241)+(-v & (1080)+v(1104))*v(243)+(-v(1077)+v(1101))*v(245)+v(1098)*v(247 & )+v(1095)*v(249)+v(1092)*v(250)+v(1089)*v(251) v(264)=v(1133)*v(202)+v(1130)*v(207)+v(1127)*v(212)+v(1124)*v & (217)+v(1121)*v(222)+v(1118)*v(227)+v(1115)*v(232)+v(1112)*v & (237)+(-v(1085)+v(1109))*v(239)+(-v(1082)+v(1106))*v(241)+(-v & (1079)+v(1103))*v(243)+(-v(1076)+v(1100))*v(245)+v(1097)*v(247 & )+v(1094)*v(249)+v(1091)*v(250)+v(1088)*v(251) v(265)=v(1147)*v(262)+v(1148)*v(263)+v(1149)*v(264) v(266)=(-(v(261)*v(263))+v(260)*v(264))/v(265) v(1155)=v(16)*v(266) v(1152)=-(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(1156)=v(16)*v(267) v(1153)=-(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(1157)=v(16)*v(268) v(1154)=-(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(1147)/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(1148)/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(1149)/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(1150)*v(287)+v(293) v(291)=v(1150)*v(288)+v(295) v(292)=v(1150)*v(289)+v(297) v(294)=v(1151)*v(287)-v(293) v(296)=v(1151)*v(288)-v(295) v(298)=v(1151)*v(289)-v(297) 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(1137)*v(284)+v(320) v(303)=v(1137)*v(285)+v(322) v(304)=v(1137)*v(286)+v(324) v(305)=v(1137)*v(281)+v(326) v(306)=v(1137)*v(282)+v(328) v(307)=v(1137)*v(283)+v(330) 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(1142)*v(284)-v(320) v(323)=v(1142)*v(285)-v(322) v(325)=v(1142)*v(286)-v(324) v(327)=v(1142)*v(281)-v(326) v(329)=v(1142)*v(282)-v(328) v(331)=v(1142)*v(283)-v(330) 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(198)*(v(1152)+v(17)*v(340))-v(200)*v(371) v(348)=v(198)*(v(1153)+v(17)*v(343))-v(200)*v(373) v(349)=v(198)*(v(1154)+v(17)*v(346))-v(200)*v(375) v(350)=v(203)*(-v(1152)+v(13)*v(340))+v(205)*v(377) v(351)=v(203)*(-v(1153)+v(13)*v(343))+v(205)*v(379) v(352)=v(203)*(-v(1154)+v(13)*v(346))+v(205)*v(381) v(353)=v(208)*(v(1155)+v(13)*v(333))+v(210)*v(359) v(354)=v(208)*(v(1156)+v(13)*v(335))+v(210)*v(361) v(355)=v(208)*(v(1157)+v(13)*v(337))+v(210)*v(363) v(356)=v(213)*(-v(1155)+v(17)*v(333))-v(215)*v(365) v(357)=v(213)*(-v(1156)+v(17)*v(335))-v(215)*v(367) v(358)=v(213)*(-v(1157)+v(17)*v(337))-v(215)*v(369) v(360)=v(1141)*v(315)-v(22)*(v(218)*v(266)+v(17)*v(359)) v(362)=v(1141)*v(317)-v(22)*(v(218)*v(267)+v(17)*v(361)) v(364)=v(1141)*v(319)-v(22)*(v(218)*v(268)+v(17)*v(363)) 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(1138)*v(308)-v(28)*(v(233)*v(266)+v(17)*v(377)) v(380)=v(1138)*v(309)-v(28)*(v(233)*v(267)+v(17)*v(379)) v(382)=v(1138)*v(310)-v(28)*(v(233)*v(268)+v(17)*v(381)) 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(1073)*v(347)+v(1070)*v(350)+v(1067)*v(353)+v(1064 & )*v(356)+v(1061)*v(360)+v(1058)*v(366)+v(1055)*v(372)+v(1052 & )*v(378)+v(1049)*v(383)+v(1046)*v(386)+v(1043)*v(390)+v(1040 & )*v(396)+v(1037)*v(401)+v(1034)*v(405)+v(1031)*v(410)+v(1028 & )*v(414)+v(1025)*v(419)+v(1022)*v(422)+v(1019)*v(426)+v(1016 & )*v(432) v(453)=v(1073)*v(348)+v(1070)*v(351)+v(1067)*v(354)+v(1064)*v & (357)+v(1061)*v(362)+v(1058)*v(368)+v(1055)*v(374)+v(1052)*v & (380)+v(1049)*v(384)+v(1046)*v(387)+v(1043)*v(392)+v(1040)*v & (398)+v(1037)*v(402)+v(1034)*v(407)+v(1031)*v(411)+v(1028)*v & (416)+v(1025)*v(420)+v(1022)*v(423)+v(1019)*v(428)+v(1016)*v & (434) v(454)=v(1073)*v(349)+v(1070)*v(352)+v(1067)*v(355)+v(1064)*v & (358)+v(1061)*v(364)+v(1058)*v(370)+v(1055)*v(376)+v(1052)*v & (382)+v(1049)*v(385)+v(1046)*v(388)+v(1043)*v(394)+v(1040)*v & (400)+v(1037)*v(403)+v(1034)*v(409)+v(1031)*v(412)+v(1028)*v & (418)+v(1025)*v(421)+v(1022)*v(424)+v(1019)*v(430)+v(1016)*v & (436) v(455)=v(1072)*v(347)+v(1069)*v(350)+v(1066)*v(353)+v(1063)*v & (356)+v(1060)*v(360)+v(1057)*v(366)+v(1054)*v(372)+v(1051)*v & (378)+v(1048)*v(383)+v(1045)*v(386)+v(1042)*v(390)+v(1039)*v & (396)+v(1036)*v(401)+v(1033)*v(405)+v(1030)*v(410)+v(1027)*v & (414)+v(1024)*v(419)+v(1021)*v(422)+v(1018)*v(426)+v(1015)*v & (432) v(456)=1.+v(1072)*v(348)+v(1069)*v(351)+v(1066)*v(354)+v(1063 & )*v(357)+v(1060)*v(362)+v(1057)*v(368)+v(1054)*v(374)+v(1051 & )*v(380)+v(1048)*v(384)+v(1045)*v(387)+v(1042)*v(392)+v(1039 & )*v(398)+v(1036)*v(402)+v(1033)*v(407)+v(1030)*v(411)+v(1027 & )*v(416)+v(1024)*v(420)+v(1021)*v(423)+v(1018)*v(428)+v(1015 & )*v(434) v(508)=-(v(453)*v(455))+v(452)*v(456) v(457)=v(1072)*v(349)+v(1069)*v(352)+v(1066)*v(355)+v(1063)*v & (358)+v(1060)*v(364)+v(1057)*v(370)+v(1054)*v(376)+v(1051)*v & (382)+v(1048)*v(385)+v(1045)*v(388)+v(1042)*v(394)+v(1039)*v & (400)+v(1036)*v(403)+v(1033)*v(409)+v(1030)*v(412)+v(1027)*v & (418)+v(1024)*v(421)+v(1021)*v(424)+v(1018)*v(430)+v(1015)*v & (436) v(507)=v(454)*v(455)-v(452)*v(457) v(506)=-(v(454)*v(456))+v(453)*v(457) v(458)=v(1071)*v(347)+v(1068)*v(350)+v(1065)*v(353)+v(1062)*v & (356)+v(1059)*v(360)+v(1056)*v(366)+v(1053)*v(372)+v(1050)*v & (378)+v(1047)*v(383)+v(1044)*v(386)+v(1041)*v(390)+v(1038)*v & (396)+v(1035)*v(401)+v(1032)*v(405)+v(1029)*v(410)+v(1026)*v & (414)+v(1023)*v(419)+v(1020)*v(422)+v(1017)*v(426)+v(1014)*v & (432) v(459)=v(1071)*v(348)+v(1068)*v(351)+v(1065)*v(354)+v(1062)*v & (357)+v(1059)*v(362)+v(1056)*v(368)+v(1053)*v(374)+v(1050)*v & (380)+v(1047)*v(384)+v(1044)*v(387)+v(1041)*v(392)+v(1038)*v & (398)+v(1035)*v(402)+v(1032)*v(407)+v(1029)*v(411)+v(1026)*v & (416)+v(1023)*v(420)+v(1020)*v(423)+v(1017)*v(428)+v(1014)*v & (434) v(505)=v(453)*v(458)-v(452)*v(459) v(502)=-(v(456)*v(458))+v(455)*v(459) v(460)=1.+v(1071)*v(349)+v(1068)*v(352)+v(1065)*v(355)+v(1062 & )*v(358)+v(1059)*v(364)+v(1056)*v(370)+v(1053)*v(376)+v(1050 & )*v(382)+v(1047)*v(385)+v(1044)*v(388)+v(1041)*v(394)+v(1038 & )*v(400)+v(1035)*v(403)+v(1032)*v(409)+v(1029)*v(412)+v(1026 & )*v(418)+v(1023)*v(421)+v(1020)*v(424)+v(1017)*v(430)+v(1014 & )*v(436) 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(1167)=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(1167)+v(449))*(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(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 !******************* S U B R O U T I N E ********************** SUBROUTINE SPP10(v,d,ul,ul0,xl,s,p,ht,hp,sg,sg0,sxd,gpost,npost &) IMPLICIT NONE include 'sms.h' INTEGER i522,i523 DOUBLE PRECISION v(17608),d(2),ul(3,20),ul0(3,20),xl(3 &,20),s(60,60),p(60),ht(*),hp(*),sg(*),sg0(*),sxd(60),gpost(64 &,19),npost(20,6) v(1289)=xl(1,1) v(1288)=xl(2,1) v(1287)=xl(3,1) v(1286)=xl(1,2) v(1285)=xl(2,2) v(1284)=xl(3,2) v(1283)=xl(1,3) v(1282)=xl(2,3) v(1281)=xl(3,3) v(1280)=xl(1,4) v(1279)=xl(2,4) v(1278)=xl(3,4) v(1277)=xl(1,5) v(1276)=xl(2,5) v(1275)=xl(3,5) v(1274)=xl(1,6) v(1273)=xl(2,6) v(1272)=xl(3,6) v(1271)=xl(1,7) v(1270)=xl(2,7) v(1269)=xl(3,7) v(1268)=xl(1,8) v(1267)=xl(2,8) v(1266)=xl(3,8) v(1265)=xl(1,9) v(1264)=xl(2,9) v(1263)=xl(3,9) v(1262)=xl(1,10) v(1261)=xl(2,10) v(1260)=xl(3,10) v(1259)=xl(1,11) v(1258)=xl(2,11) v(1257)=xl(3,11) v(1256)=xl(1,12) v(1255)=xl(2,12) v(1254)=xl(3,12) v(1253)=xl(1,13) v(1252)=xl(2,13) v(1251)=xl(3,13) v(1250)=xl(1,14) v(1249)=xl(2,14) v(1248)=xl(3,14) v(1247)=xl(1,15) v(1246)=xl(2,15) v(1245)=xl(3,15) v(1244)=xl(1,16) v(1243)=xl(2,16) v(1242)=xl(3,16) v(1241)=xl(1,17) v(1240)=xl(2,17) v(1239)=xl(3,17) v(1238)=xl(1,18) v(1237)=xl(2,18) v(1236)=xl(3,18) v(1235)=xl(1,19) v(1234)=xl(2,19) v(1233)=xl(3,19) v(1232)=xl(1,20) v(1231)=xl(2,20) v(1230)=xl(3,20) v(1229)=2.*d(2) v(1228)=1d0/(1.+d(2)) v(1227)=ul(1,1) v(1226)=ul(2,1) v(1225)=ul(3,1) v(1224)=ul(1,2) v(1223)=ul(2,2) v(1222)=ul(3,2) v(1221)=ul(1,3) v(1220)=ul(2,3) v(1219)=ul(3,3) v(1218)=ul(1,4) v(1217)=ul(2,4) v(1216)=ul(3,4) v(1215)=ul(1,5) v(1214)=ul(2,5) v(1213)=ul(3,5) v(1212)=ul(1,6) v(1211)=ul(2,6) v(1210)=ul(3,6) v(1209)=ul(1,7) v(1208)=ul(2,7) v(1207)=ul(3,7) v(1206)=ul(1,8) v(1205)=ul(2,8) v(1204)=ul(3,8) v(1203)=ul(1,9) v(1202)=ul(2,9) v(1201)=ul(3,9) v(1200)=ul(1,10) v(1199)=ul(2,10) v(1198)=ul(3,10) v(1197)=ul(1,11) v(1196)=ul(2,11) v(1195)=ul(3,11) v(1194)=ul(1,12) v(1193)=ul(2,12) v(1192)=ul(3,12) v(1191)=ul(1,13) v(1190)=ul(2,13) v(1189)=ul(3,13) v(1188)=ul(1,14) v(1187)=ul(2,14) v(1186)=ul(3,14) v(1185)=ul(1,15) v(1184)=ul(2,15) v(1183)=ul(3,15) v(1182)=ul(1,16) v(1181)=ul(2,16) v(1180)=ul(3,16) v(1179)=ul(1,17) v(1178)=ul(2,17) v(1177)=ul(3,17) v(1176)=ul(1,18) v(1175)=ul(2,18) v(1174)=ul(3,18) v(1173)=ul(1,19) v(1172)=ul(2,19) v(1171)=ul(3,19) v(1170)=ul(1,20) v(1169)=ul(2,20) v(1168)=ul(3,20) v(970)=(d(1)*v(1228))/2. v(968)=(v(1229)*v(970))/(1.-v(1229)) i522=int(ngpo) DO i523=1,i522 v(524)=gp(1,i523) v(701)=-2.*v(524) v(563)=1.-(v(524)*v(524)) v(710)=v(563)/4. v(536)=1.-v(524) v(532)=1.+v(524) v(525)=gp(2,i523) v(702)=-2.*v(525) v(565)=1.-(v(525)*v(525)) v(708)=v(565)/4. v(764)=-(v(536)*v(708)) v(760)=-(v(532)*v(708)) v(544)=1.+v(525) v(1304)=v(544)/4. v(762)=-(v(544)*v(710)) v(712)=v(544)/8. v(539)=1.-v(525) v(1305)=v(539)/4. v(758)=-(v(539)*v(710)) v(713)=v(539)/8. v(526)=gp(3,i523) v(703)=-2.*v(526) v(1290)=v(703)/4. v(706)=v(1290)*v(539) v(768)=v(532)*v(706) v(766)=v(536)*v(706) v(705)=v(1290)*v(544) v(770)=v(536)*v(705) v(769)=v(532)*v(705) v(559)=1.-(v(526)*v(526)) v(704)=v(559)/4. v(767)=-(v(532)*v(704)) v(765)=-(v(536)*v(704)) v(561)=v(544)*v(704) v(557)=v(539)*v(704) v(548)=-2.+v(525)+v(526) v(752)=-v(524)+v(548) v(1292)=v(536)*v(752) v(747)=v(524)+v(548) v(1293)=v(532)*v(747) v(545)=1.+v(526) v(1291)=v(545)/4. v(711)=v(545)/8. v(709)=v(1291)*v(701) v(773)=v(544)*v(709) v(771)=v(539)*v(709) v(707)=v(1291)*v(702) v(774)=v(536)*v(707) v(772)=v(532)*v(707) v(569)=v(1291)*v(565) v(567)=v(1291)*v(563) v(547)=v(544)*v(711) v(754)=-(v(536)*v(547)) v(756)=v(1292)*v(712)-v(754) v(755)=v(1292)*v(711)-v(754) v(753)=-(v(547)*v(752))+v(754) v(749)=v(532)*v(547) v(751)=v(1293)*v(712)+v(749) v(750)=v(1293)*v(711)+v(749) v(748)=v(547)*v(747)+v(749) v(542)=v(548)+v(702) v(742)=v(524)+v(542) v(1294)=v(532)*v(742) v(737)=-v(524)+v(542) v(1295)=v(536)*v(737) v(541)=v(539)*v(711) v(744)=v(532)*v(541) v(746)=v(1294)*v(713)+v(744) v(745)=-(v(1294)*v(711))-v(744) v(743)=v(541)*v(742)+v(744) v(739)=-(v(536)*v(541)) v(741)=v(1295)*v(713)-v(739) v(740)=-(v(1295)*v(711))+v(739) v(738)=-(v(541)*v(737))+v(739) v(537)=v(548)+v(703) v(732)=-v(524)+v(537) v(1297)=v(536)*v(732) v(727)=v(524)+v(537) v(1298)=v(532)*v(727) v(533)=1.-v(526) v(1296)=v(533)/4. v(716)=v(533)/8. v(715)=v(1296)*v(701) v(761)=v(544)*v(715) v(757)=v(539)*v(715) v(714)=v(1296)*v(702) v(763)=v(536)*v(714) v(759)=v(532)*v(714) v(554)=v(1296)*v(565) v(552)=v(1296)*v(563) v(535)=v(544)*v(716) v(734)=-(v(535)*v(536)) v(736)=-(v(1297)*v(712))+v(734) v(735)=v(1297)*v(716)-v(734) v(733)=-(v(535)*v(732))+v(734) v(729)=v(532)*v(535) v(731)=-(v(1298)*v(712))-v(729) v(730)=v(1298)*v(716)+v(729) v(728)=v(535)*v(727)+v(729) v(530)=-2.-v(525)-v(526) v(722)=v(524)+v(530) v(1299)=-(v(532)*v(722)) v(717)=-v(524)+v(530) v(1300)=-(v(536)*v(717)) v(529)=v(539)*v(716) v(724)=v(529)*v(532) v(726)=v(1299)*v(713)-v(724) v(725)=v(1299)*v(716)-v(724) v(723)=v(529)*v(722)+v(724) v(719)=-(v(529)*v(536)) v(721)=v(1300)*v(713)+v(719) v(720)=v(1300)*v(716)+v(719) v(718)=-(v(529)*v(717))+v(719) v(775)=(-v(1256)+v(1262))*v(554)+(v(1250)-v(1253))*v(557)+(-v & (1244)+v(1247))*v(561)+(-v(1232)+v(1238))*v(569)+v(1289)*v(718 & )+v(1286)*v(723)+v(1283)*v(728)+v(1280)*v(733)+v(1277)*v(738) & +v(1274)*v(743)+v(1271)*v(748)+v(1268)*v(753)+v(1265)*v(757)+v & (1259)*v(761)+v(1241)*v(771)+v(1235)*v(773) v(776)=(-v(1255)+v(1261))*v(554)+(v(1249)-v(1252))*v(557)+(-v & (1243)+v(1246))*v(561)+(-v(1231)+v(1237))*v(569)+v(1288)*v(718 & )+v(1285)*v(723)+v(1282)*v(728)+v(1279)*v(733)+v(1276)*v(738) & +v(1273)*v(743)+v(1270)*v(748)+v(1267)*v(753)+v(1264)*v(757)+v & (1258)*v(761)+v(1240)*v(771)+v(1234)*v(773) v(777)=(-v(1254)+v(1260))*v(554)+(v(1248)-v(1251))*v(557)+(-v & (1242)+v(1245))*v(561)+(-v(1230)+v(1236))*v(569)+v(1287)*v(718 & )+v(1284)*v(723)+v(1281)*v(728)+v(1278)*v(733)+v(1275)*v(738) & +v(1272)*v(743)+v(1269)*v(748)+v(1266)*v(753)+v(1263)*v(757)+v & (1257)*v(761)+v(1239)*v(771)+v(1233)*v(773) v(778)=(v(1259)-v(1265))*v(552)+(v(1235)-v(1241))*v(567)+v & (1289)*v(720)+v(1286)*v(725)+v(1283)*v(730)+v(1280)*v(735)+v & (1277)*v(740)+v(1274)*v(745)+v(1271)*v(750)+v(1268)*v(755)+v & (1262)*v(759)+v(1256)*v(763)+(-v(1244)+v(1253))*v(765)+(-v & (1247)+v(1250))*v(767)+v(1238)*v(772)+v(1232)*v(774) v(779)=(v(1258)-v(1264))*v(552)+(v(1234)-v(1240))*v(567)+v & (1288)*v(720)+v(1285)*v(725)+v(1282)*v(730)+v(1279)*v(735)+v & (1276)*v(740)+v(1273)*v(745)+v(1270)*v(750)+v(1267)*v(755)+v & (1261)*v(759)+v(1255)*v(763)+(-v(1243)+v(1252))*v(765)+(-v & (1246)+v(1249))*v(767)+v(1237)*v(772)+v(1231)*v(774) v(1303)=-(v(776)*v(778))+v(775)*v(779) v(780)=(v(1257)-v(1263))*v(552)+(v(1233)-v(1239))*v(567)+v & (1287)*v(720)+v(1284)*v(725)+v(1281)*v(730)+v(1278)*v(735)+v & (1275)*v(740)+v(1272)*v(745)+v(1269)*v(750)+v(1266)*v(755)+v & (1260)*v(759)+v(1254)*v(763)+(-v(1242)+v(1251))*v(765)+(-v & (1245)+v(1248))*v(767)+v(1236)*v(772)+v(1230)*v(774) v(1302)=v(777)*v(778)-v(775)*v(780) v(1301)=-(v(777)*v(779))+v(776)*v(780) v(781)=v(1289)*v(721)+v(1286)*v(726)+v(1283)*v(731)+v(1280)*v & (736)+v(1277)*v(741)+v(1274)*v(746)+v(1271)*v(751)+v(1268)*v & (756)+(-v(1241)+v(1265))*v(758)+(-v(1238)+v(1262))*v(760)+(-v & (1235)+v(1259))*v(762)+(-v(1232)+v(1256))*v(764)+v(1253)*v(766 & )+v(1250)*v(768)+v(1247)*v(769)+v(1244)*v(770) v(782)=v(1288)*v(721)+v(1285)*v(726)+v(1282)*v(731)+v(1279)*v & (736)+v(1276)*v(741)+v(1273)*v(746)+v(1270)*v(751)+v(1267)*v & (756)+(-v(1240)+v(1264))*v(758)+(-v(1237)+v(1261))*v(760)+(-v & (1234)+v(1258))*v(762)+(-v(1231)+v(1255))*v(764)+v(1252)*v(766 & )+v(1249)*v(768)+v(1246)*v(769)+v(1243)*v(770) v(783)=v(1287)*v(721)+v(1284)*v(726)+v(1281)*v(731)+v(1278)*v & (736)+v(1275)*v(741)+v(1272)*v(746)+v(1269)*v(751)+v(1266)*v & (756)+(-v(1239)+v(1263))*v(758)+(-v(1236)+v(1260))*v(760)+(-v & (1233)+v(1257))*v(762)+(-v(1230)+v(1254))*v(764)+v(1251)*v(766 & )+v(1248)*v(768)+v(1245)*v(769)+v(1242)*v(770) v(784)=v(1301)*v(781)+v(1302)*v(782)+v(1303)*v(783) v(785)=(-(v(780)*v(782))+v(779)*v(783))/v(784) v(1309)=v(535)*v(785) v(1306)=-(v(529)*v(785)) v(950)=v(569)*v(785) v(932)=v(561)*v(785) v(923)=-(v(557)*v(785)) v(914)=v(554)*v(785) v(786)=(v(780)*v(781)-v(778)*v(783))/v(784) v(1310)=v(535)*v(786) v(1307)=-(v(529)*v(786)) v(952)=v(569)*v(786) v(934)=v(561)*v(786) v(925)=-(v(557)*v(786)) v(916)=v(554)*v(786) v(787)=(-(v(779)*v(781))+v(778)*v(782))/v(784) v(1311)=v(535)*v(787) v(1308)=-(v(529)*v(787)) v(954)=v(569)*v(787) v(936)=v(561)*v(787) v(927)=-(v(557)*v(787)) v(918)=v(554)*v(787) v(788)=(v(777)*v(782)-v(776)*v(783))/v(784) v(944)=-(v(567)*v(788)) v(908)=-(v(552)*v(788)) v(857)=v(716)*v(788) v(833)=v(711)*v(788) v(812)=v(704)*v(788) v(789)=(-(v(777)*v(781))+v(775)*v(783))/v(784) v(946)=-(v(567)*v(789)) v(910)=-(v(552)*v(789)) v(860)=v(716)*v(789) v(835)=v(711)*v(789) v(814)=v(704)*v(789) v(790)=(v(776)*v(781)-v(775)*v(782))/v(784) v(948)=-(v(567)*v(790)) v(912)=-(v(552)*v(790)) v(863)=v(716)*v(790) v(837)=v(711)*v(790) v(816)=v(704)*v(790) v(791)=v(1301)/v(784) v(858)=v(713)*v(791) v(851)=v(712)*v(791) v(845)=v(710)*v(791) v(839)=v(708)*v(791) v(792)=v(1302)/v(784) v(861)=v(713)*v(792) v(853)=v(712)*v(792) v(847)=v(710)*v(792) v(841)=v(708)*v(792) v(793)=v(1303)/v(784) v(864)=v(713)*v(793) v(855)=v(712)*v(793) v(849)=v(710)*v(793) v(843)=v(708)*v(793) v(800)=v(701)*v(785) v(801)=v(701)*v(786) v(802)=v(701)*v(787) v(803)=v(702)*v(788) v(804)=v(702)*v(789) v(805)=v(702)*v(790) v(806)=v(703)*v(791) v(807)=v(703)*v(792) v(808)=v(703)*v(793) v(809)=v(1304)*v(806)+v(812) v(810)=v(1304)*v(807)+v(814) v(811)=v(1304)*v(808)+v(816) v(813)=v(1305)*v(806)-v(812) v(815)=v(1305)*v(807)-v(814) v(817)=v(1305)*v(808)-v(816) v(818)=v(788)+v(791) v(896)=v(785)-v(818) v(890)=-v(785)-v(818) v(819)=v(789)+v(792) v(898)=v(786)-v(819) v(892)=-v(786)-v(819) v(820)=v(790)+v(793) v(900)=v(787)-v(820) v(894)=-v(787)-v(820) v(821)=v(1291)*v(803)+v(839) v(822)=v(1291)*v(804)+v(841) v(823)=v(1291)*v(805)+v(843) v(824)=v(1291)*v(800)+v(845) v(825)=v(1291)*v(801)+v(847) v(826)=v(1291)*v(802)+v(849) v(827)=v(833)+v(851) v(828)=v(835)+v(853) v(829)=v(837)+v(855) v(830)=-v(788)+v(791) v(884)=-v(785)-v(830) v(878)=v(785)-v(830) v(831)=-v(789)+v(792) v(886)=-v(786)-v(831) v(880)=v(786)-v(831) v(832)=-v(790)+v(793) v(888)=-v(787)-v(832) v(882)=v(787)-v(832) v(834)=-v(833)+v(858) v(836)=-v(835)+v(861) v(838)=-v(837)+v(864) v(840)=v(1296)*v(803)-v(839) v(842)=v(1296)*v(804)-v(841) v(844)=v(1296)*v(805)-v(843) v(846)=v(1296)*v(800)-v(845) v(848)=v(1296)*v(801)-v(847) v(850)=v(1296)*v(802)-v(849) v(852)=-v(851)+v(857) v(854)=-v(853)+v(860) v(856)=-v(855)+v(863) v(859)=-v(857)-v(858) v(862)=-v(860)-v(861) v(865)=-v(863)-v(864) v(866)=v(717)*(v(1306)+v(536)*v(859))-v(719)*v(890) v(867)=v(717)*(v(1307)+v(536)*v(862))-v(719)*v(892) v(868)=v(717)*(v(1308)+v(536)*v(865))-v(719)*v(894) v(869)=v(722)*(-v(1306)+v(532)*v(859))+v(724)*v(896) v(870)=v(722)*(-v(1307)+v(532)*v(862))+v(724)*v(898) v(871)=v(722)*(-v(1308)+v(532)*v(865))+v(724)*v(900) v(872)=v(727)*(v(1309)+v(532)*v(852))+v(729)*v(878) v(873)=v(727)*(v(1310)+v(532)*v(854))+v(729)*v(880) v(874)=v(727)*(v(1311)+v(532)*v(856))+v(729)*v(882) v(875)=v(732)*(-v(1309)+v(536)*v(852))-v(734)*v(884) v(876)=v(732)*(-v(1310)+v(536)*v(854))-v(734)*v(886) v(877)=v(732)*(-v(1311)+v(536)*v(856))-v(734)*v(888) v(879)=v(1295)*v(834)-v(541)*(v(737)*v(785)+v(536)*v(878)) v(881)=v(1295)*v(836)-v(541)*(v(737)*v(786)+v(536)*v(880)) v(883)=v(1295)*v(838)-v(541)*(v(737)*v(787)+v(536)*v(882)) v(885)=v(742)*(v(541)*v(785)+v(532)*v(834))-v(744)*v(884) v(887)=v(742)*(v(541)*v(786)+v(532)*v(836))-v(744)*v(886) v(889)=v(742)*(v(541)*v(787)+v(532)*v(838))-v(744)*v(888) v(891)=v(747)*(v(547)*v(785)+v(532)*v(827))-v(749)*v(890) v(893)=v(747)*(v(547)*v(786)+v(532)*v(828))-v(749)*v(892) v(895)=v(747)*(v(547)*v(787)+v(532)*v(829))-v(749)*v(894) v(897)=v(1292)*v(827)-v(547)*(v(752)*v(785)+v(536)*v(896)) v(899)=v(1292)*v(828)-v(547)*(v(752)*v(786)+v(536)*v(898)) v(901)=v(1292)*v(829)-v(547)*(v(752)*v(787)+v(536)*v(900)) v(902)=v(539)*v(846)+v(908) v(903)=v(539)*v(848)+v(910) v(904)=v(539)*v(850)+v(912) v(905)=v(532)*v(840)+v(914) v(906)=v(532)*v(842)+v(916) v(907)=v(532)*v(844)+v(918) v(909)=v(544)*v(846)-v(908) v(911)=v(544)*v(848)-v(910) v(913)=v(544)*v(850)-v(912) v(915)=v(536)*v(840)-v(914) v(917)=v(536)*v(842)-v(916) v(919)=v(536)*v(844)-v(918) v(920)=v(536)*v(813)+v(923) v(921)=v(536)*v(815)+v(925) v(922)=v(536)*v(817)+v(927) v(924)=v(532)*v(813)-v(923) v(926)=v(532)*v(815)-v(925) v(928)=v(532)*v(817)-v(927) v(929)=v(532)*v(809)+v(932) v(930)=v(532)*v(810)+v(934) v(931)=v(532)*v(811)+v(936) v(933)=v(536)*v(809)-v(932) v(935)=v(536)*v(810)-v(934) v(937)=v(536)*v(811)-v(936) v(938)=v(539)*v(824)+v(944) v(939)=v(539)*v(825)+v(946) v(940)=v(539)*v(826)+v(948) v(941)=v(532)*v(821)+v(950) v(942)=v(532)*v(822)+v(952) v(943)=v(532)*v(823)+v(954) v(945)=v(544)*v(824)-v(944) v(947)=v(544)*v(825)-v(946) v(949)=v(544)*v(826)-v(948) v(951)=v(536)*v(821)-v(950) v(953)=v(536)*v(822)-v(952) v(955)=v(536)*v(823)-v(954) v(971)=1.+v(1227)*v(866)+v(1224)*v(869)+v(1221)*v(872)+v(1218 & )*v(875)+v(1215)*v(879)+v(1212)*v(885)+v(1209)*v(891)+v(1206 & )*v(897)+v(1203)*v(902)+v(1200)*v(905)+v(1197)*v(909)+v(1194 & )*v(915)+v(1191)*v(920)+v(1188)*v(924)+v(1185)*v(929)+v(1182 & )*v(933)+v(1179)*v(938)+v(1176)*v(941)+v(1173)*v(945)+v(1170 & )*v(951) v(972)=v(1227)*v(867)+v(1224)*v(870)+v(1221)*v(873)+v(1218)*v & (876)+v(1215)*v(881)+v(1212)*v(887)+v(1209)*v(893)+v(1206)*v & (899)+v(1203)*v(903)+v(1200)*v(906)+v(1197)*v(911)+v(1194)*v & (917)+v(1191)*v(921)+v(1188)*v(926)+v(1185)*v(930)+v(1182)*v & (935)+v(1179)*v(939)+v(1176)*v(942)+v(1173)*v(947)+v(1170)*v & (953) v(973)=v(1227)*v(868)+v(1224)*v(871)+v(1221)*v(874)+v(1218)*v & (877)+v(1215)*v(883)+v(1212)*v(889)+v(1209)*v(895)+v(1206)*v & (901)+v(1203)*v(904)+v(1200)*v(907)+v(1197)*v(913)+v(1194)*v & (919)+v(1191)*v(922)+v(1188)*v(928)+v(1185)*v(931)+v(1182)*v & (937)+v(1179)*v(940)+v(1176)*v(943)+v(1173)*v(949)+v(1170)*v & (955) v(974)=v(1226)*v(866)+v(1223)*v(869)+v(1220)*v(872)+v(1217)*v & (875)+v(1214)*v(879)+v(1211)*v(885)+v(1208)*v(891)+v(1205)*v & (897)+v(1202)*v(902)+v(1199)*v(905)+v(1196)*v(909)+v(1193)*v & (915)+v(1190)*v(920)+v(1187)*v(924)+v(1184)*v(929)+v(1181)*v & (933)+v(1178)*v(938)+v(1175)*v(941)+v(1172)*v(945)+v(1169)*v & (951) v(975)=1.+v(1226)*v(867)+v(1223)*v(870)+v(1220)*v(873)+v(1217 & )*v(876)+v(1214)*v(881)+v(1211)*v(887)+v(1208)*v(893)+v(1205 & )*v(899)+v(1202)*v(903)+v(1199)*v(906)+v(1196)*v(911)+v(1193 & )*v(917)+v(1190)*v(921)+v(1187)*v(926)+v(1184)*v(930)+v(1181 & )*v(935)+v(1178)*v(939)+v(1175)*v(942)+v(1172)*v(947)+v(1169 & )*v(953) v(976)=v(1226)*v(868)+v(1223)*v(871)+v(1220)*v(874)+v(1217)*v & (877)+v(1214)*v(883)+v(1211)*v(889)+v(1208)*v(895)+v(1205)*v & (901)+v(1202)*v(904)+v(1199)*v(907)+v(1196)*v(913)+v(1193)*v & (919)+v(1190)*v(922)+v(1187)*v(928)+v(1184)*v(931)+v(1181)*v & (937)+v(1178)*v(940)+v(1175)*v(943)+v(1172)*v(949)+v(1169)*v & (955) v(977)=v(1225)*v(866)+v(1222)*v(869)+v(1219)*v(872)+v(1216)*v & (875)+v(1213)*v(879)+v(1210)*v(885)+v(1207)*v(891)+v(1204)*v & (897)+v(1201)*v(902)+v(1198)*v(905)+v(1195)*v(909)+v(1192)*v & (915)+v(1189)*v(920)+v(1186)*v(924)+v(1183)*v(929)+v(1180)*v & (933)+v(1177)*v(938)+v(1174)*v(941)+v(1171)*v(945)+v(1168)*v & (951) v(1312)=(v(977)*v(977)) v(978)=v(1225)*v(867)+v(1222)*v(870)+v(1219)*v(873)+v(1216)*v & (876)+v(1213)*v(881)+v(1210)*v(887)+v(1207)*v(893)+v(1204)*v & (899)+v(1201)*v(903)+v(1198)*v(906)+v(1195)*v(911)+v(1192)*v & (917)+v(1189)*v(921)+v(1186)*v(926)+v(1183)*v(930)+v(1180)*v & (935)+v(1177)*v(939)+v(1174)*v(942)+v(1171)*v(947)+v(1168)*v & (953) v(1313)=(v(978)*v(978)) v(979)=1.+v(1225)*v(868)+v(1222)*v(871)+v(1219)*v(874)+v(1216 & )*v(877)+v(1213)*v(883)+v(1210)*v(889)+v(1207)*v(895)+v(1204 & )*v(901)+v(1201)*v(904)+v(1198)*v(907)+v(1195)*v(913)+v(1192 & )*v(919)+v(1189)*v(922)+v(1186)*v(928)+v(1183)*v(931)+v(1180 & )*v(937)+v(1177)*v(940)+v(1174)*v(943)+v(1171)*v(949)+v(1168 & )*v(955) v(1314)=(v(979)*v(979)) v(986)=(-(v(973)*v(975))+v(972)*v(976))*v(977)+(v(973)*v(974) & -v(971)*v(976))*v(978)+(-(v(972)*v(974))+v(971)*v(975))*v(979 & ) v(989)=(v(971)*v(972)+v(974)*v(975)+v(977)*v(978))/2. v(990)=(v(971)*v(973)+v(974)*v(976)+v(977)*v(979))/2. v(992)=(v(972)*v(973)+v(975)*v(976)+v(978)*v(979))/2. v(995)=v(968)*(-1.+v(986))-v(970)/v(986) v(996)=v(970)*v(971)+(-(v(976)*v(978))+v(975)*v(979))*v(995) v(997)=v(970)*v(972)+(v(976)*v(977)-v(974)*v(979))*v(995) v(998)=v(970)*v(973)+(-(v(975)*v(977))+v(974)*v(978))*v(995) v(999)=v(970)*v(974)+(v(973)*v(978)-v(972)*v(979))*v(995) v(1000)=v(970)*v(975)+(-(v(973)*v(977))+v(971)*v(979))*v(995) v(1001)=v(970)*v(976)+(v(972)*v(977)-v(971)*v(978))*v(995) v(1005)=(v(971)*v(996)+v(972)*v(997)+v(973)*v(998))/v(986) v(1006)=(v(974)*v(996)+v(975)*v(997)+v(976)*v(998))/v(986) v(1007)=(v(977)*v(996)+v(978)*v(997)+v(979)*v(998))/v(986) v(1008)=(v(1000)*v(975)+v(1001)*v(976)+v(974)*v(999))/v(986) v(1009)=(v(1000)*v(978)+v(1001)*v(979)+v(977)*v(999))/v(986) v(1010)=((v(1312)+v(1313)+v(1314))*v(970)+v(986)*v(995))/v(986 & ) v(1011)=(-v(1005)-v(1008)-v(1010))/3. gpost(i523,1)=v(1005) gpost(i523,2)=v(1006) gpost(i523,3)=v(1007) gpost(i523,4)=v(1006) gpost(i523,5)=v(1008) gpost(i523,6)=v(1009) gpost(i523,7)=v(1007) gpost(i523,8)=v(1009) gpost(i523,9)=v(1010) gpost(i523,10)=(-1.+v(1312)+(v(971)*v(971))+(v(974)*v(974))) & /2. gpost(i523,11)=v(989) gpost(i523,12)=v(990) gpost(i523,13)=v(989) gpost(i523,14)=(-1.+v(1313)+(v(972)*v(972))+(v(975)*v(975))) & /2. gpost(i523,15)=v(992) gpost(i523,16)=v(990) gpost(i523,17)=v(992) gpost(i523,18)=(-1.+v(1314)+(v(973)*v(973))+(v(976)*v(976))) & /2. gpost(i523,19)=sqrt(1.5*(2.*v(1006)**2+2.*v(1007)**2+2.*v(1009 & )**2+(v(1005)+v(1011))**2+(v(1008)+v(1011))**2+(v(1010)+v(1011 & ))**2)) ENDDO npost(1,1)=v(1227) npost(2,1)=v(1224) npost(3,1)=v(1221) npost(4,1)=v(1218) npost(5,1)=v(1215) npost(6,1)=v(1212) npost(7,1)=v(1209) npost(8,1)=v(1206) npost(9,1)=v(1203) npost(10,1)=v(1200) npost(11,1)=v(1197) npost(12,1)=v(1194) npost(13,1)=v(1191) npost(14,1)=v(1188) npost(15,1)=v(1185) npost(16,1)=v(1182) npost(17,1)=v(1179) npost(18,1)=v(1176) npost(19,1)=v(1173) npost(20,1)=v(1170) npost(1,2)=v(1226) npost(2,2)=v(1223) npost(3,2)=v(1220) npost(4,2)=v(1217) npost(5,2)=v(1214) npost(6,2)=v(1211) npost(7,2)=v(1208) npost(8,2)=v(1205) npost(9,2)=v(1202) npost(10,2)=v(1199) npost(11,2)=v(1196) npost(12,2)=v(1193) npost(13,2)=v(1190) npost(14,2)=v(1187) npost(15,2)=v(1184) npost(16,2)=v(1181) npost(17,2)=v(1178) npost(18,2)=v(1175) npost(19,2)=v(1172) npost(20,2)=v(1169) npost(1,3)=v(1225) npost(2,3)=v(1222) npost(3,3)=v(1219) npost(4,3)=v(1216) npost(5,3)=v(1213) npost(6,3)=v(1210) npost(7,3)=v(1207) npost(8,3)=v(1204) npost(9,3)=v(1201) npost(10,3)=v(1198) npost(11,3)=v(1195) npost(12,3)=v(1192) npost(13,3)=v(1189) npost(14,3)=v(1186) npost(15,3)=v(1183) npost(16,3)=v(1180) npost(17,3)=v(1177) npost(18,3)=v(1174) npost(19,3)=v(1171) npost(20,3)=v(1168) npost(1,4)=v(1227) npost(2,4)=v(1224) npost(3,4)=v(1221) npost(4,4)=v(1218) npost(5,4)=v(1215) npost(6,4)=v(1212) npost(7,4)=v(1209) npost(8,4)=v(1206) npost(9,4)=v(1203) npost(10,4)=v(1200) npost(11,4)=v(1197) npost(12,4)=v(1194) npost(13,4)=v(1191) npost(14,4)=v(1188) npost(15,4)=v(1185) npost(16,4)=v(1182) npost(17,4)=v(1179) npost(18,4)=v(1176) npost(19,4)=v(1173) npost(20,4)=v(1170) npost(1,5)=v(1226) npost(2,5)=v(1223) npost(3,5)=v(1220) npost(4,5)=v(1217) npost(5,5)=v(1214) npost(6,5)=v(1211) npost(7,5)=v(1208) npost(8,5)=v(1205) npost(9,5)=v(1202) npost(10,5)=v(1199) npost(11,5)=v(1196) npost(12,5)=v(1193) npost(13,5)=v(1190) npost(14,5)=v(1187) npost(15,5)=v(1184) npost(16,5)=v(1181) npost(17,5)=v(1178) npost(18,5)=v(1175) npost(19,5)=v(1172) npost(20,5)=v(1169) npost(1,6)=v(1225) npost(2,6)=v(1222) npost(3,6)=v(1219) npost(4,6)=v(1216) npost(5,6)=v(1213) npost(6,6)=v(1210) npost(7,6)=v(1207) npost(8,6)=v(1204) npost(9,6)=v(1201) npost(10,6)=v(1198) npost(11,6)=v(1195) npost(12,6)=v(1192) npost(13,6)=v(1189) npost(14,6)=v(1186) npost(15,6)=v(1183) npost(16,6)=v(1180) npost(17,6)=v(1177) npost(18,6)=v(1174) npost(19,6)=v(1171) npost(20,6)=v(1168) END