1if (abs(tt - 1.0) < eps)
13RealT
ww = 1.0 / (1.0 - tt);
15basis[0] = 0.25 * (-rr - ss - 1.0) * ((1.0 - rr) * (1.0 - ss) - tt + rr * ss * tt *
ww);
16basis[1] = 0.25 * (rr - ss - 1.0) * ((1.0 + rr) * (1.0 - ss) - tt - rr * ss * tt *
ww);
17basis[2] = 0.25 * (rr + ss - 1.0) * ((1.0 + rr) * (1.0 + ss) - tt + rr * ss * tt *
ww);
18basis[3] = 0.25 * (-rr + ss - 1.0) * ((1.0 - rr) * (1.0 + ss) - tt - rr * ss * tt *
ww);
20basis[4] = tt * (2.0 * tt - 1.0);
22basis[5] = 0.5 * (1.0 + rr - tt) * (1.0 - rr - tt) * (1.0 - ss - tt) *
ww;
23basis[6] = 0.5 * (1.0 + ss - tt) * (1.0 - ss - tt) * (1.0 + rr - tt) *
ww;
24basis[7] = 0.5 * (1.0 + rr - tt) * (1.0 - rr - tt) * (1.0 + ss - tt) *
ww;
25basis[8] = 0.5 * (1.0 + ss - tt) * (1.0 - ss - tt) * (1.0 - rr - tt) *
ww;
27basis[9] = tt * (1.0 - rr - tt) * (1.0 - ss - tt) *
ww;
28basis[10] = tt * (1.0 + rr - tt) * (1.0 - ss - tt) *
ww;
29basis[11] = tt * (1.0 + rr - tt) * (1.0 + ss - tt) *
ww;
30basis[12] = tt * (1.0 - rr - tt) * (1.0 + ss - tt) *
ww;