#include <OCCCurve.hpp>

Inheritance diagram for OCCCurve:

Public Member Functions
	OCCCurve (TopoDS_Edge *theEdge)
virtual	~OCCCurve ()
void	set_myMarked (CubitBoolean marked)
void	add_loop (OCCLoop *loop)
DLIList< OCCLoop * >	loops ()
void	remove_loop (OCCLoop *loop)
void	clean_loops ()
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 double	measure ()
virtual double	length_from_u (double parameter1, double parameter2)
virtual CubitBoolean	is_periodic (double &period)
virtual CubitBoolean	get_param_range (double &lower_bound, double &upper_bound)
virtual CubitStatus	get_interior_extrema (DLIList< CubitVector * > &interior_points, CubitSense &return_sense)
CubitStatus	get_interior_extrema_in_direction (DLIList< CubitVector * > &, CubitVector)
virtual CubitStatus	closest_point (CubitVector const &location, CubitVector &closest_location, CubitVector tangent_ptr=NULL, CubitVector curvature_ptr=NULL, double *param=NULL)
void	get_tangent (CubitVector const &location, CubitVector &tangent)
void	get_curvature (CubitVector const &location, CubitVector &curvature)
virtual CubitStatus	position_from_u (double u_value, CubitVector &output_position)
double	u_from_position (const CubitVector &input_position)
virtual double	u_from_arc_length (double root_param, double arc_length)
virtual CubitBoolean	is_position_on (const CubitVector &test_position)
GeometryType	geometry_type ()
CubitStatus	get_point_direction (CubitVector &origin, CubitVector &direction)
CubitStatus	get_center_radius (CubitVector &center, double &radius)
virtual double	start_param ()
virtual double	end_param ()
virtual CubitBoolean	G1_discontinuous (double param, CubitVector minus_tangent=NULL, CubitVector plus_tangent=NULL)
virtual CubitPointContainment	point_containment (const CubitVector &point)
virtual void	get_parents_virt (DLIList< TopologyBridge * > &parents)
virtual void	get_children_virt (DLIList< TopologyBridge * > &children)
virtual CubitStatus	get_spline_params (bool &rational, int &degree, DLIList< CubitVector > &cntrl_pts, DLIList< double > &cntrl_pt_weights, DLIList< double > &knots) const
virtual CubitStatus	get_ellipse_params (CubitVector &center, CubitVector &normal, CubitVector &major_axis, double &radius_ratio) const
void	get_points (DLIList< OCCPoint * > &point_list)
TopoDS_Edge *	get_TopoDS_Edge ()
void	set_TopoDS_Edge (TopoDS_Edge edge)
void	update_OCC_entity (BRepBuilderAPI_ModifyShape aBRepTrsf, BRepAlgoAPI_BooleanOperation op=NULL)
Curve *	project_curve (Surface face_ptr, DLIList< TBPoint > &normal_proj_points, CubitBoolean closed, const CubitVector *third_point)
virtual CubitStatus	get_spline_params (bool &rational, int &degree, DLIList< CubitVector > &cntrl_pts, DLIList< double > &cntrl_pt_weights, DLIList< double > &knots, bool &spline_is_reversed) const
Private Member Functions
void	adjust_periodic_parameter (double &param)
Private Attributes
TopoDS_Edge *	myTopoDSEdge
DLIList< OCCLoop * >	myLoopList
CubitBoolean	myMarked

Detailed Description

Definition at line 41 of file OCCCurve.hpp.

Constructor & Destructor Documentation

OCCCurve::OCCCurve ( TopoDS_Edge * theEdge )

Definition at line 99 of file OCCCurve.cpp.

{
  myTopoDSEdge = theEdge;
  myMarked = CUBIT_FALSE;
  assert (myTopoDSEdge->ShapeType() == TopAbs_EDGE);
}

OCCCurve::~OCCCurve ( ) [virtual]

Definition at line 115 of file OCCCurve.cpp.

{
  if (myTopoDSEdge)
  {
    myTopoDSEdge->Nullify();
    delete (TopoDS_Edge *)myTopoDSEdge;
    myTopoDSEdge = NULL;
  }
}

Member Function Documentation

void OCCCurve::add_loop ( OCCLoop * loop )

Definition at line 134 of file OCCCurve.cpp.

{ 
  //before add the loop into the looplist, check to make sure the 
  //looplist is up-to-date
  //It should be done somewhere else, however, in large mcnp2cad test, there's
  //no way to check where the code forget to remove the out-dated loop.
  for (int i = 0; i < myLoopList.size(); i++)
  {
    OCCLoop* myLoop = CAST_TO(myLoopList.get(), OCCLoop);

    if(!myLoop) 
      this->remove_loop(myLoopList.get());
    else if(myLoop->get_TopoDS_Wire() < (void*) 0x1000)
      this->remove_loop(myLoop);
    else if(myLoop->coedges().size() == 0)
      this->remove_loop(myLoop);
    myLoopList.step();
  }
  if(loop != NULL)
    myLoopList.append_unique(loop);
}

