Lines Matching refs:xi
115 double Hex_eB(double xi[3], int sign2, int sign3) in Hex_eB()
120 return 0.25*(1+sign2*xi[1])*(1+sign3*xi[2]); in Hex_eB()
123 double dHEBdxi1(double xi[3], int sign2, int sign3) in dHEBdxi1()
131 double dHEBdxi2(double xi[3], int sign2, int sign3) in dHEBdxi2()
138 return 0.25*sign2*(1+sign3*xi[2]); in dHEBdxi2()
141 double dHEBdxi3(double xi[3], int sign2, int sign3) in dHEBdxi3()
145 return 0.25*sign3*(1+sign2*xi[1]); in dHEBdxi3()
151 double Hex_fB(double xi[3], int sign3) in Hex_fB()
155 return 0.5*(1+sign3*xi[2]); in Hex_fB()
159 double dHFBdxi1(double xi[3], int sign3) in dHFBdxi1()
163 double dHFBdxi2(double xi[3], int sign3) in dHFBdxi2()
167 double dHFBdxi3(double xi[3], int sign3) in dHFBdxi3()
175 double Pyr_eB(double xi[3], int sign[3], int k, int m, int along) in Pyr_eB()
178 double tmp0 =1/(1-xi[2]); // xi3->xi[2] in Pyr_eB()
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()
193 double dPeBdxi(double xi[3], int sign[3], int k, int m, int along, int byWhich) in dPeBdxi()
196 double tmp0 =1/(1-xi[2]); // xi3->xi[2] in dPeBdxi()
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()
224 double Pyr_fB (double xi[3], int sign[3], int k, int m, int faceType) in Pyr_fB()
226 double tmp0 =1/(1-xi[1]); // xi2->xi[1] in Pyr_fB()
231 double tmp1 =2*xi[0]*tmp0; in Pyr_fB()
232 double tmp2 =2*xi[2]*tmp0; in Pyr_fB()
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()
240 psi =0.125*(1+sign[k]*2*xi[k]*tmp0)*(1-tmp1*tmp1)*(1-xi[1]*xi[1]); in Pyr_fB()
245 double dPfBdxi(double xi[3], int sign[3], int k, int m, int faceType, int byWhich) in dPfBdxi()
248 double tmp0 =1/(1-xi[2]); // xi3->xi[2] in dPfBdxi()
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()
270 dPsidxi =0.25*sign[k]*(1-tmp1*tmp1)*(1+xi[1]); 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()
277 dPsidxi =(1+xi[1])*(0.25*sign[k]*tmp1*(1-4*tmp2*tmp2)- in dPfBdxi()
279 0.25*(1+2*sign[k]*tmp1)*(1-4*tmp2*tmp2)*xi[1]; in dPfBdxi()
323 int quad_face(double xi[3], int gOrd[3], int p, double* entfn, double** edrv) in quad_face()
330 double xi1 = xi[0]; in quad_face()
331 double xi2 = xi[1]; in quad_face()
377 int hex_regn(double xi[3],int p, double* entfn, double** edrv) in hex_regn()
382 double xi1 = xi[0]; in hex_regn()
383 double xi2 = xi[1]; in hex_regn()
384 double xi3 = xi[2]; in hex_regn()
444 double xi = par[0]; in HexShapeAndDrv() local
448 double xim=1-xi; in HexShapeAndDrv()
452 double xip=1+xi; in HexShapeAndDrv()
919 double xi = par[0]; in PyrShapeAndDrv() local
932 double xipp=1+xi*zetamap; in PyrShapeAndDrv()
935 double ximp=1-xi*zetamap; in PyrShapeAndDrv()
950 dN[0][2] = 0.125 * (xi*eta*zetamap*zetamap-1); in PyrShapeAndDrv()
954 dN[1][2] =-0.125 * (xi*eta*zetamap*zetamap+1); in PyrShapeAndDrv()