non-constant Jacobians for non-simplicials

This commit is contained in:
Joachim Schöberl 2016-02-16 07:33:58 +01:00
parent fcf5fd5fbb
commit d0ba57fd3a
2 changed files with 3 additions and 3 deletions

View File

@ -906,7 +906,7 @@ namespace netgen
el.SetIndex(hi); el.SetIndex(hi);
infile >> nep; infile >> nep;
el.SetNP(nep); el.SetNP(nep);
el.SetCurved (nep != 4);
for (int j = 0; j < nep; j++) for (int j = 0; j < nep; j++)
infile >> (int&)(el[j]); infile >> (int&)(el[j]);

View File

@ -993,7 +993,7 @@ namespace netgen
default: cerr << "Element::Element: unknown element with " << np << " points" << endl; default: cerr << "Element::Element: unknown element with " << np << " points" << endl;
} }
orderx = ordery = orderz = 1; orderx = ordery = orderz = 1;
is_curved = false; is_curved = typ != TET; // false;
} }
void Element :: SetOrder (const int aorder) void Element :: SetOrder (const int aorder)
@ -1031,7 +1031,7 @@ namespace netgen
flags.deleted = 0; flags.deleted = 0;
flags.fixed = 0; flags.fixed = 0;
orderx = ordery = orderz = 1; orderx = ordery = orderz = 1;
is_curved = false; is_curved = typ != TET; // false;
} }