differentiate 1 correctly

This commit is contained in:
Joachim Schöberl 2020-10-16 10:54:34 +02:00
parent b841b1c57b
commit 4cdaa6e3df

View File

@ -281,7 +281,7 @@ namespace netgen
*testout << "tang in plane" << endl;
double dlam1 = face.w1 * v;
double dlam2 = face.w2 * v;
double dlam3 = 1-dlam1-dlam2;
double dlam3 = -dlam1-dlam2;
if (printing)
*testout << "dlam = " << dlam1 << " " << dlam2 << " " << dlam3 << endl;
bool in1 = lam1 > eps_base1 || dlam1 > -eps_base1;