-- File: [Cherry]<Thyme>System>C03>spOutput.mesa -- Last editted: -- Wilhelm April 6, 1982 11:36 AM, reformated by Barth and stored under -- [Cherry]<Barth>Thyme>1.97> . DIRECTORY spGlobals, AltoDefs, CWF, Real, RealFns, plotDefs, DisplayDefs, StreamDefs, TimeDefs, StringDefs, AbsAllocDefs; spOutput: PROGRAM IMPORTS spGlobals, CWF, P: plotDefs, DisplayDefs, Str: StreamDefs, TD: TimeDefs, StringDefs, AbsAllocDefs EXPORTS spGlobals = BEGIN OPEN spGlobals; t0: REAL; argumentPtr: TYPE = LONG POINTER TO argument; argument: TYPE = RECORD[node: nodePtr, branch: branchPtr, getComVal: BOOLEAN, dvdt: BOOLEAN, current: BOOLEAN]; printBlkPtr: TYPE = LONG POINTER TO printBlk; printBlk: TYPE = RECORD[nextBlk: printBlkPtr, relativeTo: REAL, arg: argument]; plotBlkPtr: TYPE = LONG POINTER TO plotBlk; plotBlk: TYPE = RECORD[nextBlk: plotBlkPtr, graph: P.GraphHandle, relativeTo: REAL, arg: argument]; maxPrints: CARDINAL = 10; numPrints: CARDINAL ← 0; prints: ARRAY [0..maxPrints) OF RECORD[printList: printBlkPtr, printTimeRelTo: REAL, printStream: Str.StreamHandle]; maxPlots: CARDINAL = 10; numPlots: CARDINAL ← 0; plots: ARRAY [0..maxPlots) OF RECORD[plotList: plotBlkPtr, title: LONG STRING, plotTimeRelTo: REAL, yMin, yMax: REAL, plotCount: CARDINAL]; displayBlanked: BOOLEAN; systemWindow: P.TWHandle ← NIL; thymeLog: Str.StreamHandle ← NIL; time: STRING = [40]; bwDisplay: BOOLEAN ← FALSE; makePrintBlk: PROCEDURE RETURNS[pb: printBlkPtr] = BEGIN pb ← AbsAllocDefs.Allocate[SIZE[printBlk]]; pb↑ ← [NIL, 1.0, [NIL, NIL, FALSE, FALSE, FALSE]] END; makePlotBlk: PROCEDURE RETURNS[pb: plotBlkPtr] = BEGIN pb ← AbsAllocDefs.Allocate[SIZE[plotBlk]]; pb↑ ← [NIL, NIL, 1.0, [NIL, NIL, FALSE, FALSE, FALSE]] END; restoreDisplay: PUBLIC PROCEDURE = BEGIN IF displayBlanked THEN DisplayDefs.DisplayOn[] END; killPlot: PUBLIC PROCEDURE = BEGIN i: CARDINAL; time: STRING = [40]; time.length ← 0; TD.AppendDayTime[time, TD.UnpackDT[TD.CurrentDayTime[]]]; CWF.FWF1[thymeLog, "*n*n%s", time]; CWF.WF1["*n*n%s", time]; IF thymeLog # NIL THEN thymeLog.destroy[thymeLog]; FOR i IN [0..numPrints) DO prints[i].printStream.destroy[prints[i].printStream] ENDLOOP; P.CloseWindow[] END; openSysWindow: PUBLIC PROCEDURE[herald: STRING] = BEGIN logName: STRING = [40]; CWF.SWF1[logName, "%s.log", inputFileName]; thymeLog ← Str.NewByteStream[logName, Str.WriteAppend]; thymeLog.reset[thymeLog]; time.length ← 0; TD.AppendDayTime[time, TD.UnpackDT[TD.CurrentDayTime[]]]; CWF.FWF2[thymeLog, "%s*n*n%s*n*n", herald, time]; CWF.WF2["%s*n*n%s*n*n", herald, time] END; printSysWindow: PUBLIC PROCEDURE[s: STRING] = BEGIN IF systemWindow # NIL THEN P.WriteLine[s, systemWindow]; CWF.FWF1[thymeLog, "%s*n", s] END; initPlot: PUBLIC PROCEDURE[tMin, tMax: REAL, blank: BOOLEAN] = BEGIN window: P.WindowHandle; plotWindow: P.PWHandle; plotColor: P.ColorHandle; i, color: CARDINAL; pl: plotBlkPtr; cy, cyTitle, h, w: INTEGER; w2, wTitle: P.TWHandle; s: STRING = [20]; displayBlanked ← blank AND ~bwDisplay; IF displayBlanked THEN DisplayDefs.DisplayOff[white]; t0 ← tMin; P.SetTempFileName[inputFileName]; FOR i IN [0..numPlots) DO window ← P.OpenWindow[NIL, TRUE, ~bwDisplay, FALSE]; [w, h] ← P.GetWindowSize[window]; cy ← P.HeightOfLines[MAX[plots[i].plotCount, 6]]; cyTitle ← P.HeightOfLines[3]; w2 ← P.GetTextWindow[1500, h - cy, w/2, h - 15,,,window]; wTitle ← P.GetTextWindow[w/2, h - cyTitle, w, h - 15,,,window]; P.WriteString[inputFileName, wTitle]; P.WriteString[" ", wTitle]; P.WriteLine[time, wTitle]; IF plots[i].title # NIL THEN P.WriteLine[plots[i].title, wTitle]; plotWindow ← P.GetPlotWindow[0, 0, w, h - cy - 50, tMin/plots[i].plotTimeRelTo, tMax/plots[i].plotTimeRelTo, plots[i].yMin, plots[i].yMax, , window]; IF i = 0 THEN systemWindow ← P.GetTextWindow[w/2, h - cy, w, h - cyTitle,, FALSE]; color ← 0; pl ← plots[i].plotList; UNTIL pl = NIL DO plotColor ← P.GetStandardColor[color]; P.SetWindowColor[plotColor, w2]; pl↑.graph ← SELECT color/8 FROM 0 => P.GetGraphHandle[plotWindow, plotColor], 1 => P.GetGraphHandle[plotWindow, plotColor,,, dashed], ENDCASE => P.GetGraphHandle[plotWindow, plotColor,,, dotted]; color ← color + 1; IF pl↑.arg.getComVal THEN P.WriteString["Value of ", w2] ELSE IF pl↑.arg.current THEN P.WriteString["Current through ", w2] ELSE IF pl↑.arg.dvdt THEN P.WriteString["Derivative at ", w2] ELSE P.WriteString["Voltage at ", w2]; P.WriteString[makeStringNB[pl↑.arg.node, pl↑.arg.branch], w2]; IF pl↑.relativeTo # 1.0 THEN BEGIN CWF.SWF1[s, "/%-9.2f", @pl↑.relativeTo]; P.WriteString[s, w2] END; P.WriteLine["", w2]; pl ← pl↑.nextBlk ENDLOOP ENDLOOP END; evaluateArgument: PROCEDURE[a: argumentPtr] RETURNS[iv: REAL] = BEGIN b: branchPtr; IF a↑.getComVal THEN RETURN[a↑.branch↑.comVal]; IF a↑.node # NIL THEN IF a↑.dvdt THEN iv ← a↑.node↑.nHist.f0 ELSE iv ← a↑.node↑.nHist.y ELSE BEGIN b ← a.branch; IF a↑.current THEN WITH b↑ SELECT FROM resistor => iv ← (posNode↑.nHist.y - negNode↑.nHist.y)/comVal; capacitor => iv ← (posNode↑.nHist.f0 - negNode↑.nHist.f0)*comVal; inductor => iv ← iHist.y/comVal; vSource => iv ← vsCurrent; iSource => iv ← comVal ENDCASE ELSE iv ← b↑.posNode↑.nHist.y - b↑.negNode↑.nHist.y END END; plotFromList: PUBLIC PROCEDURE[t: REAL] = BEGIN i: CARDINAL; p: plotBlkPtr; r: REAL; FOR i IN [0..numPlots) DO p ← plots[i].plotList; UNTIL p = NIL DO r ← evaluateArgument[@p↑.arg]/p↑.relativeTo; IF t > t0 THEN P.DrawGraph[p↑.graph, t/plots[i].plotTimeRelTo, r] ELSE P.PositionPen[p↑.graph, t0/plots[i].plotTimeRelTo, r]; p ← p↑.nextBlk ENDLOOP ENDLOOP END; printFromList: PUBLIC PROCEDURE[ni: CARDINAL, t: REAL, printStep: BOOLEAN] = BEGIN i: CARDINAL; p: printBlkPtr; r, trel: REAL; FOR i IN [0..numPrints) DO p ← prints[i].printList; trel ← t/prints[i].printTimeRelTo; CWF.FWF1[prints[i].printStream, "%9.3f", @trel]; IF printStep THEN CWF.FWF1[prints[i].printStream, "(%2d)", @ni]; UNTIL p = NIL DO r ← evaluateArgument[@p↑.arg]/p↑.relativeTo; CWF.FWF1[prints[i].printStream, " %9.3f", @r]; p ← p↑.nextBlk ENDLOOP; CWF.FWFCR[prints[i].printStream] ENDLOOP END; scale: PROCEDURE RETURNS[s: REAL ← 1.0] = BEGIN IF item = colon THEN BEGIN next[]; s ← getSignedNumber[] END END; makePrintList: PUBLIC PROCEDURE = BEGIN I, D, cv: BOOLEAN; pb, pb2, pb3: printBlkPtr; n: nodePtr; b: branchPtr; name: STRING = [40]; IF numPrints >= maxPrints THEN BEGIN error[631]; numPrints ← numPrints - 1 END; IF item = leftB THEN next[] ELSE error[600]; IF item = colon THEN BEGIN prints[numPrints].printTimeRelTo ← scale[]; IF item = comma THEN next[] ELSE error[603,, FALSE] END ELSE prints[numPrints].printTimeRelTo ← 1.0; prints[numPrints].printList ← NIL; UNTIL item # name DO [n, b] ← findNodeOrBranch[]; IF n = NIL AND b = NIL THEN error[620] ELSE BEGIN I ← b # NIL AND item = upArrow; D ← n # NIL AND item = quote; cv ← b # NIL AND item = atSign; IF I OR D OR cv THEN next[]; pb ← makePrintBlk[]; pb↑ ← [prints[numPrints].printList, scale[], [n, b, cv, D, I]]; prints[numPrints].printList ← pb END; IF item = comma THEN next[] ELSE EXIT ENDLOOP; pb ← NIL; pb3 ← prints[numPrints].printList; UNTIL pb3 = NIL DO pb2 ← pb3; pb3 ← pb3↑.nextBlk; pb2↑.nextBlk ← pb; pb ← pb2 ENDLOOP; prints[numPrints].printList ← pb; CWF.SWF2[name, "%s.out%d", inputFileName, @numPrints]; prints[numPrints].printStream ← Str.NewByteStream[name, Str.WriteAppend]; prints[numPrints].printStream.reset[prints[numPrints].printStream]; numPrints ← numPrints + 1; IF item = rightB THEN next[] ELSE error[601, TRUE] END; getPlotNum: PROCEDURE RETURNS[r: REAL ← 1.0] = BEGIN r ← getSignedNumber[]; IF item = comma THEN next[] ELSE error[603,, FALSE] END; makePlotList: PUBLIC PROCEDURE[bw: BOOLEAN] = BEGIN I, D, cv: BOOLEAN; pb: plotBlkPtr; n: nodePtr; b: branchPtr; bwDisplay ← bw; IF numPlots >= maxPlots THEN BEGIN error[632, FALSE]; numPlots ← numPlots - 1 END; IF item = leftB THEN next[] ELSE error[600]; IF item = string THEN BEGIN plots[numPlots].title ← newString; next[]; IF item = comma THEN next[] ELSE error[603,, FALSE] END ELSE plots[numPlots].title ← NIL; IF item = colon THEN BEGIN plots[numPlots].plotTimeRelTo ← scale[]; IF item = comma THEN next[] ELSE error[603,, FALSE] END ELSE plots[numPlots].plotTimeRelTo ← 1.0; plots[numPlots].yMin ← getPlotNum[]; plots[numPlots].yMax ← getPlotNum[]; IF plots[numPlots].yMin >= plots[numPlots].yMax THEN error[630, FALSE]; plots[numPlots].plotList ← NIL; plots[numPlots].plotCount ← 0; UNTIL item # name DO plots[numPlots].plotCount ← plots[numPlots].plotCount + 1; [n, b] ← findNodeOrBranch[]; IF n = NIL AND b = NIL THEN error[620] ELSE BEGIN I ← b # NIL AND item = upArrow; D ← n # NIL AND item = quote; cv ← b # NIL AND item = atSign; IF I OR D OR cv THEN next[]; pb ← makePlotBlk[]; pb↑ ← [plots[numPlots].plotList, NIL, scale[], [n, b, cv, D, I]]; plots[numPlots].plotList ← pb END; IF item = comma THEN next[] ELSE EXIT ENDLOOP; IF item = rightB THEN next[] ELSE error[601, TRUE]; numPlots ← numPlots + 1 END; dumpAll: PUBLIC PROCEDURE[t: REAL] = BEGIN dump: Str.StreamHandle; dname: STRING = [40]; nodes: nodePtr ← nodeList; inds: inductorPtr ← inductorList; dname.length ← 0; StringDefs.AppendString[dname, "Dumped "]; TD.AppendDayTime[dname, TD.UnpackDT[TD.CurrentDayTime[]]]; printSysWindow[dname]; CWF.SWF1[dname, "%s.dump", inputFileName]; dump ← Str.NewByteStream[dname, Str.WriteAppend]; dump.reset[dump]; CWF.FWF1[dump, "ic[%f,*n", @t]; UNTIL nodes = NIL DO CWF.FWF2[dump, "%s←%f", makeStringNB[nodes, NIL, FALSE], @nodes↑.nHist.y]; nodes ← nodes↑.nextNode; IF nodes # NIL THEN CWF.FWF[dump, ",*n"] ELSE CWF.FWFCR[dump] ENDLOOP; UNTIL inds = NIL DO CWF.FWF2[dump, "%s←%f", makeStringNB[NIL, inds, FALSE], @inds↑.iHist.y]; inds ← inds↑.nextInductor; IF inds # NIL THEN CWF.FWF[dump, ",*n"] ELSE CWF.FWFCR[dump] ENDLOOP; CWF.FWF[dump, "];*n"]; dump.destroy[dump] END; END.