def calc(axis1_vec, axis2_vec, M_vec):
axis1 = VectorField(mesh); axis1.fill(axis1_vec)
axis2 = VectorField(mesh); axis2.fill(axis2_vec)
M = VectorField(mesh); M.fill(M_vec)
H = VectorField(mesh); H.fill((0,0,0))
E = magneto.cubic_anisotropy(axis1, axis2, k, Ms, M, H) * mesh.cell_volume
return H.average(), E
def calc(axis_vec, M_vec):
axis = VectorField(mesh)
axis.fill(axis_vec)
M = VectorField(mesh)
M.fill(M_vec)
H = VectorField(mesh)
H.fill((0, 0, 0))
E = magneto.uniaxial_anisotropy(axis, k, Ms, M,
H) * mesh.cell_volume
return H.average(), E
def calc(axis1_vec, axis2_vec, M_vec):
axis1 = VectorField(mesh)
axis1.fill(axis1_vec)
axis2 = VectorField(mesh)
axis2.fill(axis2_vec)
M = VectorField(mesh)
M.fill(M_vec)
H = VectorField(mesh)
H.fill((0, 0, 0))
E = magneto.cubic_anisotropy(axis1, axis2, k, Ms, M,
H) * mesh.cell_volume
return H.average(), E
def calc(axis_vec, M_vec):
axis = VectorField(mesh); axis.fill(axis_vec)
M = VectorField(mesh); M.fill(M_vec)
H = VectorField(mesh); H.fill((0, 0, 0))
E = magneto.uniaxial_anisotropy(axis, k, Ms, M, H) * mesh.cell_volume
return H.average(), E
