Fix translational periodic BCs for general planes

This commit is contained in:
Paul Romano 2021-04-19 07:37:00 -05:00
parent f5a0caa51a
commit 261557e596
2 changed files with 4 additions and 0 deletions

View file

@ -1168,6 +1168,8 @@ void read_surfaces(pugi::xml_node node)
// Compute the dot product of the surface normals
Direction norm1 = surf1.normal({0, 0, 0});
Direction norm2 = surf2.normal({0, 0, 0});
norm1 /= norm1.norm();
norm2 /= norm2.norm();
double dot_prod = norm1.dot(norm2);
// If the dot product is 1 (to within floating point precision) then the

View file

@ -0,0 +1,2 @@
k-combined:
2.276564E+00 7.905769E-04