@@ -19,7 +19,7 @@ template <size_t order> void runTest(Omega_h::Mesh mesh) {
1919 auto y = coords[vtx * dim + 1 ];
2020 shape_field (vtx, 0 , 0 , MeshField::Mesh_Topology::Vertex) = f (x, y);
2121 });
22- if constexpr (order == 2 ) {
22+ if (order == 2 ) {
2323 Kokkos::parallel_for (
2424 mesh.nedges (), KOKKOS_LAMBDA (int edge) {
2525 const auto left = edge2vtx[edge * 2 ];
@@ -30,11 +30,11 @@ template <size_t order> void runTest(Omega_h::Mesh mesh) {
3030 shape_field (edge, 0 , 0 , MeshField::Mesh_Topology::Edge) = f (x, y);
3131 });
3232 }
33- const auto numPtsPerElem = order == 2 ? 6 : 3 ;
34- Kokkos::View<double **> local_coords (" " , mesh.nelems () * numPtsPerElem , 3 );
33+ const auto numNodesPerElem = order == 2 ? 6 : 3 ;
34+ Kokkos::View<double **> local_coords (" " , mesh.nelems () * numNodesPerElem , 3 );
3535 Kokkos::parallel_for (
36- " set" , mesh.nelems () * numPtsPerElem , KOKKOS_LAMBDA (const int &i) {
37- const auto val = i % numPtsPerElem ;
36+ " set" , mesh.nelems () * numNodesPerElem , KOKKOS_LAMBDA (const int &i) {
37+ const auto val = i % numNodesPerElem ;
3838 local_coords (i, 0 ) = (val == 0 );
3939 local_coords (i, 1 ) = (val == 1 );
4040 local_coords (i, 2 ) = (val == 2 );
@@ -53,7 +53,7 @@ template <size_t order> void runTest(Omega_h::Mesh mesh) {
5353 });
5454
5555 auto eval_results = mesh_field.triangleLocalPointEval (
56- local_coords, numPtsPerElem , shape_field);
56+ local_coords, numNodesPerElem , shape_field);
5757
5858 int errors = 0 ;
5959 const auto triVerts = mesh.ask_elem_verts ();
@@ -72,15 +72,15 @@ template <size_t order> void runTest(Omega_h::Mesh mesh) {
7272 auto x = coords[vtx * dim];
7373 auto y = coords[vtx * dim + 1 ];
7474 auto expected = f (x, y);
75- auto actual = eval_results (tri * numPtsPerElem + node, 0 );
75+ auto actual = eval_results (tri * numNodesPerElem + node, 0 );
7676 if (Kokkos::fabs (expected - actual) > MeshField::MachinePrecision) {
7777 ++errors;
7878 Kokkos::printf (
7979 " expected: %lf, actual: %lf, element: %d, node(vtx): %d\n " ,
8080 expected, actual, tri, node);
8181 }
8282 }
83- for (int node = 3 ; node < numPtsPerElem ; ++node) {
83+ for (int node = 3 ; node < numNodesPerElem ; ++node) {
8484 const auto triDim = 2 ;
8585 const auto edgeDim = 1 ;
8686 const auto triToEdgeDegree = Omega_h::simplex_degree (triDim, edgeDim);
@@ -92,7 +92,7 @@ template <size_t order> void runTest(Omega_h::Mesh mesh) {
9292 auto x = (coords[left * dim] + coords[right * dim]) / 2.0 ;
9393 auto y = (coords[left * dim + 1 ] + coords[right * dim + 1 ]) / 2.0 ;
9494 auto expected = f (x, y);
95- auto actual = eval_results (tri * numPtsPerElem + node, 0 );
95+ auto actual = eval_results (tri * numNodesPerElem + node, 0 );
9696 if (Kokkos::fabs (expected - actual) > MeshField::MachinePrecision) {
9797 ++errors;
9898 Kokkos::printf (
0 commit comments