33 typedef typename LB::Traits::RangeFieldType Scalar;
35 for (
size_t i=0; i<5; i++)
36 sign[i] = (s[i]) ? -1.0 : 1.0;
38 Scalar r = 1/std::sqrt(2);
40 n[0] = { 0.0, 0.0, -1.0};
41 n[1] = {-1.0, 0.0, 0.0};
43 n[3] = { 0.0, -1.0, 0.0};
48 c[2] = 1/2.0 * std::sqrt(2);
50 c[4] = 1/2.0 * std::sqrt(2);
52 m[0] = { 0.5, 0.5, 0.0};
53 m[1] = { 0.0, 1/3.0, 1/3.0};
54 m[2] = { 2/3.0, 1/3.0, 1/3.0};
55 m[3] = { 1/3.0, 0.0, 1/3.0};
56 m[4] = { 1/3.0, 2/3.0, 1/3.0};