ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-11-14 10:08:32 +05:00
Code Issues Packages Projects Releases Wiki Activity

Introduced fac2 to fix issue with face splits

Browse Source
This commit is contained in:
Umberto Zerbinati 2024-01-15 23:54:44 +00:00
parent d2bba6cb3c
commit 87c4e543ad
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
libsrc/meshing/hprefinement.cpp
View File

@ -677,8 +677,9 @@ namespace netgen
INDEX_3_HASHTABLE<int> newfacepts(elements.Size()+1); INDEX_3_HASHTABLE<int> newfacepts(elements.Size()+1);
// prepare new points // prepare new points
// fac1 = max(0.001,min(0.33,fac1)); // fac1 = max(0.001,min(0.33,fac1));
double fac2;
fac2 = max(0.001,min(0.33,fac1));
PrintMessage(3, " in HP-REFINEMENT with fac1 ", fac1); PrintMessage(3, " in HP-REFINEMENT with fac1 ", fac1);
*testout << " in HP-REFINEMENT with fac1 " << fac1 << endl; *testout << " in HP-REFINEMENT with fac1 " << fac1 << endl;
@ -728,8 +729,8 @@ namespace netgen
{ {
Point<3> np; Point<3> np;
for( int l=0;l<3;l++) for( int l=0;l<3;l++)
np(l) = (1-2*fac1)*mesh.Point(i3.I1())(l) np(l) = (1-2*fac2)*mesh.Point(i3.I1())(l)
+ fac1*mesh.Point(i3.I2())(l) + fac1*mesh.Point(i3.I3())(l); + fac2*mesh.Point(i3.I2())(l) + fac2*mesh.Point(i3.I3())(l);
int npi = mesh.AddPoint (np); int npi = mesh.AddPoint (np);
newfacepts.Set (i3, npi); newfacepts.Set (i3, npi);
} }
@ -817,9 +818,9 @@ namespace netgen
for (int l = 0; l < 3; l++) for (int l = 0; l < 3; l++)
newparam[hprs->splitfaces[j][3]-1][l] = newparam[hprs->splitfaces[j][3]-1][l] =
(1-2*fac1) * el.param[hprs->splitfaces[j][0]-1][l] + (1-2*fac2) * el.param[hprs->splitfaces[j][0]-1][l] +
fac1 * el.param[hprs->splitfaces[j][1]-1][l] + fac2 * el.param[hprs->splitfaces[j][1]-1][l] +
fac1 * el.param[hprs->splitfaces[j][2]-1][l]; fac2 * el.param[hprs->splitfaces[j][2]-1][l];
j++; j++;
} }
// split elements // split elements