/*************************************************************
* AceGen    2.004 Windows (6 Nov 07)                         *
*           Co. J. Korelc  2007            6 Nov 07 22:09:26 *
**************************************************************
User : USER
Evaluation time                 : 28 s    Mode  : Optimal
Number of formulae              : 537     Method: Automatic
Subroutine                      : SKR size :6448
Subroutine                      : SPP size :4067
Total size of Mathematica  code : 10515 subexpressions
Total size of C code            : 29634 bytes*/
#include "sms.h"
void SKR(double v[5409],ElementSpec *es,ElementData *ed,NodeSpec **ns,NodeData **nd,double *rdata,int *idata,double *p,double **s);
void SPP(double v[5409],ElementSpec *es,ElementData *ed,NodeSpec **ns,NodeData **nd,double *rdata,int *idata,double **gpost,double **npost);
FILE *SMTFile;
int MMAInitialisationCode[]={
0,0
};


void CondensationData(ElementSpec *es,ElementData *ed, NodeSpec **ns,NodeData **nd, double *rdata,int *idata,double **h, 
                      double **ha, double **hb,double **dh){
 *h=&ed->ht[0];
 *hb=&ed->ht[9];
 *ha=&ed->ht[18];
 if(idata[ID_NoSensParameters]>0){
  *dh=&ed->ht[234];
 }else{
  *dh=NULL;
 };
};

DLLEXPORT int SMTSetElSpec(ElementSpec *es,int *idata,int ic,int ng)
{ int intc,nd,i;
  static int pn[31]={1, 4, 3, 2, 0, 2, 6, 5, 1, 0, 1, 5, 8, 4, 0, 4, 8, 7, 3, 0, 3, 7, 6, 2, 0, 7, 8, 5, 6, 0, 0};
  static int dof[8]={3, 3, 3, 3, 3, 3, 3, 3};
  static int nsto[8]={0, 0, 0, 0, 0, 0, 0, 0};

  static int ndat[8];
  static char *nid[]={"D","D","D","D","D","D",
                       "D","D"};
  static char *gdcs[]={"Elastic modulus","Poisson ratio"};
  static double defd[]={21000e0,0.3e0,0e0};
  static char *gpcs[]={"Sxx","Sxy","Sxz","Syx","Syy","Syz",
                       "Szx","Szy","Szz","Exx","Exy","Exz",
                       "Eyx","Eyy","Eyz","Ezx","Ezy","Ezz"};
  static char *npcs[]={"DeformedMeshX","DeformedMeshY","DeformedMeshZ","u","v","w"};
  static char *sname[]={""};
  static char *idname[]={""};
  static int idindex[1];
  static char *rdname[]={""};
  static char *cswitch[]={""};
  static int iswitch[1]={0};
  static double dswitch[1]={0e0};
  static int rdindex[1];
  static int nspecs[8];
  static double version[3]={2.004,2.004,6.};
  static double pnweights[8]={1e0,1e0,1e0,1e0,1e0,1e0,
  1e0,1e0};
  static double rnodes[24]={-1e0,-1e0,-1e0,1e0,-1e0,-1e0,
  1e0,1e0,-1e0,-1e0,1e0,-1e0,
  -1e0,-1e0,1e0,1e0,-1e0,1e0,
  1e0,1e0,1e0,-1e0,1e0,1e0};
  es->ReferenceNodes=rnodes;
  es->Data=defd;
  es->id.NoGroupData=2;
  es->Code="ExamplesHypersolid3D";es->Version=version;
  es->MainTitle="";
  es->SubTitle="";
  es->SubSubTitle="";
  es->Bibliography="";
  es->id.NoDimensions=3;es->id.NoDOFGlobal=24;es->id.NoDOFCondense=9;es->id.NoNodes=8;
  es->id.NoSegmentPoints=30;es->Segments=pn;es->PostNodeWeights=pnweights;
  es->id.NoIntSwitch=0;es->IntSwitch=iswitch;es->id.DemoLimitation=0;
  es->id.NoDoubleSwitch=0;es->DoubleSwitch=dswitch;
  es->id.NoCharSwitch=0;es->CharSwitch=cswitch;
  es->DOFGlobal=dof;es->NodeID=nid;es->id.NoGPostData=18;es->id.NoNPostData=6;
  es->id.SymmetricTangent=1;es->id.CreateDummyNodes=0;es->id.PostIterationCall=0;
  es->Topology="H1";es->GroupDataNames=gdcs;es->GPostNames=gpcs;es->NPostNames=npcs;
  es->AdditionalNodes="{}&";
  es->AdditionalGraphics="{}&";
  es->MMAInitialisation=MMAInitialisationCode;
  es->MMANextStep="";
  es->MMAStepBack="";
  es->MMAPreIteration="";
  es->IDataNames=idname;es->IDataIndex=idindex;es->RDataNames=rdname;es->RDataIndex=rdindex;
  es->id.NoIData=0;es->id.NoRData=0;
  es->id.ShapeSensitivity=0;es->id.NoSensNames=0;es->SensitivityNames=sname;es->NodeSpecs=nspecs;
  es->user.SPP=SPP;es->user.SKR=SKR;
  es->user.CondensationData=CondensationData;
  es->id.DefaultIntegrationCode=7;
  if(ic==-1){intc=7;} else {intc=ic;};
  es->id.IntCode=intc;
  SMTMultiIntPoints(&intc,idata,&es->id.NoIntPoints,
      &es->id.NoIntPointsA,&es->id.NoIntPointsB,&es->id.NoIntPointsC,0);
  es->id.NoAdditionalData=0;
  es->id.NoTimeStorage=234+9*idata[ID_NoSensParameters];
  es->id.NoElementData=0;
  nd=es->id.NoDimensions*idata[ID_NoShapeParameters];for(i=0;i<8;i++)ndat[i]=nd;
  es->NoNodeStorage=nsto;es->NoNodeData=ndat;
 if(1){
  return 0;
 }else{
  return 1;
 };
};

