ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-25 05:20:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

debug messages in hp refinement

Browse Source
This commit is contained in:
Matthias Hochsteger 2020-10-19 19:30:12 +02:00
parent 4b5054d9fd
commit a2a2da13dc
1 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
libsrc/meshing/hprefinement.cpp
View File

@ -626,7 +626,7 @@ namespace netgen
// prepare new points // prepare new points
fac1 = max(0.001,min(0.33,fac1)); fac1 = max(0.001,min(0.33,fac1));
cout << " in HP-REFINEMENT with fac1 " << fac1 << endl; PrintMessage(3, " in HP-REFINEMENT with fac1 ", fac1);
*testout << " in HP-REFINEMENT with fac1 " << fac1 << endl; *testout << " in HP-REFINEMENT with fac1 " << fac1 << endl;
@ -1283,7 +1283,7 @@ namespace netgen
// cout << nwrong << " wrong prisms, " << nright << " right prisms" << endl; // cout << nwrong << " wrong prisms, " << nright << " right prisms" << endl;
} }
cout << nwrong << " wrong prisms, " << nright << " right prisms" << endl; PrintMessage(3, nwrong, " wrong prisms, ", nright, " right prisms");
NgArray<MeshPoint, 1> hpts(mesh.GetNP()); NgArray<MeshPoint, 1> hpts(mesh.GetNP());
@ -1338,7 +1338,7 @@ namespace netgen
sing = true; // iterate at least once sing = true; // iterate at least once
while(sing) while(sing)
{ {
cout << " Start new hp-refinement: step " << act_ref << endl; PrintMessage(3, " Start new hp-refinement: step ", act_ref);
DoRefinement (mesh, hpelements, ref, fac1); DoRefinement (mesh, hpelements, ref, fac1);
DoRefineDummies (mesh, hpelements, ref); DoRefineDummies (mesh, hpelements, ref);
@ -1435,16 +1435,16 @@ namespace netgen
int(Get_HPRef_Struct (hpel.type) -> geom)); int(Get_HPRef_Struct (hpel.type) -> geom));
} }
} }
cout << " Start with Update Topology " << endl; PrintMessage(5, " Start with Update Topology ");
mesh.UpdateTopology(); mesh.UpdateTopology();
cout << " Mesh Update Topology done " << endl; PrintMessage(5, " Mesh Update Topology done ");
act_ref++; act_ref++;
sing = ClassifyHPElements(mesh,hpelements, act_ref, levels); sing = ClassifyHPElements(mesh,hpelements, act_ref, levels);
} }
cout << " HP-Refinement done with " << --act_ref << " refinement steps." << endl; PrintMessage(3, " HP-Refinement done with ", --act_ref, " refinement steps.");
if(act_ref>=1) if(act_ref>=1)
{ {
@ -1813,7 +1813,7 @@ bool CheckSingularities(Mesh & mesh, INDEX_2_HASHTABLE<int> & edges, INDEX_2_HAS
} }
if (!sing) if (!sing)
cout << "PrepareElements no more to do for actual refinement " << act_ref << endl; PrintMessage(3, "PrepareElements no more to do for actual refinement ", act_ref);
return(sing); return(sing);
} }
@ -1960,8 +1960,8 @@ bool CheckSingularities(Mesh & mesh, INDEX_2_HASHTABLE<int> & edges, INDEX_2_HAS
} }
cout << "undefined elements update classification: " << cnt_undef << endl; PrintMessage(3, "undefined elements update classification: ", cnt_undef);
cout << "non-implemented in update classification: " << cnt_nonimplement << endl; PrintMessage(3, "non-implemented in update classification: ", cnt_nonimplement);
for (int i = 0; i < misses.Size(); i++) for (int i = 0; i < misses.Size(); i++)
if (misses[i]) if (misses[i])