From f99db8c784a378980d1d101f8554a8e31408e2a7 Mon Sep 17 00:00:00 2001
From: Christoph Lehrenfeld <christoph.lehrenfeld@gmail.com>
Date: Sun, 17 Oct 2010 10:44:39 +0000
Subject: [PATCH] redo second order-refinement after uniform refinement

---
 libsrc/meshing/meshtype.cpp | 6 ++++++
 ng/ngpkg.cpp                | 4 ++++
 2 files changed, 10 insertions(+)

diff --git a/libsrc/meshing/meshtype.cpp b/libsrc/meshing/meshtype.cpp
index c1bf9592..739f6406 100644
--- a/libsrc/meshing/meshtype.cpp
+++ b/libsrc/meshing/meshtype.cpp
@@ -302,6 +302,12 @@ namespace netgen
           Swap (pnum[1], pnum[2]);
           break;
         }
+      case TRIG6:
+        {
+          Swap (pnum[1], pnum[2]);
+          Swap (pnum[4], pnum[5]);
+          break;
+        }
       case QUAD:
         {
           Swap (pnum[0], pnum[3]);
diff --git a/ng/ngpkg.cpp b/ng/ngpkg.cpp
index be4f5c45..82d0f7f0 100644
--- a/ng/ngpkg.cpp
+++ b/ng/ngpkg.cpp
@@ -1802,6 +1802,10 @@ namespace netgen
 	ng_geometry -> GetRefinement().Refine(*mesh);
       }
 
+//redo second order refinement if desired
+    if (mparam.secondorder)
+      const_cast<Refinement&> (ng_geometry->GetRefinement()).MakeSecondOrder(*mesh);
+
     return TCL_OK;
   }