cgma
Public Member Functions | Static Public Member Functions | Private Attributes
OCCSurface Class Reference

#include <OCCSurface.hpp>

Inheritance diagram for OCCSurface:
Surface GeometryEntity TopologyBridge

List of all members.

Public Member Functions

 OCCSurface (TopoDS_Face *theFace)
virtual ~OCCSurface ()
virtual CubitStatus closest_point_along_vector (CubitVector &from_point, CubitVector &along_vector, CubitVector &point_on_surface)
void add_hardpoint (OCCPoint *HardPoint)
void remove_hardpoint (OCCPoint *HardPoint)
DLIList< OCCPoint * > get_hardpoints ()
virtual void append_simple_attribute_virt (const CubitSimpleAttrib &)
virtual void remove_simple_attribute_virt (const CubitSimpleAttrib &)
virtual void remove_all_simple_attribute_virt ()
virtual CubitStatus get_simple_attribute (DLIList< CubitSimpleAttrib > &)
virtual CubitStatus get_simple_attribute (const CubitString &name, DLIList< CubitSimpleAttrib > &)
virtual CubitBox bounding_box () const
virtual GeometryQueryEngine * get_geometry_query_engine () const
virtual CubitStatus get_point_normal (CubitVector &, CubitVector &)
virtual void closest_point_trimmed (CubitVector from_point, CubitVector &point_on_surface)
virtual CubitStatus closest_point_uv_guess (CubitVector const &location, double &u, double &v, CubitVector *closest_location=NULL, CubitVector *unit_normal=NULL)
virtual CubitSense get_shell_sense (ShellSM *) const
virtual CubitStatus closest_point (CubitVector const &location, CubitVector *closest_location=NULL, CubitVector *unit_normal_ptr=NULL, CubitVector *curvature1_ptr=NULL, CubitVector *curvature2_ptr=NULL)
virtual CubitStatus principal_curvatures (CubitVector const &location, double &curvature_1, double &curvature_2, CubitVector *closest_location=NULL)
virtual CubitStatus evaluate (double u, double v, CubitVector *position, CubitVector *normal, CubitVector *curvature1, CubitVector *curvature2)
virtual CubitVector position_from_u_v (double u, double v)
virtual CubitStatus u_v_from_position (CubitVector const &location, double &u, double &v, CubitVector *closest_location=NULL)
virtual CubitBoolean is_periodic ()
virtual CubitBoolean is_periodic_in_U (double &period)
virtual CubitBoolean is_periodic_in_V (double &period)
virtual CubitBoolean is_singular_in_U (double u_param)
virtual CubitBoolean is_singular_in_V (double v_param)
virtual CubitBoolean is_closed_in_U ()
virtual CubitBoolean is_closed_in_V ()
virtual CubitStatus uv_derivitives (double u_param, double v_param, CubitVector &du, CubitVector &dv)
TopoDS_Face * get_TopoDS_Face ()
void set_TopoDS_Face (TopoDS_Face &face)
int get_loops (DLIList< OCCLoop * > &)
int get_coedges (DLIList< OCCCoEdge * > &)
int get_curves (DLIList< OCCCurve * > &)
int get_points (DLIList< OCCPoint * > &points)
void set_shell (OCCShell *shell)
void set_lump (OCCLump *lump)
void set_body (OCCBody *body)
OCCShell * my_shell ()
OCCLump * my_lump ()
OCCBody * my_body ()
CubitStatus get_bodies (DLIList< OCCBody * > &bodies)
virtual CubitBoolean is_parametric ()
virtual CubitBoolean get_param_range_U (double &lower_bound, double &upper_bound)
virtual CubitBoolean get_param_range_V (double &lower_bound, double &upper_bound)
virtual CubitBoolean is_position_on (CubitVector &test_position)
virtual CubitPointContainment point_containment (const CubitVector &point)
virtual CubitPointContainment point_containment (double u, double v)
virtual CubitStatus get_projected_distance_on_surface (CubitVector *pos1, CubitVector *pos2, double &distance)
virtual CubitStatus get_nurb_params (bool &rational, int &degree_u, int &degree_v, int &num_cntrl_pts_u, int &num_cntrl_pts_v, DLIList< CubitVector > &cntrl_pts, DLIList< double > &weights, DLIList< double > &u_knots, DLIList< double > &v_knots) const
virtual CubitStatus get_sphere_params (CubitVector &center, double &radius) const
virtual CubitStatus get_cone_params (CubitVector &center, CubitVector &normal, CubitVector &major_axis, double &radius_ratio, double &sine_angle, double &cos_angle) const
virtual CubitStatus get_torus_params (CubitVector &center, CubitVector &normal, double &major_radius, double &minor_radius) const
GeometryType geometry_type ()
virtual double measure ()
CubitVector center_point ()
virtual CubitSense get_geometry_sense ()
virtual void get_parents_virt (DLIList< TopologyBridge * > &parents)
virtual void get_children_virt (DLIList< TopologyBridge * > &children)
CubitBoolean is_flat ()
CubitBoolean is_spherical ()
CubitBoolean is_conical ()
CubitStatus update_OCC_entity (BRepBuilderAPI_ModifyShape *aBRepTrsf, BRepAlgoAPI_BooleanOperation *op=NULL)

Static Public Member Functions

static CubitStatus update_OCC_entity (TopoDS_Face &old_surface, TopoDS_Shape &new_surface, BRepBuilderAPI_MakeShape *op, TopoDS_Vertex *removed_vertex=NULL, LocOpe_SplitShape *sp=NULL)

Private Attributes

TopoDS_Face * myTopoDSFace
OCCShell * myShell
OCCLump * myLump
OCCBody * myBody
DLIList< OCCPoint * > myHardPoints

Detailed Description

Definition at line 56 of file OCCSurface.hpp.

Constructor & Destructor Documentation

OCCSurface::OCCSurface ( TopoDS_Face *  theFace)

Definition at line 101 of file OCCSurface.cpp.

{
  myTopoDSFace = theFace;
  myShell = NULL;
  myLump = NULL;
  myBody = NULL;
  if(myTopoDSFace && !myTopoDSFace->IsNull())
    assert(myTopoDSFace->ShapeType() == TopAbs_FACE);
}
OCCSurface::~OCCSurface ( ) [virtual]

Definition at line 112 of file OCCSurface.cpp.

{
  if(myTopoDSFace)
  {
    myTopoDSFace->Nullify();
    delete (TopoDS_Face*)myTopoDSFace;
    myTopoDSFace = NULL;
  }
}

Member Function Documentation

void OCCSurface::add_hardpoint ( OCCPoint *  HardPoint) [inline]

Definition at line 79 of file OCCSurface.hpp.