void OCCCurve::adjust_periodic_parameter ( double & param ) [private]

Definition at line 932 of file OCCCurve.cpp.

{
    // Adjustment only legal if this is a periodic curve.
  double period;
  if ( this->is_periodic(period) && (fabs(period) > CUBIT_RESABS))
  {
    double upper_bound, lower_bound;
    this->get_param_range( lower_bound, upper_bound );
    assert((upper_bound - lower_bound) > CUBIT_RESABS * 100);

    lower_bound -= CUBIT_RESABS;
    upper_bound += CUBIT_RESABS;
    
      // Make sure period is positive
    if (period < 0.)
      period = -period;

      // Move the parameter above the low param
    while (param < lower_bound)
      param += period;
      // Move the parameter below the high param
    while (param > upper_bound)
      param -= period;
  }
}

void OCCCurve::append_simple_attribute_virt ( const CubitSimpleAttrib & csa ) [virtual]

Implements TopologyBridge.

Definition at line 202 of file OCCCurve.cpp.

  { OCCAttribSet::append_attribute(csa, *myTopoDSEdge); }

CubitBox OCCCurve::bounding_box ( void ) const [virtual]

Implements GeometryEntity.

Definition at line 278 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  Bnd_Box aBox;
  BndLib_Add3dCurve::Add(acurve, 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);
}

void OCCCurve::clean_loops ( ) [inline]

Definition at line 55 of file OCCCurve.hpp.

{myLoopList.clean_out();}

CubitStatus OCCCurve::closest_point	(	CubitVector const &	location,
		CubitVector &	closest_location,
		CubitVector *	tangent_ptr = `NULL`,
		CubitVector *	curvature_ptr = `NULL`,
		double *	param = `NULL`
	)		`[virtual]`

Implements Curve.

Definition at line 483 of file OCCCurve.cpp.

{  
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  gp_Pnt p(location.x(), location.y(), location.z()), newP(0.0, 0.0, 0.0);
  Extrema_ExtPC ext(p, acurve, Precision::Approximation());
  if (!ext.IsDone())
    return CUBIT_FAILURE;

  double sqr_dist = CUBIT_DBL_MAX;
  double pparam = 0.0;
  for (int i = 1; i <= ext.NbExt(); ++i) {
      double new_sqr_dist = p.SquareDistance( ext.Point(i).Value() );
      if (new_sqr_dist < sqr_dist) {
        sqr_dist = new_sqr_dist;
        newP = ext.Point(i).Value();
        pparam = ext.Point(i).Parameter();
      }
  }

    // if we didn't find any minimum...
  if (sqr_dist == CUBIT_DBL_MAX)
    return CUBIT_FAILURE;

    // pass back closest point
  closest_location.set( newP.X(), newP.Y(), newP.Z() );
  if (param != NULL)
    *param = pparam;

    // pass back tangent
  if (tangent_ptr != NULL) {
    BRepLProp_CLProps CLP(acurve, 2, Precision::PConfusion());
    CLP.SetParameter( pparam );
    if (!CLP.IsTangentDefined())
      return CUBIT_FAILURE;

    gp_Dir tangent;
    CLP.Tangent(tangent);
    tangent_ptr->set( tangent.X(), tangent.Y(), tangent.Z() );
  }

    // if requested, pass back curvature
  if (curvature_ptr != NULL)
    get_curvature( closest_location, *curvature_ptr );

  return CUBIT_SUCCESS;
}

double OCCCurve::end_param ( ) [virtual]

Implements Curve.

Definition at line 781 of file OCCCurve.cpp.

{
   double start = 0.0, end = 0.0;
   
   get_param_range( start, end );
   return end;
}

CubitBoolean OCCCurve::G1_discontinuous	(	double	param,
		CubitVector *	minus_tangent = `NULL`,
		CubitVector *	plus_tangent = `NULL`
	)		`[virtual]`

Reimplemented from Curve.

Definition at line 846 of file OCCCurve.cpp.

{ 
  CubitBoolean is_discon = CUBIT_TRUE;
  double first, last;
  Handle(Geom_Curve) gCurve = BRep_Tool::Curve(*myTopoDSEdge, first, last);

  if (gCurve->Continuity() < GeomAbs_G1)
     return is_discon;

  assert(first <= param && param <= last );
  
  gp_Pnt P;
  gp_Vec V1;
  gCurve->D1(param, P, V1);
  
  mtan = new CubitVector(V1.X(), V1.Y(),V1.Z());
  ptan = new CubitVector(*mtan);
     
  return CUBIT_FALSE;
}

GeometryType OCCCurve::geometry_type ( ) [virtual]

Reimplemented from Curve.

