diff options
author | Karl Berry <karl@freefriends.org> | 2010-06-13 21:40:28 +0000 |
---|---|---|
committer | Karl Berry <karl@freefriends.org> | 2010-06-13 21:40:28 +0000 |
commit | 53e163508344d593ba867c3f16c157c07fba52c6 (patch) | |
tree | 2b4e9b737c872f809d47e24085051d1e9f42b0df /Master/texmf-dist/tex | |
parent | 89005db14615e4cbb9fbee050bab6809fa3e2d97 (diff) |
pst-bspline update (13jun10)
git-svn-id: svn://tug.org/texlive/trunk@18944 c570f23f-e606-0410-a88d-b1316a301751
Diffstat (limited to 'Master/texmf-dist/tex')
-rw-r--r-- | Master/texmf-dist/tex/generic/pst-bspline/pst-bspline.tex | 332 |
1 files changed, 142 insertions, 190 deletions
diff --git a/Master/texmf-dist/tex/generic/pst-bspline/pst-bspline.tex b/Master/texmf-dist/tex/generic/pst-bspline/pst-bspline.tex index 398a42827d1..633a463f474 100644 --- a/Master/texmf-dist/tex/generic/pst-bspline/pst-bspline.tex +++ b/Master/texmf-dist/tex/generic/pst-bspline/pst-bspline.tex @@ -1,63 +1,24 @@ %% BEGIN pst-bspline.tex %% Author: Michael Sharpe (msharpe at ucsd.edu) -\def\fileversion{1.13} -\def\filedate{2008/12/13} +\def\fileversion{1.2} +\def\filedate{2010/06/12} %% -\message{ v\fileversion, \filedate} +%\message{ v\fileversion, \filedate} \csname PSTBsplineLoaded\endcsname \let\PSTBsplineLoaded\endinput +\ProvidesFile{pst-bspline.tex}[\filedate\space v\fileversion\space Bspline routines for pstricks] \ifx\PSTnodesLoaded\endinput\else\input pst-node.tex \fi\relax \edef\TheAtCode{\the\catcode`\@}\catcode`\@=11 % -% Utility macros -% Build the linear combination #3[1]]*#1+#3[[2]]*#2=#4 -% #4 can be same as #1, #2 or #3--eg, \psLCNodeVar(A)(B)(3,4){A}. -% #3 may use Postscript code -\def\psLCNodeVar(#1)(#2)(#3)#4{% -\pst@getcoor{#1}\pst@tempA% -\pst@getcoor{#2}\pst@tempB% -\pnode(#3){tmpLCn@de}% -\pnode(!% - \pst@tempA /YA exch \pst@number\psyunit div def - /XA exch \pst@number\psxunit div def - \pst@tempB /YB exch \pst@number\psyunit div def - /XB exch \pst@number\psxunit div def %stack now empty - \psGetNodeCenter{tmpLCn@de}\space - XA tmpLCn@de.x mul XB tmpLCn@de.y mul add - YA tmpLCn@de.x mul YB tmpLCn@de.y mul add){tmpLCn@deA}% -\pnode(tmpLCn@deA){#4}% -} -% Useful for debugging a node location -\def\shownode(#1){%display node user coords in console -\pnode(! /nodecoors 2 array def \psGetNodeCenter{#1} - nodecoors 1 #1.y put nodecoors 0 #1.x put end nodecoors == #1.x #1.y ){temp@} -} -% Use to construct a sequence of nodes -% Eg, \pnodes{P}(0,1)(2;5)(3,4) defines nodes P0, P1, P2 with respective locations -\def\pnodes#1{% -\newcount\bsp@args\bsp@args=0\relax% -\def\bsp@root{#1}% -\next@Node}% -% -\def\next@Node{\@ifnextchar({\next@Node@i}{\advance\bsp@args by -1\typeout{Defined nodes \bsp@root0 .. \bsp@root\the\bsp@args.}}}% -% -\def\next@Node@i(#1){% -\pnode(#1){\bsp@root\the\bsp@args}% -\advance\bsp@args by 1\relax% -\next@Node} -% -%\def\nextpt#1#2#3{\typeout{Next points:#1:#2:#3:}} -\def\nextfour#1#2#3#4{\newwrite\foo \immediate\openout\foo=foo.txt \immediate\write\foo{Next chars:#1:#2:#3:#4:}\closeout\foo} -%end utility macros -\newcount\PST@args%used in several macros +%\newcount\pst@args%used in several macros--now defined in pst-node.tex \def\PST@root{} \newcount\bsp@args \newcount\bsp@cntA \newcount\bsp@cntB \newif\ifbsp@omitends \newif\ifbsp@closed -\newif\ifbsp@showpoints -\newif\ifbsp@doubleline +%\newif\ifbsp@showpoints +%\newif\ifbsp@doubleline \newif\ifbsp@seq % \def\bsp@init{%these items are used only within the plotting macros @@ -65,217 +26,207 @@ \def\bsp@arrowB{}% }% \def\refreshbspopts{% -\bsp@doublelinefalse -\bsp@closedfalse -\bsp@showpointsfalse -\bsp@omitendsfalse -\bsp@seqfalse +\bsp@closedfalse% +\bsp@omitendsfalse% +\bsp@seqfalse% }% -\def\saveparams{ -\let\save@linecolor\pslinecolor -% \let\save@arrowA\psk@arrowA -% \let\save@arrowB\psk@arrowB - \let\save@linestyle\pslinestyle} -\def\restoreparams{% settings for linestyle, arrows, linecolor may have changed - \let\pslinecolor\save@linecolor -% \let\psk@arrowA\save@arrowA -% \let\psk@arrowB\save@arrowB - \let\pslinestyle\save@linestyle} % showframe key and its default are global -\newif\ifshowframe -\def\psset@showframe#1{\@nameuse{showframe#1}} -\psset@showframe{false} +\newif\ifshowframe% +\def\psset@showframe#1{\@nameuse{showframe#1}}% +\psset@showframe{false}% % %Bspline drawing macros -\def\psBsplineNodes{\pst@object{psBsplineNodes}}% -\def\psBsplineNodesC{\pst@object{psBsplineCNodes}}% -\def\psBsplineNodesE{\pst@object{psBsplineENodes}}% -\def\psBsplineC{\pst@object{psBsplineC}}% -\def\psBsplineE{\pst@object{psBsplineE}}% -\def\psBspline{\pst@object{psBspline}}% +\def\psBsplineNodes{\def\pst@par{}\pst@object{psBsplineNodes}}% +\def\psBsplineNodesC{\def\pst@par{}\pst@object{psBsplineCNodes}}% +\def\psBsplineNodesE{\def\pst@par{}\pst@object{psBsplineENodes}}% +\def\psBsplineC{\def\pst@par{}\pst@object{psBsplineC}}% +\def\psBsplineE{\def\pst@par{}\pst@object{psBsplineE}}% +\def\psBspline{\def\pst@par{}\pst@object{psBspline}}% %Bspline internal macros, level 1 -\def\psBsplineC@i{\bsp@init\refreshbspopts\bsp@closedtrue\ps@bspline} -\def\psBsplineE@i{\bsp@init\refreshbspopts\bsp@omitendstrue\ps@bspline} -\def\psBspline@i{\bsp@init\refreshbspopts\ps@bspline} -\def\psBsplineNodes@i{\bsp@init\refreshbspopts\bsp@seqtrue\ps@bsplineNodes} -\def\psBsplineCNodes@i{\bsp@init\refreshbspopts\bsp@seqtrue\bsp@closedtrue\ps@bsplineNodes} -\def\psBsplineENodes@i{\bsp@init\refreshbspopts\bsp@seqtrue\bsp@omitendstrue\ps@bsplineNodes} +\def\psBsplineC@i{\bsp@init\refreshbspopts\bsp@closedtrue\ps@bspline}% +\def\psBsplineE@i{\bsp@init\refreshbspopts\bsp@omitendstrue\ps@bspline}% +\def\psBspline@i{\bsp@init\refreshbspopts\ps@bspline}% +\def\psBsplineNodes@i{\bsp@init\refreshbspopts\bsp@seqtrue\ps@bsplineNodes}% +\def\psBsplineCNodes@i{\bsp@init\refreshbspopts\bsp@seqtrue\bsp@closedtrue% +\ps@bsplineNodes}% +\def\psBsplineENodes@i{\bsp@init\refreshbspopts\bsp@seqtrue\bsp@omitendstrue% +\ps@bsplineNodes}% %Bspline internal macros, level 2 \def\ps@bsplineNodes#1#2{%noderoot, endindex -\def\bsp@root{#1} -\bsp@args=#2 -\psBsplineMain} +\def\bsp@root{#1}% +\bsp@args=#2\relax% +\next@Point% +}% % \def\ps@bspline{\bsp@args=0\relax% -\@ifnextchar({\psBspline@read{pn@de}}{\psBspline@read}} +\@ifnextchar({\psBspline@read{pn@de}}{\psBspline@read}}% % \def\psBspline@read#1{% \def\bsp@root{#1}% -%\show\bsp@root -\psBsplineMain}% +\next@Point% +}% % -\def\psBsplineMain{% - \ifshowpoints \bsp@showpointstrue \fi% +\def\next@Point{\@ifnextchar({\next@Point@i}{\psBsplineMain}}% +\def\next@Point@i(#1){% +\pnode(#1){\bsp@root\the\bsp@args}% +\advance\bsp@args by \@ne% +\next@Point}% +% The main drawing routine +\def\psBsplineMain{%done reading nodes +\pst@killglue% +\begingroup% +\use@par% \ifbsp@closed\begin@ClosedObj \else \begin@OpenObj \fi% -% \ifshowpoints \bsp@showpointstrue \showpointsfalse \fi -\ifpsdoubleline \bsp@doublelinetrue \fi% \let\bsp@linecolor\pslinecolor% \let\bsp@arrowA\psk@arrowA% \let\bsp@arrowB\psk@arrowB% \let\bsp@linestyle\pslinestyle% -%\tracingcommands=2\tracingmacros=2 -%\ifbsp@seq\next@Point@ii \else -\next@Point%omit the \fi or it becomes the nextchar -}% -% -\def\next@Point{\@ifnextchar({\next@Point@i}{\next@Point@ii}}% -\def\next@Point@i(#1){% -\pnode(#1){\bsp@root\the\bsp@args} -\advance\bsp@args by 1 \relax -\next@Point} -\def\next@Point@ii{%done reading nodes -%\tracingcommands=0\tracingmacros=0 -\ifbsp@seq\advance\bsp@args by 1\fi -\ifbsp@closed - \psset{doubleline=false} - \pnode(\bsp@root0){\bsp@root\the\bsp@args} - \advance\bsp@args by 1 \relax - \pnode(\bsp@root1){\bsp@root\the\bsp@args} - \advance\bsp@args by 1 \relax -\fi -\advance\bsp@args by -1 +\ifbsp@seq\advance\bsp@args by \@ne\fi% +\ifbsp@closed% + \psset{doubleline=false}% + \pnode(\bsp@root0){\bsp@root\the\bsp@args}% + \advance\bsp@args by \@ne% + \pnode(\bsp@root1){\bsp@root\the\bsp@args}% + \advance\bsp@args by \@ne% +\fi% +\advance\bsp@args by \m@ne% % Define the left and right bezier control points in each interval. These are % denoted (root+) R0, L1, R1, L2, R2, etc -\ifshowframe - \psset{linestyle=dashed,linecolor=gray,arrows=*-*} -\else - \psset{linestyle=none,linecolor=white,arrows=-} -\fi -%2008/9/20---change \ncline{*-*} to just \ncline -\multido{\ia=0+1,\ib=1+1}{\bsp@args}{\ncline{\bsp@root\ia}{\bsp@root\ib}\ncput[npos=.3333]{\pnode{\bsp@root R\ia}}\ncput[npos=.6667]{\pnode{\bsp@root L\ib}}} +\multido{\ia=0+1,\ib=1+1}{\bsp@args}{% +\ifshowframe% + \ncline[linestyle=dashed,linecolor=gray,arrows=*-*]{\bsp@root\ia}{\bsp@root\ib}% +\fi% + \nodexn{.667(\bsp@root\ia)+.333(\bsp@root\ib)}{\bsp@root R\ia}% + \nodexn{.333(\bsp@root\ia)+.667(\bsp@root\ib)}{\bsp@root L\ib}% + }% %Finally, define the bezier endpoints for each interval -\advance\bsp@args by -1 -\ifshowframe - \psset{linestyle=solid,linecolor=red} -\else - \psset{linestyle=none,linecolor=white} -\fi -%2008/9/20---change \ncline{*-*} to just \ncline -\mmultido{\ia=0+1}{\bsp@args}{\ncline{\bsp@root L\ia}{\bsp@root R\ia}\ncput[npos=.5]{\pnode{\bsp@root S\ia}}\ifshowframe \psdot[linecolor=red](\bsp@root S\ia)\fi -} -% -\ifbsp@closed - \pnode(\bsp@root S\the\bsp@args){\bsp@root S0} - \else - \advance\bsp@args by 1 \relax - \pnode(\bsp@root0){\bsp@root S0} - \pnode(\bsp@root\the\bsp@args){\bsp@root S\the\bsp@args} - \fi -\advance\bsp@args by -1 \relax +\advance\bsp@args by \m@ne% +\mmultido{\ia=0+1}{\bsp@args}{% +\ifshowframe% + \ncline[linestyle=solid,linecolor=red]{\bsp@root L\ia}{\bsp@root R\ia}% +\fi% +\midAB(\bsp@root L\ia)(\bsp@root R\ia){\bsp@root S\ia}% +\ifshowframe% + \psdot[linecolor=red](\bsp@root S\ia)% +\fi% +}% +\ifbsp@closed% + \pnode(\bsp@root S\the\bsp@args){\bsp@root S0}% + \else% + \advance\bsp@args by \@ne% + \pnode(\bsp@root0){\bsp@root S0}% + \pnode(\bsp@root\the\bsp@args){\bsp@root S\the\bsp@args}% + \fi% +\advance\bsp@args by \m@ne% %the pieces of bezier curve -\let\pslinecolor\bsp@linecolor -\let\pslinestyle\bsp@linestyle -\ifbsp@closed +\ifbsp@closed% \pscustom[arrows=-]{% -\psbezier[showpoints=false](\bsp@root S0)(\bsp@root R0)(\bsp@root L1)(\bsp@root S1) -\multido{\ia=1+1,\ib=2+1}{\bsp@args}{\psbezier(\bsp@root R\ia)(\bsp@root L\ib)(\bsp@root S\ib)}%within \pscustom, \psbezier takes 3 args, using last point of previous +\psbezier[showpoints=false](\bsp@root S0)(\bsp@root R0)(\bsp@root L1)(\bsp@root S1)% +\multido{\ia=1+1,\ib=2+1}{\bsp@args}{\psbezier(\bsp@root R\ia)(\bsp@root L\ib)% +(\bsp@root S\ib)}% +%within \pscustom, \psbezier takes 3 args, using last point of previous }% \else %open curve---make arrows, don't use \pscustom - \ifbsp@doubleline \psset{doubleline=true} \fi - \advance\bsp@args by -1 - \ifbsp@omitends \bsp@cntA=2 \bsp@cntB=1 \advance\bsp@args by -2 \else \bsp@cntA=1\bsp@cntB=0 \fi - \psbezier[arrows=-\bsp@arrowA,showpoints=flse](\bsp@root S\the\bsp@cntA)(\bsp@root L\the\bsp@cntA)(\bsp@root R\the\bsp@cntB)(\bsp@root S\the\bsp@cntB) - \mmultido{\ia=\the\bsp@cntA+1,\ib=\the\bsp@cntB+1}{\bsp@args}{\psbezier[arrows=-,showpoints=false](\bsp@root S\ib)(\bsp@root R\ib)(\bsp@root L\ia)(\bsp@root S\ia)} + \advance\bsp@args by \m@ne% + \ifbsp@omitends \bsp@cntA=2\relax \bsp@cntB=\@ne \advance\bsp@args by -2\relax% + \else \bsp@cntA=\@ne\bsp@cntB=\z@ \fi% + \psbezier[arrows=-\bsp@arrowA,showpoints=false](\bsp@root S\the\bsp@cntA)% + (\bsp@root L\the\bsp@cntA)(\bsp@root R\the\bsp@cntB)(\bsp@root S\the\bsp@cntB)% + \mmultido{\ia=\the\bsp@cntA+1,\ib=\the\bsp@cntB+1}{\bsp@args}% + {\psbezier[arrows=-,showpoints=false](\bsp@root S\ib)(\bsp@root R\ib)% + (\bsp@root L\ia)(\bsp@root S\ia)}% % - \ifbsp@omitends \advance\bsp@args by 1 \fi - \advance\bsp@args by 1 \relax - \bsp@cntA=\bsp@args - \bsp@cntB=\bsp@cntA - \advance\bsp@cntB by 1 \relax - \psbezier[arrows=-\bsp@arrowB,showpoints=false](\bsp@root S\the\bsp@cntA)(\bsp@root R\the\bsp@cntA)(\bsp@root L\the\bsp@cntB)(\bsp@root S\the\bsp@cntB) - \ifshowframe \psdot[linecolor=red](\bsp@root L\the\bsp@cntB) \fi -\fi % end \ifbsp@closed -\ifshowframe - \psdot[linecolor=red](\bsp@root R0) - \psdot[linecolor=red](\bsp@root L\the\bsp@args) -\fi -\ifshowpoints - \multido{\i=0+1}{\bsp@args}{\psdot(\bsp@root S\i)} -\fi - \ifbsp@closed + \ifbsp@omitends \advance\bsp@args by \@ne \fi% + \advance\bsp@args by \@ne% + \bsp@cntA=\bsp@args% + \bsp@cntB=\bsp@cntA% + \advance\bsp@cntB by \@ne% + \psbezier[arrows=-\bsp@arrowB,showpoints=false](\bsp@root S\the\bsp@cntA)% + (\bsp@root R\the\bsp@cntA)(\bsp@root L\the\bsp@cntB)(\bsp@root S\the\bsp@cntB)% + \ifshowframe \psdot[linecolor=red](\bsp@root L\the\bsp@cntB) \fi% +\fi% end \ifbsp@closed +\ifshowframe% + \psdot[linecolor=red](\bsp@root R0)% + \psdot[linecolor=red](\bsp@root L\the\bsp@args)% +\fi% +\ifshowpoints% + \multido{\i=0+1}{\bsp@args}{\psdot(\bsp@root S\i)}% +\fi% + \ifbsp@closed% \end@ClosedObj% - \else + \else% \end@OpenObj% - \fi -\ifbsp@showpoints \showpointstrue \fi -%\restoreparams -\ignorespaces + \fi% +\endgroup% +\ignorespaces% } % -\def\psbspline{\pst@object{psbspline}}% +\def\psbspline{\def\pst@par{}\pst@object{psbspline}}% \def\psbspline@i{\pst@getarrows\psbspline@ii}% \def\psbspline@ii{% \pnodesX%make the following sequences into nodes named S@0,... -%last index=\PST@args +%last index=\pst@args }% \def\psbspline@iii{% -\psBsplineInterp{S@}{\the\PST@args}%creates S@B0... +\pst@killglue% +\psBsplineInterp{S@}{\the\pst@args}%creates S@B0... \begingroup% \use@par% -\psBsplineNodes{S@B}{\the\PST@args}%draw curve -\endgroup}% +\psBsplineNodes{S@B}{\the\pst@args}%draw curve +\endgroup\ignorespaces}% % \def\pnodesX{%version of \pnodes where the last step moves to \psbspline@iii -\PST@args=0\relax% +\pst@args=\z@% \def\PST@root{S@}% \next@NodeX}% % -\def\next@NodeX{\@ifnextchar({\next@NodeX@i}{\advance\PST@args by -1\typeout{Defined nodes \PST@root0 .. \PST@root\the\PST@args.}\psbspline@iii}}% +\def\next@NodeX{\@ifnextchar({\next@NodeX@i}{\advance\pst@args by \m@ne\typeout{Defined nodes \PST@root0 .. \PST@root\the\pst@args.}\psbspline@iii}}% % \def\next@NodeX@i(#1){% -\pnode(#1){\PST@root\the\PST@args}% -\advance\PST@args by 1\relax% -\next@NodeX} +\pnode(#1){\PST@root\the\pst@args}% +\advance\pst@args\@ne% +\next@NodeX}% % \def\make@mcoeff#1{% -\pstVerb{ /mcoeff #1\space 1 add array def - mcoeff 1 0.25 put - 1 1 #1\space 1 sub +\pstVerb{ /mcoeff #1\space 1 add array def % + mcoeff 1 0.25 put % + 1 1 #1\space 1 sub % {/j exch def %j is counter - mcoeff j 1 add 1 4 mcoeff j get sub div put } for }} + mcoeff j 1 add 1 4 mcoeff j get sub div put } for }}% % - \def\psBsplineInterp#1#2{% - % #1=node name root to interpolate, #2=end index, - % eg, if we want to interpolate S0,..., S100, #1=S, #2=100, +\def\psBsplineInterp#1#2{% +% #1=node name root to interpolate, #2=end index, +% eg, if we want to interpolate S0,..., S100, #1=S, #2=100, \newcount\top@mone \top@mone=#2\relax%100 \newcount\top@ndx \top@ndx=\top@mone% \advance\top@ndx by -1\relax%99 \advance\top@mone by -2\relax%98 -\newcount\top@mtwo \top@mtwo=\top@mone \advance\top@mtwo by -1\relax%97 +\newcount\top@mtwo \top@mtwo=\top@mone \advance\top@mtwo by \m@ne%97 \make@mcoeff{#2}% \def\noderoot{#1}% \pnode(\noderoot 0){\noderoot B0}% \pnode(\noderoot #2){\noderoot B#2}% \psLCNode(\noderoot 1){6}(\noderoot 0){-1}{\noderoot B1}% 6*S1-S0->SB1 \multido{\iA=2+1,\iAA=1+1}{\top@mone}{\psLCNodeVar(\noderoot\iA)(\noderoot B\iAA)(! 6 mcoeff \iAA\space get neg ){\noderoot B\iA}}%2 to n-2, 6*S_k-m_{k-1}*SB_{k-1}->SB_k -\psLCNodeVar(\noderoot \the\top@ndx)(\noderoot #2)(! mcoeff \the\top@ndx\space get dup 6 mul exch neg ){tmp@node} %m_{99}*(6*S99-S100)->tmp@node +\psLCNodeVar(\noderoot \the\top@ndx)(\noderoot #2)(! mcoeff \the\top@ndx\space get dup 6 mul exch neg ){tmp@node}%m_{99}*(6*S99-S100)->tmp@node \psLCNodeVar(tmp@node)(\noderoot B\the\top@mone)(! 1 mcoeff \the\top@ndx\space get mcoeff \the\top@mone\space get mul neg ){\noderoot B\the\top@ndx}% temp@node -m_{98}*m_{99}*SB98->SB99 \multido{\iA=\the\top@mone+-1,\iAA=\the\top@ndx+-1}{\the\top@mone}{\psLCNodeVar(\noderoot B\iA)(\noderoot B\iAA)(! mcoeff \iA\space get dup neg ){\noderoot B\iA}% m_k*SB_k-m_k*SB_{k+1}->SB_k, k=98 downto 1 }% end multido }% end def - \def\psBsplineInterpC#1#2{% +\def\psBsplineInterpC#1#2{% % #1=node name root, #2=top index, % eg, if we want to interpolate S0,..., S100, #1=100, #2=S % interpolation in this case is periodic (closed curve) -\newcount\top@mone \top@mone=#2 % +\pst@killglue +\newcount\top@mone \top@mone=#2\relax% \newcount\top@ndx \top@ndx=\top@mone% -\advance\top@ndx by 1\relax%number of nodes to interpolate -\newcount\top@mtwo \top@mtwo=\top@mone \advance\top@mtwo by -1\relax% +\advance\top@ndx by \@ne%number of nodes to interpolate +\newcount\top@mtwo \top@mtwo=\top@mone \advance\top@mtwo\m@ne% % After copying S0 to S101, \top@ndx, \top@mone, \top@mtwo are top index, top index-1, top index-two---101, 100, 99 in the example \make@mcoeff{\the\top@ndx}% \def\noderoot{#1}% \pstVerb{ /bcoeff #2\space 2 add array def /ccoeff #2\space 2 add array def /xcoeff #2\space 2 add array def /ycoeff #2\space 2 add array def -% initialize b, c, x ,y values +%% initialize b, c, x ,y values 1 1 \the\top@ndx\space %1 add {/j exch def %j is counter ccoeff j 0 put bcoeff j 0 put } for @@ -283,7 +234,7 @@ bcoeff 1 1 put bcoeff #2\space 1 put bcoeff #2\space 1 add 4 put ccoeff 1 0.25 put ccoeff #2\space 1 put }% %end \pstVerb \pnode(\noderoot 0){\noderoot\the\top@ndx}%copy node 0 to top end -%\pstVerb{ xcoeff 1 1 put ycoeff 1 1.5 put } +\pstVerb{ xcoeff 1 1 put ycoeff 1 1.5 put } \pnode(! \psGetNodeCenter{\noderoot 1} xcoeff 1 \noderoot 1.x 1.5 mul put ycoeff 1 \noderoot 1.y 1.5 mul put 1 1){temp@}% 1.5=6*0.25 \multido{\iA=2+1}{\the\top@mone}{\pnode(! ycoeff \iA\space xcoeff \iA\space \psGetNodeCenter{#1\iA} #1\iA.x 6 mul put #1\iA.y 6 mul put 1 1){temp@}}% move x, y coords times 6 to arrays \multido{\iA=2+1,\iB=1+1}{\the\top@mtwo}{\pstVerb{% @@ -308,8 +259,9 @@ ycoeff \the\top@mone\space ycoeff \the\top@mone\space get ycoeff \the\top@ndx\sp \pstVerb{ xcoeff \iA\space xcoeff \iA\space get mcoeff \iA\space get xcoeff \iB\space get mul sub ccoeff \iA\space get xcoeff \the\top@ndx\space get mul sub put % x(k)=x(k)-m(k)*x(k+1)-c(k)*x(n) ycoeff \iA\space ycoeff \iA\space get mcoeff \iA\space get ycoeff \iB\space get mul sub ccoeff \iA\space get ycoeff \the\top@ndx\space get mul sub put }}%y(k)=... -\multido{\iA=1+1}{\the\top@ndx}{\pnode(! xcoeff \iA\space get ycoeff \iA\space get ){\noderoot B\iA}}% +\multido{\iA=1+1}{\the\top@ndx}{\pnode(! xcoeff \iA\space get ycoeff \iA\space get )% +{\noderoot B\iA}}% \pnode(\noderoot B\the\top@ndx){\noderoot B0}% -} +\ignorespaces}% \catcode`\@=\TheAtCode\relax \endinput
\ No newline at end of file |