{myHardPoints.append(HardPoint);}
void OCCSurface::append_simple_attribute_virt ( const CubitSimpleAttrib &  csa) [virtual]

Implements TopologyBridge.

Definition at line 142 of file OCCSurface.cpp.

  { OCCAttribSet::append_attribute(csa, *myTopoDSFace); }
CubitBox OCCSurface::bounding_box ( void  ) const [virtual]

Implements GeometryEntity.

Definition at line 289 of file OCCSurface.cpp.

{
  TopoDS_Face face = *myTopoDSFace;
  BRepAdaptor_Surface asurface(face);
  Bnd_Box aBox;
  BndLib_AddSurface::Add(asurface, Precision::Approximation(), aBox);
  double min[3], max[3];
  aBox.Get( min[0], min[1], min[2], max[0], max[1], max[2]);
  return CubitBox(min, max);
}
CubitVector OCCSurface::center_point ( )

Definition at line 783 of file OCCSurface.cpp.

{
  GProp_GProps myProps;
  BRepGProp::SurfaceProperties(*myTopoDSFace, myProps);
  gp_Pnt pt = myProps.CentreOfMass();
  CubitVector v(pt.X(),pt.Y(), pt.Z());
  return v; 
}
CubitStatus OCCSurface::closest_point ( CubitVector const &  location,
CubitVector *  closest_location = NULL,
CubitVector *  unit_normal_ptr = NULL,
CubitVector *  curvature1_ptr = NULL,
CubitVector *  curvature2_ptr = NULL 
) [virtual]

Implements Surface.

Definition at line 325 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  gp_Pnt p(location.x(), location.y(), location.z()), newP(0.0, 0.0, 0.0);
  double minDist=0.0, u, v;
  int i;
  BRepLProp_SLProps SLP(asurface, 2, Precision::PConfusion());
  Extrema_ExtPS ext(p, asurface, Precision::Approximation(), Precision::Approximation());
  if (ext.IsDone() && (ext.NbExt() > 0)) {
      for ( i = 1 ; i <= ext.NbExt() ; i++ ) {
        if ( (i==1) || (p.Distance(ext.Point(i).Value()) < minDist) ) {
              minDist = p.Distance(ext.Point(i).Value());
              newP = ext.Point(i).Value();
              ext.Point(i).Parameter(u, v);
              SLP.SetParameters(u, v);
            }
      }
  
    if (closest_location != NULL)
        *closest_location = CubitVector(newP.X(), newP.Y(), newP.Z());
    if (unit_normal_ptr != NULL) {
      gp_Dir normal;
      if (SLP.IsNormalDefined()) {
        normal = SLP.Normal();
        *unit_normal_ptr = CubitVector(normal.X(), normal.Y(), normal.Z()); 
      }
    }
  
        gp_Dir MaxD, MinD;
        if (SLP.IsCurvatureDefined())
        {
       SLP.CurvatureDirections(MaxD, MinD);
           if (curvature_1 != NULL)
              *curvature_1 = CubitVector(MinD.X(), MinD.Y(), MinD.Z());
           if (curvature_2 != NULL)
              *curvature_2 = CubitVector(MaxD.X(), MaxD.Y(), MaxD.Z());
        }
  return CUBIT_SUCCESS;
  }
  //return as OCC did.
  return CUBIT_SUCCESS;
}
CubitStatus OCCSurface::closest_point_along_vector ( CubitVector &  from_point,
CubitVector &  along_vector,
CubitVector &  point_on_surface 
) [virtual]

Implements Surface.

Definition at line 203 of file OCCSurface.cpp.

{
  // define the search ray in a way that OCC will understand
  gp_Pnt vectorOrigin(from_point.x(), from_point.y(), from_point.z());
  gp_Dir vectorDir(along_vector.x(), along_vector.y(), along_vector.z());
  gp_Lin occVector(vectorOrigin, vectorDir);

  // perform the OCC intersection algorithm
  IntCurvesFace_ShapeIntersector csIntersector;
  csIntersector.Load(*get_TopoDS_Face(), Precision::Intersection());
  csIntersector.Perform(occVector, 0, 1e63);

  // fail if intersection could not be computed or there is no intersection
  if (!csIntersector.IsDone() || csIntersector.NbPnt() == 0)
  {
    return CUBIT_FAILURE;
  }

  // identify which of the intersection points along the vector is the closest
  double minWParam = -1;
  int minIndex = 0;
  for (int iPntIndx = 1; iPntIndx <= csIntersector.NbPnt(); ++iPntIndx)
  {
    double wParam = csIntersector.WParameter(iPntIndx);
    if (minIndex == 0 || wParam < minWParam)
    {
      minIndex = iPntIndx;
      minWParam = wParam;
    }
  }

  // package and return the closest intersection point
  gp_Pnt closestIntPnt = csIntersector.Pnt(minIndex);
  point_on_surface.set(closestIntPnt.X(),
      closestIntPnt.Y(), closestIntPnt.Z());
  return CUBIT_SUCCESS;
}
void OCCSurface::closest_point_trimmed ( CubitVector  from_point,
CubitVector &  point_on_surface 
) [virtual]

Implements Surface.

Definition at line 379 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  gp_Pnt p(from_point.x(), from_point.y(), from_point.z()), newP(0.0, 0.0, 0.0);
  double minDist=0.0;
  int i;
  Extrema_ExtPS ext(p, asurface, Precision::Approximation(), Precision::Approximation());
  if (ext.IsDone() && (ext.NbExt() > 0)) {
      for ( i = 1 ; i <= ext.NbExt() ; i++ ) {
         if ( (i==1) || (p.Distance(ext.Point(i).Value()) < minDist) ) {
              minDist = p.Distance(ext.Point(i).Value());
              newP = ext.Point(i).Value();
          }
      }
          point_on_surface = CubitVector(newP.X(), newP.Y(), newP.Z());
  }
  else
    return;

  CubitPointContainment pos = point_containment(point_on_surface);
  if(pos == CUBIT_PNT_OUTSIDE)
  {
    DLIList<OCCCurve*> curves;
    int num_curve = get_curves(curves);
    double d_min = 0., d;
    OCCCurve* the_curve = NULL;
    gp_Pnt pt = gp_Pnt(point_on_surface.x(), point_on_surface.y(), 
                       point_on_surface.z());
    TopoDS_Vertex theVertex = BRepBuilderAPI_MakeVertex(pt);
    CubitVector closest_location;
    do
    {
      for (i = 0; i < num_curve; i++)
      {
        OCCCurve* curve = curves.get_and_step();
        TopoDS_Edge* theEdge = curve->get_TopoDS_Edge( );
        BRepExtrema_DistShapeShape distShapeShape(*theEdge, theVertex);
        d = distShapeShape.Value();
        if ( i == 0 || d_min > d)
        {
          d_min = d;
          the_curve = curve;
        }
     } 
     the_curve->closest_point(point_on_surface, closest_location);       
     pos = the_curve->point_containment(closest_location);
     if(pos == CUBIT_PNT_OUTSIDE)
       curves.remove(the_curve);
     else
       break;
    }while(curves.size() > 0);
    point_on_surface = closest_location;
  }
}
CubitStatus OCCSurface::closest_point_uv_guess ( CubitVector const &  location,
double &  u,
double &  v,
CubitVector *  closest_location = NULL,
CubitVector *  unit_normal = NULL 
) [virtual]

