
DIRECTORY Commander, Controls, Draw2d, G3dBasic, G3dControl, G3dDraw, G3dFunction, G3dMatrix, G3dSpline, G3dTool, Imager, IO, Real, RealFns, Rope, ViewerOps;

G3dGraphCmdImpl: CEDAR PROGRAM
IMPORTS Controls, G3dControl, G3dDraw, G3dFunction, Draw2d, G3dSpline, G3dTool, Imager, IO, Real, RealFns, Rope, ViewerOps
~ BEGIN
ROPE:						TYPE ~ Rope.ROPE;
Viewer:					TYPE ~ Controls.Viewer;
ClickProc:				TYPE ~ Controls.ClickProc;
OuterData:				TYPE ~ Controls.OuterData;
Control:					TYPE ~ Controls.Control;
Request:					TYPE ~ Controls.Request;
Triple:					TYPE ~ G3dBasic.Triple;
Camera:					TYPE ~ G3dControl.CameraControl;
Matrix:					TYPE ~ G3dMatrix.Matrix;
arraySize:			NAT ~ 40;

GraphData:		TYPE ~ REF GraphDataRep;
GraphDataRep:	TYPE ~ RECORD [
a:							Control м NIL,
xMin:						Control м NIL,
xMax:						Control м NIL,
yMin:						Control м NIL,
yMax:						Control м NIL,
zScale:					Control м NIL,
grain:						Control м NIL,
values:					ARRAY [0..arraySize) OF ARRAY [0..arraySize) OF REAL,
graphics:					Viewer м NIL,
outer:						Viewer м NIL,
outerData:				OuterData м NIL,
camera:					Camera м NIL,
view:						Matrix м NIL,
axesOn:					BOOL м TRUE,
graphProc:				GraphProc м Gauss
];

GraphProc:		TYPE ~ PROC [x: REAL, g: GraphData] RETURNS [v: REAL];
RealFunction:	TYPE ~ PROC [v: REAL] RETURNS [REAL];
Function:			TYPE ~ RECORD [name: ROPE, use: ROPE, proc: GraphProc];
Graph: Commander.CommandProc ~ {
g: GraphData м NEW[GraphDataRep];
g.camera м G3dControl.InitCameraControl[proc: CameraControl, clientData: g];
Controls.SetSliderDialValue[g.camera.par.zMov, -0.5];
g.a м Controls.NewControl["A", , g, , , , Controller];
g.xMin м Controls.NewControl["Xmin", , g, -1.0, 1.0, 0.0, Controller];
g.yMin м Controls.NewControl["Ymin", , g, -1.0, 1.0, 0.0, Controller];
g.xMax м Controls.NewControl["Xmax", , g, -1.0, 1.0, 1.0, Controller];
g.yMax м Controls.NewControl["Ymax", , g, -1.0, 1.0, 1.0, Controller];
g.zScale м Controls.NewControl["Zscale", , g, -5.0, 5.0, 1.0, Controller];
g.grain м Controls.NewControl["Grain", , g, 0.05, 1.0, 0.1, Controller];
ComputeValues[g];
g.outerData м Controls.OuterViewer[
name: "3dGraph Gauss",
buttons: LIST[
Controls.ClickButton["Interpress-Out", InterpressOut, g],
Controls.ClickButton["Function", FunctionSelect, g],
Controls.ClickButton["Axes-On", ToggleAxes, g]
],
controls: LIST[
g.camera.proxySelect,
g.camera.proxy.xMov, g.camera.proxy.yMov, g.camera.proxy.zMov,
g.camera.proxy.xRot, g.camera.proxy.yRot, g.camera.proxy.zRot,
g.xMin, g.xMax, g.yMin, g.yMax, g.zScale, g.grain],
controlSizes: [17, 200, 60, 20, 60, 150, 150],
graphicsHeight: 300,
drawProc: DrawProc,
typescriptHeight: 18,
clientData: g];
g.outer м g.outerData.parent;
g.graphics м g.outerData.graphics;
};

