diff options
Diffstat (limited to 'modules/output_stream/macros')
| -rwxr-xr-x | modules/output_stream/macros/buildmacros.bat | 1 | ||||
| -rwxr-xr-x | modules/output_stream/macros/buildmacros.sce | 9 | ||||
| -rwxr-xr-x | modules/output_stream/macros/cleanmacros.bat | 3 | ||||
| -rwxr-xr-x | modules/output_stream/macros/lib | bin | 0 -> 336 bytes | |||
| -rwxr-xr-x | modules/output_stream/macros/names | 2 | ||||
| -rwxr-xr-x | modules/output_stream/macros/prettyprint.bin | bin | 0 -> 78016 bytes | |||
| -rwxr-xr-x | modules/output_stream/macros/prettyprint.sci | 596 | ||||
| -rwxr-xr-x | modules/output_stream/macros/ssprint.bin | bin | 0 -> 23164 bytes | |||
| -rwxr-xr-x | modules/output_stream/macros/ssprint.sci | 184 |
9 files changed, 795 insertions, 0 deletions
diff --git a/modules/output_stream/macros/buildmacros.bat b/modules/output_stream/macros/buildmacros.bat new file mode 100755 index 000000000..c4e35ec40 --- /dev/null +++ b/modules/output_stream/macros/buildmacros.bat @@ -0,0 +1 @@ +@..\..\..\bin\scilex -nwni -ns -e exec('buildmacros.sce');quit;
\ No newline at end of file diff --git a/modules/output_stream/macros/buildmacros.sce b/modules/output_stream/macros/buildmacros.sce new file mode 100755 index 000000000..a6aad99f4 --- /dev/null +++ b/modules/output_stream/macros/buildmacros.sce @@ -0,0 +1,9 @@ +//------------------------------------ +// Allan CORNET INRIA 2007 +//------------------------------------ +if (isdef("genlib") == %f) then + exec(SCI+"/modules/functions/scripts/buildmacros/loadgenlib.sce"); +end +//------------------------------------ +genlib("output_streamlib","SCI/modules/output_stream/macros",%f,%t); +//------------------------------------ diff --git a/modules/output_stream/macros/cleanmacros.bat b/modules/output_stream/macros/cleanmacros.bat new file mode 100755 index 000000000..5079dfd71 --- /dev/null +++ b/modules/output_stream/macros/cleanmacros.bat @@ -0,0 +1,3 @@ +@del *.bin 2>NUL +@del lib 2>NUL +@del names 2>NUL
\ No newline at end of file diff --git a/modules/output_stream/macros/lib b/modules/output_stream/macros/lib Binary files differnew file mode 100755 index 000000000..85699a3dd --- /dev/null +++ b/modules/output_stream/macros/lib diff --git a/modules/output_stream/macros/names b/modules/output_stream/macros/names new file mode 100755 index 000000000..d875460fc --- /dev/null +++ b/modules/output_stream/macros/names @@ -0,0 +1,2 @@ +prettyprint +ssprint diff --git a/modules/output_stream/macros/prettyprint.bin b/modules/output_stream/macros/prettyprint.bin Binary files differnew file mode 100755 index 000000000..f20dbe6d9 --- /dev/null +++ b/modules/output_stream/macros/prettyprint.bin diff --git a/modules/output_stream/macros/prettyprint.sci b/modules/output_stream/macros/prettyprint.sci new file mode 100755 index 000000000..d35009d2b --- /dev/null +++ b/modules/output_stream/macros/prettyprint.sci @@ -0,0 +1,596 @@ +// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab +// Copyright (C) 2009-2010 - Calixte Denizet +// +// This file must be used under the terms of the CeCILL. +// This source file is licensed as described in the file COPYING, which +// you should have received as part of this distribution. The terms +// are also available at +// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt + +function str = prettyprint(a, exportFormat, delimiter, processByElement, isWrapped) + // From any Scilab datatype and provide a representation to the TeX, LaTeX or MathML formats + // + // Calling Sequence + // str = prettyprint(a) // Show the variable a with the default format (LaTeX) + // str = prettyprint(a,exportFormat) // Show the variable a with the specified format + // str = prettyprint(a,exportFormat, delim) // As above but change the delimiter + // str = prettyprint(a,exportFormat, delim, processByElement) // As above but process each element independently + // str = prettyprint(a,exportFormat, delim, processByElement, isWrapped) // As above Add the special keyword of Tex or LaTeX export + // + // Parameters + // a: is a Scilab variable + //exportFormat: is the format, if omitted 'latex' is used by default, it can be 'latex', 'tex' or 'mathml'. + // delimiter: is a string indicating the delimiter to use for the resulting matrix, it's only used if isWrapped is true. The delimiter can be '(', '{', '[', '|', '||' or ')' + //processByElement: is a boolean to indicate if the resulting matrix must be converted into a single string. + //isWrapped: is a boolean to indicate if the result must be wrapped inside delimiters ('$' for latex and tex or nothing for mathml) to be used with xstring or xtitle + //str: the representation of the variable a + //If the type of a is not handled, then an error is returned. If a user wants to handle the type foo when exporting with latex, he must define the function foo2latex. + // + // Description + // Taking a variable, the prettyprint function will provide a formated representation of it. + // Formats can be TeX, LaTeX or MathML. + // They can be used in third party applications but also within Scilab with the most of the + // <link linkend="math_rendering_features_in_graphic">Scilab graphic features</link>. + // The following types are handled by this function: + // <itemizedlist> + // <listitem><para>Real / Complex matrices</para></listitem> + // <listitem><para>Polynomial types</para></listitem> + // <listitem><para>Boolean</para></listitem> + // <listitem><para>Integer</para></listitem> + // <listitem><para>String</para></listitem> + // <listitem><para>Tlist</para></listitem> + // <listitem><para>Rationnal</para></listitem> + // <listitem><para>Cell</para></listitem> + // </itemizedlist> + // + // Examples + // str = prettyprint(rand(3,3)) // Return the LaTeX representation of a 3,3 matrix + // xstring(0.2,0.2,str) // Show the representation in a graphic Windows + // + // prettyprint(rand(3,4),"mathml") // Return the MathML representation of a 3,4 matrix + // prettyprint(rand(3,4),"mathml","[") // Return the MathML representation of a 3,4 matrix with '[' as delimiter + // + // s=poly(0,'s'); G=[1,s;1+s^2,3*s^3]; + // xstring(0.2,0.2,prettyprint(G*s-1)); // Show a polynomial through a LaTeX representation + // + // See also + // math_rendering_features_in_graphic + // xtitle + // axes_properties + // label_properties + // legend_properties + // text_properties + // xstringb + // xstringl + // xstring + // + // Authors + // Calixte Denizet + + nargs = argn(2); + + select nargs + case 0 then + error(msprintf(gettext("%s: Wrong number of input argument(s): %d to %d expected."),"prettyprint",1,4)); + case 1 then + exportFormat="latex"; + delimiter = "("; + processByElement = %F; + isWrapped = %T; + case 2 then + if type(exportFormat) <> 10 then + error(msprintf(gettext("%s: Wrong type for input argument #%d: String expected.\n"),"prettyprint",2)); + end + delimiter = "("; + processByElement = %F; + isWrapped = %T; + case 3 then + if type(exportFormat) <> 10 then + error(msprintf(gettext("%s: Wrong type for input argument #%d: String expected.\n"),"prettyprint",2)); + end + if type(delimiter) <> 10 | (delimiter <> "(" & delimiter <> "[" & delimiter <> "|" & delimiter <> "||" & delimiter <> "{" & delimiter <> "") then + error(msprintf(gettext("%s: Wrong value for input argument #%d: ''%s'', ''%s'', ''%s'', ''%s'', ''%s'' or ''%s'' expected.\n"),"prettyprint",3,"(","[","|","||","{","")); + end + processByElement = %F; + isWrapped = %T; + case 4 then + if type(exportFormat) <> 10 then + error(msprintf(gettext("%s: Wrong type for input argument #%d: String expected.\n"),"prettyprint",2)); + end + if type(delimiter) <> 10 | (delimiter <> "(" & delimiter <> "[" & delimiter <> "|" & delimiter <> "||" & delimiter <> "{" & delimiter <> "") then + error(msprintf(gettext("%s: Wrong value for input argument #%d: ''%s'', ''%s'', ''%s'', ''%s'', ''%s'' or ''%s'' expected.\n"),"prettyprint",3,"(","[","|","||","{","")); + end + if type(processByElement) <> 4 then + error(msprintf(gettext("%s: Wrong type for argument #%d: Boolean expected.\n"),"prettyprint",4)); + end + isWrapped = %T; + case 5 + if type(exportFormat) <> 10 then + error(msprintf(gettext("%s: Wrong type for input argument #%d: String expected.\n"),"prettyprint",2)); + end + if type(delimiter) <> 10 | (delimiter <> "(" & delimiter <> "[" & delimiter <> "|" & delimiter <> "||" & delimiter <> "{" & delimiter <> "") then + error(msprintf(gettext("%s: Wrong value for input argument #%d: ''%s'', ''%s'', ''%s'', ''%s'', ''%s'' or ''%s'' expected.\n"),"prettyprint",3,"(","[","|","||","{","")); + end + if type(processByElement) <> 4 then + error(msprintf(gettext("%s: Wrong type for argument #%d: Boolean expected.\n"),"prettyprint",4)); + end + if type(isWrapped) <> 4 then + error(msprintf(gettext("%s: Wrong type for argument #%d: Boolean expected.\n"),"prettyprint",5)); + end + else + error(msprintf(gettext("%s: Wrong number of input arguments: %d to %d expected"),"prettyprint",1,4)); + end + try + execstr("[plus,minus,img,op,cp,ow,cw,d2s]=" + exportFormat + "conf()"); + catch + error(msprintf(gettext("%s: Wrong export format: %s"),"prettyprint",exportFormat)); + end + + typ = type(a); + + select typ + //Scalar type : real or complex + case 1 then + if or(isinf(a) | isnan(a)) | norm(imag(a)) > %eps * norm(real(a)) then + str = comp2str(a,plus,minus,img,d2s); + else + str = comp2str(real(a),plus,minus,img,d2s); + end + //Polynomial type + case 2 then + x = varn(a); + C = coeff(a); + [m,n] = size(a); + maxd = max(degree(a)); + str = emptystr(a); + for i = 0:maxd do + execstr("expo=" + exportFormat + "exp(x,i)"); + A = C(1:m,(1 + i * n):((i + 1) * n)); + if or(isinf(a) | isnan(a)) | norm(imag(A)) > %eps * norm(real(A)) then + str = str + comp2coef(A,expo,plus,minus,img,op,cp,d2s); + else + str = str + comp2coef(real(A),expo,plus,minus,img,op,cp,d2s); + end + end + L = length(plus); + for i = 1:m do + for j = 1:n do + if part(str(i,j),1:L) == plus then + str(i,j) = part(str(i,j),L + 1:length(str(i,j))); + end + end + end + //The null polynomial is represented by '' so we must replace it by '0' + str(str == "") = d2s(0); + //Boolean type + case 4 then + str = "F" + emptystr(a); + str(a) = "T"; + execstr("str=" + exportFormat + "exp(str,1)"); + //Int type + case 8 then + str = d2s(a); + //String type + case 10 then + str = a; + //Tlist type + case 16 then + t = a(1); + select t(1) + //Rationnal type + case "r" then + num = prettyprint(a("num"),exportFormat,"(",%T,%F); + den = prettyprint(a("den"),exportFormat,"(",%T,%F); + execstr("str=rational2" + exportFormat + "(num,den)"); + //Linear state space type + case "lss" then + execstr("str=lss2" + exportFormat + "(a)"); + else + str = unknown_type(t(1),a,exportFormat); + return; + end + case 17 then + a1=getfield(1,a); + select a1(1) + //Cell type + case "ce" then + dim = double(a.dims); + L = length(dim); + if L >= 3 then + str = unknown_type("ce",a,exportFormat); + return; + end + str = emptystr(dim(1),dim(2)); + for i = 1:dim(1) do + for j = 1:dim(2) do + str(i,j) = prettyprint(a(i,j).entries,exportFormat,delimiter,%F,%F); + end + end + else + str = unknown_type(a1(1),a,exportFormat); + return; + end + else + str = unknown_type(typeof(a),a,exportFormat); + return; + end + + [m,n] = size(a); + if m*n <> 1 & ~processByElement then + execstr("str=" + exportFormat + "matrix(str,''" + delimiter + "'')"); + end + + if isWrapped then + str = ow + str + cw; + end +endfunction + +//This function gives different elements of configuration +//plus is '+', minus is '-', img is sqrt(-1), op and cp are opening ang closing +//parenthesis (useful to generate complex coeff in polynomials), ow and cw are +//opening and closing delimiters used to wrap the expression to be displayed +//with xstring or xtitle, d2s is the function used to convert double into mathml. +function [plus,minus,img,op,cp,ow,cw,d2s] = mathmlconf() + plus = "<mo>+</mo>"; + minus = "<mo>-</mo>"; + img = "<mi>i</mi>"; + op = "<mfenced separator=""""><mrow>"; + cp = "</mrow></mfenced>"; + ow = ""; + cw = ""; + d2s = mathmldbl2str; +endfunction + +//This function generates var^n (useful for polynomials) +function str = mathmlexp(var,n) + if n >= 2 then + str = "<msup><mi>" + var + "</mi><mn>" + string(n) + "</mn></msup>"; + elseif n == 1 then + str = "<mi>" + var + "</mi>"; + else + str = ""; + end +endfunction + +//This function creates a fraction with the given numerator and denominator +function str = rational2mathml(num,den) + str = "<mfrac><mrow>" + num + "</mrow><mrow>" + den + "</mrow></mfrac>"; +endfunction + +//This function converts a double into a mathml string +//Negative and positive are handled and the scientific notation too +function str = mathmldbl2str(x) + if x == %inf then + str = "<infinity/>"; + elseif x == -%inf then + str = "<mo>-</mo><infinity/>"; + elseif isnan(x) then + str = "<mi>NaN</mi>"; + else + str= strsubst(string(x),"-","<mo>-</mo><mn>"); + if strindex(str,"D") <> [] then + str = strsubst(str,"+","<mo>+</mo><mn>"); + str = strsubst(str,"D","</mn><mo>×</mo><msup><mn>10</mn><mrow>"); + if x > 0 then + str = "<mn>" + str; + end + str = str + "</mn></mrow></msup>"; + elseif x >= 0 then + str = "<mn>" + str + "</mn>"; + else + str = str + "</mn>"; + end + end +endfunction + +//This function generates a matrix with the given delimiter +function str = mathmlmatrix(mat,delimiter) + if argn(2) == 1 then + delimiter = "("; + end + select delimiter + case "[" then + com = "<mfenced open=""["" close=""]"">"; + case "(" then + com = "<mfenced open=""("" close="")"">"; + case "|" then + com = "<mfenced open=""∣"" close=""∣"">"; + case "" then + com = "<mfenced>"; + case "{" then + com = "<mfenced open=""{"" close=""}"">"; + case "||" then + com = "<mfenced open=""∥"" close=""∥"">"; + else + com = "<mfenced open=""("" close="")"">"; + end + [m,n] = size(mat); + str = com + "<mtable>"; + for i=1:m do + str = str + "<mtr><mtd>" + strcat(mat(i,:),"</mtd><mtd>") + "</mtd></mtr>"; + end + str = str + "</mtable></mfenced>"; +endfunction + +//This function handles the syslin +function str = lss2mathml(sys) + ptp = "<mo>(</mo><mi>t</mi><mo>)</mo>"; + if sys(7) == "c" then + der = "<mrow><mover><mi>X</mi><mo>.</mo></mover>" + ptp + "</mrow>"; + else + der = "<mrow><mover><mi>X</mi><mo>+</mo></mover>" + ptp + "</mrow>" + end + str = "<mfenced close="""" open=""{""><mtable columnalign=""right center left""><mtr><mtd>" + der + "</mtd><mtd><mo>=</mo></mtd><mtd>" + prettyprint(sys(2),"mathml","(",%F,%F) + "<mrow><mi>X</mi>" + ptp + "</mrow><mo>+</mo>" + prettyprint(sys(3),"mathml","(",%F,%F) + "<mrow><mi>U</mi>" + ptp + "</mrow></mtd></mtr><mtr><mtd><mrow><mi>Y</mi>" + ptp + "</mrow></mtd><mtd><mo>=</mo></mtd><mtd>" + prettyprint(sys(4),"mathml","(",%F,%F) + "<mrow><mi>X</mi>" + ptp + "</mrow>"; + try + if norm(sys(5),1) == 0 then + str = str + "</mtd></mtr></mtable></mfenced>"; + else + str = str + "<mo>+</mo>" + prettyprint(sys(5),"mathml","(",%F,%F) + "<mrow><mi>U</mi>" + ptp + "</mrow></mtd></mtr></mtable></mfenced>"; + end + catch + str = str + "<mo>+</mo>" + prettyprint(sys(5),"mathml","(",%F,%F) + "<mrow><mi>U</mi>" + ptp + "</mrow></mtd></mtr></mtable></mfenced>"; + end +endfunction + +function [plus,minus,img,op,cp,ow,cw,d2s] = latexconf() + plus = "+"; + minus = "-"; + img = "i"; + op = "("; + cp = ")"; + ow = "$"; + cw = "$"; + d2s = latexdbl2str; +endfunction + +function str = latexexp(var,n) + if n >= 2 then + str = var + "^{" + string(n) + "}"; + elseif n == 1 then + str = var; + else + str = ""; + end +endfunction + +function str = rational2latex(num,den) + str = "{\frac{" + num + "}{" + den + "}}" +endfunction + +function str = latexdbl2str(x) + if x == %inf then + str = "{\infty}"; + elseif x == -%inf then + str = "{-\infty}"; + elseif isnan(x) then + str = "{\mathrm{NaN}}"; + else + str = string(x); + if strindex(str,"D") <> [] then + str = strsubst(str,"D","\times10^{"); + str = str + "}"; + end + end +endfunction + +function str = lss2latex(sys) + if sys(7) == "c" then + der = "\dot{X}(t)"; + else + der = "\stackrel{+}{X}(t)" + end + str = "{\left\{\begin{array}{rcl}" + der + "&=&" + prettyprint(sys(2),"latex","(",%F,%F) + " X(t)+" + prettyprint(sys(3),"latex","(",%F,%F) + "U(t)\cr Y(t) &=& " + prettyprint(sys(4),"latex","(",%F,%F) + " X(t) "; + try + if norm(sys(5),1) == 0 then + str = str + "\end{array}\right.}"; + else + str = str + " + " + prettyprint(sys(5),"latex","(",%F,%F) + " U(t)\end{array}\right.}"; + end + catch + str = str + " + " + prettyprint(sys(5),"latex","(",%F,%F) + " U(t)\end{array}\right.}"; + end +endfunction + +function str = latexmatrix(mat,delimiter) + if argn(2) == 1 then + delimiter = "("; + end + select delimiter + case "[" then + com = "bmatrix"; + case "(" then + com = "pmatrix"; + case "|" then + com = "vmatrix"; + case "" then + com = "matrix"; + case "{" then + com = "Bmatrix"; + case "||" then + com = "Vmatrix"; + else + com = "pmatrix"; + end + [m,n] = size(mat); + str = "{\begin{" + com + "}"; + if m*n == 0 then + str = str + "\ "; + else + for i=1:m do + str = str + strcat(mat(i,:),"&") + "\cr "; + end + end + str = str + "\end{" + com + "}}"; +endfunction + +function [plus,minus,img,op,cp,ow,cw,d2s] = texconf() + [plus,minus,img,op,cp,ow,cw,d2s] = latexconf() +endfunction + +function str = texexp(var,n) + str = latexexp(var,n); +endfunction + +function str = rational2tex(num,den) + str = "{{" + num + "}\over{" + den + "}}" +endfunction + +function str = texdbl2str(x) + if x == %inf then + str = "{\infty}"; + elseif x == -%inf then + str = "{-\infty}"; + elseif isnan(x) then + str = "{\rm{NaN}}"; + else + str = string(x); + if strindex(str,"D") <> [] then + str = strsubst(str,"D","\times10^{"); + str = str + "}"; + end + end +endfunction + +function str = lss2tex(sys) + if sys(7) == "c" then + der = "\dot{X}(t)"; + else + der = "{\buildrel + \over X(t)}" + end + str = "{\left\{\eqalign{" + der + " &= " + prettyprint(sys(2),"tex","(",%F,%F) + " X(t)+" + prettyprint(sys(3),"tex","(",%F,%F) + "U(t)\cr Y(t) &= " + prettyprint(sys(4),"tex","(",%F,%F) + " X(t) "; + try + if norm(sys(5),1) == 0 then + str = str + "}\right.}"; + else + str = str + " + " + prettyprint(sys(5),"tex","(",%F,%F) + " U(t)}\right.}"; + end + catch + str = str + " + " + prettyprint(sys(5),"tex","(",%F,%F) + " U(t)}\right.}"; + end +endfunction + +function str = texmatrix(mat,delimiter) + if argn(2) == 1 then + delimiter = "("; + end + select delimiter + case "[" then + com = "bmatrix"; + case "(" then + com = "pmatrix"; + case "|" then + com = "vmatrix"; + case "" then + com = "matrix"; + case "{" then + com = "Bmatrix"; + case "||" then + com = "Vmatrix"; + else + com = "pmatrix"; + end + [m,n] = size(mat); + str = "{\" + com + "{"; + if m*n == 0 then + str = str + "\ "; + else + for i=1:m do + str = str + strcat(mat(i,:),"&") + "\cr "; + end + end + str = str + "}}"; +endfunction + +//This function converts a complex (or a double) into a string in using the +//function d2s (double2string) fixed by mathmlconf or latexconf. +function str = comp2str(z,plus,minus,img,d2s) + [m,n] = size(z); + if m*n <> 1 then + str = emptystr(z); + for i=1:m do + for j=1:n do + str(i,j) = comp2str(z(i,j),plus,minus,img,d2s); + end + end + return; + end + re = real(z); + im = imag(z); + if re == 0 & im == 0 then + str = d2s(0); + return; + end; + if re == 0 then + if im == 1 then + str = img; + elseif im == -1 then + str = minus + img; + else + str = d2s(im) + img; + end + elseif im == 0 then + if re > 0 then + str = d2s(re); + else + str = d2s(re); + end + elseif im > 0 then + if im == 1 then + str = d2s(re) + plus + img; + else + str = d2s(re) + plus + d2s(im) + img; + end + else + if im == -1 then + str = d2s(re) + minus + img; + else + str = d2s(re) + d2s(im) + img; + end + end +endfunction + +//This function convert a complex into a coefficient (parenthesis are eventually +// added around the coeff, if var=='' (i.e. var^0), no parenthesis are added. +function str = comp2coef(z,var,plus,minus,img,op,cp,d2s) + [m,n] = size(z); + if m*n <> 1 then + str = emptystr(z); + for i=1:m do + for j=1:n do + str(i,j) = comp2coef(z(i,j),var,plus,minus,img,op,cp,d2s); + end + end + return; + end + re = real(z); + im = imag(z); + if var <> "" then + if re <> 0 & im <> 0 then + str = plus + op + comp2str(z,plus,minus,img,d2s) + cp + var; + elseif re == 0 & im == 0 then + str = ""; + elseif re == 1 & im == 0 then + str = plus + var; + elseif re == -1 & im == 0 then + str = minus + var; + elseif re > 0 | im > 0 then + str = plus + comp2str(z,plus,minus,img,d2s) + var; + else + str = comp2str(z,plus,minus,img,d2s) + var; + end + else + if re == 0 & im == 0 then + str = ""; + else + str = comp2str(z,plus,minus,img,d2s); + end + end +endfunction + +function str = unknown_type(typ,a,exportFormat) + try + execstr("str=" + typ + "2" + exportFormat + "(a)"); + catch + error(msprintf(gettext("%s: Type %s is not handled : Define the function %s2%s."),"prettyprint",typ,typ,exportFormat)) + end +endfunction
\ No newline at end of file diff --git a/modules/output_stream/macros/ssprint.bin b/modules/output_stream/macros/ssprint.bin Binary files differnew file mode 100755 index 000000000..eb8112620 --- /dev/null +++ b/modules/output_stream/macros/ssprint.bin diff --git a/modules/output_stream/macros/ssprint.sci b/modules/output_stream/macros/ssprint.sci new file mode 100755 index 000000000..93d0385f2 --- /dev/null +++ b/modules/output_stream/macros/ssprint.sci @@ -0,0 +1,184 @@ +// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab +// Copyright (C) XXXX-2008 INRIA +// Copyright (C) 2012 - DIGITE0 - Allan CORNET +// +// This file must be used under the terms of the CeCILL. +// This source file is licensed as described in the file COPYING, which +// you should have received as part of this distribution. The terms +// are also available at +// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt + +function [] = ssprint(sl, out) + + //ssprint(sl [,out]) pretty print of a linear system in state-space form + // sl = (A, B, C, D) syslin list. + // out = output (default out=%io(2)) + // See also texprint. + //Examples: + // a = [1 1;0 1]; b = [0 1;1 0]; c = [1, 1]; d = [3, 2]; + // ssprint(syslin('c', a, b, c, d)) + // ssprint(syslin('d', a, b, c, d)) + //! + + [lhs, rhs] = argn(0); + if rhs == 0 then + error(msprintf(gettext("%s: Wrong number of input argument(s).\n"), "ssprint")); + end + fil = %f; + if rhs == 1 then + out = %io(2); + else + if type(out) == 10 then + if size(out, "*") <> 1 then + error(msprintf(gettext("%s: Wrong size for input argument #%d: A character string expected.\n"), "ssprint", 2)); + end + out = file("open", out, "unknown"); + fil = %t + else + error(msprintf(gettext("%s: Wrong type for input argument #%d: A character string expected.\n"), "ssprint", 2)); + end + end + + if typeof(sl) <> "state-space" then + error(msprintf(gettext("%s: Wrong type for input argument #%d: A state-space model expected.\n"), "ssprint", 1)); + end + + function [ta] = %cv(x) + [m, n] = size(x); + if m * n == 0 then ta = " ", return, end + frmt = format(); + frmt = 10**frmt(2)/max([1, norm(coeff(x))]); + x = round(frmt * x) / frmt; + t = []; + for k = 1:m, t = [t; "|"], end; + ta = t; + for k = 1:n + aa = string(x(:, k)); + for l = 1:m + if part(aa(l), 1) <> "-" then + aa(l) = " " + aa(l); + end + end + n = max(length(aa)); + aa = part(aa + blank, 1:n); + ta = ta + aa + part(blank, 1); + end; + ta = ta + t; + endfunction + + // + // d(x)=ax + bu + //------------- + write(out, " ") + sgn = "."; + if sl(7) <> "c" then sgn = "+",end + ll = lines(); + ll = ll(1); + [a, b, c, d] = sl(2:5); + na = size(a, "r"); + [nc, nb] = size(d); + if na > 0 then + blank = []; + for k = 1:na, blank = [blank; " "],end + ta = %cv(a); + tb = %cv(b); + // + blank = part(blank, 1:4); + if na == 1 then + t = [sgn + " "; "x = "] + [" "; ta] + [" "; "x"] + else + blank([na/2, na/2+1]) = [sgn + " "; "x = "] + t = blank + ta; + end + if nb > 0 then + if na == 1 then + t = t + [" ";" + "]+[" ";tb]+[" ";"u "] + else + blank([na/2, na/2+1]) = [" "; "x + "] + t = t + blank + tb + t(na/2+1) = t(na/2 + 1) + "u " + end + end + n1 = max(length(t)) + 1 + // + n2 = max(length(t)) + if n2 < ll then + write(out, t) + else + if n1 < ll, + if n2 - n1 < ll then + write(out, part(t, 1:n1)), + write(out, " ") + write(out, part(t, n1+1:n2), "(3x,a)") + else + error(msprintf(gettext("%s: system cannot be displayed in this page.\n"), "ssprint")); + end; + else + error(msprintf(gettext("%s: system cannot be displayed in this page.\n"), "ssprint")); + end; + end; + end + // + //y = cx + du + //----------- + if nc == 0 then if fil then file("close",out);end;return,end + write(out, " "); + blank=[]; + for k = 1:nc, blank = [blank; " "], end + if na == 0 then + td = %cv(d) + if nc == 1 then + t = "y = " + td + "u" + else + blank(nc/2+1)="y = " + t = blank + td + t(nc/2 + 1) = t(nc/2+1) + "u " + end + n1 = max(length(t))+1 + else + tc = %cv(c); + td = %cv(d); + blank = part(blank, 1:4); + if nc == 1 then + t = "y = " + tc; + else + blank(nc/2 + 1) = "y = "; + t = blank + tc; + end + withd = d <> []; + if withd then withd = norm(d, 1) > 0, end + if withd then + if nc == 1 then + t = t + "x + " + td; + else + blank(nc/2 + 1) = "x + "; + t = t + blank + td; + t(nc/2 + 1) = t(nc/2 + 1) + "u "; + end + n1 = max(length(t)) + 1; + else + if nc == 1 then + t = t + "x "; + else + t(nc/2 + 1) = t(nc/2 + 1) + "x "; + end + end; + end + n2 = max(length(t)); + if n2 < ll then + write(out, t); + else + if n1 < ll, + if n2 - n1 < ll then + write(out, part(t, 1:n1)), + write(out, " "); + write(out, part(t, n1+1:n2), "(3x,a)") + else + error(msprintf(gettext("%s: system cannot be displayed in this page.\n"), "ssprint")); + end + else + error(msprintf(gettext("%s: system cannot be displayed in this page.\n"), "ssprint")); + end + end + if fil then file("close", out),end +endfunction |