def load_mesh(file_name = ""): """read(file_name:string). Function to load a mesh from the file file_name. """ raw_mesh = ocaml.mesh_readfile(file_name) # read mesh mesh_info = ocaml.mesh_plotinfo(raw_mesh) # extract data from the mesh print "read mesh: ok" return mesh_info
def load_mesh(self, mesh, file_name = ""): """load_mesh(file_name:string). Function to load a mesh from the file file_name. """ new_mesh = mesh_module.mesh() new_mesh.raw_mesh = ocaml.mesh_readfile(file_name) # read mesh print "read mesh: ok" return new_mesh
# # (C) 2006 Dr. Thomas Fischbacher # Checking bulk and surface contribs for a homogeneously magnetized octahedron. import nfem import sys,math,time import ocaml #need this as long as we use ocaml.probe_field mesh = ocaml.mesh_readfile("./debug-octa.mesh") print "MESH: ",mesh elem_m=ocaml.make_element("m_X",[3],3,1) elem_H_demag=ocaml.make_element("H_demag",[3],3,1) mwe_m=ocaml.make_mwe("mwe_m",mesh,[(1,elem_m)]) mwe_h=ocaml.make_mwe("mwe_h",mesh,[(1,elem_H_demag)]) def initial_m(dof,pos): # radially outward, zero at center ix=dof[1][0] return pos[ix] field_m=ocaml.raw_make_field(mwe_m,[initial_m],"") make_field_h=ocaml.ddd_demag_fun_3d("<S||d/dxj m_X(j)>, j:3",mwe_h,mwe_m) field_h=make_field_h(field_m) print field_h
# # (C) 2006 Dr. Thomas Fischbacher # Checking bulk and surface contribs for a homogeneously magnetized octahedron. import nfem import sys, math, time import ocaml #need this as long as we use ocaml.probe_field mesh = ocaml.mesh_readfile("./debug-octa.mesh") print "MESH: ", mesh elem_m = ocaml.make_element("m_X", [3], 3, 1) elem_H_demag = ocaml.make_element("H_demag", [3], 3, 1) mwe_m = ocaml.make_mwe("mwe_m", mesh, [(1, elem_m)]) mwe_h = ocaml.make_mwe("mwe_h", mesh, [(1, elem_H_demag)]) def initial_m(dof, pos): # radially outward, zero at center ix = dof[1][0] return pos[ix] field_m = ocaml.raw_make_field(mwe_m, [initial_m], "") make_field_h = ocaml.ddd_demag_fun_3d("<S||d/dxj m_X(j)>, j:3", mwe_h, mwe_m) field_h = make_field_h(field_m) print field_h