/******************* S U B R O U T I N E *********************/
void SKR(double v[5409],ElementSpec *es,ElementData *ed,NodeSpec **ns
     ,NodeData **nd,double *rdata,int *idata,double *p,double **s)
{
int i1,i250,i275;
v[200]=ed->ht[8];
v[199]=ed->ht[7];
v[198]=ed->ht[6];
v[197]=ed->ht[5];
v[196]=ed->ht[4];
v[195]=ed->ht[3];
v[194]=ed->ht[2];
v[193]=ed->ht[1];
v[192]=ed->ht[0];
v[51]=es->Data[1];
v[583]=2e0*v[51];
v[247]=1e0/(1e0+v[51]);
v[50]=es->Data[0];
v[248]=(v[247]*v[50])/2e0;
v[246]=(v[248]*v[583])/(1e0-v[583]);
v[49]=nd[7]->at[2];
v[48]=nd[6]->at[2];
v[47]=nd[5]->at[2];
v[46]=nd[4]->at[2];
v[45]=nd[3]->at[2];
v[44]=nd[2]->at[2];
v[43]=nd[1]->at[2];
v[42]=nd[0]->at[2];
v[41]=nd[7]->at[1];
v[40]=nd[6]->at[1];
v[39]=nd[5]->at[1];
v[38]=nd[4]->at[1];
v[37]=nd[3]->at[1];
v[36]=nd[2]->at[1];
v[35]=nd[1]->at[1];
v[34]=nd[0]->at[1];
v[33]=nd[7]->at[0];
v[32]=nd[6]->at[0];
v[31]=nd[5]->at[0];
v[30]=nd[4]->at[0];
v[29]=nd[3]->at[0];
v[28]=nd[2]->at[0];
v[27]=nd[1]->at[0];
v[26]=nd[0]->at[0];
v[25]=nd[7]->X[2];
v[24]=nd[6]->X[2];
v[23]=nd[5]->X[2];
v[22]=nd[4]->X[2];
v[21]=nd[3]->X[2];
v[591]=-v[21]+v[23];
v[20]=nd[2]->X[2];
v[19]=nd[1]->X[2];
v[590]=-v[19]+v[25];
v[18]=nd[0]->X[2];
v[592]=-v[18]+v[24];
v[187]=v[20]-v[22]+v[592];
v[17]=nd[7]->X[1];
v[16]=nd[6]->X[1];
v[15]=nd[5]->X[1];
v[14]=nd[4]->X[1];
v[13]=nd[3]->X[1];
v[588]=-v[13]+v[15];
v[12]=nd[2]->X[1];
v[11]=nd[1]->X[1];
v[587]=-v[11]+v[17];
v[10]=nd[0]->X[1];
v[589]=-v[10]+v[16];
v[185]=v[12]-v[14]+v[589];
v[9]=nd[7]->X[0];
v[8]=nd[6]->X[0];
v[7]=nd[5]->X[0];
v[6]=nd[4]->X[0];
v[5]=nd[3]->X[0];
v[585]=-v[5]+v[7];
v[4]=nd[2]->X[0];
v[3]=nd[1]->X[0];
v[584]=-v[3]+v[9];
v[2]=nd[0]->X[0];
v[586]=-v[2]+v[8];
v[183]=v[4]+v[586]-v[6];
for(i1=1;i1<=(int)(es->id.NoIntPoints);i1++){
 v[52]=es->IntPoints[4*(-1+i1)];
 v[63]=1e0-v[52];
 v[59]=1e0+v[52];
 v[53]=es->IntPoints[1+4*(-1+i1)];
 v[69]=1e0+v[53];
 v[78]=v[69]/8e0;
 v[86]=-(v[63]*v[78]);
 v[85]=-(v[59]*v[78]);
 v[65]=1e0-v[53];
 v[79]=v[65]/8e0;
 v[84]=-(v[59]*v[79]);
 v[82]=-(v[63]*v[79]);
 v[54]=es->IntPoints[2+4*(-1+i1)];
 v[70]=1e0+v[54];
 v[77]=v[70]/8e0;
 v[88]=-(v[59]*v[77]);
 v[87]=-(v[63]*v[77]);
 v[72]=v[69]*v[77];
 v[67]=v[65]*v[77];
 v[60]=1e0-v[54];
 v[80]=v[60]/8e0;
 v[83]=-(v[59]*v[80]);
 v[81]=-(v[63]*v[80]);
 v[62]=v[69]*v[80];
 v[57]=v[65]*v[80];
 v[55]=es->IntPoints[3+4*(-1+i1)];
 v[89]=(-v[2]+v[3])*v[57]+(v[4]-v[5])*v[62]+v[67]*(-v[6]+v[7])+v[72]*(v[8]-v[9]);
 v[180]=(v[183]+v[3]+v[585]-v[9])/8e0;
 v[90]=(-v[10]+v[11])*v[57]+(v[12]-v[13])*v[62]+(-v[14]+v[15])*v[67]+(v[16]-v[17])*v[72];
 v[181]=(v[11]-v[17]+v[185]+v[588])/8e0;
 v[91]=(-v[18]+v[19])*v[57]+(v[20]-v[21])*v[62]+(-v[22]+v[23])*v[67]+(v[24]-v[25])*v[72];
 v[182]=(v[187]+v[19]-v[25]+v[591])/8e0;
 v[92]=(v[2]-v[5])*v[81]+(v[3]-v[4])*v[83]+(v[7]-v[8])*v[88]+v[87]*(v[6]-v[9]);
 v[184]=(v[183]+v[5]+v[584]-v[7])/8e0;
 v[93]=(v[10]-v[13])*v[81]+(v[11]-v[12])*v[83]+(v[14]-v[17])*v[87]+(v[15]-v[16])*v[88];
 v[596]=-(v[90]*v[92])+v[89]*v[93];
 v[186]=(v[13]-v[15]+v[185]+v[587])/8e0;
 v[222]=-(v[181]*v[184])+v[180]*v[186];
 v[94]=(v[18]-v[21])*v[81]+(v[19]-v[20])*v[83]+(v[22]-v[25])*v[87]+(v[23]-v[24])*v[88];
 v[595]=v[91]*v[92]-v[89]*v[94];
 v[594]=-(v[91]*v[93])+v[90]*v[94];
 v[188]=(v[187]+v[21]-v[23]+v[590])/8e0;
 v[218]=v[182]*v[184]-v[180]*v[188];
 v[211]=-(v[182]*v[186])+v[181]*v[188];
 v[95]=(v[2]-v[6])*v[82]+(v[3]-v[7])*v[84]+(v[4]-v[8])*v[85]+v[86]*(v[5]-v[9]);
 v[189]=(-v[4]+v[584]+v[585]+v[586]+v[6])/8e0;
 v[96]=(v[10]-v[14])*v[82]+(v[11]-v[15])*v[84]+(v[12]-v[16])*v[85]+(v[13]-v[17])*v[86];
 v[190]=(-v[12]+v[14]+v[587]+v[588]+v[589])/8e0;
 v[221]=v[181]*v[189]-v[180]*v[190];
 v[220]=-(v[186]*v[189])+v[184]*v[190];
 v[97]=(v[18]-v[22])*v[82]+(v[19]-v[23])*v[84]+(v[20]-v[24])*v[85]+(v[21]-v[25])*v[86];
 v[191]=(-v[20]+v[22]+v[590]+v[591]+v[592])/8e0;
 v[216]=-(v[182]*v[189])+v[180]*v[191];
 v[214]=v[188]*v[189]-v[184]*v[191];
 v[210]=v[182]*v[190]-v[181]*v[191];
 v[209]=-(v[188]*v[190])+v[186]*v[191];
 v[202]=v[189]*v[211]+v[190]*v[218]+v[191]*v[222];
 v[599]=v[52]/v[202];
 v[598]=v[53]/v[202];
 v[597]=v[54]/v[202];
 v[98]=v[594]*v[95]+v[595]*v[96]+v[596]*v[97];
 v[601]=v[55]*v[98];
 v[593]=v[202]/v[98];
 v[271]=v[211]*v[593];
 v[270]=v[218]*v[593];
 v[269]=v[222]*v[593];
 v[268]=v[210]*v[593];
 v[267]=v[216]*v[593];
 v[266]=v[221]*v[593];
 v[265]=v[209]*v[593];
 v[264]=v[214]*v[593];
 v[263]=v[220]*v[593];
 v[102]=(-(v[94]*v[96])+v[93]*v[97])/v[98];
 v[165]=v[102]*v[72];
 v[156]=-(v[102]*v[67]);
 v[147]=v[102]*v[62];
 v[138]=-(v[102]*v[57]);
 v[103]=(v[94]*v[95]-v[92]*v[97])/v[98];
 v[167]=v[103]*v[72];
 v[158]=-(v[103]*v[67]);
 v[149]=v[103]*v[62];
 v[140]=-(v[103]*v[57]);
 v[104]=(-(v[93]*v[95])+v[92]*v[96])/v[98];
 v[169]=v[104]*v[72];
 v[160]=-(v[104]*v[67]);
 v[151]=v[104]*v[62];
 v[142]=-(v[104]*v[57]);
 v[105]=(v[91]*v[96]-v[90]*v[97])/v[98];
 v[126]=v[105]*v[80];
 v[114]=v[105]*v[77];
 v[106]=(-(v[91]*v[95])+v[89]*v[97])/v[98];
 v[129]=v[106]*v[80];
 v[116]=v[106]*v[77];
 v[107]=(v[90]*v[95]-v[89]*v[96])/v[98];
 v[132]=v[107]*v[80];
 v[118]=v[107]*v[77];
 v[108]=v[594]/v[98];
 v[127]=v[108]*v[79];
 v[120]=v[108]*v[78];
 v[109]=v[595]/v[98];
 v[130]=v[109]*v[79];
 v[122]=v[109]*v[78];
 v[110]=v[596]/v[98];
 v[133]=v[110]*v[79];
 v[124]=v[110]*v[78];
 v[111]=v[114]+v[120];
 v[112]=v[116]+v[122];
 v[113]=v[118]+v[124];
 v[115]=-v[114]+v[127];
 v[117]=-v[116]+v[130];
 v[119]=-v[118]+v[133];
 v[121]=-v[120]+v[126];
 v[123]=-v[122]+v[129];
 v[125]=-v[124]+v[132];
 v[128]=-v[126]-v[127];
 v[131]=-v[129]-v[130];
 v[134]=-v[132]-v[133];
 v[135]=v[138]+v[128]*v[63];
 v[136]=v[140]+v[131]*v[63];
 v[137]=v[142]+v[134]*v[63];
 v[139]=-v[138]+v[128]*v[59];
 v[141]=-v[140]+v[131]*v[59];
 v[143]=-v[142]+v[134]*v[59];
 v[144]=v[147]+v[121]*v[59];
 v[145]=v[149]+v[123]*v[59];
 v[146]=v[151]+v[125]*v[59];
 v[148]=-v[147]+v[121]*v[63];
 v[150]=-v[149]+v[123]*v[63];
 v[152]=-v[151]+v[125]*v[63];
 v[153]=v[156]+v[115]*v[63];
 v[154]=v[158]+v[117]*v[63];
 v[155]=v[160]+v[119]*v[63];
 v[157]=-v[156]+v[115]*v[59];
 v[159]=-v[158]+v[117]*v[59];
 v[161]=-v[160]+v[119]*v[59];
 v[162]=v[165]+v[111]*v[59];
 v[163]=v[167]+v[112]*v[59];
 v[164]=v[169]+v[113]*v[59];
 v[166]=-v[165]+v[111]*v[63];
 v[168]=-v[167]+v[112]*v[63];
 v[170]=-v[169]+v[113]*v[63];
 v[291]=v[269]*v[597];
 v[290]=v[266]*v[598];
 v[289]=v[263]*v[599];
 v[5342]=0e0;
 v[5343]=0e0;
 v[5344]=v[137];
 v[5345]=0e0;
 v[5346]=0e0;
 v[5347]=v[143];
 v[5348]=0e0;
 v[5349]=0e0;
 v[5350]=v[146];
 v[5351]=0e0;
 v[5352]=0e0;
 v[5353]=v[152];
 v[5354]=0e0;
 v[5355]=0e0;
 v[5356]=v[155];
 v[5357]=0e0;
 v[5358]=0e0;
 v[5359]=v[161];
 v[5360]=0e0;
 v[5361]=0e0;
 v[5362]=v[164];
 v[5363]=0e0;
 v[5364]=0e0;
 v[5365]=v[170];
 v[5366]=0e0;
 v[5367]=0e0;
 v[5368]=0e0;
 v[5369]=0e0;
 v[5370]=0e0;
 v[5371]=0e0;
 v[5372]=v[289];
 v[5373]=v[290];
 v[5374]=v[291];
 v[5243]=0e0;
 v[5244]=v[137];
 v[5245]=0e0;
 v[5246]=0e0;
 v[5247]=v[143];
 v[5248]=0e0;
 v[5249]=0e0;
 v[5250]=v[146];
 v[5251]=0e0;
 v[5252]=0e0;
 v[5253]=v[152];
 v[5254]=0e0;
 v[5255]=0e0;
 v[5256]=v[155];
 v[5257]=0e0;
 v[5258]=0e0;
 v[5259]=v[161];
 v[5260]=0e0;
 v[5261]=0e0;
 v[5262]=v[164];
 v[5263]=0e0;
 v[5264]=0e0;
 v[5265]=v[170];
 v[5266]=0e0;
 v[5267]=0e0;
 v[5268]=0e0;
 v[5269]=0e0;
 v[5270]=v[289];
 v[5271]=v[290];
 v[5272]=v[291];
 v[5273]=0e0;
 v[5274]=0e0;
 v[5275]=0e0;
 v[5144]=v[137];
 v[5145]=0e0;
 v[5146]=0e0;
 v[5147]=v[143];
 v[5148]=0e0;
 v[5149]=0e0;
 v[5150]=v[146];
 v[5151]=0e0;
 v[5152]=0e0;
 v[5153]=v[152];
 v[5154]=0e0;
 v[5155]=0e0;
 v[5156]=v[155];
 v[5157]=0e0;
 v[5158]=0e0;
 v[5159]=v[161];
 v[5160]=0e0;
 v[5161]=0e0;
 v[5162]=v[164];
 v[5163]=0e0;
 v[5164]=0e0;
 v[5165]=v[170];
 v[5166]=0e0;
 v[5167]=0e0;
 v[5168]=v[289];
 v[5169]=v[290];
 v[5170]=v[291];
 v[5171]=0e0;
 v[5172]=0e0;
 v[5173]=0e0;
 v[5174]=0e0;
 v[5175]=0e0;
 v[5176]=0e0;
 v[287]=v[270]*v[597];
 v[286]=v[267]*v[598];
 v[285]=v[264]*v[599];
 v[5309]=0e0;
 v[5310]=0e0;
 v[5311]=v[136];
 v[5312]=0e0;
 v[5313]=0e0;
 v[5314]=v[141];
 v[5315]=0e0;
 v[5316]=0e0;
 v[5317]=v[145];
 v[5318]=0e0;
 v[5319]=0e0;
 v[5320]=v[150];
 v[5321]=0e0;
 v[5322]=0e0;
 v[5323]=v[154];
 v[5324]=0e0;
 v[5325]=0e0;
 v[5326]=v[159];
 v[5327]=0e0;
 v[5328]=0e0;
 v[5329]=v[163];
 v[5330]=0e0;
 v[5331]=0e0;
 v[5332]=v[168];
 v[5333]=0e0;
 v[5334]=0e0;
 v[5335]=0e0;
 v[5336]=0e0;
 v[5337]=0e0;
 v[5338]=0e0;
 v[5339]=v[285];
 v[5340]=v[286];
 v[5341]=v[287];
 v[5210]=0e0;
 v[5211]=v[136];
 v[5212]=0e0;
 v[5213]=0e0;
 v[5214]=v[141];
 v[5215]=0e0;
 v[5216]=0e0;
 v[5217]=v[145];
 v[5218]=0e0;
 v[5219]=0e0;
 v[5220]=v[150];
 v[5221]=0e0;
 v[5222]=0e0;
 v[5223]=v[154];
 v[5224]=0e0;
 v[5225]=0e0;
 v[5226]=v[159];
 v[5227]=0e0;
 v[5228]=0e0;
 v[5229]=v[163];
 v[5230]=0e0;
 v[5231]=0e0;
 v[5232]=v[168];
 v[5233]=0e0;
 v[5234]=0e0;
 v[5235]=0e0;
 v[5236]=0e0;
 v[5237]=v[285];
 v[5238]=v[286];
 v[5239]=v[287];
 v[5240]=0e0;
 v[5241]=0e0;
 v[5242]=0e0;
 v[5111]=v[136];
 v[5112]=0e0;
 v[5113]=0e0;
 v[5114]=v[141];
 v[5115]=0e0;
 v[5116]=0e0;
 v[5117]=v[145];
 v[5118]=0e0;
 v[5119]=0e0;
 v[5120]=v[150];
 v[5121]=0e0;
 v[5122]=0e0;
 v[5123]=v[154];
 v[5124]=0e0;
 v[5125]=0e0;
 v[5126]=v[159];
 v[5127]=0e0;
 v[5128]=0e0;
 v[5129]=v[163];
 v[5130]=0e0;
 v[5131]=0e0;
 v[5132]=v[168];
 v[5133]=0e0;
 v[5134]=0e0;
 v[5135]=v[285];
 v[5136]=v[286];
 v[5137]=v[287];
 v[5138]=0e0;
 v[5139]=0e0;
 v[5140]=0e0;
 v[5141]=0e0;
 v[5142]=0e0;
 v[5143]=0e0;
 v[283]=v[271]*v[597];
 v[282]=v[268]*v[598];
 v[281]=v[265]*v[599];
 v[5276]=0e0;
 v[5277]=0e0;
 v[5278]=v[135];
 v[5279]=0e0;
 v[5280]=0e0;
 v[5281]=v[139];
 v[5282]=0e0;
 v[5283]=0e0;
 v[5284]=v[144];
 v[5285]=0e0;
 v[5286]=0e0;
 v[5287]=v[148];
 v[5288]=0e0;
 v[5289]=0e0;
 v[5290]=v[153];
 v[5291]=0e0;
 v[5292]=0e0;
 v[5293]=v[157];
 v[5294]=0e0;
 v[5295]=0e0;
 v[5296]=v[162];
 v[5297]=0e0;
 v[5298]=0e0;
 v[5299]=v[166];
 v[5300]=0e0;
 v[5301]=0e0;
 v[5302]=0e0;
 v[5303]=0e0;
 v[5304]=0e0;
 v[5305]=0e0;
 v[5306]=v[281];
 v[5307]=v[282];
 v[5308]=v[283];
 v[5177]=0e0;
 v[5178]=v[135];
 v[5179]=0e0;
 v[5180]=0e0;
 v[5181]=v[139];
 v[5182]=0e0;
 v[5183]=0e0;
 v[5184]=v[144];
 v[5185]=0e0;
 v[5186]=0e0;
 v[5187]=v[148];
 v[5188]=0e0;
 v[5189]=0e0;
 v[5190]=v[153];
 v[5191]=0e0;
 v[5192]=0e0;
 v[5193]=v[157];
 v[5194]=0e0;
 v[5195]=0e0;
 v[5196]=v[162];
 v[5197]=0e0;
 v[5198]=0e0;
 v[5199]=v[166];
 v[5200]=0e0;
 v[5201]=0e0;
 v[5202]=0e0;
 v[5203]=0e0;
 v[5204]=v[281];
 v[5205]=v[282];
 v[5206]=v[283];
 v[5207]=0e0;
 v[5208]=0e0;
 v[5209]=0e0;
 v[5078]=v[135];
 v[5079]=0e0;
 v[5080]=0e0;
 v[5081]=v[139];
 v[5082]=0e0;
 v[5083]=0e0;
 v[5084]=v[144];
 v[5085]=0e0;
 v[5086]=0e0;
 v[5087]=v[148];
 v[5088]=0e0;
 v[5089]=0e0;
 v[5090]=v[153];
 v[5091]=0e0;
 v[5092]=0e0;
 v[5093]=v[157];
 v[5094]=0e0;
 v[5095]=0e0;
 v[5096]=v[162];
 v[5097]=0e0;
 v[5098]=0e0;
 v[5099]=v[166];
 v[5100]=0e0;
 v[5101]=0e0;
 v[5102]=v[281];
 v[5103]=v[282];
 v[5104]=v[283];
 v[5105]=0e0;
 v[5106]=0e0;
 v[5107]=0e0;
 v[5108]=0e0;
 v[5109]=0e0;
 v[5110]=0e0;
 v[227]=v[200]*v[597];
 v[226]=v[199]*v[598];
 v[225]=v[198]*v[599];
 v[217]=v[197]*v[597];
 v[215]=v[196]*v[598];
 v[213]=v[195]*v[599];
 v[206]=v[194]*v[597];
 v[205]=v[193]*v[598];
 v[204]=v[192]*v[599];
 v[230]=1e0+v[135]*v[26]+v[139]*v[27]+v[144]*v[28]+v[148]*v[29]+v[153]*v[30]+v[157]*v[31]
  +v[162]*v[32]+v[166]*v[33]+(v[204]*v[209]+v[205]*v[210]+v[206]*v[211])*v[593];
 v[231]=v[136]*v[26]+v[141]*v[27]+v[145]*v[28]+v[150]*v[29]+v[154]*v[30]+v[159]*v[31]+v[163]*v[32]
  +v[168]*v[33]+(v[204]*v[214]+v[205]*v[216]+v[206]*v[218])*v[593];
 v[232]=v[137]*v[26]+v[143]*v[27]+v[146]*v[28]+v[152]*v[29]+v[155]*v[30]+v[161]*v[31]+v[164]*v[32]
  +v[170]*v[33]+(v[204]*v[220]+v[205]*v[221]+v[206]*v[222])*v[593];
 v[233]=v[135]*v[34]+v[139]*v[35]+v[144]*v[36]+v[148]*v[37]+v[153]*v[38]+v[157]*v[39]+v[162]*v[40]
  +v[166]*v[41]+(v[209]*v[213]+v[210]*v[215]+v[211]*v[217])*v[593];
 v[234]=1e0+v[136]*v[34]+v[141]*v[35]+v[145]*v[36]+v[150]*v[37]+v[154]*v[38]+v[159]*v[39]
  +v[163]*v[40]+v[168]*v[41]+(v[213]*v[214]+v[215]*v[216]+v[217]*v[218])*v[593];
 v[304]=-(v[231]*v[233])+v[230]*v[234];
 v[235]=v[137]*v[34]+v[143]*v[35]+v[146]*v[36]+v[152]*v[37]+v[155]*v[38]+v[161]*v[39]+v[164]*v[40]
  +v[170]*v[41]+(v[213]*v[220]+v[215]*v[221]+v[217]*v[222])*v[593];
 v[303]=v[232]*v[233]-v[230]*v[235];
 v[302]=-(v[232]*v[234])+v[231]*v[235];
 v[236]=v[135]*v[42]+v[139]*v[43]+v[144]*v[44]+v[148]*v[45]+v[153]*v[46]+v[157]*v[47]+v[162]*v[48]
  +v[166]*v[49]+(v[209]*v[225]+v[210]*v[226]+v[211]*v[227])*v[593];
 v[237]=v[136]*v[42]+v[141]*v[43]+v[145]*v[44]+v[150]*v[45]+v[154]*v[46]+v[159]*v[47]+v[163]*v[48]
  +v[168]*v[49]+(v[214]*v[225]+v[216]*v[226]+v[218]*v[227])*v[593];
 v[301]=v[231]*v[236]-v[230]*v[237];
 v[298]=-(v[234]*v[236])+v[233]*v[237];
 v[238]=1e0+v[137]*v[42]+v[143]*v[43]+v[146]*v[44]+v[152]*v[45]+v[155]*v[46]+v[161]*v[47]
  +v[164]*v[48]+v[170]*v[49]+(v[220]*v[225]+v[221]*v[226]+v[222]*v[227])*v[593];
 v[300]=-(v[232]*v[236])+v[230]*v[238];
 v[299]=v[232]*v[237]-v[231]*v[238];
 v[297]=v[235]*v[236]-v[233]*v[238];
 v[296]=-(v[235]*v[237])+v[234]*v[238];
 v[245]=v[236]*v[302]+v[237]*v[303]+v[238]*v[304];
 v[253]=(-1e0+v[245])*v[246]-v[248]/v[245];
 v[254]=v[238]*v[248]+v[253]*v[304];
 v[255]=v[237]*v[248]+v[253]*v[303];
 v[256]=v[236]*v[248]+v[253]*v[302];
 v[257]=v[235]*v[248]+v[253]*v[301];
 v[258]=v[234]*v[248]+v[253]*v[300];
 v[259]=v[233]*v[248]+v[253]*v[299];
 v[260]=v[232]*v[248]+v[253]*v[298];
 v[261]=v[231]*v[248]+v[253]*v[297];
 v[262]=v[230]*v[248]+v[253]*v[296];
 v[5041]=v[137]*v[260]+v[136]*v[261]+v[135]*v[262];
 v[5042]=v[137]*v[257]+v[136]*v[258]+v[135]*v[259];
 v[5043]=v[137]*v[254]+v[136]*v[255]+v[135]*v[256];
 v[5044]=v[143]*v[260]+v[141]*v[261]+v[139]*v[262];
 v[5045]=v[143]*v[257]+v[141]*v[258]+v[139]*v[259];
 v[5046]=v[143]*v[254]+v[141]*v[255]+v[139]*v[256];
 v[5047]=v[146]*v[260]+v[145]*v[261]+v[144]*v[262];
 v[5048]=v[146]*v[257]+v[145]*v[258]+v[144]*v[259];
 v[5049]=v[146]*v[254]+v[145]*v[255]+v[144]*v[256];
 v[5050]=v[152]*v[260]+v[150]*v[261]+v[148]*v[262];
 v[5051]=v[152]*v[257]+v[150]*v[258]+v[148]*v[259];
 v[5052]=v[152]*v[254]+v[150]*v[255]+v[148]*v[256];
 v[5053]=v[155]*v[260]+v[154]*v[261]+v[153]*v[262];
 v[5054]=v[155]*v[257]+v[154]*v[258]+v[153]*v[259];
 v[5055]=v[155]*v[254]+v[154]*v[255]+v[153]*v[256];
 v[5056]=v[161]*v[260]+v[159]*v[261]+v[157]*v[262];
 v[5057]=v[161]*v[257]+v[159]*v[258]+v[157]*v[259];
 v[5058]=v[161]*v[254]+v[159]*v[255]+v[157]*v[256];
 v[5059]=v[164]*v[260]+v[163]*v[261]+v[162]*v[262];
 v[5060]=v[164]*v[257]+v[163]*v[258]+v[162]*v[259];
 v[5061]=v[164]*v[254]+v[163]*v[255]+v[162]*v[256];
 v[5062]=v[170]*v[260]+v[168]*v[261]+v[166]*v[262];
 v[5063]=v[170]*v[257]+v[168]*v[258]+v[166]*v[259];
 v[5064]=v[170]*v[254]+v[168]*v[255]+v[166]*v[256];
 v[5065]=(v[260]*v[263]+v[261]*v[264]+v[262]*v[265])*v[599];
 v[5066]=(v[260]*v[266]+v[261]*v[267]+v[262]*v[268])*v[598];
 v[5067]=(v[260]*v[269]+v[261]*v[270]+v[262]*v[271])*v[597];
 v[5068]=(v[257]*v[263]+v[258]*v[264]+v[259]*v[265])*v[599];
 v[5069]=(v[257]*v[266]+v[258]*v[267]+v[259]*v[268])*v[598];
 v[5070]=(v[257]*v[269]+v[258]*v[270]+v[259]*v[271])*v[597];
 v[5071]=(v[254]*v[263]+v[255]*v[264]+v[256]*v[265])*v[599];
 v[5072]=(v[254]*v[266]+v[255]*v[267]+v[256]*v[268])*v[598];
 v[5073]=(v[254]*v[269]+v[255]*v[270]+v[256]*v[271])*v[597];
 for(i250=1;i250<=33;i250++){
  v[278]=v[5077+i250]*v[601];
  v[279]=v[5110+i250]*v[601];
  v[280]=v[5143+i250]*v[601];
  v[284]=v[5176+i250]*v[601];
  v[288]=v[5209+i250]*v[601];
  v[292]=v[5242+i250]*v[601];
  v[293]=v[5275+i250]*v[601];
  v[294]=v[5308+i250]*v[601];
  v[295]=v[5341+i250]*v[601];
  v[305]=(v[246]+v[248]/Power(v[245],2))*(v[278]*v[296]+v[279]*v[297]+v[280]*v[298]+v[284]*v[299]
   +v[288]*v[300]+v[292]*v[301]+v[293]*v[302]+v[294]*v[303]+v[295]*v[304]);
  v[306]=v[253]*(v[234]*v[278]-v[233]*v[279]-v[231]*v[284]+v[230]*v[288])+v[248]*v[295]
   +v[304]*v[305];
  v[307]=v[253]*(-(v[235]*v[278])+v[233]*v[280]+v[232]*v[284]-v[230]*v[292])+v[248]*v[294]
   +v[303]*v[305];
  v[308]=v[253]*(v[235]*v[279]-v[234]*v[280]-v[232]*v[288]+v[231]*v[292])+v[248]*v[293]
   +v[302]*v[305];
  v[309]=v[248]*v[292]+v[253]*(-(v[237]*v[278])+v[236]*v[279]+v[231]*v[293]-v[230]*v[294])
   +v[301]*v[305];
  v[310]=v[248]*v[288]+v[253]*(v[238]*v[278]-v[236]*v[280]-v[232]*v[293]+v[230]*v[295])
   +v[300]*v[305];
  v[311]=v[248]*v[284]+v[253]*(-(v[238]*v[279])+v[237]*v[280]+v[232]*v[294]-v[231]*v[295])
   +v[299]*v[305];
  v[312]=v[248]*v[280]+v[253]*(v[237]*v[284]-v[236]*v[288]-v[234]*v[293]+v[233]*v[294])
   +v[298]*v[305];
  v[313]=v[248]*v[279]+v[253]*(-(v[238]*v[284])+v[236]*v[292]+v[235]*v[293]-v[233]*v[295])
   +v[297]*v[305];
  v[314]=v[248]*v[278]+v[253]*(v[238]*v[288]-v[237]*v[292]-v[235]*v[294]+v[234]*v[295])
   +v[296]*v[305];
  v[5375]=v[137]*v[312]+v[136]*v[313]+v[135]*v[314];
  v[5376]=v[137]*v[309]+v[136]*v[310]+v[135]*v[311];
  v[5377]=v[137]*v[306]+v[136]*v[307]+v[135]*v[308];
  v[5378]=v[143]*v[312]+v[141]*v[313]+v[139]*v[314];
  v[5379]=v[143]*v[309]+v[141]*v[310]+v[139]*v[311];
  v[5380]=v[143]*v[306]+v[141]*v[307]+v[139]*v[308];
  v[5381]=v[146]*v[312]+v[145]*v[313]+v[144]*v[314];
  v[5382]=v[146]*v[309]+v[145]*v[310]+v[144]*v[311];
  v[5383]=v[146]*v[306]+v[145]*v[307]+v[144]*v[308];
  v[5384]=v[152]*v[312]+v[150]*v[313]+v[148]*v[314];
  v[5385]=v[152]*v[309]+v[150]*v[310]+v[148]*v[311];
  v[5386]=v[152]*v[306]+v[150]*v[307]+v[148]*v[308];
  v[5387]=v[155]*v[312]+v[154]*v[313]+v[153]*v[314];
  v[5388]=v[155]*v[309]+v[154]*v[310]+v[153]*v[311];
  v[5389]=v[155]*v[306]+v[154]*v[307]+v[153]*v[308];
  v[5390]=v[161]*v[312]+v[159]*v[313]+v[157]*v[314];
  v[5391]=v[161]*v[309]+v[159]*v[310]+v[157]*v[311];
  v[5392]=v[161]*v[306]+v[159]*v[307]+v[157]*v[308];
  v[5393]=v[164]*v[312]+v[163]*v[313]+v[162]*v[314];
  v[5394]=v[164]*v[309]+v[163]*v[310]+v[162]*v[311];
  v[5395]=v[164]*v[306]+v[163]*v[307]+v[162]*v[308];
  v[5396]=v[170]*v[312]+v[168]*v[313]+v[166]*v[314];
  v[5397]=v[170]*v[309]+v[168]*v[310]+v[166]*v[311];
  v[5398]=v[170]*v[306]+v[168]*v[307]+v[166]*v[308];
  v[5399]=(v[263]*v[312]+v[264]*v[313]+v[265]*v[314])*v[599];
  v[5400]=(v[266]*v[312]+v[267]*v[313]+v[268]*v[314])*v[598];
  v[5401]=(v[269]*v[312]+v[270]*v[313]+v[271]*v[314])*v[597];
  v[5402]=(v[263]*v[309]+v[264]*v[310]+v[265]*v[311])*v[599];
  v[5403]=(v[266]*v[309]+v[267]*v[310]+v[268]*v[311])*v[598];
  v[5404]=(v[269]*v[309]+v[270]*v[310]+v[271]*v[311])*v[597];
  v[5405]=(v[263]*v[306]+v[264]*v[307]+v[265]*v[308])*v[599];
  v[5406]=(v[266]*v[306]+v[267]*v[307]+v[268]*v[308])*v[598];
  v[5407]=(v[269]*v[306]+v[270]*v[307]+v[271]*v[308])*v[597];
  p[i250-1]+=v[5040+i250]*v[601];
  for(i275=i250;i275<=33;i275++){
   s[i250-1][i275-1]+=v[5374+i275];
  };/* end for */
 };/* end for */
};/* end for */
};

