
DIRECTORY Commander, Controls, Draw2d, FileNames, G3dDraw, G3dMatrix, G3dTool, G3dTube, ImplicitDefs, ImplicitSlice, ImplicitTube, IO, Menus, Rope;

ImplicitSliceCmdsImpl: CEDAR PROGRAM

IMPORTS Controls, FileNames, ImplicitSlice, ImplicitTube, IO, G3dDraw, G3dTool, G3dTube, Rope

~ BEGIN
Button:				TYPE ~ Controls.Button;
DrawProc:			TYPE ~ Draw2d.DrawProc;
ClickProc:			TYPE ~ Menus.ClickProc;
DiagramProc:			TYPE ~ ImplicitDefs.DiagramProc;
ValueProc:			TYPE ~ ImplicitDefs.ValueProc;
Tool:					TYPE ~ ImplicitSlice.Tool;
Function:				TYPE ~ ImplicitSlice.Function;
Tube:					TYPE ~ G3dTube.Tube;
ROPE:					TYPE ~ Rope.ROPE;
TubeData:			TYPE ~ REF TubeDataRep;
TubeDataRep:		TYPE ~ RECORD [
tool:							Tool м NIL,
tube:							Tube м NIL,
drawLines:					BOOL м FALSE
];

TubeSlice: Commander.CommandProc ~ {
t: TubeData м NEW[TubeDataRep];
f: LIST OF Function м LIST[["segment", ValueSegment, DrawSegment]];
b: LIST OF Button м LIST[
Controls.ClickButton["GetTube", GetTube, t],
Controls.ClickButton["Lines: Off", ToggleDrawLines, t]];
t.tool м ImplicitSlice.MakeTool["ImplicitSlice Tube", f, DrawTube, t,, b];
};

DrawTube: DrawProc ~ {
t: TubeData м NARROW[data];
v: G3dMTool.View м G3dTool.GetView[viewer, t.tool.renderTool];
G3dTube.DrawSplines[context, t.tube, xform];
IF t.drawLines THEN {
G3dTube.DrawContours[context, t.tube, v.camera.matrix, v.viewport];
G3dTube.DrawLines[context, t.tube, v.camera.matrix, v.viewport];
};
};

DrawSegment: DiagramProc ~ {
t: TubeData м NARROW[clientData];
xform: Matrix м G3dTool.GetView[viewer, t.tool.renderTool].camera.matrix;
vp: G3dMatrix.Viewport м G3dView.GetViewport[viewer];
ImplicitTube.DiagramTube[context, xform, vp, point, t.tube, $Segment, inverseSqrd, 0.0, 0.75];
G3dDraw.Mark[context, t.tool.point, xform, vp, IO.PutFR["  %6.3f", IO.real[t.tool.value]]];
};

ValueSegment: ValueProc ~ {
t: TubeData м NARROW[clientData];
value м ImplicitTube.SampleTube[point, t.tube, $Segment, inverseSqrd, 0.0, 0.75].value;
};

ToggleDrawLines: ClickProc ~ {
t: TubeData м NARROW[clientData];
t.drawLines м NOT t.drawLines;
Controls.ButtonToggle[t.tool.renderTool.outerData, t.drawLines, "Lines: On", "Lines: Off"];
ImplicitSlice.Repaint[t.tool];
};

GetTube: ClickProc ~ {
t: TubeData м NARROW[clientData];
t.tube м G3dTube.LastTubePicked[TRUE];
IF t.tube = NIL THEN {
r: ROPE мFileNames.ResolveRelativePath[Controls.TypescriptReadFileName[t.tool.typescript]];
IF (t.tube м G3dTube.TubeFromFile[r]) = NIL THEN {
Controls.TypescriptWrite[t.tool.typescript, Rope.Cat["Bad file: ", r, "\n"]];
RETURN;
};
};
G3dTube.SelectAll[t.tube];
IF t.tube.xSections = NIL THEN G3dTube.MakeXSections[t.tube, TRUE, TRUE];
ImplicitSlice.Repaint[t.tool];
};
ImplicitSlice.Register["Tube", TubeSlice, "Explore various tube designs."];