Definition at line 662 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  if (acurve.GetType() == GeomAbs_BezierCurve)
     return SPLINE_CURVE_TYPE;
  if (acurve.GetType() == GeomAbs_BSplineCurve)
     return SPLINE_CURVE_TYPE;
  if (acurve.GetType() == GeomAbs_Line)
     return STRAIGHT_CURVE_TYPE;
  if (acurve.GetType() == GeomAbs_Parabola)
     return PARABOLA_CURVE_TYPE;
  if (acurve.GetType() == GeomAbs_Hyperbola)
     return HYPERBOLA_CURVE_TYPE;
  if (acurve.GetType() == GeomAbs_Circle)
      return ARC_CURVE_TYPE;
  if (acurve.GetType() == GeomAbs_Ellipse)
     return ELLIPSE_CURVE_TYPE;
  return UNDEFINED_CURVE_TYPE;
}

CubitStatus OCCCurve::get_center_radius	(	CubitVector &	center,
		double &	radius
	)		`[virtual]`

Implements Curve.

Definition at line 725 of file OCCCurve.cpp.

{
  if( geometry_type() != ELLIPSE_CURVE_TYPE &&
      geometry_type() != ARC_CURVE_TYPE )
    return CUBIT_FAILURE;
 
  //get the Geom_Curve of the OCCCurve
  Standard_Real first;
  Standard_Real last;
  Handle(Geom_Curve) myCurve = BRep_Tool::Curve(*myTopoDSEdge, first, last); 

  if (Handle(Geom_Circle) gCircle = Handle(Geom_Circle)::DownCast(myCurve))
  {
     radius = gCircle->Radius();
     gp_Circ gp_circ = gCircle->Circ();
     gp_Pnt  gp_p = gp_circ.Location();
     center.set(gp_p.X(), gp_p.Y(), gp_p.Z());
  }

  else //ellipse
  {
     Handle(Geom_Ellipse) gEllipse = Handle(Geom_Ellipse)::DownCast(myCurve);
     radius = gEllipse->MajorRadius();
     gp_Elips gp_ellip = gEllipse->Elips();
     gp_Pnt  gp_p = gp_ellip.Location();
     center.set(gp_p.X(), gp_p.Y(), gp_p.Z());
  }
  return CUBIT_SUCCESS;
}

void OCCCurve::get_children_virt ( DLIList< TopologyBridge * > & children ) [virtual]

Implements TopologyBridge.

Definition at line 810 of file OCCCurve.cpp.

{
    TopTools_IndexedMapOfShape M;
    TopExp::MapShapes(*myTopoDSEdge, TopAbs_VERTEX, M);
        int ii ;
        for (ii = M.Extent();  ii > 0; ii--) {
             TopologyBridge *point = OCCQueryEngine::instance()->occ_to_cgm(M(ii));
              if (point)
                children.append_unique(point);
       }
       if(children.size() ==0)
         return;
       //make sure the points are in the start-end order
       double start , end;     
       get_param_range( start, end );
       CubitVector v;
       position_from_u(start, v);
       children.reset();
       if(!v.about_equal(CAST_TO(children.get(), TBPoint)->coordinates())) 
         children.reverse();
}

void OCCCurve::get_curvature	(	CubitVector const &	location,
		CubitVector &	curvature
	)

Definition at line 905 of file OCCCurve.cpp.

{  
    double u = u_from_position(location);
    Standard_Real first;
    Standard_Real last;
    Handle(Geom_Curve) myCurve = BRep_Tool::Curve(*myTopoDSEdge, first, last);

    gp_Pnt p;
    gp_Vec tan, cur;
    GeomLProp_CurveTool::D2(myCurve, u , p, tan, cur) ;
    curvature.set(cur.X(), cur.Y(), cur.Z()); 
}

CubitStatus OCCCurve::get_ellipse_params	(	CubitVector &	center,
		CubitVector &	normal,
		CubitVector &	major_axis,
		double &	radius_ratio
	)		const `[virtual]`

Implements Curve.

Definition at line 1021 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  gp_Elips ellipse;
  if (acurve.GetType() == GeomAbs_Ellipse)
    ellipse = acurve.Ellipse();
  else
    return CUBIT_FAILURE;

  gp_Pnt center = ellipse.Location();
  center_pt = CubitVector(center.X(), center.Y(), center.Z());
  
  gp_Ax1 normal_Axis = ellipse.Axis();
  gp_Dir normal_dir = normal_Axis.Direction();
  normal = CubitVector(normal_dir.X(), normal_dir.Y(), normal_dir.Z());

  gp_Ax1 major_Axis = ellipse.Directrix1();
  gp_Dir major_dir = major_Axis.Direction();
  major_axis = CubitVector(major_dir.X(), major_dir.Y(), major_dir.Z());

  double major = ellipse.MajorRadius();
  double minor = ellipse.MinorRadius();
  radius_ratio = major/minor;

  return CUBIT_SUCCESS;
}

GeometryQueryEngine * OCCCurve::get_geometry_query_engine ( ) const [virtual]

Implements TopologyBridge.

Definition at line 264 of file OCCCurve.cpp.

