| /libCEED/tests/ |
| H A D | t232-elemrestriction.c | 10 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
11 CeedInt ind[(num_elem + 1) + num_points]; in main()
27 point_index = (point_index + 1) % num_points; in main()
33 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_U… in main()
39 CeedScalar array[num_points]; in main()
46 point_index = (point_index + 1) % num_points; in main()
68 point_index = (point_index + 1) % num_points; in main()
|
| H A D | t595-operator.c | 14 …* (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem = 4, num_points = num_elem * num… in main() local
26 CeedVectorCreate(ceed, dim * num_points, &x_points); in main()
28 CeedScalar x_array[dim * num_points]; in main()
41 CeedInt ind_x[num_elem + 1 + num_points]; in main()
44 for (CeedInt i = 0; i < num_points; i++) ind_x[num_elem + 1 + i] = i; in main()
45 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim, num_points * dim, CEED_MEM_HOST… in main()
47 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
51 CeedVectorCreate(ceed, num_points, &q_data); in main()
|
| H A D | t230-elemrestriction.c | 8 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
9 CeedInt ind[(num_elem + 1) + num_points]; in main()
24 point_index = (point_index + 1) % num_points; in main()
30 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_U… in main()
|
| H A D | t233-elemrestriction.c | 9 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
10 CeedInt ind[(num_elem + 1) + num_points]; in main()
26 point_index = (point_index + 1) % num_points; in main()
32 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_U… in main()
37 CeedScalar array[num_points]; in main()
44 point_index = (point_index + 1) % num_points; in main()
|
| H A D | t231-elemrestriction.c | 10 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
11 CeedInt ind[(num_elem + 1) + num_points]; in main()
27 point_index = (point_index + 1) % num_points; in main()
33 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_U… in main()
39 CeedScalar array[num_points]; in main()
46 point_index = (point_index + 1) % num_points; in main()
|
| H A D | t234-elemrestriction.c | 10 CeedInt num_elem = 3, num_points = num_elem * 2; in main() local
11 CeedInt ind[(num_elem + 1) + num_points]; in main()
27 point_index = (point_index + 1) % num_points; in main()
33 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
69 point_index = (point_index + 1) % num_points; in main()
|
| H A D | t365-basis.c | 34 …const CeedInt p = 9, q = 9, num_points = 4, x_dim = CeedIntPow(2, dim), p_dim = CeedIntPow(p, dim); in main() local
39 CeedVectorCreate(ceed, num_points * dim, &x_points); in main()
41 CeedVectorCreate(ceed, num_points * dim, &u_points); in main()
43 CeedVectorCreate(ceed, num_points * dim, &ones); in main()
85 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_GRAD, x_points, u, u_p… in main()
86 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_TRANSPOSE, CEED_EVAL_GRAD, x_points, ones, v); in main()
88 …CeedBasisApplyAddAtPoints(basis_u, 1, &num_points, CEED_TRANSPOSE, CEED_EVAL_GRAD, x_points, ones,… in main()
96 for (CeedInt i = 0; i < num_points * dim; i++) sum_2 += u_points_array[i]; in main()
|
| H A D | t356-basis.c | 35 …const CeedInt p = 9, q = 9, num_points = 4, x_dim = CeedIntPow(2, dim), p_dim = CeedIntPow(p, dim); in main() local
39 CeedVectorCreate(ceed, num_points * dim, &x_points); in main()
41 CeedVectorCreate(ceed, num_points * dim, &v); in main()
78 … CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_GRAD, x_points, u, v); in main()
85 for (CeedInt i = 0; i < num_points; i++) { in main()
88 for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[d * num_points + i]; in main()
92 if (fabs(v_array[d * num_points + i] - dfx) > 1E-3) { in main()
94 … printf("[%" CeedInt_FMT "] %f != %f = df(%f", dim, v_array[d * num_points + i], dfx, coord[0]); in main()
|
| H A D | t599-operator.c | 13 …* (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem = 4, num_points = num_elem * num… in main() local
23 CeedVectorCreate(ceed, dim * num_points, &x_points); in main()
25 CeedScalar x_array[dim * num_points]; in main()
38 CeedInt ind_x[num_elem + 1 + num_points]; in main()
41 for (CeedInt i = 0; i < num_points; i++) ind_x[num_elem + 1 + i] = i; in main()
42 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim, num_points * dim, CEED_MEM_HOST… in main()
44 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
|
| H A D | t357-basis.c | 34 …const CeedInt p = 9, q = 9, num_points = 4, x_dim = CeedIntPow(2, dim), p_dim = CeedIntPow(p, dim); in main() local
39 CeedVectorCreate(ceed, num_points * dim, &x_points); in main()
41 CeedVectorCreate(ceed, num_points * dim, &u_points); in main()
43 CeedVectorCreate(ceed, num_points * dim, &ones); in main()
85 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_GRAD, x_points, u, u_p… in main()
86 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_TRANSPOSE, CEED_EVAL_GRAD, x_points, ones, v); in main()
94 for (CeedInt i = 0; i < num_points * dim; i++) sum_2 += u_points_array[i]; in main()
|
| H A D | t352-basis.c | 25 …const CeedInt p = 9, q = 9, num_comp = 3, num_points = 4, x_dim = CeedIntPow(2, dim), p_dim = Ceed… in main() local
29 CeedVectorCreate(ceed, num_points * dim, &x_points); in main()
31 CeedVectorCreate(ceed, num_comp * num_points, &v); in main()
68 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, x_points, u, v… in main()
75 for (CeedInt i = 0; i < num_points; i++) { in main()
78 for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[d * num_points + i]; in main()
81 if (fabs(v_array[c * num_points + i] - fx) > 1E-4) { in main()
|
| H A D | t590-operator.c | 14 …* (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem = 4, num_points = num_elem * num… in main() local
24 CeedVectorCreate(ceed, dim * num_points, &x_points); in main()
26 CeedScalar x_array[dim * num_points]; in main()
39 CeedInt ind_x[num_elem + 1 + num_points]; in main()
42 for (CeedInt i = 0; i < num_points; i++) ind_x[num_elem + 1 + i] = i; in main()
43 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim, num_points * dim, CEED_MEM_HOST… in main()
|
| H A D | t591-operator.c | 14 …* (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem = 4, num_points = num_elem * num… in main() local
25 CeedVectorCreate(ceed, dim * num_points, &x_points); in main()
27 CeedScalar x_array[dim * num_points]; in main()
40 CeedInt ind_x[num_elem + 1 + num_points]; in main()
43 for (CeedInt i = 0; i < num_points; i++) ind_x[num_elem + 1 + i] = i; in main()
44 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim, num_points * dim, CEED_MEM_HOST… in main()
46 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
50 CeedVectorCreate(ceed, num_points, &q_data); in main()
|
| H A D | t593-operator.c | 14 …em + 1, num_nodes_u = num_elem * (p - 1) + 1, num_points_per_elem = 4, num_points = num_elem * num… in main() local
15 CeedInt ind_x[num_elem * 2], ind_u[num_elem * p], ind_x_points[num_elem + 1 + num_points]; in main()
16 CeedScalar x_array_mesh[num_nodes_x], x_array_points[num_points]; in main()
73 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
75 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
79 CeedVectorCreate(ceed, num_points, &q_data); in main()
|
| H A D | t351-basis.c | 25 …const CeedInt p = 9, q = 10, num_points = 4, x_dim = CeedIntPow(2, dim), p_dim = CeedIntPow(p, dim… in main() local
29 CeedVectorCreate(ceed, num_points * dim, &x_points); in main()
31 CeedVectorCreate(ceed, num_points, &v); in main()
68 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, x_points, u, v… in main()
75 for (CeedInt i = 0; i < num_points; i++) { in main()
78 for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[d * num_points + i]; in main()
|
| H A D | t592-operator.c | 14 …* (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem = 4, num_points = num_elem * num… in main() local
25 CeedVectorCreate(ceed, dim * num_points, &x_points); in main()
27 CeedScalar x_array[dim * num_points]; in main()
40 CeedInt ind_x[num_elem + 1 + num_points]; in main()
43 for (CeedInt i = 0; i < num_points; i++) ind_x[num_elem + 1 + i] = i; in main()
44 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim, num_points * dim, CEED_MEM_HOST… in main()
46 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
50 CeedVectorCreate(ceed, num_points, &q_data); in main()
176 for (CeedInt i = 0; i < num_points; i++) { in main()
|
| H A D | t597-operator.c | 25 … num_dofs = (n_x * 2 + 1) * (n_y * 2 + 1), num_points_per_elem = 4, num_points = num_elem * num… in main() local
31 CeedVectorCreate(ceed, dim * num_points, &x_points); in main()
33 CeedScalar x_array[dim * num_points]; in main()
46 CeedInt ind_x[num_elem + 1 + num_points]; in main()
49 for (CeedInt i = 0; i < num_points; i++) ind_x[num_elem + 1 + i] = i; in main()
50 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim, num_points * dim, CEED_MEM_HOST… in main()
52 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim * (dim + 1) / 2, num_points * di… in main()
71 CeedVectorCreate(ceed, num_points * dim * (dim + 1) / 2, &q_data); in main()
|
| H A D | t596-operator.c | 25 … num_dofs = (n_x * 2 + 1) * (n_y * 2 + 1), num_points_per_elem = 4, num_points = num_elem * num… in main() local
31 CeedVectorCreate(ceed, dim * num_points, &x_points); in main()
33 CeedScalar x_array[dim * num_points]; in main()
46 CeedInt ind_x[num_elem + 1 + num_points]; in main()
49 for (CeedInt i = 0; i < num_points; i++) ind_x[num_elem + 1 + i] = i; in main()
50 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim, num_points * dim, CEED_MEM_HOST… in main()
52 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
70 CeedVectorCreate(ceed, num_points, &q_data); in main()
|
| H A D | t350-basis.c | 20 const CeedInt p = 5, q = 5, num_points = 4; in main() local
27 CeedVectorCreate(ceed, num_points, &x_points); in main()
29 CeedVectorCreate(ceed, num_points, &v); in main()
59 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, x_points, u, v… in main()
66 for (CeedInt i = 0; i < num_points; i++) { in main()
|
| H A D | t355-basis.c | 26 const CeedInt p = 5, q = 5, num_points = 4; in main() local
33 CeedVectorCreate(ceed, num_points, &x_points); in main()
35 CeedVectorCreate(ceed, num_points, &v); in main()
65 CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_GRAD, x_points, u, v); in main()
72 for (CeedInt i = 0; i < num_points; i++) { in main()
|
| H A D | t594-operator.c | 14 …em + 1, num_nodes_u = num_elem * (p - 1) + 1, num_points_per_elem = 4, num_points = num_elem * num… in main() local
15 CeedInt ind_x[num_elem * 2], ind_u[num_elem * p], ind_x_points[num_elem + 1 + num_points]; in main()
16 CeedScalar x_array_mesh[num_nodes_x], x_array_points[num_points], assembled_true[num_nodes_u]; in main()
73 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
75 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C… in main()
79 CeedVectorCreate(ceed, num_points, &q_data); in main()
|
| H A D | t354-basis.c | 25 …const CeedInt p = 9, q = 9, num_points = 4, x_dim = CeedIntPow(2, dim), p_dim = CeedIntPow(p, dim); in main() local
29 CeedVectorCreate(ceed, num_points * dim, &x_points); in main()
32 CeedVectorCreate(ceed, num_points, &v); in main()
72 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, x_points, u, v… in main()
74 for (CeedInt i = 0; i < num_points; i++) { in main()
83 for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[d * num_points + i]; in main()
|
| H A D | t353-basis.c | 20 const CeedInt p = 5, q = 5, num_points = 4; in main() local
27 CeedVectorCreate(ceed, num_points, &x_points); in main()
30 CeedVectorCreate(ceed, num_points, &v); in main()
63 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, x_points, u, v… in main()
65 for (CeedInt i = 0; i < num_points; i++) { in main()
|
| H A D | t364-basis.c | 20 const CeedInt p = 5, q = 5, num_points = 4; in main() local
27 CeedVectorCreate(ceed, num_points, &x_points); in main()
30 CeedVectorCreate(ceed, num_points, &v); in main()
63 …CeedBasisApplyAtPoints(basis_u, 1, &num_points, CEED_NOTRANSPOSE, CEED_EVAL_INTERP, x_points, u, v… in main()
65 for (CeedInt i = 0; i < num_points; i++) { in main()
|
| /libCEED/examples/petsc/src/ |
| H A D | swarmutils.c | 55 PetscInt num_points = cell_points[num_elem + 1] - num_elem - 2; in DMSwarmCeedContextCreate() local
56 CeedInt offsets[num_elem + 1 + num_points]; in DMSwarmCeedContextCreate()
59 …for (PetscInt i = num_elem + 1; i < num_points + num_elem + 1; i++) offsets[i] = cell_points[i + 1… in DMSwarmCeedContextCreate()
63 …CeedElemRestrictionCreateAtPoints((*ctx)->ceed, num_elem, num_points, num_comp, num_points * num_c… in DMSwarmCeedContextCreate()
65 …CeedElemRestrictionCreateAtPoints((*ctx)->ceed, num_elem, num_points, dim, num_points * dim, CEED_… in DMSwarmCeedContextCreate()
67 …CeedElemRestrictionCreateAtPoints((*ctx)->ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST… in DMSwarmCeedContextCreate()
79 CeedVectorCreate((*ctx)->ceed, num_points * dim, &x_ref_points); in DMSwarmCeedContextCreate()
256 …ePointLocations(DM dm_swarm, PointSwarmType point_swarm_type, PetscInt num_points, PetscInt num_po… in DMSwarmInitalizePointLocations() argument
302 for (PetscInt p = 0; p < num_points; p++) { in DMSwarmInitalizePointLocations()
303 …point_coords[p * dim + 0] = -PetscCosReal((PetscReal)(p + 1) / (PetscReal)(num_points + 1) * PETSC… in DMSwarmInitalizePointLocations()
[all …]
|