END.
..

SqDistToSegment: PROC [a, b, q: Pair] RETURNS [REAL] ~ {
n: Pair ~ [b.x-a.x, b.y-a.y];
nSq: REAL ~ n.x*n.x+n.y*n.y;
qSq: REAL ~ SqDist[q, a];
IF nSq < 0.1*qSq
THEN RETURN[qSq]
ELSE {
ua: Pair ~ [q.x-a.x, q.y-a.y];
alpha: REAL ~ (ua.x*n.x+ua.y*n.y)/nSq;
RETURN[SqDist[q, SELECT alpha FROM
<= 0.0 => a,
>= 1.0 => b,
ENDCASE => [a.x+alpha*n.x, a.y+alpha*n.y]]];
};
};

SqDist: PROC [p0, p1: Pair] RETURNS [REAL] ~ {
d: Pair м [p1.x-p0.x, p1.y-p0.y];
RETURN[d.x*d.x+d.y*d.y];
}; 

SqDistToContour: PROC [contour: Contour, q: Pair] RETURNS [min: REAL м 10000.0] ~ {
pairs: PairSequence ~ contour.pairs;
n, n0, n1: NAT;
FOR i: NAT IN [0..pairs.length) DO
d: REAL ~ SqDist[q, pairs[i]];
IF d < min THEN {min м d; n м i};
ENDLOOP;
n0 м IF n > 0 THEN n-1 ELSE IF contour.closed THEN pairs.length-1 ELSE n;
n1 м IF n < pairs.length-1 THEN n+1 ELSE IF contour.closed THEN 0 ELSE n;
RETURN[MIN[
SqDistToSegment[pairs[n0], pairs[n], q],
SqDistToSegment[pairs[n], pairs[n1], q]]];
};

SqDistToArc: PROC [q: Pair, arc: REAL] RETURNS [REAL] ~ {
IF q.y > arc THEN RETURN[
MIN[SqDist[q, [arc, arc]], SqDist[q, [-arc, arc]]]]
ELSE {
yDif: REAL ~ q.y-arc;
d: REAL ~ arc-RealFns.SqRt[q.x*q.x+yDif*yDif];
RETURN[d*d];
};
};

FOR s: NAT IN [pm.sMin..pm.sMin+pm.sSize) DO
IF p.abort THEN RETURN;
FOR f: NAT IN [pm.fMin..pm.fMin+pm.fSize) DO
q: Pair ~ PairCoords[[f, s]];
potential: REAL м 0.0;
IF p.options[radius] AND p.rad.value # 0.0 THEN {
d: REAL ~ q.x*q.x+q.y*q.y;
potential м IF d # 0.0 THEN 10.0*p.rad.value/d ELSE 255.0;
};
IF p.options[arc] THEN {
d: REAL ~ SqDistToArc[q, p.arc.value];
potential м IF d # 0.0 THEN potential+1.0/d ELSE 255.0;
};
IF p.options[line] AND ABS[q.x] < p.line.value THEN {
potential м IF q.y = 0.0
THEN 255.0
ELSE potential+10.0*(p.line.value-ABS[q.x])*(p.line.value-ABS[q.x])/ABS[q.y];
};
IF p.options[contour] THEN {
d: REAL ~ SqDistToContour[c, q];
potential м IF d # 0.0 THEN potential+1.0/d ELSE 255.0;
};
CtBasic.PutPixel[pm, f, s, Real.RoundI[MIN[255.0, MAX[0.0, potential]]]];
ENDLOOP;
ENDLOOP;
  °  
