OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added ASMu3D::getCorner()

Browse Source
This commit is contained in:
Kjetil Andre Johannessen 2016-11-18 13:48:29 +01:00 committed by Arne Morten Kvarving
parent 44781a3ee7
commit ba25c130b8
1 changed files with 21 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
src/ASM/LR/ASMu3D.C
View File

@ -2007,6 +2007,25 @@ bool ASMu3D::getOrder (int& p1, int& p2, int& p3) const
int ASMu3D::getCorner(int I, int J, int K, int basis) const
{
std::cerr << "ASMu3D::getCorner not implemented properly yet" << std::endl;
exit(776654);
int edge = LR::NONE;
if (I > 0) edge |= LR::EAST;
else edge |= LR::WEST;
if (J > 0) edge |= LR::NORTH;
else edge |= LR::SOUTH;
if (K > 0) edge |= LR::TOP;
else edge |= LR::BOTTOM;
const LR::LRSplineVolume* vol = this->getBasis(basis);
std::vector<LR::Basisfunction*> corner; // vector of one function for corner-input
vol->getEdgeFunctions(corner, (LR::parameterEdge) edge);
if( corner.empty() )
return -1;
size_t ofs = 1;
for (int i = 1; i < basis; i++)
ofs += this->getNoNodes(i);
return corner.front()->getId()+ofs;
}