Controller: Controls.ControlProc ~ {
g: GraphData м NARROW[control.clientData];
Set: PROC [control: Control, value: REAL] ~ {Controls.SetSliderDialValue[control, value]};
SELECT control FROM
g.xMin	=> IF g.xMin.value > g.xMax.value THEN Set[g.xMax, g.xMin.value];
g.yMin	=> IF g.yMin.value > g.yMax.value THEN Set[g.yMax, g.yMin.value];
g.xMax	=> IF g.xMax.value < g.xMin.value THEN Set[g.xMin, g.xMax.value];
g.yMax	=> IF g.xMax.value < g.xMin.value THEN Set[g.yMin, g.yMax.value];
ENDCASE;
Set[g.grain, MAX[g.grain.value, (g.xMax.value-g.xMin.value)/(arraySize-1), (g.yMax.value-g.yMin.value)/(arraySize-1)]];
IF control.whatChanged = $TypedIn OR control.mouse.button = right
THEN Repaint[g, $Controller];
};
InterpressOut: ClickProc ~ {
g: GraphData м NARROW[clientData];
Draw2d.IPOut[Controls.TypescriptReadFileName[g.outerData.typescript], DrawAction, g];
};

ToggleAxes: ClickProc ~ {
g: GraphData м NARROW[clientData];
Controls.ButtonToggle[g.outerData, g.axesOn м NOT g.axesOn, "Axes-On", "Axes-Off"];
IF mouseButton = blue THEN Repaint[g, $Axes];
};

FunctionSelect: ClickProc ~ {
Reverse: PROC [in: LIST OF Request] RETURNS [out: LIST OF Request] ~ {
FOR r: LIST OF Request м in, r.rest WHILE r # NIL DO out м CONS[r.first, out]; ENDLOOP;
};
g: GraphData м NARROW[clientData];
n, choice: INTEGER м 1;
requests: LIST OF Request м NIL;
FOR l: LIST OF Function м functions, l.rest WHILE l # NIL DO
requests м CONS[[l.first.name, l.first.use], requests];
ENDLOOP;
choice м Controls.PopUpRequest[["Function Select"], Reverse[requests]];
FOR l: LIST OF Function м functions, l.rest WHILE l # NIL DO
IF n = choice THEN {
g.graphProc м l.first.proc;
g.outer.name м Rope.Concat["3dGraph ", l.first.name];
ViewerOps.PaintViewer[g.outer, caption];
EXIT;
};
n м n+1;
ENDLOOP;
Repaint[g, $FunctionSelect];
};
CameraControl: Controls.ControlProc ~ {
g: GraphData м NARROW[clientData];
IF control.mouse.button # right THEN RETURN;
Repaint[g, $Camera];
};
Repaint: PROC [g: GraphData, whatChanged: REF ANY м NIL] ~ {
IF whatChanged # $Camera THEN ComputeValues[g];
ViewerOps.PaintViewer[g.graphics, client, FALSE, whatChanged];
};

DrawProc: Controls.DrawProc ~ {
g: GraphData ~ NARROW[clientData];
Action: PROC ~ {
g.view м G3dControl.InitContext[context, g.camera, viewer, , g.view];
DrawAction[context, g, whatChanged, viewer];
};
Draw2d.DoWithBuffer[context, Action];
};

ComputeValues: PROC [g: GraphData] ~ {
ix, iy: NAT м 0;
FOR x: REAL м g.xMin.value, x+g.grain.value WHILE x < g.xMax.value DO
xValue: REAL м 0.0;
xValue м g.graphProc[x, g ! Real.RealException => CONTINUE];
iy м 0;
FOR y: REAL м g.yMin.value, y+g.grain.value WHILE y < g.yMax.value DO
yValue: REAL м 0.0;
yValue м g.graphProc[y, g ! Real.RealException => CONTINUE];
g.values[ix][iy] м g.zScale.value*xValue*yValue;
iy м iy+1;
ENDLOOP;
ix м ix+1;
ENDLOOP;
};