ImplicitSliceCmdsImpl.mesa
Copyright ╙ 1985, 1990 by Xerox Corporation.  All rights reserved.
Bloomenthal, August 11, 1992 4:06 pm PDT
Types
Tube Testing
Start Code
Fancy Tube Display
DrawTube: PROC [context: Context, tube: Tube] ~ {
tubeProc: G3dTube.TubeProc ~ {
i: G3dSpline.SplinePlaneIntersect м
G3dSpline.IntersectSplineAndPlane[tube.coeffs, t.slice.plane];
FOR n: NAT IN [0..i.nRoots) DO
Draw3d.Mark[context, G3dSpline.Position[s.tube.coeffs, i.roots[n]], t.view, , x];
ENDLOOP;
DrawSpline[context, tube.coeffs];
};
Imager.SetColor[context, Imager.black];
Imager.SetStrokeWidth[context, 2.0];
G3dTube.ApplyToTube[tube, tubeProc];
Imager.SetStrokeWidth[context, 4.0];
};

DrawSpline: PROC [context: Context, coeffs: Coeffs] ~ {
PatternProc: PROC [i: NAT] RETURNS [REAL] ~ {RETURN[2]};
behind: BOOL м Plane3d.DistanceToPt[Spline3d.Position[coeffs, 0.0], s.slice.plane] < 0;
was: BOOL м NOT behind;
solid, dash: Imager.Trajectory м NIL;
FOR t: REAL м 0.0, t+0.005 WHILE t <= 1.004 DO
q: Triple м G3dSpline.Position[coeffs, t];
v: Pair м G3dMatrix.TransformD[q, s.view];
behind м G3dPlane.DistanceToPoint[q, s.slice.plane] < 0;
SELECT TRUE FROM
was AND NOT behind => solid м ImagerPath.MoveTo[v];
NOT was AND behind => dash м ImagerPath.MoveTo[v];
behind => dash м ImagerPath.LineTo[dash, v];
ENDCASE => solid м ImagerPath.LineTo[solid, v];
was м behind;
ENDLOOP;
IF dash # NIL THEN
Imager.MaskDashedStrokeTrajectory[context, dash, 2, PatternProc, 0, 20];
IF solid # NIL THEN Imager.MaskStrokeTrajectory[context, solid];
};
Discarded Value Procs
Line: InfoProc ~ {
FOR n: NAT IN [0..info.length) DO
i: Info м info[n];
FOR nn: NAT IN (n..info.length) DO
ii: Info м info[nn];
IF i.point # ii.point THEN {
near: NearSegment3d м G3dVector.NearestTo3dSegment[i.point, ii.point, point];
distance: REAL м SELECT near.w0 FROM
0.0			=> ii.distance,
1.0			=> i.distance,
ENDCASE	=> G3dVector.Distance[near.point, point];
aveRadius: REAL м 0.5*(i.radius+ii.radius);
value м value+Strength[Ratio[distance, aveRadius, spread]];
}
ELSE value м value+i.strength;
ENDLOOP;
ENDLOOP;
IF info.length > 2 THEN value м 2.0*value/(info.length*(info.length-1));
};

MaxInfo: PROC [info: REF Infos] RETURNS [max: NAT м 0] ~ {
strength: REAL м 0.0;
FOR n: NAT IN [0..info.length) DO
IF info[n].strength > strength THEN {strength м info[n].strength; max м n};
ENDLOOP;
};

TotalValue: PROC [info: REF Infos, proc: PROC [i, ii: Info] RETURNS [REAL]] RETURNS [REAL]
~ {
max: NAT м MaxInfo[info];
i: Info м info[max];
value: REAL м i.strength;
FOR n: NAT IN [0..info.length) DO
IF n # max THEN value м value+proc[i, info[n]];
ENDLOOP;
RETURN[value];
};

CosineAndSynergy: InfoProc ~ {
Proc: PROC [i, ii: Info] RETURNS [REAL] ~ {
tot: REAL м i.strength+ii.strength;
acos: REAL м G3dVector.Dot[i.vector, ii.vector];
RETURN[(1.0-acos)* (IF tot > 1.0 THEN tot*tot*ii.strength ELSE ii.strength)];
};
RETURN[TotalValue[info, Proc]];
};

