187         double tmp1 =2*xi[k]*tmp0;  in Pyr_eB()  local
188         psi =0.125*(tmp1*tmp1-1)*(1+sign[m]*(2*xi[m]*tmp0))*(1-xi[2]);  in Pyr_eB()
211         double tmp1 =2*xi[k]*tmp0;  in dPeBdxi()  local
213             dPsidxi =0.5*tmp1*(1+sign[m]*(2*xi[m]*tmp0));  in dPeBdxi()
215             dPsidxi =0.25*(tmp1*tmp1-1)*sign[m];  in dPeBdxi()
217             dPsidxi =0.25*(tmp1*tmp1)*(1+sign[m]*(2*xi[m]*tmp0))-0.125*(tmp1*tmp1-1);  in dPeBdxi()
231         double tmp1 =2*xi[0]*tmp0;  in Pyr_fB()  local
234         psi =0.125*(1-tmp1*tmp1)*(1-tmp2*tmp2)*(1-xi[1]);  in Pyr_fB()
238         double tmp1 =2*xi[m]*tmp0;  in Pyr_fB()  local
240         psi =0.125*(1+sign[k]*2*xi[k]*tmp0)*(1-tmp1*tmp1)*(1-xi[1]*xi[1]);  in Pyr_fB()
252         double tmp1 =2*xi[0]*tmp0;  in dPfBdxi()  local
257             dPsidxi =-0.5*tmp1*(1-tmp2*tmp2);  in dPfBdxi()
260             dPsidxi = -0.125* (1+tmp1*tmp1+tmp2*tmp2-3*tmp1*tmp1*tmp2*tmp2);  in dPfBdxi()
263             dPsidxi =-0.5*(1-tmp1*tmp1)*tmp2;  in dPfBdxi()
269             double tmp1 =2*xi[m]*tmp0;  in dPfBdxi()  local
270             dPsidxi =0.25*sign[k]*(1-tmp1*tmp1)*(1+xi[1]);  in dPfBdxi()
275             double tmp1 =xi[k]*tmp0;  in dPfBdxi()  local
277             dPsidxi =(1+xi[1])*(0.25*sign[k]*tmp1*(1-4*tmp2*tmp2)-  in dPfBdxi()
278                                 (1+2*sign[k]*tmp1)*tmp2*tmp2)-  in dPfBdxi()
279                 0.25*(1+2*sign[k]*tmp1)*(1-4*tmp2*tmp2)*xi[1];  in dPfBdxi()
423     int tmp1[4];  in HexShapeAndDrv()  local
544                     tmp1[3]=z;  in HexShapeAndDrv()
545                 tmp1[z]=edge[e][1][z];  in HexShapeAndDrv()
548             if((arg2[0] = tmp1[3]) == 0) { arg2[1]=1;arg2[2]=2 ;}  in HexShapeAndDrv()
549             else if(tmp1[3]==1) { arg2[1]=2;arg2[2]=0;}  in HexShapeAndDrv()
560             EdgeBlend = Hex_eB(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()
562             if( tmp1[3] == 0){  in HexShapeAndDrv()
563                 dEBdxi    = dHEBdxi1(arg,tmp1[a],tmp1[b]);   in HexShapeAndDrv()
564                 dEBdeta   = dHEBdxi2(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()
565                 dEBdzeta  = dHEBdxi3(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()
566             } else if( tmp1[3] == 1 ){  in HexShapeAndDrv()
567                 dEBdeta   = dHEBdxi1(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()
568                 dEBdzeta  = dHEBdxi2(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()
569                 dEBdxi    = dHEBdxi3(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()
571                 dEBdzeta  = dHEBdxi1(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()
572                 dEBdxi    = dHEBdxi2(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()
573                 dEBdeta   = dHEBdxi3(arg,tmp1[a],tmp1[b]);  in HexShapeAndDrv()