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, -1.0, 0.0};
41 n[1] = {-1.0, 0.0, 0.0};
43 n[3] = { 0.0, 0.0, -1.0};
44 n[4] = { 0.0, 0.0, 1.0};
52 m[0] = { 0.5, 0.0, 0.5};
53 m[1] = { 0.0, 0.5, 0.5};
54 m[2] = { 0.5, 0.5, 0.5};
55 m[3] = { 1/3.0, 1/3.0, 0.0};
56 m[4] = { 1/3.0, 1/3.0, 1.0};