/******************* S U B R O U T I N E *********************/
void SPP(double v[5409],ElementSpec *es,ElementData *ed,NodeSpec **ns
     ,NodeData **nd,double *rdata,int *idata,double **gpost,double **npost)
{
int i318;
v[517]=ed->ht[8];
v[516]=ed->ht[7];
v[515]=ed->ht[6];
v[514]=ed->ht[5];
v[513]=ed->ht[4];
v[512]=ed->ht[3];
v[511]=ed->ht[2];
v[510]=ed->ht[1];
v[509]=ed->ht[0];
v[368]=es->Data[1];
v[610]=2e0*v[368];
v[564]=1e0/(1e0+v[368]);
v[367]=es->Data[0];
v[565]=(v[367]*v[564])/2e0;
v[563]=(v[565]*v[610])/(1e0-v[610]);
v[366]=nd[7]->at[2];
v[365]=nd[6]->at[2];
v[364]=nd[5]->at[2];
v[363]=nd[4]->at[2];
v[362]=nd[3]->at[2];
v[361]=nd[2]->at[2];
v[360]=nd[1]->at[2];
v[359]=nd[0]->at[2];
v[358]=nd[7]->at[1];
v[357]=nd[6]->at[1];
v[356]=nd[5]->at[1];
v[355]=nd[4]->at[1];
v[354]=nd[3]->at[1];
v[353]=nd[2]->at[1];
v[352]=nd[1]->at[1];
v[351]=nd[0]->at[1];
v[350]=nd[7]->at[0];
v[349]=nd[6]->at[0];
v[348]=nd[5]->at[0];
v[347]=nd[4]->at[0];
v[346]=nd[3]->at[0];
v[345]=nd[2]->at[0];
v[344]=nd[1]->at[0];
v[343]=nd[0]->at[0];
v[342]=nd[7]->X[2];
v[341]=nd[6]->X[2];
v[340]=nd[5]->X[2];
v[339]=nd[4]->X[2];
v[338]=nd[3]->X[2];
v[618]=-v[338]+v[340];
v[337]=nd[2]->X[2];
v[336]=nd[1]->X[2];
v[617]=-v[336]+v[342];
v[335]=nd[0]->X[2];
v[619]=-v[335]+v[341];
v[504]=v[337]-v[339]+v[619];
v[334]=nd[7]->X[1];
v[333]=nd[6]->X[1];
v[332]=nd[5]->X[1];
v[331]=nd[4]->X[1];
v[330]=nd[3]->X[1];
v[615]=-v[330]+v[332];
v[329]=nd[2]->X[1];
v[328]=nd[1]->X[1];
v[614]=-v[328]+v[334];
v[327]=nd[0]->X[1];
v[616]=-v[327]+v[333];
v[502]=v[329]-v[331]+v[616];
v[326]=nd[7]->X[0];
v[325]=nd[6]->X[0];
v[324]=nd[5]->X[0];
v[323]=nd[4]->X[0];
v[322]=nd[3]->X[0];
v[612]=-v[322]+v[324];
v[321]=nd[2]->X[0];
v[320]=nd[1]->X[0];
v[611]=-v[320]+v[326];
v[319]=nd[0]->X[0];
v[613]=-v[319]+v[325];
v[500]=v[321]-v[323]+v[613];
for(i318=1;i318<=(int)(es->id.NoIntPoints);i318++){
 v[369]=es->IntPoints[4*(-1+i318)];
 v[380]=1e0-v[369];
 v[376]=1e0+v[369];
 v[370]=es->IntPoints[1+4*(-1+i318)];
 v[386]=1e0+v[370];
 v[395]=v[386]/8e0;
 v[403]=-(v[380]*v[395]);
 v[402]=-(v[376]*v[395]);
 v[382]=1e0-v[370];
 v[396]=v[382]/8e0;
 v[401]=-(v[376]*v[396]);
 v[399]=-(v[380]*v[396]);
 v[371]=es->IntPoints[2+4*(-1+i318)];
 v[387]=1e0+v[371];
 v[394]=v[387]/8e0;
 v[405]=-(v[376]*v[394]);
 v[404]=-(v[380]*v[394]);
 v[389]=v[386]*v[394];
 v[384]=v[382]*v[394];
 v[377]=1e0-v[371];
 v[397]=v[377]/8e0;
 v[400]=-(v[376]*v[397]);
 v[398]=-(v[380]*v[397]);
 v[379]=v[386]*v[397];
 v[374]=v[382]*v[397];
 v[406]=(-v[319]+v[320])*v[374]+(v[321]-v[322])*v[379]+(-v[323]+v[324])*v[384]+(v[325]-v[326]
  )*v[389];
 v[497]=(v[320]-v[326]+v[500]+v[612])/8e0;
 v[407]=(-v[327]+v[328])*v[374]+(v[329]-v[330])*v[379]+(-v[331]+v[332])*v[384]+(v[333]-v[334]
  )*v[389];
 v[498]=(v[328]-v[334]+v[502]+v[615])/8e0;
 v[408]=(-v[335]+v[336])*v[374]+(v[337]-v[338])*v[379]+(-v[339]+v[340])*v[384]+(v[341]-v[342]
  )*v[389];
 v[499]=(v[336]-v[342]+v[504]+v[618])/8e0;
 v[409]=(v[319]-v[322])*v[398]+(v[320]-v[321])*v[400]+(v[323]-v[326])*v[404]+(v[324]-v[325]
  )*v[405];
 v[501]=(v[322]-v[324]+v[500]+v[611])/8e0;
 v[410]=(v[327]-v[330])*v[398]+(v[328]-v[329])*v[400]+(v[331]-v[334])*v[404]+(v[332]-v[333]
  )*v[405];
 v[622]=-(v[407]*v[409])+v[406]*v[410];
 v[503]=(v[330]-v[332]+v[502]+v[614])/8e0;
 v[539]=-(v[498]*v[501])+v[497]*v[503];
 v[411]=(v[335]-v[338])*v[398]+(v[336]-v[337])*v[400]+(v[339]-v[342])*v[404]+(v[340]-v[341]
  )*v[405];
 v[621]=v[408]*v[409]-v[406]*v[411];
 v[620]=-(v[408]*v[410])+v[407]*v[411];
 v[505]=(v[338]-v[340]+v[504]+v[617])/8e0;
 v[535]=v[499]*v[501]-v[497]*v[505];
 v[528]=-(v[499]*v[503])+v[498]*v[505];
 v[412]=(v[319]-v[323])*v[399]+(v[320]-v[324])*v[401]+(v[321]-v[325])*v[402]+(v[322]-v[326]
  )*v[403];
 v[506]=(-v[321]+v[323]+v[611]+v[612]+v[613])/8e0;
 v[413]=(v[327]-v[331])*v[399]+(v[328]-v[332])*v[401]+(v[329]-v[333])*v[402]+(v[330]-v[334]
  )*v[403];
 v[507]=(-v[329]+v[331]+v[614]+v[615]+v[616])/8e0;
 v[538]=v[498]*v[506]-v[497]*v[507];
 v[537]=-(v[503]*v[506])+v[501]*v[507];
 v[414]=(v[335]-v[339])*v[399]+(v[336]-v[340])*v[401]+(v[337]-v[341])*v[402]+(v[338]-v[342]
  )*v[403];
 v[508]=(-v[337]+v[339]+v[617]+v[618]+v[619])/8e0;
 v[533]=-(v[499]*v[506])+v[497]*v[508];
 v[531]=v[505]*v[506]-v[501]*v[508];
 v[527]=v[499]*v[507]-v[498]*v[508];
 v[526]=-(v[505]*v[507])+v[503]*v[508];
 v[519]=v[506]*v[528]+v[507]*v[535]+v[508]*v[539];
 v[625]=v[369]/v[519];
 v[624]=v[370]/v[519];
 v[623]=v[371]/v[519];
 v[415]=v[412]*v[620]+v[413]*v[621]+v[414]*v[622];
 v[626]=v[519]/v[415];
 v[419]=(-(v[411]*v[413])+v[410]*v[414])/v[415];
 v[482]=v[389]*v[419];
 v[473]=-(v[384]*v[419]);
 v[464]=v[379]*v[419];
 v[455]=-(v[374]*v[419]);
 v[420]=(v[411]*v[412]-v[409]*v[414])/v[415];
 v[484]=v[389]*v[420];
 v[475]=-(v[384]*v[420]);
 v[466]=v[379]*v[420];
 v[457]=-(v[374]*v[420]);
 v[421]=(-(v[410]*v[412])+v[409]*v[413])/v[415];
 v[486]=v[389]*v[421];
 v[477]=-(v[384]*v[421]);
 v[468]=v[379]*v[421];
 v[459]=-(v[374]*v[421]);
 v[422]=(v[408]*v[413]-v[407]*v[414])/v[415];
 v[443]=v[397]*v[422];
 v[431]=v[394]*v[422];
 v[423]=(-(v[408]*v[412])+v[406]*v[414])/v[415];
 v[446]=v[397]*v[423];
 v[433]=v[394]*v[423];
 v[424]=(v[407]*v[412]-v[406]*v[413])/v[415];
 v[449]=v[397]*v[424];
 v[435]=v[394]*v[424];
 v[425]=v[620]/v[415];
 v[444]=v[396]*v[425];
 v[437]=v[395]*v[425];
 v[426]=v[621]/v[415];
 v[447]=v[396]*v[426];
 v[439]=v[395]*v[426];
 v[427]=v[622]/v[415];
 v[450]=v[396]*v[427];
 v[441]=v[395]*v[427];
 v[428]=v[431]+v[437];
 v[429]=v[433]+v[439];
 v[430]=v[435]+v[441];
 v[432]=-v[431]+v[444];
 v[434]=-v[433]+v[447];
 v[436]=-v[435]+v[450];
 v[438]=-v[437]+v[443];
 v[440]=-v[439]+v[446];
 v[442]=-v[441]+v[449];
 v[445]=-v[443]-v[444];
 v[448]=-v[446]-v[447];
 v[451]=-v[449]-v[450];
 v[452]=v[380]*v[445]+v[455];
 v[453]=v[380]*v[448]+v[457];
 v[454]=v[380]*v[451]+v[459];
 v[456]=v[376]*v[445]-v[455];
 v[458]=v[376]*v[448]-v[457];
 v[460]=v[376]*v[451]-v[459];
 v[461]=v[376]*v[438]+v[464];
 v[462]=v[376]*v[440]+v[466];
 v[463]=v[376]*v[442]+v[468];
 v[465]=v[380]*v[438]-v[464];
 v[467]=v[380]*v[440]-v[466];
 v[469]=v[380]*v[442]-v[468];
 v[470]=v[380]*v[432]+v[473];
 v[471]=v[380]*v[434]+v[475];
 v[472]=v[380]*v[436]+v[477];
 v[474]=v[376]*v[432]-v[473];
 v[476]=v[376]*v[434]-v[475];
 v[478]=v[376]*v[436]-v[477];
 v[479]=v[376]*v[428]+v[482];
 v[480]=v[376]*v[429]+v[484];
 v[481]=v[376]*v[430]+v[486];
 v[483]=v[380]*v[428]-v[482];
 v[485]=v[380]*v[429]-v[484];
 v[487]=v[380]*v[430]-v[486];
 v[544]=v[517]*v[623];
 v[543]=v[516]*v[624];
 v[542]=v[515]*v[625];
 v[534]=v[514]*v[623];
 v[532]=v[513]*v[624];
 v[530]=v[512]*v[625];
 v[523]=v[511]*v[623];
 v[522]=v[510]*v[624];
 v[521]=v[509]*v[625];
 v[547]=1e0+v[343]*v[452]+v[344]*v[456]+v[345]*v[461]+v[346]*v[465]+v[347]*v[470]+v[348]*v[474]
  +v[349]*v[479]+v[350]*v[483]+(v[521]*v[526]+v[522]*v[527]+v[523]*v[528])*v[626];
 v[548]=v[343]*v[453]+v[344]*v[458]+v[345]*v[462]+v[346]*v[467]+v[347]*v[471]+v[348]*v[476]
  +v[349]*v[480]+v[350]*v[485]+(v[521]*v[531]+v[522]*v[533]+v[523]*v[535])*v[626];
 v[549]=v[343]*v[454]+v[344]*v[460]+v[345]*v[463]+v[346]*v[469]+v[347]*v[472]+v[348]*v[478]
  +v[349]*v[481]+v[350]*v[487]+(v[521]*v[537]+v[522]*v[538]+v[523]*v[539])*v[626];
 v[550]=v[351]*v[452]+v[352]*v[456]+v[353]*v[461]+v[354]*v[465]+v[355]*v[470]+v[356]*v[474]
  +v[357]*v[479]+v[358]*v[483]+(v[526]*v[530]+v[527]*v[532]+v[528]*v[534])*v[626];
 v[551]=1e0+v[351]*v[453]+v[352]*v[458]+v[353]*v[462]+v[354]*v[467]+v[355]*v[471]+v[356]*v[476]
  +v[357]*v[480]+v[358]*v[485]+(v[530]*v[531]+v[532]*v[533]+v[534]*v[535])*v[626];
 v[627]=-(v[548]*v[550])+v[547]*v[551];
 v[552]=v[351]*v[454]+v[352]*v[460]+v[353]*v[463]+v[354]*v[469]+v[355]*v[472]+v[356]*v[478]
  +v[357]*v[481]+v[358]*v[487]+(v[530]*v[537]+v[532]*v[538]+v[534]*v[539])*v[626];
 v[629]=-(v[549]*v[551])+v[548]*v[552];
 v[628]=v[549]*v[550]-v[547]*v[552];
 v[553]=v[359]*v[452]+v[360]*v[456]+v[361]*v[461]+v[362]*v[465]+v[363]*v[470]+v[364]*v[474]
  +v[365]*v[479]+v[366]*v[483]+(v[526]*v[542]+v[527]*v[543]+v[528]*v[544])*v[626];
 v[554]=v[359]*v[453]+v[360]*v[458]+v[361]*v[462]+v[362]*v[467]+v[363]*v[471]+v[364]*v[476]
  +v[365]*v[480]+v[366]*v[485]+(v[531]*v[542]+v[533]*v[543]+v[535]*v[544])*v[626];
 v[555]=1e0+v[359]*v[454]+v[360]*v[460]+v[361]*v[463]+v[362]*v[469]+v[363]*v[472]+v[364]*v[478]
  +v[365]*v[481]+v[366]*v[487]+(v[537]*v[542]+v[538]*v[543]+v[539]*v[544])*v[626];
 v[562]=v[555]*v[627]+v[554]*v[628]+v[553]*v[629];
 npost[0][0]=v[343];
 npost[0][1]=v[351];
 npost[0][2]=v[359];
 npost[0][3]=v[343];
 npost[0][4]=v[351];
 npost[0][5]=v[359];
 npost[1][0]=v[344];
 npost[1][1]=v[352];
 npost[1][2]=v[360];
 npost[1][3]=v[344];
 npost[1][4]=v[352];
 npost[1][5]=v[360];
 npost[2][0]=v[345];
 npost[2][1]=v[353];
 npost[2][2]=v[361];
 npost[2][3]=v[345];
 npost[2][4]=v[353];
 npost[2][5]=v[361];
 npost[3][0]=v[346];
 npost[3][1]=v[354];
 npost[3][2]=v[362];
 npost[3][3]=v[346];
 npost[3][4]=v[354];
 npost[3][5]=v[362];
 npost[4][0]=v[347];
 npost[4][1]=v[355];
 npost[4][2]=v[363];
 npost[4][3]=v[347];
 npost[4][4]=v[355];
 npost[4][5]=v[363];
 npost[5][0]=v[348];
 npost[5][1]=v[356];
 npost[5][2]=v[364];
 npost[5][3]=v[348];
 npost[5][4]=v[356];
 npost[5][5]=v[364];
 npost[6][0]=v[349];
 npost[6][1]=v[357];
 npost[6][2]=v[365];
 npost[6][3]=v[349];
 npost[6][4]=v[357];
 npost[6][5]=v[365];
 npost[7][0]=v[350];
 npost[7][1]=v[358];
 npost[7][2]=v[366];
 npost[7][3]=v[350];
 npost[7][4]=v[358];
 npost[7][5]=v[366];
 v[569]=(-1e0+v[562])*v[563]-v[565]/v[562];
 v[570]=v[549]*v[565]+(-(v[551]*v[553])+v[550]*v[554])*v[569];
 v[571]=v[548]*v[565]+(v[552]*v[553]-v[550]*v[555])*v[569];
 v[572]=v[547]*v[565]+(-(v[552]*v[554])+v[551]*v[555])*v[569];
 v[573]=(v[552]*v[570]+v[551]*v[571]+v[550]*v[572])/v[562];
 v[574]=(v[555]*v[570]+v[554]*v[571]+v[553]*v[572])/v[562];
 v[575]=v[552]*v[565]+(v[548]*v[553]-v[547]*v[554])*v[569];
 v[576]=v[551]*v[565]+(-(v[549]*v[553])+v[547]*v[555])*v[569];
 v[577]=v[550]*v[565]+(v[549]*v[554]-v[548]*v[555])*v[569];
 v[578]=(v[555]*v[575]+v[554]*v[576]+v[553]*v[577])/v[562];
 v[579]=(v[547]*v[548]+v[550]*v[551]+v[553]*v[554])/2e0;
 v[580]=(v[547]*v[549]+v[550]*v[552]+v[553]*v[555])/2e0;
 v[581]=(v[548]*v[549]+v[551]*v[552]+v[554]*v[555])/2e0;
 gpost[i318-1][0]=(v[549]*v[570]+v[548]*v[571]+v[547]*v[572])/v[562];
 gpost[i318-1][1]=v[573];
 gpost[i318-1][2]=v[574];
 gpost[i318-1][3]=v[573];
 gpost[i318-1][4]=(v[552]*v[575]+v[551]*v[576]+v[550]*v[577])/v[562];
 gpost[i318-1][5]=v[578];
 gpost[i318-1][6]=v[574];
 gpost[i318-1][7]=v[578];
 gpost[i318-1][8]=(v[555]*(v[555]*v[565]+v[569]*v[627])+v[554]*(v[554]*v[565]+v[569]*v[628])
  +v[553]*(v[553]*v[565]+v[569]*v[629]))/v[562];
 gpost[i318-1][9]=(-1e0+(v[547]*v[547])+(v[550]*v[550])+(v[553]*v[553]))/2e0;
 gpost[i318-1][10]=v[579];
 gpost[i318-1][11]=v[580];
 gpost[i318-1][12]=v[579];
 gpost[i318-1][13]=(-1e0+(v[548]*v[548])+(v[551]*v[551])+(v[554]*v[554]))/2e0;
 gpost[i318-1][14]=v[581];
 gpost[i318-1][15]=v[580];
 gpost[i318-1][16]=v[581];
 gpost[i318-1][17]=(-1e0+(v[549]*v[549])+(v[552]*v[552])+(v[555]*v[555]))/2e0;
};/* end for */
};