178     double tmp0 =1/(1-xi[2]); // xi3->xi[2]  in Pyr_eB()  local
184         psi =0.125*(xi[2]*xi[2]-1)*(1+sign[k]*(2*xi[k]*tmp0))*(1+sign[m]*(2*xi[m]*tmp0));  in Pyr_eB()
187         double tmp1 =2*xi[k]*tmp0;  in Pyr_eB()
188         psi =0.125*(tmp1*tmp1-1)*(1+sign[m]*(2*xi[m]*tmp0))*(1-xi[2]);  in Pyr_eB()
196     double tmp0 =1/(1-xi[2]); // xi3->xi[2]  in dPeBdxi()  local
203             dPsidxi =-0.25*(1+xi[2])*sign[k]*(1+sign[m]*(2*xi[m]*tmp0));  in dPeBdxi()
205             dPsidxi =-0.25*(1+xi[2])*sign[m]*(1+sign[k]*(2*xi[k]*tmp0));  in dPeBdxi()
207             dPsidxi =0.25*(xi[2]*(1+sign[k]*2*xi[k]*tmp0)*(1+sign[m]*2*xi[m]*tmp0)-  in dPeBdxi()
208                            (1+xi[2])*tmp0*(sign[k]*xi[k]*(1+sign[m]*2*xi[m]*tmp0)+  in dPeBdxi()
209                                            sign[m]*xi[m]*(1+sign[k]*2*xi[k]*tmp0)));  in dPeBdxi()
211         double tmp1 =2*xi[k]*tmp0;  in dPeBdxi()
213             dPsidxi =0.5*tmp1*(1+sign[m]*(2*xi[m]*tmp0));  in dPeBdxi()
217             dPsidxi =0.25*(tmp1*tmp1)*(1+sign[m]*(2*xi[m]*tmp0))-0.125*(tmp1*tmp1-1);  in dPeBdxi()
226     double tmp0 =1/(1-xi[1]); // xi2->xi[1]  in Pyr_fB()  local
231         double tmp1 =2*xi[0]*tmp0;  in Pyr_fB()
232         double tmp2 =2*xi[2]*tmp0;  in Pyr_fB()
238         double tmp1 =2*xi[m]*tmp0;  in Pyr_fB()
240         psi =0.125*(1+sign[k]*2*xi[k]*tmp0)*(1-tmp1*tmp1)*(1-xi[1]*xi[1]);  in Pyr_fB()
248     double tmp0 =1/(1-xi[2]); // xi3->xi[2]  in dPfBdxi()  local
252         double tmp1 =2*xi[0]*tmp0;  in dPfBdxi()
253         double tmp2 =2*xi[2]*tmp0;  in dPfBdxi()
269             double tmp1 =2*xi[m]*tmp0;  in dPfBdxi()
273             dPsidxi =-(1+sign[k]*2*xi[k]*tmp0)*(xi[m]*tmp0)*(1+xi[1]);  in dPfBdxi()
275             double tmp1 =xi[k]*tmp0;  in dPfBdxi()
276             double tmp2 =xi[m]*tmp0;  in dPfBdxi()