From d0ba57fd3a2a977294ff2f35e1d9c1e1fbd5b544 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Joachim=20Sch=C3=B6berl?= <joachim.schoeberl@tuwien.ac.at>
Date: Tue, 16 Feb 2016 07:33:58 +0100
Subject: [PATCH] non-constant Jacobians for non-simplicials

---
 libsrc/meshing/meshclass.cpp | 2 +-
 libsrc/meshing/meshtype.cpp  | 4 ++--
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/libsrc/meshing/meshclass.cpp b/libsrc/meshing/meshclass.cpp
index dba48a4d..c4326b35 100644
--- a/libsrc/meshing/meshclass.cpp
+++ b/libsrc/meshing/meshclass.cpp
@@ -906,7 +906,7 @@ namespace netgen
                 el.SetIndex(hi);
                 infile >> nep;
                 el.SetNP(nep);
-
+                el.SetCurved (nep != 4);
                 for (int j = 0; j < nep; j++)
                   infile >> (int&)(el[j]);
 
diff --git a/libsrc/meshing/meshtype.cpp b/libsrc/meshing/meshtype.cpp
index 120efffc..2cc74fad 100644
--- a/libsrc/meshing/meshtype.cpp
+++ b/libsrc/meshing/meshtype.cpp
@@ -993,7 +993,7 @@ namespace netgen
       default: cerr << "Element::Element: unknown element with " << np << " points" << endl;
       }
     orderx = ordery = orderz = 1;
-    is_curved = false;
+    is_curved = typ != TET; // false;
   }
 
   void Element :: SetOrder (const int aorder) 
@@ -1031,7 +1031,7 @@ namespace netgen
     flags.deleted = 0;
     flags.fixed = 0;
     orderx = ordery = orderz = 1;
-    is_curved = false;
+    is_curved =  typ != TET; // false;
   }