MultiLatticeSearch: Works

2026-06-05 15:05:11 +02:00
parent b49c7a766d
commit 3522568d49
2 changed files with 7 additions and 9 deletions
@@ -44,8 +44,6 @@ TEST_CASE("MultiLatticeSearch_Empty") {
    CHECK(result.empty());
 }

-#include <iostream>
-
 TEST_CASE("MultiLatticeSearch_EP") {
    // Real EP case
    CrystalLattice cell1(Coord(-13.2, -30.0, -29.6),
@@ -58,8 +56,12 @@ TEST_CASE("MultiLatticeSearch_EP") {

    auto result = MultiLatticeSearch({cell1, cell2}, 0.1, 3);

-
    REQUIRE(result.size() == 1);
-    std::cout << result[0].rotation_vector << std::endl;
-    std::cout << result[0].rotation_vector.Length() * 180.0 / M_PI << std::endl;
+
+    RotMatrix matrix(result[0].rotation_vector.Length(), result[0].rotation_vector.Normalize());
+    auto cell2_rot = cell1.Multiply(matrix);
+
+    CHECK((cell2_rot.Vec0() - cell2.Vec0()).Length() < cell2.Vec0().Length() / 100);
+    CHECK((cell2_rot.Vec1() - cell2.Vec1()).Length() < cell2.Vec1().Length() / 100);
+    CHECK((cell2_rot.Vec2() - cell2.Vec2()).Length() < cell2.Vec2().Length() / 100);
 }