62   const Real dy = (y - yc)/disty;  in shape()  local
85       return dy;  in shape()
95       return dx*dy;  in shape()
98       return dy*dy;  in shape()
104       return dz*dy;  in shape()
114       return dx*dx*dy;  in shape()
117       return dx*dy*dy;  in shape()
120       return dy*dy*dy;  in shape()
126       return dx*dy*dz;  in shape()
129       return dy*dy*dz;  in shape()
135       return dy*dz*dz;  in shape()
145       return dx*dx*dx*dy;  in shape()
148       return dx*dx*dy*dy;  in shape()
151       return dx*dy*dy*dy;  in shape()
154       return dy*dy*dy*dy;  in shape()
160       return dx*dx*dy*dz;  in shape()
163       return dx*dy*dy*dz;  in shape()
166       return dy*dy*dy*dz;  in shape()
172       return dx*dy*dz*dz;  in shape()
175       return dy*dy*dz*dz;  in shape()
181       return dy*dz*dz*dz;  in shape()
198         val *= dy;  in shape()
271   const Real dy = (y - yc)/disty;  in shape_deriv()  local
308             return dy/distx;  in shape_deriv()
327             return 2.*dx*dy/distx;  in shape_deriv()
330             return dy*dy/distx;  in shape_deriv()
339             return dy*dz/distx;  in shape_deriv()
358             return 3.*dx*dx*dy/distx;  in shape_deriv()
361             return 2.*dx*dy*dy/distx;  in shape_deriv()
364             return dy*dy*dy/distx;  in shape_deriv()
373             return 2.*dx*dy*dz/distx;  in shape_deriv()
376             return dy*dy*dz/distx;  in shape_deriv()
385             return dy*dz*dz/distx;  in shape_deriv()
411               val *= dy;  in shape_deriv()
446             return 2.*dy/disty;  in shape_deriv()
465             return 2.*dx*dy/disty;  in shape_deriv()
468             return 3.*dy*dy/disty;  in shape_deriv()
477             return 2.*dy*dz/disty;  in shape_deriv()
496             return 2.*dx*dx*dy/disty;  in shape_deriv()
499             return 3.*dx*dy*dy/disty;  in shape_deriv()
502             return 4.*dy*dy*dy/disty;  in shape_deriv()
511             return 2.*dx*dy*dz/disty;  in shape_deriv()
514             return 3.*dy*dy*dz/disty;  in shape_deriv()
523             return 2.*dy*dz*dz/disty;  in shape_deriv()
546               val *= dy;  in shape_deriv()
587             return dy/distz;  in shape_deriv()
609             return dx*dy/distz;  in shape_deriv()
612             return dy*dy/distz;  in shape_deriv()
618             return 2.*dy*dz/distz;  in shape_deriv()
643             return dx*dx*dy/distz;  in shape_deriv()
646             return dx*dy*dy/distz;  in shape_deriv()
649             return dy*dy*dy/distz;  in shape_deriv()
655             return 2.*dx*dy*dz/distz;  in shape_deriv()
658             return 2.*dy*dy*dz/distz;  in shape_deriv()
664             return 3.*dy*dz*dz/distz;  in shape_deriv()
681               val *= dy;  in shape_deriv()
761   const Real dy = (y - yc)/disty;  in shape_second_deriv()  local
811             return 2.*dy/dist2x;  in shape_second_deriv()
832             return 6.*dx*dy/dist2x;  in shape_second_deriv()
835             return 2.*dy*dy/dist2x;  in shape_second_deriv()
845             return 2.*dy*dz/dist2x;  in shape_second_deriv()
873               val *= dy;  in shape_second_deriv()
916             return 2.*dy/distxy;  in shape_second_deriv()
939             return 4.*dx*dy/distxy;  in shape_second_deriv()
942             return 3.*dy*dy/distxy;  in shape_second_deriv()
952             return 2.*dy*dz/distxy;  in shape_second_deriv()
979               val *= dy;  in shape_second_deriv()
1022             return 6.*dy/dist2y;  in shape_second_deriv()
1045             return 6.*dx*dy/dist2y;  in shape_second_deriv()
1048             return 12.*dy*dy/dist2y;  in shape_second_deriv()
1058             return 6.*dy*dz/dist2y;  in shape_second_deriv()
1084               val *= dy;  in shape_second_deriv()
1130             return dy/distxz;  in shape_second_deriv()
1154             return 2.*dx*dy/distxz;  in shape_second_deriv()
1157             return dy*dy/distxz;  in shape_second_deriv()
1166             return 2.*dy*dz/distxz;  in shape_second_deriv()
1190               val *= dy;  in shape_second_deriv()
1236             return 2.*dy/distyz;  in shape_second_deriv()
1260             return 2.*dx*dy/distyz;  in shape_second_deriv()
1263             return 3.*dy*dy/distyz;  in shape_second_deriv()
1272             return 4.*dy*dz/distyz;  in shape_second_deriv()
1295               val *= dy;  in shape_second_deriv()
1342             return 2.*dy/dist2z;  in shape_second_deriv()
1363             return 2.*dx*dy/dist2z;  in shape_second_deriv()
1366             return 2.*dy*dy/dist2z;  in shape_second_deriv()
1372             return 6.*dy*dz/dist2z;  in shape_second_deriv()
1389               val *= dy;  in shape_second_deriv()