{
  return OCCQueryEngine::instance();
}

CubitStatus OCCCurve::get_interior_extrema	(	DLIList< CubitVector * > &	interior_points,
		CubitSense &	return_sense
	)		`[virtual]`

Implements Curve.

Definition at line 377 of file OCCCurve.cpp.

{
  DLIList<CubitVector*> point_list;
  CubitVector x(1.0, 0.0, 0.0);
  get_interior_extrema_in_direction(point_list, x);
  CubitVector y(0.0, 1.0, 0.0);
  get_interior_extrema_in_direction(point_list, y);
  CubitVector z(0.0, 0.0, 1.0);
  get_interior_extrema_in_direction(point_list, z );

  //like , return only points that aren't at an endpoint and are not
  //close to previous point
  const double epsilon = 30.* GEOMETRY_RESABS;
  const double epsilon_squared = epsilon*epsilon;

  //get both vertices' coordinates.
  CubitVector endpoints[2];
  int i = 0;
  TopExp_Explorer aVertexExp(*myTopoDSEdge, TopAbs_VERTEX);
  while(aVertexExp.More())
  {
     TopoDS_Vertex v = TopoDS::Vertex(aVertexExp.Current());
     gp_Pnt p = BRep_Tool::Pnt(v);
     endpoints[i].x(p.X());
     endpoints[i].y(p.Y());
     endpoints[i].z(p.Z());  
     i++;
     aVertexExp.Next();
  } 

  //compare to see if the Points in point_list are interior and far apart
  int j;
  CubitVector* cubit_position = NULL;
  CubitVector * temp_position = NULL;
  point_list.sort();
  point_list.reset();
  for (j = point_list.size(); j--; )
  {
     temp_position = point_list.get_and_step();
     // save if not equal to an endpoint, or prior point
     if (temp_position->distance_between_squared(endpoints[0]) > epsilon_squared
        &&
        temp_position->distance_between_squared(endpoints[1]) > epsilon_squared)
     {
        if (!cubit_position ||
            temp_position->distance_between_squared(*cubit_position) > epsilon_squared)
        {
          cubit_position = temp_position ;
          interior_points.append( cubit_position );
        } // If point isn't close to previous point
     } // If point isn't at an endpoint
  } // for each point

  // Return sense is whatever the sense of this curve is.
  TopAbs_Orientation sense = myTopoDSEdge->Orientation();
  return_sense = (sense == TopAbs_FORWARD ? CUBIT_FORWARD : CUBIT_REVERSED);
 
  return CUBIT_SUCCESS;
}

CubitStatus OCCCurve::get_interior_extrema_in_direction	(	DLIList< CubitVector * > &	interior_points,
		CubitVector	dir
	)

Definition at line 439 of file OCCCurve.cpp.

{
  //Create a straight line.
  gp_Pnt origin(0.0, 0.0, 0.0);
  gp_Dir adir(dir.x(), dir.y(), dir.z());
  Handle(Geom_Line) line = GC_MakeLine(origin, adir);

  //get the Geom_Curve of the OCCCurve
  Standard_Real first;
  Standard_Real last;
  Handle(Geom_Curve) myCurve = BRep_Tool::Curve(*myTopoDSEdge, first, last);
  
  GeomAPI_ExtremaCurveCurve extrema(myCurve, line);
  int nPnt = extrema.NbExtrema();
  for (int i = 1; i <= nPnt ; i++)
  {
    gp_Pnt P1, P2;
    extrema.Points(i, P1, P2);
    CubitVector* v = new CubitVector(P1.X(), P1.Y(), P1.Z());
    interior_points.append(v);
  }
  return CUBIT_SUCCESS;
}

CubitBoolean OCCCurve::get_param_range	(	double &	lower_bound,
		double &	upper_bound
	)		`[virtual]`

Implements Curve.

Definition at line 356 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  lower_bound = acurve.FirstParameter();
  upper_bound = acurve.LastParameter();
  return CUBIT_TRUE;
}

void OCCCurve::get_parents_virt ( DLIList< TopologyBridge * > & parents ) [virtual]

Implements TopologyBridge.

Definition at line 790 of file OCCCurve.cpp.

{ 
   //check to see all myLoops are up-to-date.
   add_loop(NULL);

   for(int i = 0; i < myLoopList.size(); i++) 
   {
      DLIList<OCCCoEdge*> coedges = myLoopList.get_and_step()->coedges();
    
      for(int j = 0; j < coedges.size(); j++)
      {
        OCCCoEdge * coedge = coedges.get_and_step();
    if(coedge->curve() == this)
    {
      parents.append(coedge);
          break;
    }
      }
   }
}

CubitStatus OCCCurve::get_point_direction	(	CubitVector &	origin,
		CubitVector &	direction
	)		`[virtual]`

Reimplemented from Curve.

Definition at line 691 of file OCCCurve.cpp.

