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