DrawAction: Controls.DrawProc ~ {
g: GraphData ~ NARROW[clientData];
ix, iy: NAT м 0;
FOR x: REAL м g.xMin.value, x+g.grain.value WHILE x < g.xMax.value DO
p0: Triple м [x, g.yMin.value, g.values[ix][iy м 0]];
FOR y: REAL м g.yMin.value+g.grain.value, y+g.grain.value WHILE y < g.yMax.value DO
p1: Triple м [x, y, g.values[ix][iy м iy+1]];
G3dDraw.Segment[context, p0, p1, g.view, []];
p0 м p1;
ENDLOOP;
ix м ix+1;
ENDLOOP;
iy м 0;
FOR y: REAL м g.yMin.value, y+g.grain.value WHILE y < g.yMax.value DO
p0: Triple м [g.xMin.value, y, g.values[ix м 0][iy]];
FOR x: REAL м g.xMin.value+g.grain.value, x+g.grain.value WHILE x < g.xMax.value DO
p1: Triple м [x, y, g.values[ix м ix+1][iy]];
G3dDraw.Segment[context, p0, p1, g.view, []];
p0 м p1;
ENDLOOP;
iy м iy+1;
ENDLOOP;
IF g.axesOn THEN G3dDraw.Axes[context, g.view, []];
};
Spike:		GraphProc ~ {v м G3dFunction.Spike[x]};
Gauss:		GraphProc ~ {v м G3dFunction.Gauss[x]};
Poisson:		GraphProc ~ {v м G3dFunction.Poisson[x, g.a.value]};
Power:		GraphProc ~ {v м IF x>0 THEN RealFns.Power[x, g.a.value] ELSE 0};
Sin:			GraphProc ~ {v м RealFns.Power[0.5+0.5*RealFns.Sin[x], g.a.value]};
Ln:			GraphProc ~ {v м RealFns.Ln[x]};
Log:			GraphProc ~ {v м RealFns.Log[g.a.value, x]};
Perlin:		GraphProc ~ {v м G3dFunction.Perlin[x, g.a.value]};
Wyvill:		GraphProc ~ {v м IF x>0 THEN G3dFunction.Wyvill[x/g.a.value] ELSE 0};
XSquared:	GraphProc ~ {v м IF x>0 THEN 1-x*x ELSE 0};
SlowInOut:	GraphProc ~ {v м IF x>0 THEN G3dSpline.SlowInOut[0.0, 1.0, x] ELSE 0};
functions:		LIST OF Function ~ LIST[
["Spike",					"a second order curve", 					Spike],
["Gauss",					"the normal distribution curve",			Gauss],
["Poisson",				"the poisson curve, with a the scalar",	Poisson],
["Power",					"a power curve, with a the exponent",	Power],
["Sin",					"a sin curve raised to the a power",		Sin],
["Ln",						"natural logarithm",							Ln],
["Log",					"logarithm to the base a",					Log],
["Perlin",					"Ken's curve, with a the scalar",			Perlin],
["Wyvill",				"Wyvill's soft function",					Wyvill],
["xSquared",				"simplifiied Wyvill",						XSquared],
["SlowInOut",			"Slow in and out curve",					SlowInOut]];

use: ROPE м "graph <function> [-option]\n\tfunctions are:\n";

FOR l: LIST OF Function м functions, l.rest WHILE l # NIL DO
use м IO.PutFR["%g\t\t%g:\t%g\n", IO.rope[use], IO.rope[l.first.name], IO.rope[l.first.use]];
ENDLOOP;