Implements Surface.

Definition at line 307 of file OCCSurface.cpp.

{
  // don't use u and v guesses
 return closest_point(location, closest_location, unit_normal);
}
CubitStatus OCCSurface::evaluate ( double  u,
double  v,
CubitVector *  position,
CubitVector *  normal,
CubitVector *  curvature1,
CubitVector *  curvature2 
) [virtual]

Implements Surface.

Definition at line 483 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);

  gp_Pnt p = asurface.Value(u, v);
  if(position!=NULL)
      position->set(p.X(), p.Y(), p.Z());


  BRepLProp_SLProps SLP(asurface, 2, Precision::PConfusion());
  SLP.SetParameters(u, v);


  if(normal!=NULL)
  {
    gp_Dir occ_normal;
    //normal of a RefFace point to outside of the material
    if (SLP.IsNormalDefined()) 
    {
      occ_normal = SLP.Normal();
      normal->set(occ_normal.X(), occ_normal.Y(), occ_normal.Z()); 
    }
  }


  gp_Dir MaxD, MinD;
  if (curvature1 && curvature2 && SLP.IsCurvatureDefined())
  {
    SLP.CurvatureDirections(MaxD, MinD);
    if (curvature1 != NULL)
      *curvature1 = CubitVector(MinD.X(), MinD.Y(), MinD.Z());
    if (curvature2 != NULL)
      *curvature2 = CubitVector(MaxD.X(), MaxD.Y(), MaxD.Z());
  }


  return CUBIT_SUCCESS;
}
GeometryType OCCSurface::geometry_type ( ) [virtual]

Reimplemented from Surface.

Definition at line 259 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  if (asurface.GetType() == GeomAbs_BezierSurface)
     return SPLINE_SURFACE_TYPE;
  if (asurface.GetType() == GeomAbs_BSplineSurface)
     return SPLINE_SURFACE_TYPE;
  if (asurface.GetType() == GeomAbs_Plane)      
     return PLANE_SURFACE_TYPE;
  if (asurface.GetType() == GeomAbs_Cylinder ||
      asurface.GetType() == GeomAbs_Cone)
     return CONE_SURFACE_TYPE;
  if (asurface.GetType() == GeomAbs_Sphere)
     return SPHERE_SURFACE_TYPE;
  if (asurface.GetType() == GeomAbs_Torus)
      return TORUS_SURFACE_TYPE;
  if (asurface.GetType() == GeomAbs_SurfaceOfRevolution)
     return UNDEFINED_SURFACE_TYPE;
  if (asurface.GetType() == GeomAbs_SurfaceOfExtrusion)
     return UNDEFINED_SURFACE_TYPE;
  if (asurface.GetType() == GeomAbs_OffsetSurface)
     return UNDEFINED_SURFACE_TYPE;
  return UNDEFINED_SURFACE_TYPE;  
}
CubitStatus OCCSurface::get_bodies ( DLIList< OCCBody * > &  bodies)

Definition at line 1308 of file OCCSurface.cpp.

{
   TopoDS_Face* topo_face = this->get_TopoDS_Face();
   OCCQueryEngine* oqe = OCCQueryEngine::instance();
   DLIList <OCCBody* > *all_bodies = oqe->BodyList;
   TopTools_IndexedDataMapOfShapeListOfShape M;
   OCCBody * body = NULL;
   for(int j = 0; j <  all_bodies->size(); j++)
   {
     body = all_bodies->get_and_step();
     TopExp_Explorer Ex;
     TopoDS_Face the_face;
     TopoDS_Shape* pshape;
     body->get_TopoDS_Shape(pshape);
     TopoDS_Shape ashape;
     if (pshape && !pshape->IsNull() && 
         pshape->ShapeType() <= TopAbs_FACE &&
         OCCQueryEngine::instance()->OCCMap->IsBound(*pshape) == Standard_True)
       ashape = *pshape;
     M.Clear();

     TopExp::MapShapesAndAncestors(ashape, TopAbs_FACE, TopAbs_COMPOUND, M);
     if(!M.Contains(*topo_face))
       continue;
     bodies.append_unique(body);
  }
  return CUBIT_SUCCESS;
}
void OCCSurface::get_children_virt ( DLIList< TopologyBridge * > &  children) [virtual]

Implements TopologyBridge.

Definition at line 884 of file OCCSurface.cpp.

{
  TopTools_IndexedMapOfShape M;
  TopExp::MapShapes(*myTopoDSFace, TopAbs_WIRE, M);
  int ii;
  for (ii=1; ii<=M.Extent(); ii++) {
     TopologyBridge *loop = OCCQueryEngine::instance()->occ_to_cgm(M(ii));
     if(loop)
       children.append_unique(loop);
  }
}
int OCCSurface::get_coedges ( DLIList< OCCCoEdge * > &  result_list)

Definition at line 959 of file OCCSurface.cpp.

{
  DLIList<OCCLoop*> loop_list;
  get_loops( loop_list );
  loop_list.reset();
  for ( int i = 0; i < loop_list.size(); i++ )
  {
    result_list += loop_list.next(i)->coedges( );
  }
  return result_list.size();
}
CubitStatus OCCSurface::get_cone_params ( CubitVector &  center,
CubitVector &  normal,
CubitVector &  major_axis,
double &  radius_ratio,
double &  sine_angle,
double &  cos_angle 
) const [virtual]

Implements Surface.

Definition at line 1447 of file OCCSurface.cpp.

{
  if(const_cast<OCCSurface*> (this)->geometry_type() != CONE_SURFACE_TYPE)
    return CUBIT_FAILURE;

  //all cone type surface in OCC have circular base.
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  gp_Cone cone = asurface.Cone();
  double half_angle = cone.SemiAngle();
  sine_angle = sin(half_angle);
  cos_angle = cos(half_angle);

  gp_Ax1 axis = cone.Axis();
  gp_Dir dir = axis.Direction();
  normal.set(dir.X(), dir.Y(), dir.Z());

  gp_Pnt center_pt = cone.Location();
  center = CubitVector(center_pt.X(), center_pt.Y(), center_pt.Z());

  gp_Ax1 x_axis = cone.XAxis();
  dir = x_axis.Direction();   
  major_axis.set(dir.X(), dir.Y(), dir.Z());

  radius_ratio = 1;
  return CUBIT_SUCCESS;
}
int OCCSurface::get_curves ( DLIList< OCCCurve * > &  result_list)