CosineAndClose: InfoProc ~ {
Proc: PROC [i, ii: Info] RETURNS [REAL] ~ {
totalStrength: REAL м i.strength+ii.strength;
IF totalStrength # 0.0 THEN {
acos: REAL м G3dVector.Dot[i.vector, ii.vector];
closeness: REAL м 1.0-ABS[i.strength-ii.strength]/totalStrength;
RETURN[closeness*(1.0-acos)*ii.strength];
};
RETURN[0.0];
};
RETURN[TotalValue[info, Proc]];
};

Summation: InfoProc ~ {
FOR n: NAT IN [0..info.length) DO value м value+info[n].strength; ENDLOOP;
};

Strength: PROC [distance, radius, spread: REAL] RETURNS [REAL] ~ {
Attempt 2:
dOverR: REAL м distance/radius;
IF dOverR > 4.0
THEN RETURN[0.0]
ELSE {
dif: REAL м dOverR-4.0;
RETURN[4.0-Real.SqRt[16.0-dif*dif]];
};
};

Strength: PROC [distance, radius, spread: REAL] RETURNS [REAL] ~ {
Attempt 1:
dOverR: REAL м distance/radius;
RETURN[SELECT dOverR FROM
IN[0.00..0.25)	=> 4.0-4.00*dOverR,
IN[0.25..0.50)	=> 4.0-4.00*dOverR,
IN[0.50..0.75)	=> 3.0-2.00*dOverR,
IN[0.75..1.00)	=> 3.0-2.00*dOverR,
IN[1.00..1.50)	=> 1.5-0.50*dOverR,
IN[1.50..2.00)	=> 1.5-0.50*dOverR,
IN[2.00..3.00)	=> 1.0-0.25*dOverR,
IN[3.00..4.00)	=> 1.0-0.25*dOverR,
ENDCASE		=> 0.0];
};
ValueOfTube: PROC [tube: Tube, point: Triple, spread: REAL] RETURNS [value: REAL м 0.]~ {
max, nBranches: NAT м G3dTube.NBranches[tube]+1;
FOR n: NAT IN [0..nBranches) DO
branch: Tube м IF n = 0 THEN tube.next ELSE tube.branches[n-1];
nearest: NearSpline3d м G3dSpline.PreciseNearestPoint[point, branch.xCoeffs];
vector: Triple м G3dVector.Normalize[G3dVector.Sub[nearest.point, point]];
radius: REAL м G3dTube.Radius[branch, nearest.t];
strength: REAL м Strength[nearest.distance, radius, spread];
IF strength > value THEN {value м strength; max м n};
info[n] м [nearest.point, vector, nearest.distance, radius, strength];
ENDLOOP;
FOR n: NAT IN [0..nBranches) DO
IF n # max THEN {
acos: REAL м G3dVector.Dot[info[max].vec, info[n].vec];
value м value+0.5*(1.0-acos)*info[n].strength;
};
ENDLOOP;
value м (4.0/3.0)*value;
};

ELSE potential+(p.line.value-ABS[q.x])/(q.y*q.y);
d: REAL ~ SELECT q.x FROM
< -p.line.value => SqDist[q, [-p.line.value, 0.0]],
> p.line.value => SqDist[q, [p.line.value, 0.0]],
ENDCASE => q.y*q.y;
potential м IF d # 0.0 THEN potential+1.0/d ELSE 255.0;
 Э╩ж  ХNewlineDelimiter
Щ ЩJЪЬBЩBJЩ(JШ JЪ╧k	ЬКШУJШ ЧЪ╨blЬЭЬЭШ$JШ JЪЭЬVШ]JШ ЧJЪЬЭШheadЪ╧lЩJЪЬ
ЭЬШ"JЪЬЭЬШ#JЪЬЭЬШ$JЪЬЭЬШ/JЪЬЭЬШ+JЪЬ
ЭЬШ$JЪЬЭЬШ+JЪЬЭЬШJЪЭЬЭЬЭЬШЧЪЯЩJЪЬЭЬЭЬШ#ЪЬЭЬЭЬШJЪЬЭЬШJЪЬЭЬШJЪЬЭЬЭШJШJШ ЧЪ╨bn	ЬШ$JЪЬЭЬШJЪЬЭЬЭЬЭЬ)ШCЪЬЭЬЭЬ
ЭЬШJШ,JШ8ЧJШJJШJШ ЧЪ╧nЬШJЪЬЭЬШJШ>JШ,ЪЭЬЭЬШJШCJШ@JШЧJШJШ ЧЪбЬШJЪЬЭЬШ!JШIJШ5JШ^JЪЬ/ЭЬЭЬШ[JШJШ ЧЪбЬШJЪЬЭЬШ!JШWJШJШ ЧЪ╧bЬШJЪЬЭЬШ!JЪЬЭЬШJШ[JШJШJШ ЧЪвЬШJЪЬЭЬШ!JЪЬ ЭЬШ&ЪЭЬ
ЭЬЭЬШJЪЬ╧sЭгЬSШ[ЪЭЬ&ЭЬЭЬШ2JЪЬMШMJЪЭЬШJШЧJШЧJШJЪ
ЭЬЭЬЭЬЭЬЭЬШIJШJШЧЧЪЯ
Щ
JШKЧJШ ЪЭЬШJШЧЪЯЩЪбЬЭЬ#Щ1ЪбЬЩcodeЩ#LЩ>ЧЪЭЬЭЬЭЬЭЩLЩQLЪЭЬЩЧJЩ!JЩЧJЩ'JЩ$JЩ$JЩ$JЩJЩ ЧЪб
ЬЭЬ'Щ7JЪбЬЭЬЭЬЭЬЭЬЭЬЩ8JЪЬЭЬHЩWJЪЬЭЬЭЬЩJЪЬ!ЭЬЩ%ЪЭЬЭЬЭЬЭЩ.LЩ*LЩ*LЩ8ЪЭЬЭЬЭЩJЪЬЭЬЭЬ(Щ3JЪЭЬЭЬ'Щ2JЩ,JЪЭЬ(Щ/ЧJЩLЪЭЬЩЧЪЭЬЭЬЭЩLЩHЧLЪЭЬ	ЭЬЭЬ-Щ@JЩЧЧЪЯЩЪбЬЩЪЭЬЭЬЭЬЭЩ!JЩЪЭЬЭЬЭЬЭЩ"JЩЪЭЬЭЬЩJЩMЪЬ
ЭЬЭЬ	ЭЩ$JЩJЩJЪЭЬ*Щ1ЧJЪЬЭЬЩ+JЩ;JЩЧJЪЭЬЩJЪЭЬЩЧJЪЭЬЩЧJЪЭЬЭЬ1ЩHJЩJЩ ЧЪ
бЬЭЬЭЬЭЬЭЬ	Щ:JЪЬ
ЭЬЩЪЭЬЭЬЭЬЭЩ!JЪЭЬЭЬ(ЩKJЪЭЬЩЧJЩJЩ ЧЪб
ЬЭЬЭЬЭЬЭЬЭЬЭЬЭЬЩZJЩJЪЬЭЬЩJЩJЪЬЭЬЩЪЭЬЭЬЭЬЭЩ!JЪЭЬ	ЭЬ Щ/JЪЭЬЩЧJЪЭЬЩJЩJЩ ЧЪбЬЩЪбЬЭЬЭЬЭЬЩ+JЪЬЭЬЩ#JЪЬЭЬ&Щ0JЪЭЬЭЬЭЬЭЬЩMJЩЧJЪЭЬЩJЩJЩ ЧЪбЬЩЪбЬЭЬЭЬЭЬЩ+JЪЬЭЬЩ-ЪЭЬЭЬЩJЪЬЭЬ&Щ0JЪЬЭЬЭЬ'Щ@JЪЭЬ#Щ)JЩЧJЪЭЬЩJЩЧJЪЭЬЩJЩJЩ ЧЪб	ЬЩJЪ
ЭЬЭЬЭЬЭЬ ЭЬЩJJЩJЩ ЧЪ
бЬЭЬЭЬЭЬЭЬЩBJЩ
JЪЬгЬЩЪЭЬЩJЪЭЬЭЬЩЪЭЬЩJЪЬЭЬЩJЪЭЬЩ$JЩЧЧЩJЩ ЧЧЪ
бЬЭЬЭЬЭЬЭЬЩBJЩ
JЪЬгЬЩЪЭЬЭЬЭЩJЪЭЬ Щ"JЪЭЬ Щ"JЪЭЬ Щ"JЪЭЬ Щ"JЪЭЬ Щ"JЪЭЬ Щ"JЪЭЬ Щ"JЪЭЬ Щ"JЪЭЬ
ЩЧJЩJШ ЧЪ
бЬЭЬ%ЭЬЭЬ	ЭЬ	ЩYJЪЬЭЬЩ0ЪЭЬЭЬЭЬЭЩJЪЬЭЬЭЬЭЬЩ?JЩMJЩJJЪЬЭЬ%Щ1JЪЬ
ЭЬ.Щ<JЪЭЬЭЬЩ5JЩFJЪЭЬЩЧЪЭЬЭЬЭЬЭЩЪЭЬ	ЭЬЩJЪЬЭЬ-Щ7JЩ.JЩЧJЪЭЬЩЧJЩJЩJЩ ЧЪбЬЭЬЭЬЭЬШ8JШJЪЬЭЬШJЪЬЭЬШЪЭЬШJЪЭЬЭЬШЪЭЬШJШJЪЬЭЬШ&ЪЭЬЭЬЭШ"JШJШJЪЭЬ%Ш,ЧJШЧЧJШJШ ЧЪбЬЭЬЭЬЭЬШ.JШ!JЪЭЬШJШJШ ЧЪбЬЭЬЭЬЭЬШSJШ$JЪЬЭЬШЪЭЬЭЬЭЬЭШ"JЪЬЭЬШJЪЭЬ	ЭЬШ!JЪЭЬШЧJЪЬЭЬЭЬЭЬЭЬЭЬЭЬШIJЪЬЭЬЭЬЭЬЭЬЭЬЭЬШIЪЭЬЭЬШJШ(JШ*ЧШLШ ЧЧЪ
бЬЭЬЭЬЭЬЭЬШ9ЪЭЬ
ЭЬЭЬШJЪЭЬ0Ш3ЧЪЭЬШJЪЬЭЬШJЪЬЭЬ'Ш.JЪЭЬШJШЧJШJШ ЧЪЭЬЭЬЭЬЭШ,JЪЭЬ	ЭЬЭЬШЪЭЬЭЬЭЬЭШ,JШJЪЬЭЬШЪЭЬЭЬЭЬШ1JЪЬЭЬШJЪЬЭЬ	ЭЬЭЬШ:JШЧЪЭЬЭЬШJЪЬЭЬШ&JЪЬЭЬ	ЭЬЭЬШ7JШЧЪЭЬЭЬЭЬЭЬШ5ЪЬЭЬ
ШJЪЭЬШ
JЪЭЬЭЬЭЬЭЬШMJЪЭЬЭЬЩ1ЧЪЬЭЬЭЬЭЩJЩ3JЩ1JЪЭЬЩЧJЪЬЭЬ	ЭЬЭЬЩ7JШЧЪЭЬЭЬШJЪЬЭЬШ JЪЬЭЬ	ЭЬЭЬШ7JШЧJЪЬ'ЭЬЭЬШIJЪЭЬШЧJЪЭЬШЧЧЧ ЕЧ    ╕  2V  