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

removed more output

Browse Source
This commit is contained in:
Lukas 2018-07-20 13:56:27 +02:00
parent 7f424a543a
commit b2ae6210fc
1 changed files with 0 additions and 38 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
libsrc/meshing/parallelmesh.cpp
View File

@ -484,9 +484,6 @@ namespace netgen
} }
const Element2d & sel2 = (*this)[sei]; const Element2d & sel2 = (*this)[sei];
FlatArray<const PointIndex> points2 = sel2.PNums(); FlatArray<const PointIndex> points2 = sel2.PNums();
for (int l = 0; l < points2.Size(); l++)
cout << points2[l] << " ";
cout << endl << endl;
has_ided_sels = true; has_ided_sels = true;
ided_sel[sei] = os1[0]; ided_sel[sei] = os1[0];
ided_sel[os1[0]] = sei; ided_sel[os1[0]] = sei;
@ -640,7 +637,6 @@ namespace netgen
const Segment & seg = (*this)[segi]; const Segment & seg = (*this)[segi];
dests.SetSize(0); dests.SetSize(0);
dests.Append(seg.GetPartition()); dests.Append(seg.GetPartition());
// cout << "partition of seg " << segi << " is " << dests[0] << endl;
for (int l = 0; l < per_seg_trans[segi].Size(); l++) for (int l = 0; l < per_seg_trans[segi].Size(); l++)
{ {
int dest2 = (*this)[per_seg_trans[segi][l]].GetPartition(); int dest2 = (*this)[per_seg_trans[segi][l]].GetPartition();
@ -773,8 +769,6 @@ namespace netgen
MyMPI_Recv (verts, 0, MPI_TAG_MESH+1); MyMPI_Recv (verts, 0, MPI_TAG_MESH+1);
cout << " got vertex data : " << endl << verts << endl;
int numvert = verts.Size(); int numvert = verts.Size();
paralleltop -> SetNV (numvert); paralleltop -> SetNV (numvert);
@ -807,7 +801,6 @@ namespace netgen
idents.SetType(idnr, (Identifications::ID_TYPE)pp_data[idnr]); idents.SetType(idnr, (Identifications::ID_TYPE)pp_data[idnr]);
} }
cout << "vertex identifications: " << endl;
int offset = 2*maxidentnr+1; int offset = 2*maxidentnr+1;
for(int idnr = 1; idnr < maxidentnr+1; idnr++) for(int idnr = 1; idnr < maxidentnr+1; idnr++)
{ {
@ -817,8 +810,6 @@ namespace netgen
for (int k = 0; k<npairs; k++) { for (int k = 0; k<npairs; k++) {
PointIndex loc1 = glob2loc_vert_ht.Get(pairdata[2*k]); PointIndex loc1 = glob2loc_vert_ht.Get(pairdata[2*k]);
PointIndex loc2 = glob2loc_vert_ht.Get(pairdata[2*k+1]); PointIndex loc2 = glob2loc_vert_ht.Get(pairdata[2*k+1]);
cout << loc1 << " <--> " << loc2 << " || ";
cout << pairdata[2*k] << " <--> " << pairdata[2*k+1] << endl;
idents.Add(loc1, loc2, idnr); idents.Add(loc1, loc2, idnr);
} }
} }
@ -848,19 +839,8 @@ namespace netgen
el.SetIndex(elarray[ind++]); el.SetIndex(elarray[ind++]);
el.SetNP(elarray[ind++]); el.SetNP(elarray[ind++]);
int ind0 = ind;
int ind1 = ind;
for ( int j = 0; j < el.GetNP(); j++) for ( int j = 0; j < el.GetNP(); j++)
el[j] = glob2loc_vert_ht.Get (elarray[ind++]); el[j] = glob2loc_vert_ht.Get (elarray[ind++]);
cout << "Add element, glob: " << endl;
for ( int j = 0; j < el.GetNP(); j++)
cout << elarray[ind0++] << " ";
cout << endl;
cout << "-> LOC: " << endl;
for ( int j = 0; j < el.GetNP(); j++)
cout << glob2loc_vert_ht.Get (elarray[ind1++]) << " ";
cout << endl;
AddVolumeElement (el); AddVolumeElement (el);
} }
@ -899,25 +879,11 @@ namespace netgen
int nep = selbuf[ii++]; int nep = selbuf[ii++];
Element2d tri(nep); Element2d tri(nep);
tri.SetIndex(faceind); tri.SetIndex(faceind);
int ii0 = ii;
int ii1 = ii;
cout << "Add Surfel, glob: " << endl;
for(int j = 1; j <= nep; j++) {
cout << selbuf[ii0++] << " ";
ii0++;
}
cout << endl;
for(int j = 1; j <= nep; j++) for(int j = 1; j <= nep; j++)
{ {
tri.PNum(j) = glob2loc_vert_ht.Get (selbuf[ii++]); tri.PNum(j) = glob2loc_vert_ht.Get (selbuf[ii++]);
tri.GeomInfoPi(j).trignum = selbuf[ii++]; tri.GeomInfoPi(j).trignum = selbuf[ii++];
} }
cout << "-> LOC: " << endl;
for(int j = 1; j <= nep; j++) {
cout << selbuf[ii1++] << " ";
glob2loc_vert_ht.Get (selbuf[ii1++]);
}
cout << endl;
paralleltop->SetLoc2Glob_SurfEl ( sel+1, globsel ); paralleltop->SetLoc2Glob_SurfEl ( sel+1, globsel );
AddSurfaceElement (tri); AddSurfaceElement (tri);
sel ++; sel ++;
@ -942,12 +908,8 @@ namespace netgen
globsegi = int (segmbuf[ii++]); globsegi = int (segmbuf[ii++]);
seg.si = int (segmbuf[ii++]); seg.si = int (segmbuf[ii++]);
cout << "add segm, glob: " << endl;
cout << segmbuf[ii] << " " << segmbuf[ii+1] << endl;
seg.pnums[0] = glob2loc_vert_ht.Get (int(segmbuf[ii++])); seg.pnums[0] = glob2loc_vert_ht.Get (int(segmbuf[ii++]));
seg.pnums[1] = glob2loc_vert_ht.Get (int(segmbuf[ii++])); seg.pnums[1] = glob2loc_vert_ht.Get (int(segmbuf[ii++]));
cout << "-> LOC: " << endl;
cout << seg.pnums[0] << " " << seg.pnums[1] << endl;
seg.geominfo[0].trignum = int( segmbuf[ii++] ); seg.geominfo[0].trignum = int( segmbuf[ii++] );
seg.geominfo[1].trignum = int ( segmbuf[ii++]); seg.geominfo[1].trignum = int ( segmbuf[ii++]);
seg.surfnr1 = int ( segmbuf[ii++]); seg.surfnr1 = int ( segmbuf[ii++]);