{
  if (geometry_type() != STRAIGHT_CURVE_TYPE)
    return CUBIT_FAILURE;

  //get the underlying geometry curve
  double first,last;
  Handle(Geom_Curve) gCurve = BRep_Tool::Curve(*myTopoDSEdge, first, last);

  //get the origin and direction of the underlying curve
  Handle(Geom_Line) gLine = Handle(Geom_Line)::DownCast(gCurve);
  gp_Ax1 axis = gLine->Position();
  gp_Pnt loc = axis.Location();
  gp_Dir dir = axis.Direction();
  point.set(loc.X(), loc.Y(), loc.Z());

  //Based on the TopoDS_Edge's orientation, give the unit vector.
  if (myTopoDSEdge->Orientation() == TopAbs_FORWARD)
    direction.set(dir.X(), dir.Y(), dir.Z());
  else if(myTopoDSEdge->Orientation() == TopAbs_REVERSED)
    direction.set(-dir.X(), -dir.Y(), -dir.Z());
  return CUBIT_SUCCESS;  
}

void OCCCurve::get_points ( DLIList< OCCPoint * > & point_list )

Definition at line 868 of file OCCCurve.cpp.

{
  TopTools_IndexedMapOfShape M;
  TopExp::MapShapes(*myTopoDSEdge, TopAbs_VERTEX, M);
  int ii;
  for (ii=M.Extent(); ii>0; ii--) {
      TopologyBridge *point = OCCQueryEngine::instance()->occ_to_cgm(M(ii));
          if (point)
        result_list.append_unique(dynamic_cast<OCCPoint*>(point));
  }
  if(result_list.size() ==0)
    return;
  //make sure the points are in the start-end order
  double start , end;     
  get_param_range( start, end );
  CubitVector v;
  position_from_u(start, v);
  result_list.reset();
  if(!v.about_equal(result_list.get()->coordinates()))
    result_list.reverse();

}

CubitStatus OCCCurve::get_simple_attribute ( DLIList< CubitSimpleAttrib > & csa_list ) [virtual]

Implements TopologyBridge.

Definition at line 247 of file OCCCurve.cpp.

  { return OCCAttribSet::get_attributes(*myTopoDSEdge, csa_list); }

CubitStatus OCCCurve::get_simple_attribute	(	const CubitString &	name,
		DLIList< CubitSimpleAttrib > &	csa_list
	)		`[virtual]`

Implements TopologyBridge.

Definition at line 251 of file OCCCurve.cpp.

  { return OCCAttribSet::get_attributes( name, *myTopoDSEdge, csa_list ); }

CubitStatus OCCCurve::get_spline_params	(	bool &	rational,
		int &	degree,
		DLIList< CubitVector > &	cntrl_pts,
		DLIList< double > &	cntrl_pt_weights,
		DLIList< double > &	knots
	)		const `[virtual]`

Definition at line 968 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  Handle_Geom_BSplineCurve h_S = NULL;
  if (acurve.GetType() == GeomAbs_BSplineCurve)
    h_S = acurve.BSpline();
  else
    return CUBIT_FAILURE;
  assert ( h_S != NULL);

  rational = h_S->IsRational();
  degree =   h_S->Degree();
  
  TColStd_Array1OfReal K(1, h_S->NbKnots());
  h_S->Knots(K);
  for (int i = K.Lower(); i <= K.Upper(); i++)
    knots.append(K.Value(i));

  TColgp_Array1OfPnt P(1, h_S->NbPoles());
  h_S->Poles(P);
  for (int i = P.Lower(); i <= P.Upper(); i++) 
  {
    gp_Pnt point = P.Value(i);
    CubitVector v = CubitVector(point.X(), point.Y(), point.Z());
    cntrl_pts.append(v); 
  }

  if(rational)
  {
    TColStd_Array1OfReal W(1, h_S->NbPoles());
    h_S->Weights(W);
    for (int i = W.Lower(); i <= W.Upper(); i++)
      cntrl_pt_weights.append(W.Value(i));
  }
  return CUBIT_SUCCESS;
}

CubitStatus OCCCurve::get_spline_params	(	bool &	rational,
		int &	degree,
		DLIList< CubitVector > &	cntrl_pts,
		DLIList< double > &	cntrl_pt_weights,
		DLIList< double > &	knots,
		bool &	spline_is_reversed
	)		const `[virtual]`

Implements Curve.

Definition at line 1321 of file OCCCurve.cpp.

  {
      return CUBIT_FAILURE;
  }

void OCCCurve::get_tangent	(	CubitVector const &	location,
		CubitVector &	tangent
	)

Definition at line 891 of file OCCCurve.cpp.

{
    double u = u_from_position(location);
    Standard_Real first;
    Standard_Real last;
    Handle(Geom_Curve) myCurve = BRep_Tool::Curve(*myTopoDSEdge, first, last);
 
    gp_Pnt p;
    gp_Vec tan;
    GeomLProp_CurveTool::D1(myCurve, u , p, tan) ;
    tangent.set(tan.X(), tan.Y(), tan.Z());
}

