Esempio n. 1
0
def run(mesh):
    mc = MonteCarlo(mesh, name='test1')
    mc.set_m(random_m)
    J = 50 * const.k_B
    mc.set_options(H=[0, 0, 0.0], J=J, D=0.27 * J, T=4.0)
    mc.run(steps=100000,
           save_m_steps=None,
           save_vtk_steps=50000,
           save_data_steps=50000)
def test_mc_run(Hz=6.0, T=5.0):
    np.random.seed(100)
    #mesh = CuboidMesh(nx=28*2, ny=16*3, nz=1, periodicity = (True, True, False))
    mesh = HexagonalMesh(1, 28, 16*2, periodicity=(True, True))
    mc = MonteCarlo(mesh, name='dyn')
    mc.set_m(random_m)
    mc.set_options(H=[0,0,Hz], J=50.0, D=0.5*50, T=T)
    mc.run(steps=5000, save_m_steps=None, save_vtk_steps=None, save_data_steps=1000)
    skx_number = mc.skyrmion_number()
    assert(skx_number<-3)
    assert(skx_number>-3.5)
Esempio n. 3
0
def test_mc_run(Hz=6.0, T=5.0):
    #This fast test just shows whether mc can be run or not.
    np.random.seed(100)
    #mesh = CuboidMesh(nx=28*2, ny=16*3, nz=1, periodicity = (True, True, False))
    mesh = HexagonalMesh(1, 28, 16 * 2, periodicity=(True, True))
    mc = MonteCarlo(mesh, name='test_mc')
    mc.set_m(random_m)
    J = 50 * const.k_B
    mc.set_options(H=[0, 0, Hz], J=J, D=0.5 * J, T=T)
    mc.run(steps=5000,
           save_m_steps=None,
           save_vtk_steps=None,
           save_data_steps=1000)
    skx_number = mc.skyrmion_number()
    assert (skx_number < -2.5)
    assert (skx_number > -3.5)
Esempio n. 4
0
def run(mesh):
    
    mc = MonteCarlo(mesh, name='test1')
    mc.set_m(random_m)
    mc.set_options(H=[0,0,0.0], J=50.0, D=0.27*50, T=5.0, Kc=50*0.1)
    mc.run(steps=20000, save_m_steps=None, save_vtk_steps=1000, save_data_steps=10)