UnitCell operator== and operator!=

Framework/Geometry/inc/MantidGeometry/Crystal/UnitCell.h

@@ -161,6 +161,8 @@ public:
   double volume() const;
   double recVolume() const;
   virtual void recalculateFromGstar(const Kernel::Matrix<double> &NewGstar);
+  bool operator==(const UnitCell &other) const;
+  bool operator!=(const UnitCell &other) const;
 
 protected:
   /// Lattice parameter a,b,c,alpha,beta,gamma (in \f$ \mbox{ \AA } \f$ and

Framework/Geometry/src/Crystal/UnitCell.cpp

@@ -880,6 +880,13 @@ void UnitCell::recalculateFromGstar(const DblMatrix &NewGstar) {
   calculateB();
 }
 
+bool UnitCell::operator==(const UnitCell &other) const {
+  return da == other.da; // da error not used in comparison
+}
+bool UnitCell::operator!=(const UnitCell &other) const {
+  return !this->operator==(other);
+}
+
 std::ostream &operator<<(std::ostream &out, const UnitCell &unitCell) {
   // always show the lattice constants
   out << "Lattice Parameters:" << std::fixed << std::setprecision(6)

Framework/Geometry/test/UnitCellTest.h

@@ -170,4 +170,53 @@ public:
     TS_ASSERT_DIFFERS(precisionLimit.c(), precisionLimitOther.c());
     TS_ASSERT_DELTA(precisionLimit.c(), precisionLimitOther.c(), 1e-9);
   }
+
+  void test_equals_when_unitcell_identical() {
+    const UnitCell a(2.0, 4.0, 5.0, 90.0, 100.0, 102.0);
+    const UnitCell b(a);
+    TS_ASSERT_EQUALS(a, b);
+    TS_ASSERT(!(a != b));
+  }
+
+  void test_not_equals_when_unitcell_differs_in_a() {
+    const UnitCell a(1.0, 1.0, 1.0, 90.0, 90.0, 90.0);
+    UnitCell b(a);
+    b.seta(2.0);
+    TS_ASSERT_DIFFERS(a, b);
+  }
+
+  void test_not_equals_when_unitcell_differs_in_b() {
+    const UnitCell a(1.0, 1.0, 1.0, 90.0, 90.0, 90.0);
+    UnitCell b(a);
+    b.setb(2.0);
+    TS_ASSERT_DIFFERS(a, b);
+  }
+
+  void test_not_equals_when_unitcell_differs_in_c() {
+    const UnitCell a(1.0, 1.0, 1.0, 90.0, 90.0, 90.0);
+    UnitCell b(a);
+    b.setc(2.0);
+    TS_ASSERT_DIFFERS(a, b);
+  }
+
+  void test_not_equals_when_unitcell_differs_in_alpha() {
+    const UnitCell a(1.0, 1.0, 1.0, 90.0, 90.0, 90.0);
+    UnitCell b(a);
+    b.setalpha(100);
+    TS_ASSERT_DIFFERS(a, b);
+  }
+
+  void test_not_equals_when_unitcell_differs_in_beta() {
+    const UnitCell a(1.0, 1.0, 1.0, 90.0, 90.0, 90.0);
+    UnitCell b(a);
+    b.setbeta(100);
+    TS_ASSERT_DIFFERS(a, b);
+  }
+
+  void test_not_equals_when_unitcell_differs_in_gamma() {
+    const UnitCell a(1.0, 1.0, 1.0, 90.0, 90.0, 90.0);
+    UnitCell b(a);
+    b.setgamma(100);
+    TS_ASSERT_DIFFERS(a, b);
+  }
 };