def test_init_invalid_args(self): for K, u1, u2 in self.invalid_args: with pytest.raises((TypeError, ValueError)): term = mm.CubicAnisotropy(K=K, u1=u1, u2=u2) with pytest.raises(AttributeError): term = mm.CubicAnisotropy(wrong=1)
def test_repr(self): for K1, u1, u2 in self.valid_args: anisotropy = mm.CubicAnisotropy(K1=K1, u1=u1, u2=u2) exp_str = ("CubicAnisotropy(K1={}, u1={}, u2={}, " "name=\"{}\")").format(K1, u1, u2, "cubicanisotropy") assert repr(anisotropy) == exp_str anisotropy = mm.CubicAnisotropy(1000, (0, 1, 0), (0, 0, 1), name="test_name") assert repr(anisotropy) == ("CubicAnisotropy(K1=1000, u1=(0, 1, 0), " "u2=(0, 0, 1), name=\"test_name\")")
def test_repr(self): for K1, u1, u2 in self.valid_args: anisotropy = mm.CubicAnisotropy(K1=K1, u1=u1, u2=u2) exp_str = ('CubicAnisotropy(K1={}, u1={}, u2={}, ' 'name=\'{}\')').format(K1, u1, u2, 'cubicanisotropy') assert repr(anisotropy) == exp_str anisotropy = mm.CubicAnisotropy(1000, (0, 1, 0), (0, 0, 1), name='test_name') assert repr(anisotropy) == ('CubicAnisotropy(K1=1000, u1=(0, 1, 0), ' 'u2=(0, 0, 1), name=\'test_name\')')
def setup(self): A = 1e-12 self.exchange = mm.Exchange(A=A) H = (0, 0, 1.2e6) self.zeeman = mm.Zeeman(H=H) K1 = 1e4 K2 = 3e2 u = (0, 1, 0) self.uniaxialanisotropy = mm.UniaxialAnisotropy(K1=K1, K2=K2, u=u) self.demag = mm.Demag() D = 1e-3 crystalclass = 't' self.dmi = mm.DMI(D=D, crystalclass=crystalclass) K1 = 5e6 u1 = (0, 0, 1) u2 = (0, 1, 0) self.cubicanisotropy = mm.CubicAnisotropy(K1=K1, u1=u1, u2=u2) self.terms = [ self.exchange, self.zeeman, self.uniaxialanisotropy, self.demag, self.dmi, self.cubicanisotropy ] self.invalid_terms = [ 1, 2.5, 0, 'abc', [3, 7e-12], [self.exchange, self.zeeman] ]
def test_exchange_cubicanisotropy(self): name = 'exchange_cubicanisotropy' A = {'r1': 1e-12, 'r2': 0} K = 1e5 u1 = (1, 0, 0) u2 = (0, 1, 0) Ms = 1e6 mesh = df.Mesh(region=self.region, cell=self.cell, subregions=self.subregions) system = mm.System(name=name) system.energy = mm.Exchange(A=A) + \ mm.CubicAnisotropy(K=K, u1=u1, u2=u2) system.m = df.Field(mesh, dim=3, value=(1, 0.3, 0), norm=Ms) md = self.calculator.MinDriver() md.drive(system) value = system.m(mesh.region.random_point()) assert np.linalg.norm(np.subtract(value, (Ms, 0, 0))) < 1e-3 self.calculator.delete(system)
def test_field_vector_vector(self): name = 'cubicanisotropy_field_vector_vector' mesh = df.Mesh(region=self.region, cell=self.cell) def K_fun(pos): x, y, z = pos if x <= 0: return 0 else: return 1e5 K = df.Field(mesh, dim=1, value=K_fun) u1 = (0, 0, 1) u2 = (0, 1, 0) Ms = 1e6 system = mm.System(name=name) system.energy = mm.CubicAnisotropy(K=K, u1=u1, u2=u2) system.m = df.Field(mesh, dim=3, value=(0, 0.3, 1), norm=Ms) md = self.calculator.MinDriver() md.drive(system) value = system.m((-2e-9, 1e-9, 1e-9)) assert np.linalg.norm(np.cross(value, (0, 0.3 * Ms, Ms))) < 1e-3 value = system.m((2e-9, 2e-9, 2e-9)) assert np.linalg.norm(np.subtract(value, (0, 0, Ms))) < 1e-3 self.calculator.delete(system)
def test_scalar_vector_vector(self): name = 'cubicanisotropy_scalar_vector_vector' mesh = df.Mesh(region=self.region, cell=self.cell) K = 1e5 u1 = (0, 0, 1) u2 = (0, 1, 0) Ms = 1e6 system = mm.System(name=name) system.energy = mm.CubicAnisotropy(K=K, u1=u1, u2=u2) def m_fun(pos): x, y, z = pos if x <= 0: return (0, 0.2, 1) else: return (0, 1, 0.2) system.m = df.Field(mesh, dim=3, value=m_fun, norm=Ms) md = self.calculator.MinDriver() md.drive(system) value = system.m((-1e-9, 2e-9, 2e-9)) assert np.linalg.norm(np.subtract(value, (0, 0, Ms))) < 1e-3 value = system.m((1e-9, 2e-9, 2e-9)) assert np.linalg.norm(np.subtract(value, (0, Ms, 0))) < 1e-3 self.calculator.delete(system)
def test_dict_vector_vector(self): name = 'cubicanisotropy_dict_vector_vector' mesh = df.Mesh(region=self.region, cell=self.cell, subregions=self.subregions) K = {'r1': 0, 'r2': 1e5} u1 = (0, 0, 1) u2 = (0, 1, 0) Ms = 1e6 system = mm.System(name=name) system.energy = mm.CubicAnisotropy(K=K, u1=u1, u2=u2) system.m = df.Field(mesh, dim=3, value=(0, 0.3, 1), norm=Ms) md = self.calculator.MinDriver() md.drive(system) value = system.m((-2e-9, 1e-9, 1e-9)) assert np.linalg.norm(np.cross(value, (0, 0.3 * Ms, Ms))) < 1e-3 value = system.m((2e-9, 2e-9, 2e-9)) assert np.linalg.norm(np.subtract(value, (0, 0, Ms))) < 1e-3 self.calculator.delete(system)
def test_init_valid_args(self): for K, u1, u2 in self.valid_args: term = mm.CubicAnisotropy(K=K, u1=u1, u2=u2) check_term(term) assert hasattr(term, 'K') assert hasattr(term, 'u1') assert hasattr(term, 'u2') assert term.name == 'cubicanisotropy' assert re.search(r'^CubicAnisotropy\(K=.+, u1=.+\, u2=.+\)$', repr(term))
def test_init_valid_args(self): for K1, u1, u2 in self.valid_args: anisotropy = mm.CubicAnisotropy(K1=K1, u1=u1, u2=u2) assert anisotropy.K1 == K1 assert isinstance(anisotropy.K1, numbers.Real) assert isinstance(anisotropy.u1, (tuple, list, np.ndarray)) assert isinstance(anisotropy.u2, (tuple, list, np.ndarray)) assert len(anisotropy.u1) == 3 assert len(anisotropy.u2) == 3 assert all([isinstance(i, numbers.Real) for i in anisotropy.u1]) assert all([isinstance(i, numbers.Real) for i in anisotropy.u2])
def setup(self): name = 'compute_tests' p1 = (0, 0, 0) p2 = (10e-9, 2e-9, 2e-9) cell = (2e-9, 2e-9, 2e-9) region = df.Region(p1=p1, p2=p2) mesh = df.Mesh(region=region, cell=cell) self.system = mm.System(name=name) self.system.energy = ( mm.Exchange(A=1e-12) + mm.Demag() + mm.Zeeman(H=(8e6, 0, 0)) + mm.UniaxialAnisotropy(K=1e4, u=(0, 0, 1)) + mm.CubicAnisotropy(K=1e3, u1=(1, 0, 0), u2=(0, 1, 0))) self.system.m = df.Field(mesh, dim=3, value=(0, 0, 1), norm=8e6)
def test_repr_latex(self): for K1, u1, u2 in self.valid_args: anisotropy = mm.CubicAnisotropy(K1=K1, u1=u1, u2=u2) latex = anisotropy._repr_latex_() # Assert some characteristics of LaTeX string. assert isinstance(latex, str) assert latex[0] == latex[-1] == '$' assert 'K_{1}' in latex assert '\mathbf{u}_{1}' in latex assert '\mathbf{u}_{2}' in latex assert '\mathbf{m}' in latex assert '^{2}' in latex assert '\cdot' in latex assert latex.count("\mathbf{u}_{1}") == 2 assert latex.count("\mathbf{u}_{2}") == 2 assert latex.count("\mathbf{u}_{3}") == 2
def setup(self): self.exchange = mm.Exchange(A=1e-12) self.zeeman = mm.Zeeman(H=(0, 0, 1.2e6)) self.uniaxialanisotropy = mm.UniaxialAnisotropy(K=1e4, u=(0, 1, 0)) self.demag = mm.Demag() self.dmi = mm.DMI(D=1e-3, crystalclass='T') self.cubicanisotropy = mm.CubicAnisotropy(K={ 'r1': 1e6, 'r2': 5e6 }, u1=(0, 0, 1), u2=(0, 1, 0)) self.terms = [ self.exchange, self.zeeman, self.uniaxialanisotropy, self.demag, self.dmi, self.cubicanisotropy ] self.invalid_terms = [1, 2.5, 0, 'abc', [3, 7e-12], [self.exchange, 2]]
def test_script(self): for K1, u1, u2 in self.valid_args: anisotropy = mm.CubicAnisotropy(K1=K1, u1=u1, u2=u2) with pytest.raises(NotImplementedError): script = anisotropy._script
def test_name(self): for K1, u1, u2 in self.valid_args: anisotropy = mm.CubicAnisotropy(K1=K1, u1=u1, u2=u2) assert anisotropy.name == 'cubicanisotropy'
def test_init_invalid_args(self): for K1, u1, u2 in self.invalid_args: with pytest.raises(Exception): anisotropy = mm.CubicAnisotropy(K1=K1, u1=u1, u2=u2)