With the following unit tests:

1. HomogeneousMaterial test
2. HomogeneousMagneticMaterial test

Tests/UnitTests/TestCore/HomogeneousMagneticMaterialTest.h

+#ifndef HOMOGENEOUSMAGNETICMATERIALTEST_H
+#define HOMOGENEOUSMAGNETICMATERIALTEST_H
+
+#include "HomogeneousMagneticMaterial.h"
+#include "gtest/gtest.h"
+
+class HomogeneousMagneticMaterialTest : public ::testing::Test
+{
+public:
+    HomogeneousMagneticMaterialTest() {}
+    virtual ~HomogeneousMagneticMaterialTest() {}
+};
+
+TEST_F(HomogeneousMagneticMaterialTest, HomogeneousMagneticMaterialWithRefIndex)
+{
+    complex_t refIndex = complex_t(1.0, 2.0);
+    kvector_t magnetism = kvector_t(3.0, 4.0, 5.0);
+    HomogeneousMagneticMaterial material("MagMaterial", refIndex, magnetism);
+    EXPECT_EQ("MagMaterial", material.getName());
+    EXPECT_EQ(refIndex, material.getRefractiveIndex());
+    EXPECT_EQ(magnetism, material.getMagneticField());
+
+    Eigen::Matrix2cd matrix = material.getScatteringMatrix(-2.91042993836710484e-3);
+    EXPECT_EQ(complex_t(2,4.0), matrix(0,0));
+    EXPECT_EQ(complex_t(3.0,-4.0), matrix(0,1));
+    EXPECT_EQ(complex_t(3.0, 4.0), matrix(1,0));
+    EXPECT_EQ(complex_t(-8,4.0), matrix(1,1));
+
+    complex_t refIndex2 = complex_t(2.0, 3.0);
+    material.setRefractiveIndex(refIndex2);
+    EXPECT_EQ(refIndex2, material.getRefractiveIndex());
+
+    Eigen::Matrix2cd matrix2 = material.getScatteringMatrix(-2.91042993836710484e-3);
+    EXPECT_EQ(complex_t(0.0,12.0), matrix2(0,0));
+    EXPECT_EQ(complex_t(3.0, -4.0), matrix2(0,1));
+    EXPECT_EQ(complex_t(3.0,4.0), matrix2(1,0));
+    EXPECT_EQ(complex_t(-10.0,12.0), matrix2(1,1));
+
+
+    kvector_t magnetism2 = kvector_t(5.0, 6.0, 7.0);
+    material.setMagneticField(magnetism2);
+    EXPECT_EQ(magnetism2, material.getMagneticField());
+
+}
+
+TEST_F(HomogeneousMagneticMaterialTest, HomogeneousMagneticMaterialWithRefIndexAndMagField)
+{
+
+    kvector_t magnetism = kvector_t(3.0, 4.0, 5.0);
+    HomogeneousMagneticMaterial material("MagMaterial", 2.0, 2.0, magnetism);
+    EXPECT_EQ("MagMaterial", material.getName());
+    EXPECT_EQ(complex_t(-1.0, 2.0), material.getRefractiveIndex());
+    EXPECT_EQ(magnetism, material.getMagneticField());
+
+    Eigen::Matrix2cd matrix = material.getScatteringMatrix(-2.91042993836710484e-3);
+    EXPECT_EQ(complex_t(2,-4.0), matrix(0,0));
+    EXPECT_EQ(complex_t(3.0,-4.0), matrix(0,1));
+    EXPECT_EQ(complex_t(3.0, 4.0), matrix(1,0));
+    EXPECT_EQ(complex_t(-8,-4.0), matrix(1,1));
+
+    complex_t refIndex2 = complex_t(2.0, 3.0);
+    material.setRefractiveIndex(refIndex2);
+    EXPECT_EQ(refIndex2, material.getRefractiveIndex());
+
+    kvector_t magnetism2 = kvector_t(5.0, 6.0, 7.0);
+    material.setMagneticField(magnetism2);
+    EXPECT_EQ(magnetism2, material.getMagneticField());
+
+}
+
+TEST_F(HomogeneousMagneticMaterialTest, HomogeneousMagneticMaterialTransform)
+{
+    complex_t refIndex = complex_t(0.0, 0.0);
+    kvector_t magnetism = kvector_t(0.0, 0.0, 0.0);
+    HomogeneousMagneticMaterial material("MagMaterial", refIndex, magnetism);
+
+    Geometry::Transform3D transform =
+            Geometry::Transform3D::createRotateZ(45.*Units::degree);
+    const IMaterial * tMaterial = material.createTransformedMaterial(transform);
+
+    EXPECT_EQ("MagMaterial", tMaterial->getName());
+    EXPECT_EQ(refIndex, tMaterial->getRefractiveIndex());
+
+    Eigen::Matrix2cd matrix = tMaterial->getScatteringMatrix(-2.91042993836710484e-3);
+    EXPECT_EQ(complex_t(0.0, 0.0), matrix(0,0));
+    EXPECT_EQ(complex_t(0.0, 0.0), matrix(0,1));
+    EXPECT_EQ(complex_t(0.0, 0.0), matrix(1,0));
+    EXPECT_EQ(complex_t(0.0, 0.0), matrix(1,1));
+
+
+    delete tMaterial;
+}
+
+TEST_F(HomogeneousMagneticMaterialTest, HomogeneousMagneticMaterialClone)
+{
+    complex_t refIndex = complex_t(1.0, 2.0);
+    kvector_t magnetism = kvector_t(3.0, 4.0, 5.0);
+    HomogeneousMagneticMaterial material("MagMaterial", refIndex, magnetism);
+
+    HomogeneousMagneticMaterial * clone = material.clone();
+
+    EXPECT_EQ("MagMaterial", clone->getName());
+    EXPECT_EQ(refIndex, clone->getRefractiveIndex());
+    EXPECT_EQ(magnetism, clone->getMagneticField());
+
+    Eigen::Matrix2cd matrix = clone->getScatteringMatrix(-2.91042993836710484e-3);
+    EXPECT_EQ(complex_t(2,4.0), matrix(0,0));
+    EXPECT_EQ(complex_t(3.0,-4.0), matrix(0,1));
+    EXPECT_EQ(complex_t(3.0, 4.0), matrix(1,0));
+    EXPECT_EQ(complex_t(-8,4.0), matrix(1,1));
+
+    complex_t refIndex2 = complex_t(2.0, 3.0);
+    clone->setRefractiveIndex(refIndex2);
+    EXPECT_EQ(refIndex2, clone->getRefractiveIndex());
+
+    kvector_t magnetism2 = kvector_t(5.0, 6.0, 7.0);
+    clone->setMagneticField(magnetism2);
+    EXPECT_EQ(magnetism2, clone->getMagneticField());
+
+
+    Geometry::Transform3D transform =
+            Geometry::Transform3D::createRotateZ(45.*Units::degree);
+    const IMaterial * tMaterial = clone->createTransformedMaterial(transform);
+
+    EXPECT_EQ("MagMaterial", tMaterial->getName());
+    EXPECT_EQ(refIndex2, tMaterial->getRefractiveIndex());
+
+
+    delete tMaterial;
+    delete clone;
+}
+
+
+#endif //HOMOGENEOUSMAGNETICMATERIALTEST_H