use м Rope.Concat[use, "\toptions are:
		-scale	<REAL> (default = 1)
		-xMin <REAL> (default = 0)
		-xMax <REAL> (default = 1)
		-a		<REAL> (default = 1)"];

G3dTool.Register["Graph", Graph, use];

END.

   L  
G3dGraphCmdImpl.mesa
Copyright ╙ 1988, 1992 by Xerox Corporation.  All rights reserved.
Bloomenthal, October 20, 1992 5:38 pm PDT
Glassner, September 30, 1988 12:41:52 pm PDT
Imported Types
Local Types
Viewer Setup
biScrollable: TRUE,
Options
View Control
LF Display
G3dControl.UpdateCamera[g.camera];
GraphProcs
Start Code
 Э╩:  ХNewlineDelimiter
Ц"cedarcode" styleЩ ЩJЪЬ	╧eЬ6ЩBJЩ)JЩ,JШ JЪ╧k	ЬqЮЬ!ШЭJШ ЧЪ╨bnЬЮЬЮШJЪЮЬQЮЬ ШzЧJЪЬЮШheadЪ╧lЩJЪЮЬЮЬЮЬШJЪЬЮЬШ#JЪЬЮЬШ(JЪЬЮЬШ(JЪЬЮЬШ%JЪЬЮЬШ%JЪЬЮЬШ#JЪЬЮЬШ,JЪЬЮЬШ$ЧЪаЩЪЬ
ЮЬШJШ ЧJЪЬЮЬЮЬШ$ЪЬЮЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬЮЬЮЬЮЬЮЬШAJЪЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬШJЪЬЮЬЮЬШJШJШJШ ЧJЪЬЮЬЮЬЮЬЮЬЮЬШBJЪЬЮЬЮЬЮЬЮЬЮЬШ3JЪ	ЬЮЬЮЬЮЬЮЬШCЧЪаЩЪ╧bЬШ JЪЬЮЬШ!JШLJЪЬ5Ш5JШ6JШFJШFJШFJШFJШJJШHJЪЬШШ#JШЪЬ	ЮЬШJШ9JШ4JШ.JШЧЪЬ
ЮЬШJШJШ>JШ>JШ3ЧJШ.JШJШJШJЪЬЮЬЩJШЧJШJШ"JШJШ ЧЪ╧n
ЬШ$JЪЬЮЬШ*JЪвЬЮЬЮЬ2ШZЪЮЬ	ЮШJЪЬ
ЮЬЮЬШHJЪЬ
ЮЬЮЬШHJЪЬ
ЮЬЮЬШHJЪЬ
ЮЬЮЬШHJЪЮЬШЧJЪЬЮЬgШwЪЮЬ ЮЬШAJЪЮЬШЧJШЧЧЪаЩЪвЬШJЪЬЮЬШ"JШUJШJШ ЧЪб
ЬШJЪЬЮЬШ"JЪЬ-ЮЬ"ШSJЪЮЬЮЬШ-JШJШ ЧЪбЬШЪвЬЮЬЮЬЮЬ
ЮЬЮЬЮЬШFJЪЮЬЮЬЮЬЮЬЮЬЮЬЮЬЮЬШWJШЧJЪЬЮЬШ"JЪЬЮЬШJЪЬ
ЮЬЮЬЮЬШ ЪЮЬЮЬЮЬЮЬЮЬЮШ<JЪЬЮЬ(Ш7JЪЮЬШЧJШGЪЮЬЮЬЮЬЮЬЮЬЮШ<ЪЮЬЮЬШJШJШ5JШ(JЪЮЬШJШЧJШJЪЮЬШЧJЪЬШJШЧЧЪаЩЪбЬШ'JЪЬЮЬШ"JЪЮЬЮЬЮЬШ,JШJШЧЧЪа
Щ
Ъ
вЬЮЬЮЬЮЬЮЬШ<JЪЮЬЮЬШ/JЪЬ*ЮЬШ>JШJШ ЧЪбЬШJЪЬЮЬШ"ЪвЬЮЬШJЩ"JШEJШ,JШЧJШ%JШJШ ЧЪвЬЮЬШ&JЪЬЮЬШЪЮЬЮЬ!ЮЬЮШEJЪЬЮЬШJЪЬ2ЮЬШ<JШЪЮЬЮЬ!ЮЬЮШEJЪЬЮЬШJЪЬ2ЮЬШ<JШ0JШ
JЪЮЬШЧJШ
JЪЮЬШЧJШJШ ЧЪв
ЬШ!JЪЬЮЬШ"JЪЬЮЬШЪЮЬЮЬ!ЮЬЮШEJШ5ЪЮЬЮЬ/ЮЬЮШSJШ-JШ-JШJЪЮЬШЧJШ
JЪЮЬШЧJШЪЮЬЮЬ!ЮЬЮШEJШ5ЪЮЬЮЬ/ЮЬЮШSJШ-JШ-JШJЪЮЬШЧJШ
JЪЮЬШЧJЪЮЬ
ЮЬ#Ш3JШЧЧЪа
Щ
JЪЯЬ*Ш/JЪЯЬ*Ш/JЪЯЬ7Ш>JЪЯЬЮЬЮЬЮЬШIJЪЯЬGШJJЪЯЬ$Ш&JЪЯЬ0Ш3JЪЯЬ6Ш<JЪЯЬЮЬЮЬ!ЮЬШNJЪвЬЮЬЮЬЮЬШ5JЪЯ	ЬЮЬЮЬ"ЮЬШQЧЪа
Щ
ЪЬЮЬЮЬЮЬШ$JШ2JШ8JШ@JШ;JШ7JШ+JШ0JШ;JШ4JШ5JШ9JШ ЧJЪЬЮЬ4Ш=JШ ЪЮЬЮЬЮЬЮЬЮЬЮШ<JЪ	ЬЮЬЮЬЮЬЮЬШ]JЪЮЬШЧJШ JЪ	Ь1╧sЬгЬгЬгЬШЬJШ JШ&JШ ЧЪЮЬШJШ ЧЧ ЕЧ    ╬  )T  