12   double L[4],Le[2],Lf[3],mode,blend,bdrv[3],mdrv[3],epdrv[3][2],mfdrv[2];  in TetShapeAndDrv()  local
16   L[0]=par[0];  in TetShapeAndDrv()
17   L[1]=par[1];  in TetShapeAndDrv()
18   L[2]=par[2];  in TetShapeAndDrv()
19   L[3]=1.0e0-par[0]-par[1]-par[2];  in TetShapeAndDrv()
22     N[i] = L[i];  in TetShapeAndDrv()
38       Le[0]=L[TetEMAP[i][0]];  in TetShapeAndDrv()
39       Le[1]=L[TetEMAP[i][1]];  in TetShapeAndDrv()
40       blend = R_edgeBlendTet(TetEMAP[i],L);  in TetShapeAndDrv()
41       R_edgeBlendTetDrv(TetEMAP[i],L,bdrv);  in TetShapeAndDrv()
60       Lf[0]=L[TetFMAP[i][0]];  in TetShapeAndDrv()
61       Lf[1]=L[TetFMAP[i][1]];  in TetShapeAndDrv()
62       Lf[2]=L[TetFMAP[i][2]];  in TetShapeAndDrv()
64       R_faceBlendTetDrv(TetFMAP[i],L,bdrv);  in TetShapeAndDrv()
87       blend = L[0]*L[1]*L[2]*L[3];  in TetShapeAndDrv()
88       tmp = 1.0e0-L[0]-L[1]-L[2] ;  in TetShapeAndDrv()
89       rsw = L[0]*L[1]*tmp;  in TetShapeAndDrv()
90       stw = L[1]*L[2]*tmp ;  in TetShapeAndDrv()
91       rtw = L[0]*L[2]*tmp ;  in TetShapeAndDrv()
92       rst = L[0]*L[1]*L[2] ;  in TetShapeAndDrv()
93       rstw = rst*L[3] ;  in TetShapeAndDrv()
95 	mode = R_modeShapeTet(ip,is,L);  in TetShapeAndDrv()
97 	R_modeShapeTetDrv(ip, is, L, mdrv);  in TetShapeAndDrv()
114   double L[3];  in TriShapeAndDrv()  local
118   L[0]=par[0];  in TriShapeAndDrv()
119   L[1]=par[1];  in TriShapeAndDrv()
120   L[2]=1.0e0-par[0]-par[1];  in TriShapeAndDrv()
126     N[i] = L[i];  in TriShapeAndDrv()
141     N[3] = -2.0e0*L[0]*L[1];  in TriShapeAndDrv()
142     N[4] = -2.0e0*L[1]*L[2];  in TriShapeAndDrv()
143     N[5] = -2.0e0*L[0]*L[2];  in TriShapeAndDrv()
144     dN[3][0] = -2.0e0*L[1];  in TriShapeAndDrv()
145     dN[3][1] = -2.0e0*L[0];  in TriShapeAndDrv()
146     dN[4][0] = 2.0e0*L[1];  in TriShapeAndDrv()
147     dN[4][1] = -2.0e0+2.0e0*L[0]+4.0e0*L[1];  in TriShapeAndDrv()
148     dN[5][0] = -2.0e0+4.0e0*L[0]+2.0e0*L[1];  in TriShapeAndDrv()
149     dN[5][1] = 2.0e0*L[0];  in TriShapeAndDrv()