Lines Matching refs:CeedInt
48 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[3]);
49 …CartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[3], CeedInt degree, CeedInt num_comp,…
51 int SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[3], CeedInt mesh_degree, CeedVector …
52 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords);
57 CeedInt dim = 3; // dimension of the mesh in main()
58 CeedInt num_comp_x = 3; // number of x components in main()
59 CeedInt mesh_degree = 4; // polynomial degree for the mesh in main()
60 CeedInt sol_degree = 4; // polynomial degree for the solution in main()
61 CeedInt num_qpts = sol_degree + 2; // number of 1D quadrature points in main()
62 CeedInt prob_size = -1; // approximate problem size in main()
63 CeedInt help = 0, test = 0, gallery = 0, benchmark = 0; in main()
134 CeedInt num_xyz[3]; in main()
148 CeedInt mesh_size, sol_size; in main()
205 CeedInt elem_qpts = CeedIntPow(num_qpts, dim); in main()
206 CeedInt num_elem = 1; in main()
208 for (CeedInt d = 0; d < dim; d++) num_elem *= num_xyz[d]; in main()
250 for (CeedInt i = 0; i < sol_size; i++) { in main()
252 for (CeedInt d = 0; d < dim; d++) u_array[i] += x_array[i + d * sol_size]; in main()
267 for (CeedInt i = 0; i < benchmark; i++) { in main()
279 for (CeedInt i = 0; i < sol_size; i++) surface_area += fabs(v_array[i]); in main()
314 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[3]) { in GetCartesianMeshSize()
317 CeedInt num_elem = prob_size / CeedIntPow(degree, dim); in GetCartesianMeshSize()
318 CeedInt s = 0; // find s: num_elem/2 < 2^s <= num_elem in GetCartesianMeshSize()
324 CeedInt r = s % dim; in GetCartesianMeshSize()
326 for (CeedInt d = 0; d < dim; d++) { in GetCartesianMeshSize()
327 CeedInt sd = s / dim; in GetCartesianMeshSize()
338 …CartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[3], CeedInt degree, CeedInt num_comp,… in BuildCartesianRestriction()
340 CeedInt p = degree + 1; in BuildCartesianRestriction()
341 CeedInt num_nodes = CeedIntPow(p, dim); // number of scalar nodes per element in BuildCartesianRestriction()
342 CeedInt elem_qpts = CeedIntPow(num_qpts, dim); // number of qpts per element in BuildCartesianRestriction()
343 CeedInt nd[3], num_elem = 1, scalar_size = 1; in BuildCartesianRestriction()
345 for (CeedInt d = 0; d < dim; d++) { in BuildCartesianRestriction()
354 CeedInt *el_nodes = malloc(sizeof(CeedInt) * num_elem * num_nodes); in BuildCartesianRestriction()
356 for (CeedInt e = 0; e < num_elem; e++) { in BuildCartesianRestriction()
357 CeedInt e_xyz[3] = {1, 1, 1}, re = e; in BuildCartesianRestriction()
359 for (CeedInt d = 0; d < dim; d++) { in BuildCartesianRestriction()
363 CeedInt *local_elem_nodes = el_nodes + e * num_nodes; in BuildCartesianRestriction()
365 for (CeedInt l_nodes = 0; l_nodes < num_nodes; l_nodes++) { in BuildCartesianRestriction()
366 CeedInt g_nodes = 0, g_nodes_stride = 1, r_nodes = l_nodes; in BuildCartesianRestriction()
368 for (CeedInt d = 0; d < dim; d++) { in BuildCartesianRestriction()
389 int SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[3], CeedInt mesh_degree, CeedVector mesh_co… in SetCartesianMeshCoords()
390 CeedInt p = mesh_degree + 1; in SetCartesianMeshCoords()
391 CeedInt nd[3], scalar_size = 1; in SetCartesianMeshCoords()
393 for (CeedInt d = 0; d < dim; d++) { in SetCartesianMeshCoords()
404 for (CeedInt i = 0; i < p; i++) nodes[i] = 0.5 + 0.5 * nodes[i]; in SetCartesianMeshCoords()
405 for (CeedInt gs_nodes = 0; gs_nodes < scalar_size; gs_nodes++) { in SetCartesianMeshCoords()
406 CeedInt r_nodes = gs_nodes; in SetCartesianMeshCoords()
408 for (CeedInt d = 0; d < dim; d++) { in SetCartesianMeshCoords()
409 CeedInt d1d = r_nodes % nd[d]; in SetCartesianMeshCoords()
424 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords) { in TransformMeshCoords()
429 for (CeedInt i = 0; i < mesh_size; i++) { in TransformMeshCoords()