DIRECTORYBasics USING [LowHalf],BasicTime USING [Now],Commander USING [CommandProc, Register],CommandTool USING [ParseToList],Convert USING [CardFromRope, Error],FileNames USING [CurrentWorkingDirectory],FS USING [ComponentPositions, EnumerateForNames, ExpandName, NameProc, Error, StreamOpen],IO USING [Close, EndOfStream, Error, GetCard, GetReal, GetRopeLiteral, int, PutF, PutFR, STREAM],List USING [DReverse],Plot USING [AddVector, CreateViewer, Curves, PlotSpec, PlotSpecRec, RealSequence, RopeSequence, Vector],Real USING [LargestNumber, SmallestNormalizedNumber],Rope USING [Cat, Equal, Fetch, Find, IsEmpty, Length, ROPE, Substr],ViewerClasses USING [Viewer];SimplePlotHack: CEDAR PROGRAMIMPORTS Basics, BasicTime, Commander, CommandTool, Convert, FileNames, FS, IO, List, Plot, Rope = {OPEN IO, Plot;RopeList: TYPE = LIST OF Rope.ROPE;DataError: SIGNAL[reason: Rope.ROPE _ NIL] = CODE;SPlot: Commander.CommandProc = {argList: RopeList;length: NAT;[argList, length] _ CommandTool.ParseToList[cmd];IF length <= 0 THEN msg _ "To plot data in files, type: SPlot <files>"ELSE {wDir: Rope.ROPE _ FileNames.CurrentWorkingDirectory[];nPlots: NAT _ 0;maxPlots: NAT _ 12;allVersions: BOOL _ FALSE;FOR arg: RopeList _ argList, arg.rest UNTIL arg = NIL OR msg # NIL DOIF arg.first.Fetch[0] = '- THEN {IF arg.first.Length[] >= 2 THEN SELECT arg.first.Fetch[1] FROM'A, 'a => allVersions _ TRUE;'H, 'h => allVersions _ FALSE;IN ['0..'9] => {max: NAT _ Convert.CardFromRope[arg.first.Substr[1]! Convert.Error => max _ maxPlots]; -- don't change itmaxPlots _ max;};ENDCASE; -- simply ignors illegal switches}ELSE {fileList: RopeList _ FileListFrom[arg.first, wDir, allVersions];FOR file: RopeList _ fileList, file.rest UNTIL file = NIL OR msg # NIL DOmsg _ IF nPlots >= maxPlots THENIO.PutFR["Note: max. %g plots for each command.", IO.int[maxPlots]]ELSE SimplePlot[file.first, nPlots # 0, cmd.out];nPlots _ nPlots + 1;ENDLOOP;};ENDLOOP;};}; -- SPlotFileListFrom: PROC [pattern, wDir: Rope.ROPE _ NIL, allVersions: BOOL _ FALSE] RETURNS [fileList: RopeList _ NIL] = {root, lastRoot: Rope.ROPE _ NIL;LinkIt: FS.NameProc -- PROC [fullFName] RETURNS [continue: BOOL] -- = {excl: INT _ fullFName.Find["!"];continue _ TRUE;IF fullFName.Substr[excl-6, 6].Equal[".press", FALSE] THEN RETURN;IF ~allVersions THEN {root _ fullFName.Substr[0, excl];IF root.Equal[lastRoot, FALSE] THEN {fileList.first _ fullFName;RETURN;};lastRoot _ root;};fileList _ CONS[fullFName, fileList];}; -- LinkItFS.EnumerateForNames[pattern, LinkIt, wDir];TRUSTED {fileList _ LOOPHOLE[List.DReverse[LOOPHOLE[fileList]]]};}; -- FileListFromSimplePlot: PROC[file: Rope.ROPE _ NIL, iconic: BOOL _ TRUE, out: IO.STREAM]RETURNS[msg: Rope.ROPE _ NIL] = {s: IO.STREAM _ NIL;viewer: ViewerClasses.Viewer;curves: Curves _ NIL;plotSpec: PlotSpec;ok: BOOL _ TRUE;fullName, status: Rope.ROPE _ NIL;xmin, xmax, ymin, ymax: REAL;IF file.IsEmpty[] THEN RETURN["No file name."];[fullName, ]_ FS.ExpandName[file, FileNames.CurrentWorkingDirectory[]];s _ FS. StreamOpen[fullName ! FS.Error => {status _ error.explanation; ok _ FALSE; CONTINUE} ];IF ok THEN {ENABLE {DataError => {msg _ reason;ok _ FALSE;CONTINUE;};IO.EndOfStream => {msg _ Rope.Cat["End of file reached", status];ok _ FALSE;CONTINUE;};IO.Error => {msg _ Rope.Cat[SELECT ec FROMSyntaxError => "Syntax error",Overflow => "Overflow in input conversion"ENDCASE => IO.PutFR["IO Error # %g", IO.int[LOOPHOLE[ec, CARDINAL]]],status];ok _ FALSE;CONTINUE;};ABORTED => {msg _ "SPlot aborted. ... ";ok _ FALSE;CONTINUE;};};nVector: CARDINAL;r: REAL;plotSpec _ NEW[PlotSpecRec _ [file: fullName, time: BasicTime.Now[]]];out.PutF[Rope.Cat["Reading ", fullName, "\n"]];status _ " in getting the title.";plotSpec.title _ s.GetRopeLiteral[];status _ " in getting the number of curves.";plotSpec.nCurvesMax _ Basics.LowHalf[s.GetCard[]];plotSpec.legendEntries _ NEW[RopeSequence[plotSpec.nCurvesMax]];FOR i: CARDINAL IN [0..plotSpec.nCurvesMax) DOstatus _ IO.PutFR[" in getting the name for curve #%g.", IO.int[i+1]];plotSpec.legendEntries[i] _ s.GetRopeLiteral[];ENDLOOP;status _ " in getting the number of rows.";nVector _ Basics.LowHalf[s.GetCard[]];xmin _ ymin _ Real.LargestNumber;xmax _ ymax _ Real.SmallestNormalizedNumber;FOR i: CARDINAL IN [0..nVector) DOvector: Vector _ NEW[RealSequence[plotSpec.nCurvesMax + 1]];FOR j: CARDINAL IN [0..plotSpec.nCurvesMax] DOENABLE {IO.Error => {msg _ Rope.Cat[SELECT ec FROMSyntaxError => "Syntax error",Overflow => "Overflow"ENDCASE => IO.PutFR["IO Error # %g", IO.int[LOOPHOLE[ec, CARDINAL]]],IO.PutFR[" in getting the %g-th number on the %g-th row in the table.",IO.int[j+1], IO.int[i+1]]];GOTO exit;};IO.EndOfStream => {msg _ Rope.Cat["End of file reached",IO.PutFR[" in getting the %g-th number on the %g-th row in the table.",IO.int[j+1], IO.int[i+1]]];GOTO exit;};};vector[j] _ r _ s.GetReal[];IF j # 0 THEN {ymax _ MAX[ymax, r];ymin _ MIN[ymin, r];};ENDLOOP;IF i = 0 THEN xmin _ vector[0]ELSE IF i+1 = nVector THEN xmax _ vector[0];curves _ CONS[vector, curves];out.PutF["."];REPEATexit => ok _ FALSE;ENDLOOP;};IF ok THEN {cv: Curves _ NIL;IF ymax = ymin THEN {IF ymax = 0.0 THEN { ymax _ 1.0; ymin _ -1.0 }ELSE { ymax _ ymax + ABS[ymax]; ymin _ ymin - ABS[ymin]; };};IF xmax = xmin THEN {IF xmax = 0.0 THEN { xmax _ 1.0; xmin _ -1.0 }ELSE { xmax _ xmax + ABS[xmax]; xmin _ xmin - ABS[xmin]; };};out.PutF[" plotting ... "];plotSpec.bounds _ [xmin, ymin, xmax, ymax];viewer _ CreateViewer[spec: plotSpec, iconic: iconic];FOR c: Curves _ curves, c.rest UNTIL c = NIL DOcv _ CONS[c.first, cv];ENDLOOP;FOR c: Curves _ cv, c.rest UNTIL c = NIL DOAddVector[viewer, c.first];ENDLOOP;out.PutF["done.\n"];};IF s # NIL THEN s.Close[];}; -- SimplePlotCommander.Register["SPlot", SPlot,"SPlot <data files>, use a simple minded plot routine to plot data in the data files that have the following special format:-- comment nodes are ignored.-- carriage returns are treated the same as spaces.<title> -- a text string enclosed by double quotes.<xmin> <ymin> <xmax> <ymax> -- four real numbers.<number of curves> -- Curves are functions of x.<names of the curves> -- Each name is a text string enclosed by double quotes.<number of rows>-- Each row in the table below has a leading x value followed by the corresponding y values.-- If there are n curves then the table for m different x values will be an m x n+1 matrix as follows.x1  y1(x1)  y2(x1)  y3(x1)  ...  yn(x1)x2  y1(x2)  y2(x2)  y3(x2)  ...  yn(x2)... ... ... ... ... ... ... ... ... ... ... ... ... ... xm  y1(xm)  y2(xm)  y3(xm)  ...  yn(xm)"];}.   ╪  SimplePlotHack.mesaLast Edited by: SChen, June 15, 1985 7:37:41 pm PDTSweetsun Chen, August 2, 1985 1:07:43 pm PDTset plotSpecc.first is the last one just added to the list, so we should reverse the order.Э╩o  Ш JЪЬЩЪЬ3Щ3IcodeЩ,ЧJШ Ъ╧k	Ш	JЪЬЭЬШJЪЬ
ЭЬШJЪЬ
ЭЬШ(JЪЬЭЬШ JЪЬЭЬШ$JЪЬ
ЭЬШ*JЪЭЬЭЬRШZJЪЭЬЭЬQЭЬШaJЪЬЭЬШJЪЬЭЬ^ШhJЪЬЭЬ+Ш5JЪЬЭЬ,ЭЬ
ШDJЪЬЭЬ
ШJШ ЧЪЬЭЬЭШJЪЭЬ@ЭЬЭЬШcЧJШ JЪЭЬЭЬШJШ JЪ	Ь
ЭЬЭЬЭЬЭЬШ#JШ JЪ	ЬЭЬЭЬЭЬЭЬШ2JШ ЪЬ Ш JЪЬШJЪЬЭЬШJЪЬ1Ш1JЪЭЬЭЬ3ШFЪЭЬШJЪЬЭЬ'Ш6JЪЬЭЬШJЪЬ
ЭЬШJЪЬЭЬЭЬШЪЭЬ#ЭЬЭЬЭЬЭЬЭШEЪЭЬЭЬШ!ЪЭЬЭЬЭЬЭШ>JЪЬЭЬШJЪЬЭЬШЪЭЬШЪЬЭЬ+Ш3JЪЬ$╧cШ6ЧJШJШЧJЪЭЬЮ!Ш*ЧJШЧЪЭЬШJЪЬ@Ш@ЪЭЬ&ЭЬЭЬЭЬЭЬЭШIЪЬЭЬЭШ JЪЭЬ0ЭЬШCJЪЭЬ-Ш1ЧJЪЬШJЪЭЬШЧJШЧJЪЭЬШЧJЪЬШЧJЪЬЮШJШ ЧЪ╧nЬЭЬЭЬЭЬЭЬЭЬЭЬЭЬШuJЪЬЭЬЭЬШ ЪЬЭЬ
Ю╨ckЮаЮаЮЬШGJЪЬЭЬШ JЪЬЭЬШJЪЭЬ-ЭЬЭЬЭЬШBЪЭЬЭЬШJШ!ЪЭЬЭЬЭЬШ%JЪЬШJЪЭЬШJШЧJЪЬШJШЧJЪЬЭЬШ%JЪЬЮЬШЧJШ JЪЭЬ*Ш,JЪЭЬЭЬЭЬШAJЪЬЮШЧJШ ЪЯ
ЬЭЬЭЬЭЬ
ЭЬЭЬЭЬЭЬШLJЪЭЬЭЬЭЬШ!JЪЬЭЬЭЬЭЬШJШJЪЬЭЬШJШJЪЬЭЬЭЬШJЪЬЭЬЭЬШ"JЪЬЭЬШJЪЭЬЭЬЭЬШ/JЪЬЭЬ8ШHJЪ	ЬЭЬЭЬ,ЭЬЭЬШ_ЪЭЬЭЬШЪЭЬШШJЪЬШJЪЬЭЬШJЪЭЬШ	JШЧЪЭЬШJЪЬ.Ш.JЪЬЭЬШJЪЭЬШ	JШЧЪЭЬШЪЬШЪЭЬЭШJШJШ*JЪ
ЭЬЭЬЭЬЭЬЭЬШEЧJШЧJЪЬЭЬШJЪЭЬШ	JШЧЪЭЬШJЪЬШJЪЬЭЬШJЪЭЬШ	JШЧJШЧJЪЬЩJЪЬ	ЭЬШJЪЬЭЬШJЪЬЭЬ8ШFJШ JЪЬ/Ш/JШ"JЪЬ$Ш$JШ JШ-JЪЬ2Ш2JШ JЪЬЭЬ$Ш@ЪЭЬЭЬЭЬЭШ.JЪЬ	ЭЬ.ЭЬШFJЪЬ/Ш/JЪЭЬШЧJШ JШ+JЪЬ&Ш&JШ JЪЬ!Ш!JЪЬ,Ш,ЪЭЬЭЬЭЬЭШ"JЪЬЭЬ(Ш<ЪЭЬЭЬЭЬЭШ.ЪЭЬШЪЭЬШЪЬШЪЭЬЭШJШJШJЪ
ЭЬЭЬЭЬЭЬЭЬШEЧЪЭЬEШGJЪЭЬЭЬШЧЧJЪЭЬШ
JШЧЪЭЬШЪЬШJЪЬШЪЭЬEШGJЪЭЬЭЬШЧЧJЪЭЬШ
JШЧJШЧJЪЬШЪЭЬЭЬШJЪЬЭЬ
ШJЪЬЭЬ
ШJШЧJЪЭЬШЧJЪЭЬЭЬШJЪЭЬЭЬЭЬШ,JЪЬ	ЭЬШJШЪЭШJЪЬЭЬШЧJЪЭЬШЧJШЧJШ ЪЭЬЭЬШJЪЬЭЬШЪЭЬЭЬШJЪЭЬЭЬШ.JЪЭЬЭЬЭЬ
Ш;JШЧЪЭЬЭЬШJЪЭЬЭЬШ.JЪЭЬЭЬЭЬ
Ш;JШЧJШ JШJШ+JЪЬ6Ш6JЩOЪЭЬЭЬЭЬЭШ/JЪЬЭЬШJЪЭЬШЧЪЭЬЭЬЭЬЭШ+JЪЬШJЪЭЬШЧJШJШJШ ЧJЪЭЬЭЬЭЬШJЪЬЮШЧJШ ЪЬШJЪЬШJЪЬ┐Ш┐ЧJШ JШJШ Ч ЕЧ    n  #╡  