Definition at line 971 of file OCCSurface.cpp.

{
  DLIList<OCCCoEdge*> coedge_list;
  get_coedges( coedge_list );
  coedge_list.reset();
  for ( int i = coedge_list.size(); i--; )
  {
    OCCCoEdge* coedge = coedge_list.get_and_step();
    OCCCurve* curve = dynamic_cast<OCCCurve*>(coedge->curve());
    if (curve)
      result_list.append_unique(curve);
  }
  return result_list.size();
}
GeometryQueryEngine * OCCSurface::get_geometry_query_engine ( ) const [virtual]

Implements TopologyBridge.

Definition at line 197 of file OCCSurface.cpp.

{
   return OCCQueryEngine::instance();
}
CubitSense OCCSurface::get_geometry_sense ( ) [virtual]

Implements Surface.

Definition at line 838 of file OCCSurface.cpp.

{
  return CUBIT_FORWARD;
}
DLIList<OCCPoint*> OCCSurface::get_hardpoints ( ) [inline]

Definition at line 81 of file OCCSurface.hpp.

{return myHardPoints;}
int OCCSurface::get_loops ( DLIList< OCCLoop * > &  result_list)

Definition at line 946 of file OCCSurface.cpp.

{
  TopTools_IndexedMapOfShape M;
  TopExp::MapShapes(*myTopoDSFace, TopAbs_WIRE, M);
  int ii;
  for (ii=1; ii<=M.Extent(); ii++) {
     TopologyBridge *loop = OCCQueryEngine::instance()->occ_to_cgm(M(ii));
     if(loop)
       result_list.append_unique(dynamic_cast<OCCLoop*>(loop));
  }
  return result_list.size();
}
CubitStatus OCCSurface::get_nurb_params ( bool &  rational,
int &  degree_u,
int &  degree_v,
int &  num_cntrl_pts_u,
int &  num_cntrl_pts_v,
DLIList< CubitVector > &  cntrl_pts,
DLIList< double > &  weights,
DLIList< double > &  u_knots,
DLIList< double > &  v_knots 
) const [virtual]

Implements Surface.

Definition at line 1350 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurf(*myTopoDSFace);
  Handle_Geom_BSplineSurface h_S = NULL;
  if (asurf.GetType() == GeomAbs_BSplineSurface)
    h_S = asurf.BSpline();
  else
    return CUBIT_FAILURE;
  assert ( h_S != NULL);

  if(h_S->IsURational() || h_S->IsVRational())
    rational = true;
  else 
    rational = false;

  degree_u = h_S->UDegree();
  degree_v = h_S->VDegree();

  num_cntrl_pts_u = h_S->NbUPoles();
  num_cntrl_pts_v = h_S->NbVPoles();

  TColgp_Array2OfPnt P(1, h_S->NbUPoles(), 1, h_S->NbVPoles());
  h_S->Poles(P);
  gp_Pnt cnt_p;
  for(int i = P.LowerRow(); i <= P.UpperRow(); i++)
  {
    for (int j = P.LowerCol(); j <= P.UpperCol(); j++)
    {
      cnt_p = P(i,j);
      cntrl_pts.append(CubitVector(cnt_p.X(), cnt_p.Y(), cnt_p.Z()));
    }
  } 
  TColStd_Array2OfReal W(1,h_S->NbUPoles(), 1, h_S->NbVPoles());
  h_S->Weights(W);
  for(int i = W.LowerRow(); i <= W.UpperRow(); i++)
  { 
    for (int j = W.LowerCol(); j <= W.UpperCol(); j++)
    {
      cntrl_pt_weights.append(W(i,j));
    }
  }
  TColStd_Array1OfReal Ku(1,h_S->NbUKnots()), Kv(1, h_S->NbVKnots());
  h_S->UKnots(Ku);
  h_S->VKnots(Kv);
  for(int i = Ku.Lower(); i <= Ku.Upper(); i++)
    u_knots.append(Ku.Value(i));

  for (int i = Kv.Lower(); i <= Kv.Upper(); i++)
    v_knots.append(Kv.Value(i));
   return CUBIT_SUCCESS; 
}
CubitBoolean OCCSurface::get_param_range_U ( double &  lower_bound,
double &  upper_bound 
) [virtual]

Implements Surface.

Definition at line 742 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  lower_bound = asurface.FirstUParameter();
  upper_bound = asurface.LastUParameter();
  return CUBIT_TRUE;
}
CubitBoolean OCCSurface::get_param_range_V ( double &  lower_bound,
double &  upper_bound 
) [virtual]

Implements Surface.

Definition at line 758 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  lower_bound = asurface.FirstVParameter();

  upper_bound = asurface.LastVParameter();
  return CUBIT_TRUE;
}
void OCCSurface::get_parents_virt ( DLIList< TopologyBridge * > &  parents) [virtual]

Implements TopologyBridge.

Definition at line 843 of file OCCSurface.cpp.

{ 
  if(myShell) //shell or sheet body
  {
    parents.append(myShell);
    return;
  }

  OCCQueryEngine* oqe = (OCCQueryEngine*) get_geometry_query_engine();
  OCCBody * body = NULL;
  DLIList <OCCBody* > *bodies = oqe->BodyList;
  TopTools_IndexedDataMapOfShapeListOfShape M;
  for(int i = 0; i <  bodies->size(); i++)
  {
     body = bodies->get_and_step();
     DLIList<OCCSurface*> surfaces;
     body = bodies->get_and_step();
     body->get_all_surfaces(surfaces);
     if(surfaces.move_to(this))
     {
       TopoDS_Shape* shape = oqe->instance()->
          get_TopoDS_Shape_of_entity(CAST_TO(body, TopologyBridge));
       TopExp::MapShapesAndAncestors(*shape, TopAbs_FACE, TopAbs_SHELL, M);
       if(!M.Contains(*(get_TopoDS_Face())))
         continue;
       const TopTools_ListOfShape& ListOfShapes =
                                M.FindFromKey(*(get_TopoDS_Face()));
       if (!ListOfShapes.IsEmpty())
       {
         TopTools_ListIteratorOfListOfShape it(ListOfShapes) ;
         for (;it.More(); it.Next())
         {
           TopoDS_Shell Shell = TopoDS::Shell(it.Value());
           int k = oqe->OCCMap->Find(Shell);
           parents.append((OCCShell*)(oqe->OccToCGM->find(k))->second);
         }
       }
    }
  }
}
CubitStatus OCCSurface::get_point_normal ( CubitVector &  location,
CubitVector &  normal 
) [virtual]

Implements Surface.

Definition at line 301 of file OCCSurface.cpp.

{
  return closest_point( bounding_box().center(), &location, &normal );
}
int OCCSurface::get_points ( DLIList< OCCPoint * > &  points)

Definition at line 986 of file OCCSurface.cpp.

{
  DLIList<OCCCurve*> curves;
  int num_crv = get_curves(curves);
  for(int i = 0; i < num_crv; i++)
  {
    OCCCurve* curve = curves.get_and_step();
    curve->get_points(points);
  }
  points += get_hardpoints();
  return points.size();
}
CubitStatus OCCSurface::get_projected_distance_on_surface ( CubitVector *  pos1,
CubitVector *  pos2,
double &  distance 
) [virtual]

Implements Surface.

Definition at line 1337 of file OCCSurface.cpp.

{
  CubitVector closest_point1;
  this->closest_point_trimmed(*pos1, closest_point1);  
  CubitVector closest_point2;
  this->closest_point_trimmed(*pos2, closest_point2); 
  distance = closest_point1.distance_between(closest_point2);
  return CUBIT_SUCCESS;
}
CubitSense OCCSurface::get_shell_sense ( ShellSM *  shell_ptr) const [virtual]

Implements Surface.

Definition at line 897 of file OCCSurface.cpp.

{
  OCCShell* shell = dynamic_cast<OCCShell*>(shell_ptr);
  if (!shell) // error
    return CUBIT_UNKNOWN;
    
  if (shell->is_sheet())  // relative sense is "both" for sheet
    return CUBIT_UNKNOWN;

  TopoDS_Shell* shellShapePtr = shell->get_TopoDS_Shell();
  if (!shellShapePtr) // error
    return CUBIT_UNKNOWN;

  TopExp_Explorer shellFaceExp(shellShapePtr->Oriented(TopAbs_FORWARD),
      TopAbs_FACE);
  bool isForward = false;
  bool isReversed = false;
  while (shellFaceExp.More())
  {
    const TopoDS_Shape shellFace = shellFaceExp.Current();
    if (shellFace.IsSame(*myTopoDSFace))
    {
      if (shellFace.Orientation() == TopAbs_FORWARD)
      {
        isForward = true;
      }
      else if (shellFace.Orientation() == TopAbs_REVERSED)
      {
        isReversed = true;
      }
      else
      {
        isForward = true;
        isReversed = true;
      }
    }
    shellFaceExp.Next();
  }

  if (isForward && isReversed)
    return CUBIT_UNKNOWN;
  if (isReversed)
    return CUBIT_REVERSED;
  if (isForward)
    return CUBIT_FORWARD;
  return CUBIT_UNKNOWN;
}
CubitStatus OCCSurface::get_simple_attribute ( DLIList< CubitSimpleAttrib > &  csa_list) [virtual]

Implements TopologyBridge.

Definition at line 181 of file OCCSurface.cpp.

  { return OCCAttribSet::get_attributes(*myTopoDSFace,csa_list); }
CubitStatus OCCSurface::get_simple_attribute ( const CubitString &  name,
DLIList< CubitSimpleAttrib > &  csa_list 
) [virtual]

Implements TopologyBridge.

Definition at line 185 of file OCCSurface.cpp.

  { return OCCAttribSet::get_attributes( name, *myTopoDSFace, csa_list ); }
CubitStatus OCCSurface::get_sphere_params ( CubitVector &  center,
double &  radius 
) const [virtual]

Implements Surface.

Definition at line 1421 of file OCCSurface.cpp.

{
  if(const_cast<OCCSurface*> (this)->geometry_type() != SPHERE_SURFACE_TYPE)
    return CUBIT_FAILURE;

  BRepAdaptor_Surface asurface(*myTopoDSFace);
  gp_Sphere sphere =  asurface.Sphere(); 
  gp_Pnt center_pt = sphere.Location();
  center = CubitVector(center_pt.X(), center_pt.Y(), center_pt.Z());

  radius = sphere.Radius();
  return CUBIT_SUCCESS;
}
TopoDS_Face* OCCSurface::get_TopoDS_Face ( ) [inline]

Definition at line 321 of file OCCSurface.hpp.

  {if (myTopoDSFace && !myTopoDSFace->IsNull()) assert (myTopoDSFace->ShapeType() == TopAbs_FACE); return myTopoDSFace;}
CubitStatus OCCSurface::get_torus_params ( CubitVector &  center,
CubitVector &  normal,
double &  major_radius,
double &  minor_radius 
) const [virtual]

Implements Surface.

Definition at line 1492 of file OCCSurface.cpp.

{
  if(const_cast<OCCSurface*> (this)->geometry_type() != TORUS_SURFACE_TYPE)
    return CUBIT_FAILURE;

  BRepAdaptor_Surface asurface(*myTopoDSFace);
  gp_Torus torus = asurface.Torus();
  gp_Pnt center_pt = torus.Location();
  center = CubitVector(center_pt.X(), center_pt.Y(), center_pt.Z());

  gp_Ax1 axis = torus.Axis();
  gp_Dir dir = axis.Direction();
  normal.set(dir.X(), dir.Y(), dir.Z());

  major_radius = torus.MajorRadius();
  minor_radius = torus.MinorRadius();

  return CUBIT_SUCCESS;
}
CubitBoolean OCCSurface::is_closed_in_U ( ) [virtual]

Reimplemented from Surface.

Definition at line 668 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  return (asurface.IsUClosed())?CUBIT_TRUE:CUBIT_FALSE;
}
CubitBoolean OCCSurface::is_closed_in_V ( ) [virtual]

Reimplemented from Surface.

Definition at line 677 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  return (asurface.IsVClosed())?CUBIT_TRUE:CUBIT_FALSE;
}
CubitBoolean OCCSurface::is_conical ( )
CubitBoolean OCCSurface::is_flat ( )
CubitBoolean OCCSurface::is_parametric ( ) [virtual]

Implements Surface.

Definition at line 727 of file OCCSurface.cpp.

{
  return CUBIT_TRUE;
}
CubitBoolean OCCSurface::is_periodic ( ) [virtual]

Implements Surface.

Definition at line 576 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  return (asurface.IsUPeriodic() || asurface.IsVPeriodic())?CUBIT_TRUE:CUBIT_FALSE;
}
CubitBoolean OCCSurface::is_periodic_in_U ( double &  period) [virtual]

Implements Surface.

Definition at line 587 of file OCCSurface.cpp.