TopoDS_Edge* OCCCurve::get_TopoDS_Edge ( ) [inline]

Definition at line 337 of file OCCCurve.hpp.

    { assert (myTopoDSEdge->ShapeType() == TopAbs_EDGE); return myTopoDSEdge; }

CubitBoolean OCCCurve::is_periodic ( double & period ) [virtual]

Implements Curve.

Definition at line 333 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  if (acurve.IsPeriodic())
  {
    period = acurve.Period();
    return CUBIT_TRUE;
  }
  return CUBIT_FALSE;
}

CubitBoolean OCCCurve::is_position_on ( const CubitVector & test_position ) [virtual]

Implements Curve.

Definition at line 637 of file OCCCurve.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;
}

double OCCCurve::length_from_u	(	double	parameter1,
		double	parameter2
	)		`[virtual]`

Implements Curve.

Definition at line 317 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  return GCPnts_AbscissaPoint::Length(acurve, parameter1, parameter2);
}

DLIList<OCCLoop*> OCCCurve::loops ( ) [inline]

Definition at line 53 of file OCCCurve.hpp.

{return myLoopList;}

double OCCCurve::measure ( ) [virtual]

Implements GeometryEntity.

Definition at line 298 of file OCCCurve.cpp.

{
  GProp_GProps myProps;
  BRepGProp::LinearProperties(*myTopoDSEdge, myProps);
  return myProps.Mass();
}

CubitPointContainment OCCCurve::point_containment ( const CubitVector & point ) [virtual]

Implements Curve.

Definition at line 1053 of file OCCCurve.cpp.

{
   if (is_position_on(point) == CUBIT_TRUE)
   {
     DLIList<OCCPoint*> points;
     get_points(points);
     for (int i = 0; i < points.size(); i++)
     {
    OCCPoint* pnt = points.get_and_step();
        CubitVector v = pnt->coordinates();
        double d = v.distance_between(point); 
        if (d < GEOMETRY_RESABS)
      return CUBIT_PNT_BOUNDARY; 
     }
     return CUBIT_PNT_ON;
   } 
   return CUBIT_PNT_OFF;
}

CubitStatus OCCCurve::position_from_u	(	double	u_value,
		CubitVector &	output_position
	)		`[virtual]`

Implements Curve.

Definition at line 547 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  gp_Pnt p = acurve.Value(u_value);
  output_position.x(p.X());
  output_position.y(p.Y());
  output_position.z(p.Z());
  return CUBIT_SUCCESS;
}

Curve * OCCCurve::project_curve	(	Surface *	face_ptr,
		DLIList< TBPoint * > &	normal_proj_points,
		CubitBoolean	closed,
		const CubitVector *	third_point
	)

Definition at line 1142 of file OCCCurve.cpp.

