<<>> <> <> <> 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. .. <> <> <> <> <> <> <> <> <> <<};>> <> <> <> <> <<};>> <<>> <> <> <> <> <> <> <> <> <> <> < solid ¬ ImagerPath.MoveTo[v];>> < dash ¬ ImagerPath.MoveTo[v];>> < dash ¬ ImagerPath.LineTo[dash, v];>> < solid ¬ ImagerPath.LineTo[solid, v];>> <> <> <> <> <> <<};>> <> <> <> <> <> <> <> <> <> <<0.0 => ii.distance,>> <<1.0 => i.distance,>> < G3dVector.Distance[near.point, point];>> <> <> <<}>> <> <> <> < 2 THEN value ¬ 2.0*value/(info.length*(info.length-1));>> <<};>> <<>> <> <> <> < strength THEN {strength ¬ info[n].strength; max ¬ n};>> <> <<};>> <<>> <> <<~ {>> <> <> <> <> <> <> <> <<};>> <<>> <> <> <> <> < 1.0 THEN tot*tot*ii.strength ELSE ii.strength)];>> <<};>> <> <<};>> <<>> <> <> <> <> <> <> <> <<};>> <> <<};>> <> <<};>> <<>> <> <> <<};>> <<>> <> <> <> < 4.0>> <> <> <> <> <<};>> <<};>> <<>> <> <> <> <> < 4.0-4.00*dOverR,>> < 4.0-4.00*dOverR,>> < 3.0-2.00*dOverR,>> < 3.0-2.00*dOverR,>> < 1.5-0.50*dOverR,>> < 1.5-0.50*dOverR,>> < 1.0-0.25*dOverR,>> < 1.0-0.25*dOverR,>> < 0.0];>> <<};>> <> <> <> <> <> <> <> <> < value THEN {value ¬ strength; max ¬ n};>> <> <> <> <> <> <> <<};>> <> <> <<};>> <<>> 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]; <> <> <<< -p.line.value => SqDist[q, [-p.line.value, 0.0]],>> <<> p.line.value => SqDist[q, [p.line.value, 0.0]],>> < q.y*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;