Home
last modified time | relevance | path

Searched full:degree (Results 1 – 25 of 86) sorted by relevance

1234

/libCEED/examples/petsc/
H A DREADME.md11 To run, `./bpsraw -ceed [ceed-resource] -problem bp[1-6] -degree [degree]`
23 To run, `./bps -ceed [ceed-resource] -problem bp[1-6] -degree [degree]`
35 mpiexec -n 64 ./bps -problem bp1,bp2,bp3,bp4 -degree 2,3,5,7 \
48 To run, `./multigrid -ceed [ceed-resource] -problem bp[1-6] -degree [degree]`
61 - `-degree` - Polynomial degree of tensor product basis
72 To run, `./area -problem cube -ceed [ceed-resource] -petscspace_degree [degree]`
76 `./area -problem sphere -ceed [ceed-resource] -petscspace_degree [degree]`
84 - `-petscspace_degree` - Polynomial degree of tensor product basis
H A Darea.c22 // ./area -problem cube -degree 3 -dm_refine 2
23 // ./area -problem sphere -degree 3 -dm_refine 2
27 // mpiexec -n 4 ./area -problem cube -degree 3 -dm_refine 2
28 // mpiexec -n 4 ./area -problem sphere -degree 3 -dm_refine 2
33 //TESTARGS -ceed {ceed_resource} -test -degree 3 -dm_refine 1
66 degree = 3; // default degree for finite element bases in main() local
89 …PetscCall(PetscOptionsInt("-degree", "Polynomial degree of tensor product basis", NULL, degree, &d… in main()
111 PetscCall(SetupDMByDegree(dm, degree, q_extra, num_comp_u, topo_dim, false)); in main()
148 PetscInt P = degree + 1, Q = P + q_extra; in main()
168 …PetscCall(SetupLibceedByDegree(dm, ceed, degree, topo_dim, q_extra, num_comp_x, num_comp_u, g_size… in main()
H A Dbpssphere.c21 // bpssphere -problem bp1 -degree 3
22 // bpssphere -problem bp2 -degree 3
23 // bpssphere -problem bp3 -degree 3
24 // bpssphere -problem bp4 -degree 3
25 // bpssphere -problem bp5 -degree 3 -ceed /cpu/self
26 // bpssphere -problem bp6 -degree 3 -ceed /gpu/cuda
28 //TESTARGS -ceed {ceed_resource} -test -problem bp3 -degree 3 -dm_refine 2
53 …PetscInt degree = 3, q_extra, l_size, g_size, topo_dim = 2, num_comp_x = 3, num_comp_u… in main() local
84 degree = test_mode ? 3 : 2; in main()
85 …PetscCall(PetscOptionsInt("-degree", "Polynomial degree of tensor product basis", NULL, degree, &d… in main()
[all …]
H A Dbps.c20 // ./bps -problem bp1 -degree 3
21 // ./bps -problem bp2 -degree 3
22 // ./bps -problem bp3 -degree 3
23 // ./bps -problem bp4 -degree 3
24 // ./bps -problem bp5 -degree 3 -ceed /cpu/self
25 // ./bps -problem bp6 -degree 3 -ceed /gpu/cuda
27 //TESTARGS(name="BP3, tet elements") -ceed {ceed_resource} -test -problem bp3 -degree 3 -ksp_max_it…
28 //TESTARGS(name="BP5, hex elements") -ceed {ceed_resource} -test -problem bp5 -degree 3 -ksp_max_it…
29 //TESTARGS(name="BP1+3, hex elements") -ceed {ceed_resource} -test -problem bp1_3 -degree 3 -ksp_ma…
30 //TESTARGS(name="BP2+4, hex elements") -ceed {ceed_resource} -test -problem bp2_4 -degree 3 -ksp_ma…
[all …]
H A Dmultigrid.c25 //TESTARGS(name="BP3, hex elements") -ceed {ceed_resource} -test -problem bp3 -degree 3
26 //TESTARGS(name="BP3, tet elements") -ceed {ceed_resource} -test -problem bp3 -degree 3 -simplex
52 …PetscInt degree = 3, q_extra, *l_size, *xl_size, *g_size, dim = 3, fine_level, mesh_elem[3] = {3, … in main() local
95 degree = test_mode ? 3 : 2; in main()
96 …PetscCall(PetscOptionsInt("-degree", "Polynomial degree of tensor product basis", NULL, degree, &d… in main()
97 …if (degree < 1) SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "-degree %" PetscInt_FMT " mus… in main()
149 num_levels = degree; in main()
152 num_levels = ceil(log(degree) / log(2)) + 1; in main()
164 level_degrees[fine_level] = degree; in main()
223 PetscInt P = degree + 1, Q = P + q_extra; in main()
H A Dbpsswarm.c21 // bpssphere -problem bp1 -degree 3
22 // bpssphere -problem bp2 -degree 3
23 // bpssphere -problem bp3 -degree 3
56 …PetscInt degree, q_extra, l_size, g_size, dim = 3, num_comp_u = 1, xl_size, num_points… in main() local
101 degree = 2; in main()
102 …PetscCall(PetscOptionsInt("-degree", "Polynomial degree of tensor product basis", NULL, degree, &d… in main()
189 …PetscCall(SetupDMByDegree(dm_mesh, degree, q_extra, num_comp_u, dim, bp_options[bp_choice].enforce… in main()
238 PetscInt P = degree + 1, Q = P + q_extra; in main()
H A Dbpsraw.c27 //TESTARGS -ceed {ceed_resource} -test -problem bp2 -degree 5 -q_extra 1 -ksp_max_it_clip 15,15
59 static void GlobalNodes(const PetscInt p[3], const PetscInt i_rank[3], PetscInt degree, const P… in GlobalNodes() argument
60 for (int d = 0; d < 3; d++) m_nodes[d] = degree * mesh_elem[d] + (i_rank[d] == p[d] - 1); in GlobalNodes()
62 static PetscInt GlobalStart(const PetscInt p[3], const PetscInt i_rank[3], PetscInt degree, const P… in GlobalStart() argument
70 GlobalNodes(p, ijk_rank, degree, mesh_elem, m_nodes); in GlobalStart()
342 …PetscInt degree, q_extra, local_nodes, local_elem, mesh_elem[3], m_nodes[3], p[3], i_rank[3], l_no… in main() local
378 degree = test_mode ? 3 : 1; in main()
379 …PetscCall(PetscOptionsInt("-degree", "Polynomial degree of tensor product basis", NULL, degree, &d… in main()
391 P = degree + 1; in main()
419 for (local_elem = PetscMax(1, local_nodes / (degree * degree * degree));; local_elem++) { in main()
[all …]
/libCEED/examples/python/
H A Dex_common.py25 mesh_degree: Mesh polynomial degree
26 solution_degree: Solution polynomial degree
38 parser.add_argument("-m", "--mesh-degree", type=int, default=4,
39 help="Mesh polynomial degree (default: 4)")
40 parser.add_argument("-p", "--solution-degree", type=int, default=4,
41 help="Solution polynomial degree (default: 4)")
59 def get_cartesian_mesh_size(dim, degree, prob_size): argument
64 degree: Polynomial degree
71 num_elem = prob_size // (degree ** dim)
91 def build_cartesian_restriction(ceed, dim, num_xyz, degree, num_comp, num_q_comp, num_qpts, create_… argument
[all …]
/libCEED/examples/solids/
H A DREADME.md20 ./elasticity -mesh [.exo file] -degree [degree] -nu [nu] -E [E] [boundary options] -problem [proble…
37 * - `-degree [int]`
38 - Polynomial degree of the finite element basis
65 ./elasticity -mesh [.exo file] -degree 4 -E 1e6 -nu 0.3 -bc_clamp 998,999 -bc_clamp_998_translate 0…
158 ./elasticity -mesh [mesh] -degree [degree] -nu [nu] -E [E] -forcing mms
170 The default smoother uses degree 3 Chebyshev with Jacobi preconditioning.
171 …(Lower degree is often faster, albeit less robust; try {code}`-outer_mg_levels_ksp_max_it 2`, for …
172 …Application of the linear operators for all levels with degree $p > 1$ is performed matrix-free us…
/libCEED/examples/solids/src/
H A Dcl-options.c33 app_ctx->degree = 3; in ProcessCommandLineOptions()
34 …PetscCall(PetscOptionsInt("-degree", "Polynomial degree of tensor product basis", NULL, app_ctx->d… in ProcessCommandLineOptions()
164 app_ctx->num_levels = ceil(log(app_ctx->degree) / log(2)) + 1; in ProcessCommandLineOptions()
167 app_ctx->num_levels = app_ctx->degree; in ProcessCommandLineOptions()
180 app_ctx->level_degrees[app_ctx->num_levels - 1] = app_ctx->degree; in ProcessCommandLineOptions()
186 app_ctx->level_degrees[0] = app_ctx->degree; in ProcessCommandLineOptions()
/libCEED/benchmarks/
H A Dpostprocess_plot.py93 sol_p_set = sel_runs['degree'].drop_duplicates()
97 qpts = sel_runs['quadrature_pts'].loc[pl_runs['degree'] == sol_p]
100 print('Degree: %i, quadrature points:' % sol_p, qpts[0])
102 d = [[run['degree'], run['num_elem'], 1. * run['num_unknowns'] / num_nodes / vdim,
105 pl_runs.loc[(pl_runs['degree'] == sol_p) |
126 d = [[run['degree'], run['num_elem'], 1. * run['num_unknowns'] / num_nodes / vdim,
129 pl_runs.loc[(pl_runs['degree'] == sol_p) |
H A Dpetsc-bpsraw.sh13 # -degree <1>: Polynomial degree of tensor product basis
31 local all_args=("${common_args[@]}" -degree $sol_p -local $loc_dofs -problem $bp)
H A Dpetsc-bps.sh12 # -degree <1>: Polynomial degree of tensor product basis
30 local all_args=("${common_args[@]}" -degree $sol_p -local_nodes $loc_nodes -problem $bp)
/libCEED/examples/ceed/
H A Dex3-volume.c47 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim…
48 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, …
58 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main()
59 CeedInt sol_degree = 4; // polynomial degree for the solution in main()
102 printf(" Mesh degree [-m] : %" CeedInt_FMT "\n", mesh_degree); in main()
103 printf(" Solution degree [-p] : %" CeedInt_FMT "\n", sol_degree); in main()
283 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim]) { in GetCartesianMeshSize() argument
285 // prob_size ~ num_elem * degree^dim in GetCartesianMeshSize()
286 CeedInt num_elem = prob_size / CeedIntPow(degree, dim); in GetCartesianMeshSize()
307 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt… in BuildCartesianRestriction() argument
[all …]
H A Dex1-volume.c48 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim…
49 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, …
59 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main()
60 CeedInt sol_degree = 4; // polynomial degree for the solution in main()
105 printf(" Mesh degree [-m] : %" CeedInt_FMT "\n", mesh_degree); in main()
106 printf(" Solution degree [-p] : %" CeedInt_FMT "\n", sol_degree); in main()
294 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim]) { in GetCartesianMeshSize() argument
296 // prob_size ~ num_elem * degree^dim in GetCartesianMeshSize()
297 CeedInt num_elem = prob_size / CeedIntPow(degree, dim); in GetCartesianMeshSize()
318 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt… in BuildCartesianRestriction() argument
[all …]
H A Dex2-surface.c48 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[3]);
49 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[3], CeedInt degree, Ce…
59 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main()
60 CeedInt sol_degree = 4; // polynomial degree for the solution in main()
109 printf(" Mesh degree [-m] : %" CeedInt_FMT "\n", mesh_degree); in main()
110 printf(" Solution degree [-p] : %" CeedInt_FMT "\n", sol_degree); in main()
314 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[3]) { in GetCartesianMeshSize() argument
316 // prob_size ~ num_elem * degree^dim in GetCartesianMeshSize()
317 CeedInt num_elem = prob_size / CeedIntPow(degree, dim); in GetCartesianMeshSize()
338 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[3], CeedInt degree, CeedInt n… in BuildCartesianRestriction() argument
[all …]
H A Dex1-volume-f.f9042 subroutine getcartesianmeshsize(fe_dim, degree, prob_size, num_xyz) argument
45 integer degree local
51 num_elem = prob_size/(degree**fe_dim)
55 ! prob_size ~ num_elem * degree^dim
75 subroutine buildcartesianrestriction(ceed, fe_dim, num_xyz, degree, num_comp, mesh_size, num_qpts, … argument
83 integer degree local
104 p = degree + 1
356 write (*, *) ' Mesh degree [-m] : ', mesh_degree
357 write (*, *) ' Solution degree [-p] : ', sol_degree
/libCEED/doc/papers/joss/
H A Dpaper.md71 …when the polynomial degree $p$ of the basis functions is increased, resulting in $O(p^d)$ storage …
114 Some constructors are provided for arbitrary polynomial degree $H^1$ Lagrange bases with a tensor-p…
118 …n, topology, or basis degree (see \autoref{fig:schematic}), in contrast to systems like FEniCS whe…
119 …t size and $p$ the polynomial degree of the basis functions; see @babuska1994hpfem) and $p$-multig…
127 (in this case representing a Laplacian) on element resolution, topology, and basis degree.\label{fi…
151 …sen system with NVIDIA V100 GPUs. Each curve represents fixing the basis degree $p$ and varying th…
/libCEED/examples/rust-qfunctions/
H A Dex1-volume.c49 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim…
50 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, …
60 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main()
61 CeedInt sol_degree = 4; // polynomial degree for the solution in main()
106 printf(" Mesh degree [-m] : %" CeedInt_FMT "\n", mesh_degree); in main()
107 printf(" Solution degree [-p] : %" CeedInt_FMT "\n", sol_degree); in main()
304 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim]) { in GetCartesianMeshSize() argument
306 // prob_size ~ num_elem * degree^dim in GetCartesianMeshSize()
307 CeedInt num_elem = prob_size / CeedIntPow(degree, dim); in GetCartesianMeshSize()
328 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt… in BuildCartesianRestriction() argument
[all …]
/libCEED/examples/fluids/
H A DstdoutParsing.py23 …r".*(?:^Degree of FEM Space: (\d+)).*(?:^Global FEM nodes: (\d{2,})).*(?:^dm_plex_box_faces: (\S+)…
39 Values of "dofs", "time", "error", "degree", and "box_faces"'
46 values["degree"] = match[1]
H A Dconv_test.sh21 run_flags[degree]=2
41 echo ",mesh_res,degree,rel_error" > $file_name
45 run_flags[degree]=$d
H A Dblasius.yaml14 degree: 1
20 # degree: 2
/libCEED/examples/fluids/src/
H A Ddm_utils.c402 @brief Setup `DM` with FE space of appropriate degree
408 @param[in] degree Polynomial orders of field
417 …ByOrderBegin_FEM(PetscBool setup_faces, PetscBool setup_coords, PetscInt degree, PetscInt coord_or… in DMSetupByOrderBegin_FEM() argument
419 PetscInt dim, q_order = degree + q_extra; in DMSetupByOrderBegin_FEM()
431 PetscInt q_order = degree + q_extra; in DMSetupByOrderBegin_FEM()
433 PetscCall(PetscFECreateLagrange(comm, dim, field_sizes[i], is_simplex, degree, q_order, &fe)); in DMSetupByOrderBegin_FEM()
495 @brief Setup `DM` with FE space of appropriate degree with no boundary conditions
501 @param[in] degree Polynomial orders of field
510 PetscErrorCode DMSetupByOrder_FEM(PetscBool setup_faces, PetscBool setup_coords, PetscInt degree, P… in DMSetupByOrder_FEM() argument
513 …PetscCall(DMSetupByOrderBegin_FEM(setup_faces, setup_coords, degree, coord_order, q_extra, num_fie… in DMSetupByOrder_FEM()
H A Dcloptions.c102 app_ctx->degree = 1; in ProcessCommandLineOptions()
103 …PetscCall(PetscOptionsInt("-degree", "Polynomial degree of finite elements", NULL, app_ctx->degree… in ProcessCommandLineOptions()
/libCEED/julia/LibCEED.jl/src/
H A DQuadrature.jl4 Return the Gauss-Legendre quadrature rule with `q` points (integrates polynomials of degree
23 Return the Gauss-Lobatto quadrature rule with `q` points (integrates polynomials of degree

1234