Permalink
Browse files

brute force for iga mortar mapper as default

  • Loading branch information...
Fabien Pean
Fabien Pean committed Nov 3, 2014
1 parent 5cd4e8d commit b591e2b19cb2334d169a303fb1c93cc9ba009df0
Showing with 2 additions and 2 deletions.
  1. +2 −2 Emperor/src/mapping/IGAMortarMapper.cpp
@@ -638,7 +638,7 @@ void IGAMortarMapper::computeCouplingMatrices() {
z=projectedElementFEWZ[nodeCount * 2 + 1];
polygonWZ.push_back(make_pair(w,z));
// Second point
isProjectedOnPatchBoundary=thePatch->computePointProjectionOnPatchBoundary(u, v, div, dis, P1, P2);
isProjectedOnPatchBoundary=thePatch->computePointProjectionOnPatchBoundary_Brute(u, v, div, dis, P1, P2);
if (isProjectedOnPatchBoundary && dis <= disTol) {
polygonUV.push_back(make_pair(u,v));
double P1x = projectedElementFEWZ[nodeCount * 2];
@@ -691,7 +691,7 @@ void IGAMortarMapper::computeCouplingMatrices() {
// First point
u = (*projectedCoords)[nodeIndexNext][patchIndex][0];
v = (*projectedCoords)[nodeIndexNext][patchIndex][1];
isProjectedOnPatchBoundary=thePatch->computePointProjectionOnPatchBoundary(u, v, div, dis, P2, P1);
isProjectedOnPatchBoundary=thePatch->computePointProjectionOnPatchBoundary_Brute(u, v, div, dis, P2, P1);
if (isProjectedOnPatchBoundary && dis <= disTol) {
polygonUV.push_back(make_pair(u,v));
double P1x = projectedElementFEWZ[nodeCountNext * 2];

0 comments on commit b591e2b

Please sign in to comment.