{

  BRepAdaptor_Surface asurface(*myTopoDSFace);
  if (!asurface.IsUPeriodic())
     return CUBIT_FALSE;

  period = asurface.UPeriod(); 
  return CUBIT_TRUE;
}
CubitBoolean OCCSurface::is_periodic_in_V ( double &  period) [virtual]

Implements Surface.

Definition at line 604 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  if (!asurface.IsVPeriodic())
     return CUBIT_FALSE;

  period = asurface.VPeriod();
  return CUBIT_TRUE;
}
CubitBoolean OCCSurface::is_position_on ( CubitVector &  test_position) [virtual]

Implements Surface.

Definition at line 797 of file OCCSurface.cpp.

{
  CubitVector new_point;
  CubitStatus stat = closest_point(test_position, &new_point, NULL,NULL,NULL);
  if ( !stat )
    return CUBIT_FALSE;
  CubitVector result_vec = test_position - new_point;
  if ( result_vec.length_squared() < GEOMETRY_RESABS )
    return CUBIT_TRUE;
  return CUBIT_FALSE;
}
CubitBoolean OCCSurface::is_singular_in_U ( double  u_param) [virtual]

Implements Surface.

Definition at line 621 of file OCCSurface.cpp.

{
  //from  MasterIndex.htm:
  // singular_u
  // The only singularity recognized is where every value of the 
  // nonconstant parameter generates the same object-space point,
  // and these can only occur at the ends of the parameter range
  // as returned by the functions above. A plane is nonsingular 
  // in both directions. 
  double u_lower, u_upper;
  get_param_range_U( u_lower, u_upper );

  if ( u_param < u_lower - CUBIT_RESABS ||
       u_param > u_upper + CUBIT_RESABS )
  {
    PRINT_ERROR("u parameter is outside parameter bounds.\n");
    return CUBIT_FALSE;
  }

  //Currently, haven't found any singularity check in OCC.
  return CUBIT_FALSE;
}
CubitBoolean OCCSurface::is_singular_in_V ( double  v_param) [virtual]

Implements Surface.

Definition at line 648 of file OCCSurface.cpp.

{
  double v_lower, v_upper;
  get_param_range_V( v_lower, v_upper );

  if ( v_param < v_lower - CUBIT_RESABS ||
       v_param > v_upper + CUBIT_RESABS )
  {
    PRINT_ERROR("v parameter is outside parameter bounds.\n");
    return CUBIT_FALSE;
  }

  //Currently, haven't found any singularity check in OCC.
  return CUBIT_FALSE;
}
CubitBoolean OCCSurface::is_spherical ( )
double OCCSurface::measure ( ) [virtual]

Implements GeometryEntity.

Definition at line 772 of file OCCSurface.cpp.

{
  GProp_GProps myProps;
  BRepGProp::SurfaceProperties(*myTopoDSFace, myProps);
  return myProps.Mass();
}
OCCBody* OCCSurface::my_body ( ) [inline]

Definition at line 346 of file OCCSurface.hpp.

{return myBody;}
OCCLump* OCCSurface::my_lump ( ) [inline]

Definition at line 344 of file OCCSurface.hpp.

{return myLump;}
OCCShell* OCCSurface::my_shell ( ) [inline]

Definition at line 342 of file OCCSurface.hpp.

{return myShell;}
CubitPointContainment OCCSurface::point_containment ( const CubitVector &  point) [virtual]

Implements Surface.

Definition at line 809 of file OCCSurface.cpp.

{
   TopoDS_Face *face = get_TopoDS_Face();
   gp_Pnt p(point.x(), point.y(), point.z());
   double tol = OCCQueryEngine::instance()->get_sme_resabs_tolerance();

   //It's checking the state of the projected point of THIS Point
   BRepClass_FaceClassifier face_classifier;
   face_classifier.Perform(*face, p, tol);
   TopAbs_State state = face_classifier.State();
   
   if (state == TopAbs_IN)
     return CUBIT_PNT_INSIDE;
   else if (state == TopAbs_OUT)
     return CUBIT_PNT_OUTSIDE;
   else if (state == TopAbs_ON)
     return CUBIT_PNT_BOUNDARY;

   return CUBIT_PNT_UNKNOWN;
}
CubitPointContainment OCCSurface::point_containment ( double  u,
double  v 
) [virtual]

Implements Surface.

Definition at line 830 of file OCCSurface.cpp.

{
  CubitVector point = position_from_u_v(u_param, v_param);
  return point_containment(point);
}
CubitVector OCCSurface::position_from_u_v ( double  u,
double  v 
) [virtual]

Implements Surface.

Definition at line 532 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace, Standard_False);
  gp_Pnt p = asurface.Value(u, v);
  return CubitVector (p.X(), p.Y(), p.Z());
}
CubitStatus OCCSurface::principal_curvatures ( CubitVector const &  location,
double &  curvature_1,
double &  curvature_2,
CubitVector *  closest_location = NULL 
) [virtual]

Implements Surface.

Definition at line 446 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  gp_Pnt p(location.x(), location.y(), location.z()), newP(0.0, 0.0, 0.0);
  double minDist=0.0, u, v;
  int i;
  BRepLProp_SLProps SLP(asurface, 2, Precision::PConfusion());
  Extrema_ExtPS ext(p, asurface, Precision::Approximation(), Precision::Approximation());
  if (ext.IsDone() && (ext.NbExt() > 0)) {
      for ( i = 1 ; i <= ext.NbExt() ; i++ ) {
          if ( (i==1) || (p.Distance(ext.Point(i).Value()) < minDist) ) {
              minDist = p.Distance(ext.Point(i).Value());
              newP = ext.Point(i).Value();
              ext.Point(i).Parameter(u, v);
              SLP.SetParameters(u, v);
          }
      }
  }
  if (closest_location != NULL)
    *closest_location = CubitVector(newP.X(), newP.Y(), newP.Z());

  if (SLP.IsCurvatureDefined())
  {
    curvature_1 = SLP.MinCurvature();
    curvature_2 = SLP.MaxCurvature();
  }
  return CUBIT_SUCCESS;
}
void OCCSurface::remove_all_simple_attribute_virt ( ) [virtual]

Implements TopologyBridge.

Definition at line 167 of file OCCSurface.cpp.

{
  OCCAttribSet::remove_attribute(CubitSimpleAttrib(), *myTopoDSFace);
}
void OCCSurface::remove_hardpoint ( OCCPoint *  HardPoint) [inline]

Definition at line 80 of file OCCSurface.hpp.

{myHardPoints.remove(HardPoint);}
void OCCSurface::remove_simple_attribute_virt ( const CubitSimpleAttrib &  csa) [virtual]

Implements TopologyBridge.

Definition at line 154 of file OCCSurface.cpp.

  { OCCAttribSet::remove_attribute( csa , *myTopoDSFace); }
void OCCSurface::set_body ( OCCBody *  body) [inline]

Definition at line 339 of file OCCSurface.hpp.

  { myBody = body;}
void OCCSurface::set_lump ( OCCLump *  lump) [inline]

Definition at line 336 of file OCCSurface.hpp.

  { myLump = lump;}
void OCCSurface::set_shell ( OCCShell *  shell) [inline]

Definition at line 333 of file OCCSurface.hpp.

  { myShell = shell;}
void OCCSurface::set_TopoDS_Face ( TopoDS_Face &  face)

Definition at line 122 of file OCCSurface.cpp.

{
  if(myTopoDSFace && face.IsEqual(*myTopoDSFace))
    return;

  if(myTopoDSFace)
    myTopoDSFace->Nullify();
  *myTopoDSFace = face ;
}
CubitStatus OCCSurface::u_v_from_position ( CubitVector const &  location,
double &  u,
double &  v,
CubitVector *  closest_location = NULL 
) [virtual]

Implements Surface.

Definition at line 547 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  gp_Pnt p(location.x(), location.y(), location.z()), newP(0.0, 0.0, 0.0);
  double minDist=0.0;
  int i;
  Extrema_ExtPS ext(p, asurface, Precision::Confusion(), Precision::Confusion());
  if (ext.IsDone() && (ext.NbExt() > 0)) {
      for ( i = 1 ; i <= ext.NbExt() ; i++ ) {
          if ( (i==1) || (p.Distance(ext.Point(i).Value()) < minDist) ) {
              minDist = p.Distance(ext.Point(i).Value());
              newP = ext.Point(i).Value();
              ext.Point(i).Parameter(u, v);
          }
      }
  }
  if (closest_location != NULL) *closest_location = CubitVector(newP.X(), newP.Y(), newP.Z());
  return CUBIT_SUCCESS;
}
CubitStatus OCCSurface::update_OCC_entity ( TopoDS_Face &  old_surface,
TopoDS_Shape &  new_surface,
BRepBuilderAPI_MakeShape *  op,
TopoDS_Vertex *  removed_vertex = NULL,
LocOpe_SplitShape *  sp = NULL 
) [static]

Definition at line 1059 of file OCCSurface.cpp.

{
  TopTools_IndexedMapOfShape M, M2;
  TopoDS_Shape shape_face, shape, shape2, shape_edge, shape_vertex;
  double dTOL = OCCQueryEngine::instance()->get_sme_resabs_tolerance();

  //First check if new_surface is type shell.
  TopExp::MapShapes(new_surface, TopAbs_FACE,M);
  if(M.Extent() > 1 )
  {
    OCCQueryEngine::instance()->update_OCC_map(old_surface, new_surface);
    return CUBIT_SUCCESS;
  }

  
  if(M.Extent() == 1 )
    shape_face = M(1);

      //     GfxDebug::clear();
      //    OCCDrawTool::instance()->draw_TopoDS_Shape(&shape_face,CUBIT_GREEN,false,true);
      //    GfxDebug::mouse_xforms();


  M.Clear(); 
  //set the Wires
  TopExp::MapShapes(old_surface, TopAbs_WIRE, M);

  TopTools_ListOfShape shapes;  
  BRepFilletAPI_MakeFillet2d* test_op = NULL;

  for (int ii=1; ii<=M.Extent(); ii++) 
  {
     TopoDS_Wire wire = TopoDS::Wire(M(ii));
     TopTools_ListOfShape shapes;
     if(op)
     {
       //test_op = dynamic_cast<BRepFilletAPI_MakeFillet2d*>(op);
       test_op = NULL; //casting fails on both OSX and Linux
       if(!test_op)
         shapes.Assign(op->Modified(wire));
       if(shapes.Extent() == 0)
         shapes.Assign(op->Generated(wire));
       if(!new_surface.IsNull())
         TopExp::MapShapes(new_surface,TopAbs_WIRE, M2);
     }
     else if(sp)
       shapes.Assign(sp->DescendantShapes(wire));

     if (shapes.Extent() == 1)
     {
       shape = shapes.First();
       if(M2.Extent() == 1)
       {
         shape2 = TopoDS::Wire(M2(1));
         if(!shape.IsSame(shape2))
           shape = shape2;
       }
       else if(M2.Extent() > 1)
         shape.Nullify();
     }
     else if(shapes.Extent() > 1)
       shape.Nullify();
     else if(op->IsDeleted(wire) || shapes.Extent() == 0)
     {
       TopTools_IndexedMapOfShape M_new;
       TopExp::MapShapes(new_surface, TopAbs_WIRE, M_new);
       if (M_new.Extent()== 1)
         shape = M_new(1);
       else if(!shape_face.IsNull())
       {
         M_new.Clear();
         TopExp::MapShapes(shape_face, TopAbs_WIRE, M_new);
         if (M_new.Extent()== 1)
           shape = M_new(1);
       }
       else 
         shape.Nullify();
     }
     else
     {
       shape = wire;
       continue;
     }

     //set curves
     BRepTools_WireExplorer Ex;
      
     for(Ex.Init(wire); Ex.More();Ex.Next())
     {
       TopoDS_Edge edge = Ex.Current();
       //check to see if the edge made itself into a curve.
       GProp_GProps myProps;
       BRepGProp::LinearProperties(edge, myProps);
       double length =  myProps.Mass();
       if(length < dTOL)
         continue;
       if(op && !test_op)
       {
         shapes.Assign(op->Modified(edge));
         if(shapes.Extent() == 0)
           shapes.Assign(op->Generated(edge));
       }
         
       else if(sp)
         shapes.Assign(sp->DescendantShapes(edge));

       if (shapes.Extent() == 1)
       {
        //in fillet creating mothod, one edge could generated a face, so check
        //it here.
         TopAbs_ShapeEnum type = shapes.First().TShape()->ShapeType(); 
         if(type != TopAbs_EDGE)
           shape_edge.Nullify();
         else
           shape_edge = shapes.First();
       }
       else if (shapes.Extent() > 1)
       {
         //update all attributes first.
         TopTools_ListIteratorOfListOfShape it;
         it.Initialize(shapes);
         for(; it.More(); it.Next())
         {
           shape_edge = it.Value();
           OCCQueryEngine::instance()->copy_attributes(edge, shape_edge);
         }
         shape_edge.Nullify();
       }
       else if (op->IsDeleted(edge))
         shape_edge.Nullify(); 
       else if (test_op)
       {
         if(!test_op->IsModified(edge))
           shape_edge = edge;
         else
           shape_edge = (test_op->Modified(edge)).First();
       } 
       else
         shape_edge = edge;

       //update vertex
       TopoDS_Vertex vertex = Ex.CurrentVertex();
       shapes.Clear();
       if(test_op)
         assert(removed_vertex != NULL);

       if(op && ! test_op )
       {
         shapes.Assign(op->Modified(vertex));
         if(shapes.Extent() == 0)
           shapes.Assign(op->Generated(vertex));
       }
       if(sp)
         shapes.Assign(sp->DescendantShapes(vertex));

       if (shapes.Extent() == 1)
         shape_vertex = shapes.First();

       else if(shapes.Extent() > 1)
       {
         //update all attributes first.
         TopTools_ListIteratorOfListOfShape it;
         it.Initialize(shapes);
         for(; it.More(); it.Next())
         {
           shape_vertex = it.Value();
           OCCQueryEngine::instance()->copy_attributes(vertex, shape_vertex);
         }
         shape_vertex.Nullify() ;
       }
       else if(op->IsDeleted(vertex) || (test_op && vertex.IsSame( *removed_vertex)))
       {
     if(!shape.IsNull() && !shape_edge.IsNull() && !shape_edge.Closed()) 
         //there should be a vertex corresponding to the old_vertex.
         //find the vertices within tolerance distance with old_vertex.
         {
           TopoDS_Iterator It(shape_edge);
           for(; It.More(); It.Next())
           {
             TopoDS_Vertex v = TopoDS::Vertex(It.Value());
             gp_Pnt pt1 = BRep_Tool::Pnt(v);
             gp_Pnt pt2 = BRep_Tool::Pnt(vertex);
             if(pt1.IsEqual(pt2, dTOL))
             {
               shape_vertex = v;  
               break;
             }
           }
         }
         else   
           shape_vertex.Nullify() ;
       } 
       else
         shape_vertex = vertex;
      
       if(!vertex.IsSame(shape_vertex) )
         OCCQueryEngine::instance()->update_OCC_map(vertex, shape_vertex);

       if (!edge.IsSame(shape_edge))
         OCCQueryEngine::instance()->update_OCC_map(edge, shape_edge);
     }
     if (!wire.IsSame(shape))
       OCCQueryEngine::instance()->update_OCC_map(wire, shape);
  }

  if (!old_surface.IsSame(new_surface))
  {
    TopAbs_ShapeEnum shapetype =  TopAbs_SHAPE;
    if(!new_surface.IsNull())
      shapetype = new_surface.TShape()->ShapeType();  
    if(shapetype == TopAbs_FACE || new_surface.IsNull())
      OCCQueryEngine::instance()->update_OCC_map(old_surface, new_surface);
    else 
    {
      TopTools_IndexedMapOfShape M;
      TopExp::MapShapes(new_surface, TopAbs_FACE, M);   
      TopoDS_Shape new_shape;
      if(M.Extent() == 1)
        new_shape = M(1);
      else if(M.Extent() > 1)
      {
        for(int i = 1; i <= M.Extent(); i++)
        {
          GProp_GProps myProps;
          BRepGProp::SurfaceProperties(old_surface, myProps);
          double orig_mass = myProps.Mass();
          gp_Pnt orig_pnt = myProps.CentreOfMass();
          BRepGProp::SurfaceProperties(M(i), myProps);
          double after_mass = myProps.Mass();
          gp_Pnt after_pnt = myProps.CentreOfMass();
          if(fabs(-after_mass + orig_mass) <= dTOL && 
             orig_pnt.IsEqual(after_pnt, dTOL))
          {
            new_shape = M(i);
            break;
          }
        }
      }
      OCCQueryEngine::instance()->update_OCC_map(old_surface, new_shape);
    }
  }
  return CUBIT_SUCCESS;
}
CubitStatus OCCSurface::update_OCC_entity ( BRepBuilderAPI_ModifyShape *  aBRepTrsf,
BRepAlgoAPI_BooleanOperation *  op = NULL 
)