Tests/UnitTests/TestCore/HomogeneousMaterialTest.h

+#ifndef HOMOGENEOUSMATERIALTEST_H
+#define HOMOGENEOUSMATERIALTEST_H
+
+#include "HomogeneousMaterial.h"
+#include "gtest/gtest.h"
+
+class HomogeneousMaterialTest : public ::testing :: Test
+{
+public:
+    HomogeneousMaterialTest() {}
+    virtual ~HomogeneousMaterialTest(){}
+};
+
+TEST_F(HomogeneousMaterialTest, HomogeneousMaterialWithRefIndex)
+{
+    complex_t refIndex = complex_t(1.0, 2.0);
+    HomogeneousMaterial material("Material1", refIndex);
+    EXPECT_EQ("Material1", material.getName());
+    EXPECT_EQ(refIndex, material.getRefractiveIndex());
+
+    Eigen::Matrix2cd matrix = material.getScatteringMatrix(3.0);
+    EXPECT_EQ(complex_t(-3.0,4.0), matrix(0,0));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix(0,1));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix(1,0));
+    EXPECT_EQ(complex_t(-3.0,4.0), matrix(1,1));
+
+    complex_t refIndex2 = complex_t(2.0, 3.0);
+    material.setRefractiveIndex(refIndex2);
+    EXPECT_EQ(refIndex2, material.getRefractiveIndex());
+
+    Eigen::Matrix2cd matrix2 = material.getScatteringMatrix(4.0);
+    EXPECT_EQ(complex_t(-5.0,12.0), matrix2(0,0));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix2(0,1));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix2(1,0));
+    EXPECT_EQ(complex_t(-5.0,12.0), matrix2(1,1));
+
+}
+
+TEST_F(HomogeneousMaterialTest, HomogeneousMaterialWithRefIndexParam)
+{
+    HomogeneousMaterial material("Material1", 2.0, 2.0);
+    EXPECT_EQ("Material1", material.getName());
+    EXPECT_EQ(complex_t(-1.0,2.0), material.getRefractiveIndex());
+
+    Eigen::Matrix2cd matrix = material.getScatteringMatrix(3.0);
+    EXPECT_EQ(complex_t(-3.0,-4.0), matrix(0,0));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix(0,1));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix(1,0));
+    EXPECT_EQ(complex_t(-3.0,-4.0), matrix(1,1));
+
+}
+
+TEST_F(HomogeneousMaterialTest, HomogeneousMaterialTransform)
+{
+    complex_t refIndex = complex_t(1.0, 2.0);
+    HomogeneousMaterial material("Material1", refIndex);
+
+    Geometry::Transform3D transform =
+            Geometry::Transform3D::createRotateZ(45.*Units::degree);
+    const IMaterial * tMaterial = material.createTransformedMaterial(transform);
+
+    EXPECT_EQ("Material1", tMaterial->getName());
+    EXPECT_EQ(refIndex, tMaterial->getRefractiveIndex());
+
+    Eigen::Matrix2cd matrix = tMaterial->getScatteringMatrix(3.0);
+    EXPECT_EQ(complex_t(-3.0,4.0), matrix(0,0));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix(0,1));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix(1,0));
+    EXPECT_EQ(complex_t(-3.0,4.0), matrix(1,1));
+
+    delete tMaterial;
+}
+
+TEST_F(HomogeneousMaterialTest, HomogeneousMaterialClone)
+{
+    complex_t refIndex = complex_t(1.0, 2.0);
+    HomogeneousMaterial material("Material1", refIndex);
+
+    HomogeneousMaterial * clone = material.clone();
+
+    EXPECT_EQ("Material1", clone->getName());
+    EXPECT_EQ(refIndex, clone->getRefractiveIndex());
+
+    Eigen::Matrix2cd matrix = clone->getScatteringMatrix(3.0);
+    EXPECT_EQ(complex_t(-3.0,4.0), matrix(0,0));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix(0,1));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix(1,0));
+    EXPECT_EQ(complex_t(-3.0,4.0), matrix(1,1));
+
+
+    complex_t refIndex2 = complex_t(2.0, 3.0);
+    clone->setRefractiveIndex(refIndex2);
+    EXPECT_EQ(refIndex2, clone->getRefractiveIndex());
+
+    Eigen::Matrix2cd matrix2 = clone->getScatteringMatrix(4.0);
+    EXPECT_EQ(complex_t(-5.0,12.0), matrix2(0,0));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix2(0,1));
+    EXPECT_EQ(complex_t(0.0,0.0), matrix2(1,0));
+    EXPECT_EQ(complex_t(-5.0,12.0), matrix2(1,1));
+
+
+
+    Geometry::Transform3D transform =
+            Geometry::Transform3D::createRotateZ(45.*Units::degree);
+    const IMaterial * tMaterial = clone->createTransformedMaterial(transform);
+
+    EXPECT_EQ("Material1", tMaterial->getName());
+    EXPECT_EQ(refIndex2, tMaterial->getRefractiveIndex());
+
+
+    delete tMaterial;
+    delete clone;
+
+}
+
+#endif //HOMOGENEOUSMATERIALTEST_H

Tests/UnitTests/TestCore/main.cpp

@@ -29,6 +29,8 @@
 #include "FTDistributionsTest.h"
 #include "MultiLayerTest.h"
 #include "ParticleDecorationTest.h"
+#include "HomogeneousMaterialTest.h"
+#include "HomogeneousMagneticMaterialTest.h"
 
 struct ErrorStreamRedirect {
     ErrorStreamRedirect( std::streambuf * new_buffer )