Home
last modified time | relevance | path

Searched refs:mesh (Results 1 – 25 of 40) sorted by relevance

12

/libCEED/examples/mfem/
H A Dbp1.cpp82 mfem::Mesh *mesh = new mfem::Mesh(mesh_file, 1, 1); in main() local
83 int dim = mesh->Dimension(); in main()
89 … int ref_levels = (int)floor((log(max_nnodes / mesh->GetNE()) - dim * log(order)) / log(2.) / dim); in main()
91 mesh->UniformRefinement(); in main()
94 if (mesh->GetNodalFESpace() == NULL) { in main()
95 mesh->SetCurvature(1, false, -1, mfem::Ordering::byNODES); in main()
97 if (mesh->NURBSext) { in main()
98 mesh->SetCurvature(order, false, -1, mfem::Ordering::byNODES); in main()
105 mfem::FiniteElementSpace *fespace = new mfem::FiniteElementSpace(mesh, fec); in main()
150 sol_sock << "solution\n" << *mesh << sol << std::flush; in main()
[all …]
H A Dbp3.cpp104 mfem::Mesh *mesh = new mfem::Mesh(mesh_file, 1, 1); in main() local
105 int dim = mesh->Dimension(); in main()
111 … int ref_levels = (int)floor((log(max_nnodes / mesh->GetNE()) - dim * log(order)) / log(2.) / dim); in main()
113 mesh->UniformRefinement(); in main()
116 if (mesh->GetNodalFESpace() == NULL) { in main()
117 mesh->SetCurvature(1, false, -1, mfem::Ordering::byNODES); in main()
119 if (mesh->NURBSext) { in main()
120 mesh->SetCurvature(order, false, -1, mfem::Ordering::byNODES); in main()
127 mfem::FiniteElementSpace *fespace = new mfem::FiniteElementSpace(mesh, fec); in main()
135 if (mesh->bdr_attributes.Size()) { in main()
[all …]
H A Dbp1.hpp33 mfem::Mesh *mesh = fes->GetMesh(); in FESpace2Ceed() local
37 switch (mesh->Dimension()) { in FESpace2Ceed()
77 …CeedBasisCreateTensorH1(ceed, mesh->Dimension(), fes->GetVDim(), order + 1, ir.GetNPoints(), shape… in FESpace2Ceed()
82 for (int i = 0; i < mesh->GetNE(); i++) { in FESpace2Ceed()
88 …CeedElemRestrictionCreate(ceed, mesh->GetNE(), fe->GetDof(), fes->GetVDim(), fes->GetNDofs(), (fes… in FESpace2Ceed()
95 mfem::Mesh *mesh = fes->GetMesh(); in CeedMassOperator() local
99 …CeedInt num_elem = mesh->GetNE(), dim = mesh->SpaceDimension(), ncompx = dim,… in CeedMassOperator()
103 const mfem::FiniteElementSpace *mesh_fes = mesh->GetNodalFESpace(); in CeedMassOperator()
111 CeedVectorCreate(ceed, mesh->GetNodes()->Size(), &node_coords); in CeedMassOperator()
112 CeedVectorSetArray(node_coords, CEED_MEM_HOST, CEED_USE_POINTER, mesh->GetNodes()->GetData()); in CeedMassOperator()
[all …]
H A Dbp3.hpp34 mfem::Mesh *mesh = fes->GetMesh(); in FESpace2Ceed() local
38 switch (mesh->Dimension()) { in FESpace2Ceed()
78 …CeedBasisCreateTensorH1(ceed, mesh->Dimension(), fes->GetVDim(), order + 1, ir.GetNPoints(), shape… in FESpace2Ceed()
83 for (int i = 0; i < mesh->GetNE(); i++) { in FESpace2Ceed()
89 …CeedElemRestrictionCreate(ceed, mesh->GetNE(), fe->GetDof(), fes->GetVDim(), fes->GetNDofs(), (fes… in FESpace2Ceed()
96 mfem::Mesh *mesh = fes->GetMesh(); in CeedDiffusionOperator() local
100 …CeedInt num_elem = mesh->GetNE(), dim = mesh->SpaceDimension(), ncompx = dim,… in CeedDiffusionOperator()
104 const mfem::FiniteElementSpace *mesh_fes = mesh->GetNodalFESpace(); in CeedDiffusionOperator()
112 CeedVectorCreate(ceed, mesh->GetNodes()->Size(), &node_coords); in CeedDiffusionOperator()
113 CeedVectorSetArray(node_coords, CEED_MEM_HOST, CEED_USE_POINTER, mesh->GetNodes()->GetData()); in CeedDiffusionOperator()
[all …]
/libCEED/julia/LibCEED.jl/examples/
H A Dex1-volume.jl9 # map [0,1] to [0,1] varying the mesh density
41 # Determine the mesh size based on the given approximate problem size.
45 # Build CeedElemRestriction objects describing the mesh and solution discrete
58 println("Number of mesh nodes : ", div(mesh_size, dim))
61 # Create a CeedVector with the mesh coordinates.
64 # Apply a transformation to the mesh.
131 # Compute the mesh volume using the mass operator: vol = 1^T \cdot M \cdot 1
132 print("Computing the mesh volume using the formula: vol = 1^T.M.1 ...")
141 # Compute and print the sum of the entries of 'v' giving the mesh volume.
145 @printf("Exact mesh volume : % .14g\n", exact_vol)
[all …]
H A Dex3-volume.jl9 # map [0,1] to [0,1] varying the mesh density
41 # Determine the mesh size based on the given approximate problem size.
45 # Build CeedElemRestriction objects describing the mesh and solution discrete
68 println("Number of mesh nodes : ", div(mesh_size, dim))
71 # Create a CeedVector with the mesh coordinates.
74 # Apply a transformation to the mesh.
170 # Compute the mesh volume using the mass+diffusion operator: vol = 1^T \cdot (M + K) \cdot 1
171 print("Computing the mesh volume using the formula: vol = 1^T * (M + K) * 1...")
180 # Compute and print the sum of the entries of 'v' giving the mesh volume.
184 @printf("Exact mesh volume : % .14g\n", exact_vol)
[all …]
H A Dex2-surface.jl8 # map [0,1] to [0,1] varying the mesh density
31 # Build CeedElemRestriction objects describing the mesh and solution discrete
60 println("Number of mesh nodes : ", div(mesh_size, dim))
63 # Create a CeedVector with the mesh coordinates.
67 # Apply a transformation to the mesh.
137 # Compute the mesh surface area using the diff operator:
139 print("Computing the mesh surface area using the formula: sa = 1^T.|K.x| ...")
158 @printf("Exact mesh surface area : % .14g\n", exact_sa)
159 @printf("Computed mesh surface area : % .14g\n", sa)
H A Dex1-volume-c.jl192 # map [0,1] to [0,1] varying the mesh density
245 # Determine the mesh size based on the given approximate problem size.
249 # Build CeedElemRestriction objects describing the mesh and solution discrete
262 println("Number of mesh nodes : ", div(mesh_size, dim))
265 # Create a C.CeedVector with the mesh coordinates.
269 # Apply a transformation to the mesh.
383 # Compute the mesh volume using the mass operator: vol = 1^T \cdot M \cdot 1
384 print("Computing the mesh volume using the formula: vol = 1^T.M.1 ...")
398 # Compute and print the sum of the entries of 'v' giving the mesh volume.
406 @printf("Exact mesh volume : % .14g\n", exact_vol)
[all …]
/libCEED/examples/petsc/
H A DREADME.md5 ### CEED bakeoff problems with raw mesh management - bpsraw
27 - `-mesh` - Read mesh from file
52 - `-mesh` - Read mesh from file
87 - `-mesh` - Read mesh from file
/libCEED/examples/solids/
H A DREADME.md20 ./elasticity -mesh [.exo file] -degree [degree] -nu [nu] -E [E] [boundary options] -problem [proble…
35 * - `-mesh [filename]`
36 - Path to mesh file in any format supported by PETSc.
52 This solver can use any mesh format that PETSc's `DMPlex` can read (Exodus, Gmsh, Med, etc.).
54 Note that many mesh formats require PETSc to be configured appropriately; e.g., `--download-exodusi…
57 Consider the specific example of the mesh seen below:
65 ./elasticity -mesh [.exo file] -degree 4 -E 1e6 -nu 0.3 -bc_clamp 998,999 -bc_clamp_998_translate 0…
70 As an alternative to specifying a mesh with {code}`-mesh`, the user may use a DMPlex box mesh by sp…
72 …le exploiting {code}`-dm_plex_box_faces`, we consider a {code}`4 x 4 x 4` mesh where essential (Dr…
80 If the coordinates for a particular side of a mesh are zero along the axis of rotation, it may appe…
[all …]
/libCEED/examples/rust/ex3-volume/
H A DCargo.toml12 mesh = { path = "../mesh" }
/libCEED/examples/rust/ex1-volume/
H A DCargo.toml12 mesh = { path = "../mesh" }
/libCEED/examples/rust/ex2-surface/
H A DCargo.toml12 mesh = { path = "../mesh" }
/libCEED/examples/rust/ex3-volume-vector/
H A DCargo.toml12 mesh = { path = "../mesh" }
/libCEED/examples/rust/ex1-volume-vector/
H A DCargo.toml12 mesh = { path = "../mesh" }
/libCEED/examples/rust/ex2-surface-vector/
H A DCargo.toml12 mesh = { path = "../mesh" }
/libCEED/examples/python/
H A Dtutorial-6-shell.ipynb43 "application of the mass operator. Arbitrary mesh and solution orders in 1D, 2D and 3D\n",
189 "Number of mesh nodes : 274625\n",
192 "Computing the mesh volume using the formula: vol = 1^T.M.1 ... done.\n",
193 "Exact mesh volume : 2.3561944901923\n",
194 "Computed mesh volume : 2.3561944901921\n",
220 "Arbitrary mesh and solution orders in 1D, 2D and 3D are supported from the same code.\n",
280 "Number of mesh nodes : 274625\n",
283 "Computing the mesh surface area using the formula: sa = 1^T.|K.x| ... done.\n",
284 "Exact mesh surface area : 6\n",
285 "Computed mesh surface area : 5.9773703490853\n",
/libCEED/examples/rust/ex1-volume-vector/src/
H A Dmain.rs103 let num_xyz = mesh::cartesian_mesh_size(dim, solution_degree, problem_size); in example_1_vector()
118 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_1_vector()
120 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, 1, num_qpts)?; in example_1_vector()
122 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, ncomp_u, num_qpts)?; in example_1_vector()
131 let mut mesh_coords = mesh::cartesian_mesh_coords(&ceed, dim, num_xyz, mesh_degree, mesh_size)?; in example_1_vector()
/libCEED/examples/nek/boxes/
H A Db1e.rea32 0.00000 p028 TORDER: mesh velocity (0: p28=p27)
51 0.00000 p047 vnu: mesh matieral prop.
133 F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)
/libCEED/examples/rust/ex1-volume/src/
H A Dmain.rs101 let num_xyz = mesh::cartesian_mesh_size(dim, solution_degree, problem_size); in example_1()
116 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_1()
118 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, 1, num_qpts)?; in example_1()
127 let mut mesh_coords = mesh::cartesian_mesh_coords(&ceed, dim, num_xyz, mesh_degree, mesh_size)?; in example_1()
/libCEED/examples/rust/mesh/
H A DCargo.toml2 name = "mesh"
/libCEED/examples/rust/ex3-volume/src/
H A Dmain.rs97 let num_xyz = mesh::cartesian_mesh_size(dim, solution_degree, problem_size); in example_3()
112 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_3()
113 let (_, rstr_qdata) = mesh::build_cartesian_restriction( in example_3()
122 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, 1, num_qpts)?; in example_3()
131 let mut mesh_coords = mesh::cartesian_mesh_coords(&ceed, dim, num_xyz, mesh_degree, mesh_size)?; in example_3()
/libCEED/examples/rust/ex2-surface/src/
H A Dmain.rs102 let num_xyz = mesh::cartesian_mesh_size(dim, solution_degree, problem_size); in example_2()
117 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_2()
118 let (_, rstr_qdata) = mesh::build_cartesian_restriction( in example_2()
127 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, 1, num_qpts)?; in example_2()
136 let mut mesh_coords = mesh::cartesian_mesh_coords(&ceed, dim, num_xyz, mesh_degree, mesh_size)?; in example_2()
/libCEED/examples/rust/ex3-volume-vector/src/
H A Dmain.rs99 let num_xyz = mesh::cartesian_mesh_size(dim, solution_degree, problem_size); in example_3_vector()
114 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, mesh_degree, ncomp_x, num_qpts)?; in example_3_vector()
115 let (_, rstr_qdata) = mesh::build_cartesian_restriction( in example_3_vector()
124 mesh::build_cartesian_restriction(&ceed, dim, num_xyz, solution_degree, ncomp_u, num_qpts)?; in example_3_vector()
133 let mut mesh_coords = mesh::cartesian_mesh_coords(&ceed, dim, num_xyz, mesh_degree, mesh_size)?; in example_3_vector()
/libCEED/doc/sphinx/source/api/
H A DCeedElemRestriction.rst8 Expressing element decomposition and degrees of freedom over a mesh

12