Skip to content
Snippets Groups Projects

Update waves

Merged Boman Romain requested to merge adrien into master
2 files
+ 9
6
Compare changes
  • Side-by-side
  • Inline
Files
2
+ 8
5
@@ -8,18 +8,21 @@ GaussTetra4::GaussTetra4(int n) : Gauss()
{
switch (n)
{
case 2:
case 2: {
p.resize(4);
w.resize(4); // 4 pts !
p[0] = Pt(1. / 6, 1. / 6, 1. / 6);
double a = (5. - sqrt(5.)) / 20.;
double b = (5. + 3.*sqrt(5.)) / 20.;
p[0] = Pt(a, a, a);
w[0] = 1. / 4. / 6.;
p[1] = Pt(2. / 3, 1. / 6, 1. / 6);
p[1] = Pt(a, a, b);
w[1] = 1. / 4. / 6.;
p[2] = Pt(1. / 6, 2. / 3, 1. / 6);
p[2] = Pt(a, b, a);
w[2] = 1. / 4. / 6.;
p[3] = Pt(1. / 6, 1. / 6, 2. / 3);
p[3] = Pt(b, a, a);
w[3] = 1. / 4. / 6.;
break;
}
default:
throw std::runtime_error("not implemented (GaussTetra4)!");
}
Loading