{
   TopoDS_Edge* edge = get_TopoDS_Edge();
   if (edge == NULL)
   {
        PRINT_ERROR("Cannot project the curve .\n"
                 "Possible incompatible geometry engines.\n");
        return (Curve*) NULL;
   }

   TopoDS_Face* face = CAST_TO(face_ptr, OCCSurface)->get_TopoDS_Face();
   if(face == NULL)
   {
        PRINT_ERROR("Cannot project the curve to the surface.\n"
                 "Possible incompatible geometry engines.\n");
        return (Curve*) NULL;
   }

   BRepAlgo_NormalProjection aProjection;
   aProjection.Init(*face);
   aProjection.Add(*edge);
   aProjection.Build();
   if (!aProjection.IsDone())
   {
        PRINT_ERROR("Cannot project the curve to the surface.\n"
                 "OCC engine failure.\n");
        return (Curve*) NULL;
   }

   TopoDS_Shape new_shape = aProjection.Projection();//compound shape
   int num_projection = 0;
   if (new_shape.IsNull())
   {
       PRINT_ERROR("Cannot project the curve to the surface.\n");
       return (Curve*) NULL;
   }

   else
   {
     //count how many free edges and vertices the new_shape has.
     TopExp_Explorer Ex;
     for (Ex.Init(new_shape,TopAbs_EDGE); Ex.More(); Ex.Next())
       num_projection++;
     for (Ex.Init(new_shape,TopAbs_VERTEX, TopAbs_EDGE); Ex.More(); Ex.Next())
       num_projection++;
   }

   if(num_projection == 0)
   {
      PRINT_INFO("No projection on the surface.\n");
      return (Curve*) NULL;
   }

   else if ( num_projection == 1 )
   {
      if(closed == true)
       PRINT_WARNING("Cannot project the curve to create a closed projection.\n"                 "There is only one projection segment.\n");

      TopExp_Explorer Ex;
      TopoDS_Edge new_edge;
      TopoDS_Vertex new_point;
      for (Ex.Init(new_shape,TopAbs_EDGE); Ex.More(); Ex.Next())
      {
        new_edge = TopoDS::Edge(Ex.Current());
        return OCCQueryEngine::instance()->populate_topology_bridge(new_edge,
                                                                   CUBIT_TRUE );
      }
      for(Ex.Init(new_shape,TopAbs_VERTEX);Ex.More(); Ex.Next())
      {
        new_point = TopoDS::Vertex(Ex.Current());
        normal_proj_points.append(OCCQueryEngine::instance()->populate_topology_bridge(new_point, CUBIT_TRUE));
      } 
      return (Curve*) NULL;
   }

   else if (num_projection == 2)
   {
      double d;
      double first, last;
      TopExp_Explorer Ex;
      TopoDS_Edge edge1, edge2;
      TopoDS_Vertex point;

      int count = 0;
      for (Ex.Init(new_shape,TopAbs_EDGE); Ex.More(); Ex.Next())
      {
        count++;
        if(count == 1)
          edge1 = TopoDS::Edge(Ex.Current());
        if(count == 2)
          edge2 = TopoDS::Edge(Ex.Current());
      }

      for(Ex.Init(new_shape,TopAbs_VERTEX);Ex.More(); Ex.Next())
      {
        point = TopoDS::Vertex(Ex.Current());
        normal_proj_points.append(OCCQueryEngine::instance()->populate_topology_bridge(point, CUBIT_TRUE));
      }

      if(edge1.IsNull())
        return OCCQueryEngine::instance()->populate_topology_bridge(edge2, CUBIT_TRUE);

      if(edge2.IsNull())
        return OCCQueryEngine::instance()->populate_topology_bridge(edge1, CUBIT_TRUE);

      if(edge1.IsNull() && edge2.IsNull())
        return (Curve*) NULL;

      Handle(Geom_Curve) myCurve1 =
                        BRep_Tool::Curve(edge1,first,last);

      Handle(Geom_Curve) myCurve2= BRep_Tool::Curve(edge2, first, last);
      //If the surface is periodic, so it has 2 projections, we just need to
      //find the segment to which the third_point is closer.
      if(closed == CUBIT_FALSE && third_point != NULL)
      {
        gp_Pnt P (third_point->x(), third_point->y(), third_point->z());
        GeomAPI_ProjectPointOnCurve projOncurve(P, myCurve1);
        if (projOncurve.NbPoints() == 0)
        {
          PRINT_ERROR("Cannot project the curve to the surface.\n"
                 "OCC engine failure.\n");
          return (Curve*) NULL;
        }
        d = projOncurve.LowerDistance();

        //Compare with the second solution
        GeomAPI_ProjectPointOnCurve projOncurve2(P, myCurve2);
        if (projOncurve2.NbPoints() == 0)
        {
           PRINT_ERROR("Cannot project the curve to the surface.\n"
                 "OCC engine failure.\n");
          return (Curve*) NULL;
        }

        double d2 = projOncurve2.LowerDistance();
        TopoDS_Edge new_edge =
                d > d2 ? edge2 : edge1 ;
        return OCCQueryEngine::instance()->populate_topology_bridge(new_edge, CUBIT_TRUE);
      }


      else if (closed == CUBIT_TRUE)
      {
        //connect the two segment into a closed shape. Assume both segment
        // has the same curve type, create Bezier closed curve.
        GeomAdaptor_Curve acurve1(myCurve1);
        GeomAdaptor_Curve acurve2(myCurve2);
        //get 10 points of each curve, combine them to make one Bezier curve
        int NbPoints = 10;
        GCPnts_QuasiUniformAbscissa distribution1(acurve1, NbPoints);
        GCPnts_QuasiUniformAbscissa distribution2(acurve2, NbPoints);
        TColgp_Array1OfPnt points(1, 2*NbPoints-1);
        int i;
        for (i = 1; i <= NbPoints; i++)
        {
           double u = distribution1.Parameter(i);
           gp_Pnt P = myCurve1->Value(u);
           points.SetValue(i, P);
        }

        for (int j = NbPoints-1; j >= 1; j--)
        {
           double u = distribution2.Parameter(j); 
           gp_Pnt P = myCurve2->Value(u);
           points.SetValue(++i,P); 
        }    

        Geom_BezierCurve BezierCurve(points);
        Handle(Geom_Curve) curve_ptr(&BezierCurve);
        TopoDS_Edge new_edge = BRepBuilderAPI_MakeEdge(curve_ptr);
        return OCCQueryEngine::instance()->populate_topology_bridge(new_edge);
      }
   }
   return (Curve*) NULL;
}

void OCCCurve::remove_all_simple_attribute_virt ( ) [virtual]

Implements TopologyBridge.

Definition at line 231 of file OCCCurve.cpp.

{
  OCCAttribSet::remove_attribute(CubitSimpleAttrib(), *myTopoDSEdge);
}

void OCCCurve::remove_loop ( OCCLoop * loop ) [inline]

Definition at line 54 of file OCCCurve.hpp.

{myLoopList.remove(loop);}

