1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
|
// File: ChFiKPart_ComputeData_CS.cxx
// Created: Wed Apr 13 16:57:48 1994
// Author: Isabelle GRIGNON
// <isg@phylox>
#include <ChFiKPart_ComputeData_CS.hxx>
#include <gp_Ax3.hxx>
#include <gp_Vec.hxx>
#include <gp_Pnt.hxx>
#include <gp_Dir.hxx>
#include <ElCLib.hxx>
void ChFiKPart_CornerSpine(const Handle(Adaptor3d_HSurface)& S1,
const Handle(Adaptor3d_HSurface)& S2,
const gp_Pnt2d& P1S1,
const gp_Pnt2d& P2S1,
const gp_Pnt2d& P1S2,
const gp_Pnt2d& P2S2,
const Standard_Real R,
gp_Cylinder& cyl,
gp_Circ& circ,
Standard_Real& First,
Standard_Real& Last)
{
gp_Ax3 ax = S1->Plane().Position();
gp_Vec V1(ax.XDirection());
gp_Vec V2(ax.YDirection());
gp_Pnt P;
gp_Vec du,dv;
S2->D1(P1S2.X(),P1S2.Y(),P,du,dv);
gp_Vec V(P,S1->Value(P1S1.X(),P1S1.Y()));
V = V.Dot(V1)*V1+V.Dot(V2)*V2;
V.Normalize();
gp_Pnt P2 = S2->Value(P2S2.X(),P2S2.Y());
gp_Vec Vorien(P,P2);
gp_Pnt cent;
gp_Dir dx(V);
if(V.Dot(Vorien) >= 0.){
cent.SetCoord(P.X()+R*V.X(),P.Y()+R*V.Y(),P.Z()+R*V.Z());
dx.Reverse();
}
else {
cent.SetCoord(P.X()-R*V.X(),P.Y()-R*V.Y(),P.Z()-R*V.Z());
}
gp_Dir dy(gp_Vec(cent,P2));
dy = (dx^dy)^dx;
gp_Ax2 circax2(cent,dx^dy,dx);
gp_Ax3 cylax3(circax2);
if((du^dv).Dot(dx) < 0.) cylax3.ZReverse();
First = 0.;
Last = ElCLib::CircleParameter(circax2,P2);
circ.SetPosition(circax2);
circ.SetRadius(R);
cyl.SetPosition(cylax3);
cyl.SetRadius(R);
}
|
