Add a `Field::compare` method

This adds a compare method for a instance of Field to compare itself to another one. I originally wanted to overload operator==, but implemented a normal method in order to be able to pass a tolerance. The method should return true only if the the two fields have the same number of blocks and these blocks have the same structure, and each field value are equal up to the absolute tolerance passed to the function.

This does the job for now, but we might want to look into using the boost manipulator boost::test_tools::per_element in the test assertions. That would require the field to be implemented as a container with a forward iterator: beyond my C++ skills at the minute.

tests/init_fields.hpp

@@ -32,6 +32,7 @@ class InitFields
 public:
     Field<double, stateIn> *fixt_in;
     Field<double, stateOut> *fixt_out;
+    Field<double, stateOut> *fixt_expected;
     MultiRegions::ExpListSharedPtr fixt_explist{nullptr};
     ~InitFields()
     {
@@ -67,8 +68,10 @@ public:
         // Create two Field objects with a MemoryRegionCPU backend by default
         auto f_in  = Field<double, stateIn>::create(blocks_in);
         auto f_out = Field<double, stateOut>::create(blocks_out);
+	auto f_expected = Field<double, stateOut>::create(blocks_out);
         fixt_in    = new Field<double, stateIn>(std::move(f_in));
         fixt_out   = new Field<double, stateOut>(std::move(f_out));
+	fixt_expected = new Field<double, stateOut>(std::move(f_expected));
     }
 
 protected: