ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-11 21:50:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

fix Alfeld typos

Browse Source
This commit is contained in:
Joachim Schoeberl 2023-07-12 18:24:19 -07:00
parent b603ed7cfd
commit 93b39231c4
5 changed files with 18 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
libsrc/include/nginterface_v2.hpp
View File

@ -389,7 +389,7 @@ namespace netgen
// also added from nginterface.h, still 1-based, need redesign // also added from nginterface.h, still 1-based, need redesign
void HPRefinement (int levels, double parameter = 0.125, void HPRefinement (int levels, double parameter = 0.125,
bool setorders = true,bool ref_level = false); bool setorders = true,bool ref_level = false);
void SplitAlefeld (); void SplitAlfeld ();
size_t GetNP() const; size_t GetNP() const;
int GetSurfaceElementSurfaceNumber (size_t ei) const; int GetSurfaceElementSurfaceNumber (size_t ei) const;

4
libsrc/interface/nginterface_v2.cpp
View File

@ -1229,11 +1229,11 @@ namespace netgen
::netgen::HPRefinement (*mesh, &ref, SPLIT_HP, levels, parameter, setorders, ref_level); ::netgen::HPRefinement (*mesh, &ref, SPLIT_HP, levels, parameter, setorders, ref_level);
} }
void Ngx_Mesh::SplitAlefeld () void Ngx_Mesh::SplitAlfeld ()
{ {
NgLock meshlock (mesh->MajorMutex(), true); NgLock meshlock (mesh->MajorMutex(), true);
Refinement & ref = const_cast<Refinement&> (mesh->GetGeometry()->GetRefinement()); Refinement & ref = const_cast<Refinement&> (mesh->GetGeometry()->GetRefinement());
::netgen::HPRefinement (*mesh, &ref, SPLIT_ALEFELD, 1, 0.5, true, true); ::netgen::HPRefinement (*mesh, &ref, SPLIT_ALFELD, 1, 0.5, true, true);
} }

18
libsrc/meshing/hpref_trig.hpp
View File

@ -787,34 +787,34 @@ HPRef_Struct reftrig_3singedges =
// HP_TRIG_3SINGEDGES // HP_TRIG_3SINGEDGES
int reftrig_Alefeld_splitedges[][3] = int reftrig_Alfeld_splitedges[][3] =
{ {
{ 0, 0, 0 } { 0, 0, 0 }
}; };
int reftrig_Alefeld_splitfaces[][4] = int reftrig_Alfeld_splitfaces[][4] =
{ {
{ 1, 2, 3, 4 }, { 1, 2, 3, 4 },
{ 0, 0, 0, 0 } { 0, 0, 0, 0 }
}; };
HPREF_ELEMENT_TYPE reftrig_Alefeld_newelstypes[] = HPREF_ELEMENT_TYPE reftrig_Alfeld_newelstypes[] =
{ {
HP_TRIG, HP_TRIG, HP_TRIG, HP_TRIG, HP_TRIG, HP_TRIG,
HP_NONE, HP_NONE,
}; };
int reftrig_Alefeld_newels[][8] = int reftrig_Alfeld_newels[][8] =
{ {
{ 1, 2, 4 }, { 1, 2, 4 },
{ 2, 3, 4 }, { 2, 3, 4 },
{ 3, 1, 4 }, { 3, 1, 4 },
}; };
HPRef_Struct reftrig_Alefeld = HPRef_Struct reftrig_Alfeld =
{ {
HP_TRIG, HP_TRIG,
reftrig_Alefeld_splitedges, reftrig_Alfeld_splitedges,
reftrig_Alefeld_splitfaces, reftrig_Alfeld_splitfaces,
0, 0,
reftrig_Alefeld_newelstypes, reftrig_Alfeld_newelstypes,
reftrig_Alefeld_newels reftrig_Alfeld_newels
}; };

8
libsrc/meshing/hprefinement.cpp
View File

@ -174,8 +174,8 @@ namespace netgen
case HP_TRIG_3SINGEDGES: case HP_TRIG_3SINGEDGES:
hps = &reftrig_3singedges; break; hps = &reftrig_3singedges; break;
case HP_TRIG_ALEFELD: case HP_TRIG_ALFELD:
hps = &reftrig_Alefeld; break; hps = &reftrig_Alfeld; break;
case HP_QUAD: case HP_QUAD:
@ -1846,7 +1846,7 @@ bool CheckSingularities(Mesh & mesh, INDEX_2_HASHTABLE<int> & edges, INDEX_2_HAS
cornerpoint, edgepoint, faces, face_edges, cornerpoint, edgepoint, faces, face_edges,
surf_edges, facepoint, levels, act_ref); surf_edges, facepoint, levels, act_ref);
if (act_ref == 1 && split == SPLIT_ALEFELD) if (act_ref == 1 && split == SPLIT_ALFELD)
sing = true; sing = true;
if(sing==0) return(sing); if(sing==0) return(sing);
@ -1900,7 +1900,7 @@ bool CheckSingularities(Mesh & mesh, INDEX_2_HASHTABLE<int> & edges, INDEX_2_HAS
hpel.type = ClassifyTrig(hpel, edges, edgepoint_dom, cornerpoint, edgepoint, hpel.type = ClassifyTrig(hpel, edges, edgepoint_dom, cornerpoint, edgepoint,
faces, face_edges, surf_edges, facepoint, dim, fd); faces, face_edges, surf_edges, facepoint, dim, fd);
else else
hpel.type = HP_TRIG_ALEFELD; hpel.type = HP_TRIG_ALFELD;
dd = 2; dd = 2;
break; break;

4
libsrc/meshing/hprefinement.hpp
View File

@ -41,7 +41,7 @@ enum HPREF_ELEMENT_TYPE {
HP_TRIG_SINGEDGES23, HP_TRIG_SINGEDGES23,
HP_TRIG_3SINGEDGES = 40, HP_TRIG_3SINGEDGES = 40,
HP_TRIG_ALEFELD, HP_TRIG_ALFELD,
HP_QUAD = 50, HP_QUAD = 50,
HP_QUAD_SINGCORNER, HP_QUAD_SINGCORNER,
@ -313,7 +313,7 @@ public:
}; };
enum SplittingType { SPLIT_HP, SPLIT_ALEFELD }; enum SplittingType { SPLIT_HP, SPLIT_ALFELD };
DLL_HEADER extern void HPRefinement (Mesh & mesh, Refinement * ref, SplittingType split, int levels, DLL_HEADER extern void HPRefinement (Mesh & mesh, Refinement * ref, SplittingType split, int levels,
double fac1=0.125, bool setorders=true, bool ref_level = false); double fac1=0.125, bool setorders=true, bool ref_level = false);