Definition at line 1004 of file OCCSurface.cpp.

{
  assert(aBRepTrsf != NULL || op != NULL);
  TopoDS_Shape shape;

  if(aBRepTrsf)
    shape = aBRepTrsf->ModifiedShape(*get_TopoDS_Face());

  else
  {
    TopTools_ListOfShape shapes;
    shapes.Assign(op->Modified(*get_TopoDS_Face()));
    if(shapes.Extent() == 0)
         shapes.Assign(op->Generated(*get_TopoDS_Face()));
    if (shapes.Extent() == 1)
      shape = shapes.First();
    else if(shapes.Extent() > 1)
    {
      shape = shapes.First();
    }
    else if(op->IsDeleted(*get_TopoDS_Face()))
      ;
    else
      return CUBIT_SUCCESS;
  }
 
  TopoDS_Face surface; 
  if(!shape.IsNull())
    surface = TopoDS::Face(shape);

  if (aBRepTrsf) 
  {
    //set the loops
    DLIList<OCCLoop *> loops;
    this->get_loops(loops);
    for (int i = 1; i <= loops.size(); i++)
    {
       OCCLoop *loop = loops.get_and_step();
       loop->update_OCC_entity(aBRepTrsf, op);
    }
    OCCQueryEngine::instance()->update_OCC_map(*myTopoDSFace, surface);
  }

  else if(op)
    update_OCC_entity(*myTopoDSFace, surface, op);

  return CUBIT_SUCCESS;
}
CubitStatus OCCSurface::uv_derivitives ( double  u_param,
double  v_param,
CubitVector &  du,
CubitVector &  dv 
) [virtual]

Implements Surface.

Definition at line 687 of file OCCSurface.cpp.

{
  BRepAdaptor_Surface asurface(*myTopoDSFace);
  gp_Pnt p;
  gp_Vec d1u, d1v;
  asurface.D1(u, v, p, d1u, d1v);
  du = CubitVector(d1u.X(), d1u.Y(), d1u.Z());
  dv = CubitVector(d1v.X(), d1v.Y(), d1v.Z());
  return CUBIT_SUCCESS;
}

Member Data Documentation

OCCBody* OCCSurface::myBody [private]

Definition at line 529 of file OCCSurface.hpp.

DLIList<OCCPoint*> OCCSurface::myHardPoints [private]

Definition at line 530 of file OCCSurface.hpp.

OCCLump* OCCSurface::myLump [private]

Definition at line 528 of file OCCSurface.hpp.

OCCShell* OCCSurface::myShell [private]

Definition at line 527 of file OCCSurface.hpp.

TopoDS_Face* OCCSurface::myTopoDSFace [private]

Definition at line 524 of file OCCSurface.hpp.

The documentation for this class was generated from the following files:
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines