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

curved elements: parallel and 2d

Browse Source
This commit is contained in:
Joachim Schoeberl 2011-08-09 21:23:05 +00:00
parent 8a450b506b
commit 5733e491fe
1 changed files with 41 additions and 33 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

74
libsrc/meshing/curvedelems.cpp
View File

@ -450,18 +450,20 @@ namespace netgen
Array<int> master_edgeorder; Array<int> master_edgeorder;
Array<int> master_edgecoeffsindex; Array<int> master_edgecoeffsindex;
Array<Vec<3> > master_edgecoeffs; Array<Vec<3> > master_edgecoeffs;
Array<int> master_faceorder;
Array<int> master_facecoeffsindex;
Array<Vec<3> > master_facecoeffs;
MyMPI_Bcast (master_edgeorder, 0); MyMPI_Bcast (master_edgeorder, 0);
MyMPI_Bcast (master_edgecoeffsindex, 0); MyMPI_Bcast (master_edgecoeffsindex, 0);
MyMPI_Bcast (master_edgecoeffs, 0); MyMPI_Bcast (master_edgecoeffs, 0);
Array<int> master_faceorder; if (mesh.GetDimension() == 3)
Array<int> master_facecoeffsindex; {
Array<Vec<3> > master_facecoeffs; MyMPI_Bcast (master_faceorder, 0);
MyMPI_Bcast (master_facecoeffsindex, 0);
MyMPI_Bcast (master_faceorder, 0); MyMPI_Bcast (master_facecoeffs, 0);
MyMPI_Bcast (master_facecoeffsindex, 0); }
MyMPI_Bcast (master_facecoeffs, 0);
const MeshTopology & top = mesh.GetTopology(); const MeshTopology & top = mesh.GetTopology();
@ -487,33 +489,36 @@ namespace netgen
edgecoeffs[edgecoeffsindex[i]+j] = master_edgecoeffs[master_edgecoeffsindex[glob-1]+j]; edgecoeffs[edgecoeffsindex[i]+j] = master_edgecoeffs[master_edgecoeffsindex[glob-1]+j];
} }
faceorder.SetSize (top.GetNFaces()); if (mesh.GetDimension() == 3)
facecoeffsindex.SetSize (top.GetNFaces()+1);
facecoeffsindex[0] = 0;
for (int i = 0; i < top.GetNFaces(); i++)
{ {
int glob = partop.GetDistantFaceNum (0, i+1); faceorder.SetSize (top.GetNFaces());
faceorder[i] = master_faceorder[glob-1]; facecoeffsindex.SetSize (top.GetNFaces()+1);
int ncoefs = master_facecoeffsindex[glob]-master_facecoeffsindex[glob-1]; facecoeffsindex[0] = 0;
facecoeffsindex[i+1] = facecoeffsindex[i] + ncoefs; for (int i = 0; i < top.GetNFaces(); i++)
{
int glob = partop.GetDistantFaceNum (0, i+1);
faceorder[i] = master_faceorder[glob-1];
int ncoefs = master_facecoeffsindex[glob]-master_facecoeffsindex[glob-1];
facecoeffsindex[i+1] = facecoeffsindex[i] + ncoefs;
}
facecoeffs.SetSize (facecoeffsindex[top.GetNFaces()]);
for (int i = 0; i < top.GetNFaces(); i++)
{
int glob = partop.GetDistantFaceNum (0, i+1);
int ncoefs = master_facecoeffsindex[glob]-master_facecoeffsindex[glob-1];
for (int j = 0; j < ncoefs; j++)
facecoeffs[facecoeffsindex[i]+j] = master_facecoeffs[master_facecoeffsindex[glob-1]+j];
}
} }
facecoeffs.SetSize (facecoeffsindex[top.GetNFaces()]); else
for (int i = 0; i < top.GetNFaces(); i++)
{ {
int glob = partop.GetDistantFaceNum (0, i+1); faceorder.SetSize (top.GetNFaces());
int ncoefs = master_facecoeffsindex[glob]-master_facecoeffsindex[glob-1]; faceorder = 1;
for (int j = 0; j < ncoefs; j++) facecoeffsindex.SetSize (top.GetNFaces()+1);
facecoeffs[facecoeffsindex[i]+j] = master_facecoeffs[master_facecoeffsindex[glob-1]+j]; facecoeffsindex = 0;
} }
/*
faceorder.SetSize (top.GetNFaces());
faceorder = 1;
facecoeffsindex.SetSize (top.GetNFaces()+1);
facecoeffsindex = 0;
*/
ishighorder = 1; ishighorder = 1;
return; return;
} }
@ -1040,9 +1045,12 @@ namespace netgen
MyMPI_Bcast (edgecoeffsindex, 0); MyMPI_Bcast (edgecoeffsindex, 0);
MyMPI_Bcast (edgecoeffs, 0); MyMPI_Bcast (edgecoeffs, 0);
MyMPI_Bcast (faceorder, 0); if (mesh.GetDimension() == 3)
MyMPI_Bcast (facecoeffsindex, 0); {
MyMPI_Bcast (facecoeffs, 0); MyMPI_Bcast (faceorder, 0);
MyMPI_Bcast (facecoeffsindex, 0);
MyMPI_Bcast (facecoeffs, 0);
}
} }
#endif #endif