add PointMatrix for Geometry::SEGMENT

make Element::SEGMENT one of the tensor product types
2021-08-29 10:00:01 -07:00 · 2021-08-29 09:59:27 -07:00
2 changed files with 15 additions and 2 deletions
@@ -3495,8 +3495,6 @@ void Mesh::SetMeshGen()
            mesh_geoms |= (1 << Geometry::TETRAHEDRON);
         case Element::TRIANGLE:
            mesh_geoms |= (1 << Geometry::TRIANGLE);
-         case Element::SEGMENT:
-            mesh_geoms |= (1 << Geometry::SEGMENT);
         case Element::POINT:
            mesh_geoms |= (1 << Geometry::POINT);
            meshgen |= 1;
@@ -3508,6 +3506,8 @@ void Mesh::SetMeshGen()
            mesh_geoms |= (1 << Geometry::SQUARE);
            mesh_geoms |= (1 << Geometry::SEGMENT);
            mesh_geoms |= (1 << Geometry::POINT);
+         case Element::SEGMENT:
+            mesh_geoms |= (1 << Geometry::SEGMENT);
            meshgen |= 2;
            break;

@@ -4253,6 +4253,19 @@ void NCMesh::GetPointMatrix(Geometry::Type geom, const char* ref_path,
            pm = PointMatrix(mid12, mid20, mid01);
         }
      }
+      else if (geom == Geometry::SEGMENT)
+      {
+         Point mid01(pm(0), pm(1));
+
+         if (child == 0)
+         {
+            pm = PointMatrix(pm(0), mid01);
+         }
+         else if (child == 1)
+         {
+            pm = PointMatrix(mid01, pm(1));
+         }
+      }
   }

   // write the points to the matrix
Author	SHA1	Message	Date
Robert W. Anderson	40746a2d3c	add PointMatrix for Geometry::SEGMENT	2021-08-29 10:00:01 -07:00
Robert W. Anderson	7a4b6f82f8	make Element::SEGMENT one of the tensor product types	2021-08-29 09:59:27 -07:00