void OCCCurve::remove_simple_attribute_virt ( const CubitSimpleAttrib & csa ) [virtual]

Implements TopologyBridge.

Definition at line 216 of file OCCCurve.cpp.

  { OCCAttribSet::remove_attribute(csa, *myTopoDSEdge); }

void OCCCurve::set_myMarked ( CubitBoolean marked ) [inline]

Definition at line 50 of file OCCCurve.hpp.

{myMarked = marked;}

void OCCCurve::set_TopoDS_Edge ( TopoDS_Edge edge )

Definition at line 156 of file OCCCurve.cpp.

{
  if(myTopoDSEdge && edge.IsEqual(*myTopoDSEdge))
    return;

  else
  {
    DLIList<OCCPoint*> points ;
    this->get_points(points);
    for(int i = 0; i < points.size(); i++)
    {
      OCCPoint* point = points.get_and_step();
      TopoDS_Vertex* vtx = point->get_TopoDS_Vertex();
      TopExp_Explorer Ex;
      CubitBoolean found = false;
      for (Ex.Init(edge, TopAbs_VERTEX); Ex.More(); Ex.Next())
      {
        TopoDS_Shape sh = Ex.Current();
        if(vtx->IsPartner(sh))
        {
           found = true;
           break;
         }
      }
      if (!found)
        point->remove_curve(this);
    }
  }
 
  if(myTopoDSEdge)
    myTopoDSEdge->Nullify();
  *myTopoDSEdge = edge;
}

double OCCCurve::start_param ( ) [virtual]

Implements Curve.

Definition at line 765 of file OCCCurve.cpp.

{
   double start = 0.0, end = 0.0;
   get_param_range( start, end );
   return start;
}

double OCCCurve::u_from_arc_length	(	double	root_param,
		double	arc_length
	)		`[virtual]`

Implements Curve.

Definition at line 618 of file OCCCurve.cpp.

{
  BRepAdaptor_Curve acurve(*myTopoDSEdge);
  GCPnts_AbscissaPoint abs(acurve, arc_length, root_param);
  if (abs.IsDone()) return abs.Parameter();
  else return 0.0;
}

double OCCCurve::u_from_position ( const CubitVector & input_position ) [virtual]

Implements Curve.

Definition at line 571 of file OCCCurve.cpp.

{
    // Get the closest point on the Curve to the input position
  CubitVector closest_point;
  double u_val;
  this->closest_point(input_position, closest_point,
                      NULL, NULL, &u_val);
    // closest_point already makes adjustments for sense and periodicity
  
  return u_val;
}

void OCCCurve::update_OCC_entity	(	BRepBuilderAPI_ModifyShape *	aBRepTrsf,
		BRepAlgoAPI_BooleanOperation *	op = `NULL`
	)

Definition at line 1077 of file OCCCurve.cpp.

{
  if (myMarked == 1) 
     return;

  assert(aBRepTrsf != NULL || op != NULL);
  
  TopoDS_Shape shape;
  if(aBRepTrsf)
    shape = aBRepTrsf->ModifiedShape(*get_TopoDS_Edge()); 
  else
  {
    TopTools_ListOfShape shapes;
    shapes.Assign(op->Modified(*get_TopoDS_Edge()));
    if(shapes.Extent() == 0)
      shapes.Assign(op->Generated(*get_TopoDS_Edge()));
    if(shapes.Extent() == 1)
      shape = shapes.First();
    else if(shapes.Extent() > 1)
    {
      shape = shapes.First();
    }
    else if (op->IsDeleted(*get_TopoDS_Edge()))
      ;
    else
      return ;
  }
  TopoDS_Edge curve;
  if(!shape.IsNull())
    curve = TopoDS::Edge(shape);

  //make sure the shape (edge) length is greater than 0.
  GProp_GProps myProps;
  BRepGProp::LinearProperties(curve, myProps);
  double d = myProps.Mass();
  if(d > OCCQueryEngine::instance()->get_sme_resabs_tolerance())
  {
    //set the vertices
    DLIList<TopologyBridge*> vertices;
    get_children_virt(vertices);
    for (int i = 1; i <= vertices.size(); i++)
    {
       TopologyBridge* tb = vertices.get_and_step();
       OCCPoint *point = CAST_TO(tb, OCCPoint);
       if (point)
         point->update_OCC_entity(aBRepTrsf, op);
    }
    myMarked = 1;
  }
  else
    curve.Nullify();
  OCCQueryEngine::instance()->update_OCC_map(*myTopoDSEdge, curve);
}

Member Data Documentation

DLIList<OCCLoop*> OCCCurve::myLoopList [private]

Definition at line 364 of file OCCCurve.hpp.

CubitBoolean OCCCurve::myMarked [private]

Definition at line 365 of file OCCCurve.hpp.

TopoDS_Edge* OCCCurve::myTopoDSEdge [private]

Definition at line 363 of file OCCCurve.hpp.

The documentation for this class was generated from the following files:

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation