28 for (std::size_t i=0; i<sign_.size(); i++)
29 sign_[i] = (s[i]) ? -1.0 : 1.0;
31 m_[0] = {1/3.0, 1/3.0, 0.0};
32 m_[1] = {1/3.0, 0.0, 1/3.0};
33 m_[2] = { 0.0, 1/3.0, 1/3.0};
34 m_[3] = {1/3.0, 1/3.0, 1/3.0};
35 n_[0] = { 0.0, 0.0, -1.0};
36 n_[1] = { 0.0, -1.0, 0.0};
37 n_[2] = { -1.0, 0.0, 0.0};
38 n_[3] = {1.0/sqrt(3.0), 1.0/sqrt(3.0), 1.0/sqrt(3.0)};
42 c_[3] = sqrt(2.0)/sqrt(3.0);