doc-src/isac/jrocnik/eJMT-paper/jrocnik_eJMT.tex
author Walther Neuper <neuper@ist.tugraz.at>
Fri, 07 Sep 2012 18:12:50 +0200
changeset 42468 5f8f02e1ea9f
parent 42467 1035c36360ae
child 42469 264803a0c13e
permissions -rwxr-xr-x
jrocnik: example calculation

polished in \sect 4 in accordance to prog. in \sect 3
     1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     2 % Electronic Journal of Mathematics and Technology (eJMT) %
     3 % style sheet for LaTeX.  Please do not modify sections   %
     4 % or commands marked 'eJMT'.                              %
     5 %                                                         %
     6 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     7 %                                                         %
     8 % eJMT commands                                           %
     9 %                                                         %
    10 \documentclass[12pt,a4paper]{article}%                    %
    11 \usepackage{times}                                        %
    12 \usepackage{amsfonts,amsmath,amssymb}                     %
    13 \usepackage[a4paper]{geometry}                            %
    14 \usepackage{fancyhdr}                                     %
    15 \usepackage{color}                                        %
    16 \usepackage[pdftex]{hyperref} % see note below            %
    17 \usepackage{graphicx}%                                    %
    18 \hypersetup{                                              %
    19     a4paper,                                              %
    20     breaklinks                                            %
    21 }                                                         %
    22 %                                                         %
    23 \newtheorem{theorem}{Theorem}                             %
    24 \newtheorem{acknowledgement}[theorem]{Acknowledgement}    %
    25 \newtheorem{algorithm}[theorem]{Algorithm}                %
    26 \newtheorem{axiom}[theorem]{Axiom}                        %
    27 \newtheorem{case}[theorem]{Case}                          %
    28 \newtheorem{claim}[theorem]{Claim}                        %
    29 \newtheorem{conclusion}[theorem]{Conclusion}              %
    30 \newtheorem{condition}[theorem]{Condition}                %
    31 \newtheorem{conjecture}[theorem]{Conjecture}              %
    32 \newtheorem{corollary}[theorem]{Corollary}                %
    33 \newtheorem{criterion}[theorem]{Criterion}                %
    34 \newtheorem{definition}[theorem]{Definition}              %
    35 \newtheorem{example}[theorem]{Example}                    %
    36 \newtheorem{exercise}[theorem]{Exercise}                  %
    37 \newtheorem{lemma}[theorem]{Lemma}                        %
    38 \newtheorem{notation}[theorem]{Notation}                  %
    39 \newtheorem{problem}[theorem]{Problem}                    %
    40 \newtheorem{proposition}[theorem]{Proposition}            %
    41 \newtheorem{remark}[theorem]{Remark}                      %
    42 \newtheorem{solution}[theorem]{Solution}                  %
    43 \newtheorem{summary}[theorem]{Summary}                    %
    44 \newenvironment{proof}[1][Proof]{\noindent\textbf{#1.} }  %
    45 {\ \rule{0.5em}{0.5em}}                                   %
    46 %                                                         %
    47 % eJMT page dimensions                                    %
    48 %                                                         %
    49 \geometry{left=2cm,right=2cm,top=3.2cm,bottom=4cm}        %
    50 %                                                         %
    51 % eJMT header & footer                                    %
    52 %                                                         %
    53 \newcounter{ejmtFirstpage}                                %
    54 \setcounter{ejmtFirstpage}{1}                             %
    55 \pagestyle{empty}                                         %
    56 \setlength{\headheight}{14pt}                             %
    57 \geometry{left=2cm,right=2cm,top=3.2cm,bottom=4cm}        %
    58 \pagestyle{fancyplain}                                    %
    59 \fancyhf{}                                                %
    60 \fancyhead[c]{\small The Electronic Journal of Mathematics%
    61 \ and Technology, Volume 1, Number 1, ISSN 1933-2823}     %
    62 \cfoot{%                                                  %
    63   \ifnum\value{ejmtFirstpage}=0%                          %
    64     {\vtop to\hsize{\hrule\vskip .2cm\thepage}}%          %
    65   \else\setcounter{ejmtFirstpage}{0}\fi%                  %
    66 }                                                         %
    67 %                                                         %
    68 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    69 %
    70 % Please place your own definitions here
    71 %
    72 \def\isac{${\cal I}\mkern-2mu{\cal S}\mkern-5mu{\cal AC}$}
    73 \def\sisac{\footnotesize${\cal I}\mkern-2mu{\cal S}\mkern-5mu{\cal AC}$}
    74 
    75 \usepackage{color}
    76 \definecolor{lgray}{RGB}{238,238,238}
    77 
    78 %
    79 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    80 %                                                         %
    81 % How to use hyperref                                     %
    82 % -------------------                                     %
    83 %                                                         %
    84 % Probably the only way you will need to use the hyperref %
    85 % package is as follows.  To make some text, say          %
    86 % "My Text Link", into a link to the URL                  %
    87 % http://something.somewhere.com/mystuff, use             %
    88 %                                                         %
    89 % \href{http://something.somewhere.com/mystuff}{My Text Link}
    90 %                                                         %
    91 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    92 %
    93 \begin{document}
    94 %
    95 % document title
    96 %
    97 \title{Trials with TP-based Programming
    98 \\
    99 for Interactive Course Material}%
   100 %
   101 % Single author.  Please supply at least your name,
   102 % email address, and affiliation here.
   103 %
   104 \author{\begin{tabular}{c}
   105 \textit{Jan Ro\v{c}nik} \\
   106 jan.rocnik@student.tugraz.at \\
   107 IST, SPSC\\
   108 Graz University of Technologie\\
   109 Austria\end{tabular}
   110 }%
   111 %
   112 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   113 %                                                         %
   114 % eJMT commands - do not change these                     %
   115 %                                                         %
   116 \date{}                                                   %
   117 \maketitle                                                %
   118 %                                                         %
   119 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   120 %
   121 % abstract
   122 %
   123 \begin{abstract}
   124 
   125 Traditional course material in engineering disciplines lacks an
   126 important component, interactive support for step-wise problem
   127 solving. Theorem-Proving (TP) technology is appropriate for one part
   128 of such support, in checking user-input. For the other part of such
   129 support, guiding the learner towards a solution, another kind of
   130 technology is required. %TODO ... connect to prototype ...
   131 
   132 A prototype combines TP with a programming language, the latter
   133 interpreted in a specific way: certain statements in a program, called
   134 tactics, are treated as breakpoints where control is handed over to
   135 the user. An input formula is checked by TP (using logical context
   136 built up by the interpreter); and if a learner gets stuck, a program
   137 describing the steps towards a solution of a problem ``knows the next
   138 step''. This kind of interpretation is called Lucas-Interpretation for
   139 \emph{TP-based programming languages}.
   140 
   141 This paper describes the prototype's TP-based programming language
   142 within a case study creating interactive material for an advanced
   143 course in Signal Processing: implementation of definitions and
   144 theorems in TP, formal specification of a problem and step-wise
   145 development of the program solving the problem. Experiences with the
   146 ork flow in iterative development with testing and identifying errors
   147 are described, too. The description clarifies the components missing
   148 in the prototype's language as well as deficiencies experienced during
   149 programming.
   150 \par
   151 These experiences are particularly notable, because the author is the
   152 first programmer using the language beyond the core team which
   153 developed the prototype's TP-based language interpreter.
   154 %
   155 \end{abstract}%
   156 %
   157 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   158 %                                                         %
   159 % eJMT command                                            %
   160 %                                                         %
   161 \thispagestyle{fancy}                                     %
   162 %                                                         %
   163 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   164 %
   165 % Please use the following to indicate sections, subsections,
   166 % etc.  Please also use \subsubsection{...}, \paragraph{...}
   167 % and \subparagraph{...} as necessary.
   168 %
   169 
   170 \section{Introduction}\label{intro}
   171 
   172 % \paragraph{Didactics of mathematics} 
   173 %WN: wenn man in einem high-quality paper von 'didactics' spricht, 
   174 %WN muss man am state-of-the-art ankn"upfen -- siehe
   175 %WN W.Neuper, On the Emergence of TP-based Educational Math Assistants
   176 % faces a specific issue, a gap
   177 % between (1) introduction of math concepts and skills and (2)
   178 % application of these concepts and skills, which usually are separated
   179 % into different units in curricula (for good reasons). For instance,
   180 % (1) teaching partial fraction decomposition is separated from (2)
   181 % application for inverse Z-transform in signal processing.
   182 % 
   183 % \par This gap is an obstacle for applying math as an fundamental
   184 % thinking technology in engineering: In (1) motivation is lacking
   185 % because the question ``What is this stuff good for?'' cannot be
   186 % treated sufficiently, and in (2) the ``stuff'' is not available to
   187 % students in higher semesters as widespread experience shows.
   188 % 
   189 % \paragraph{Motivation} taken by this didactic issue on the one hand,
   190 % and ongoing research and development on a novel kind of educational
   191 % mathematics assistant at Graz University of
   192 % Technology~\footnote{http://www.ist.tugraz.at/isac/} promising to
   193 % scope with this issue on the other hand, several institutes are
   194 % planning to join their expertise: the Institute for Information
   195 % Systems and Computer Media (IICM), the Institute for Software
   196 % Technology (IST), the Institutes for Mathematics, the Institute for
   197 % Signal Processing and Speech Communication (SPSC), the Institute for
   198 % Structural Analysis and the Institute of Electrical Measurement and
   199 % Measurement Signal Processing.
   200 %WN diese Information ist f"ur das Paper zu spezielle, zu aktuell 
   201 %WN und damit zu verg"anglich.
   202 % \par This thesis is the first attempt to tackle the above mentioned
   203 % issue, it focuses on Telematics, because these specific studies focus
   204 % on mathematics in \emph{STEOP}, the introductory orientation phase in
   205 % Austria. \emph{STEOP} is considered an opportunity to investigate the
   206 % impact of {\sisac}'s prototype on the issue and others.
   207 % 
   208 
   209 \paragraph{Traditional course material} in engineering disciplines lacks an
   210 important component, interactive support for step-wise problem
   211 solving. Theorem-Proving (TP) technology can provide such support by
   212 specific services. An important part of such services is called
   213 ``next-step-guidance'', generated by a specific kind of ``TP-based
   214 programming language''. In the
   215 {\sisac}-project~\footnote{http://www.ist.tugraz.at/projects/isac/} such
   216 a language is prototyped in line with~\cite{plmms10} and built upon
   217 the theorem prover
   218 Isabelle~\cite{Nipkow-Paulson-Wenzel:2002}\footnote{http://isabelle.in.tum.de/}.
   219 The TP services are coordinated by a specific interpreter for the
   220 programming language, called
   221 Lucas-Interpreter~\cite{wn:lucas-interp-12}. The language and the
   222 interpreter will be briefly re-introduced in order to make the paper
   223 self-contained.
   224 
   225 \subparagraph{The main part} of the paper is an account of first experiences
   226 with programming in this TP-based language. The experience was gained
   227 in a case study by the author. The author was considered an ideal
   228 candidate for this study for the following reasons: as a student in
   229 Telematics (computer science with focus on Signal Processing) he had
   230 general knowledge in programming as well as specific domain knowledge
   231 in Signal Processing; and he was not involved in the development of
   232 {\sisac}'s programming language and interpeter, thus a novice to the
   233 language.
   234 
   235 \subparagraph{The goal} of the case study was (1) some TP-based programs for
   236 interactive course material for a specific ``Adavanced Signal
   237 Processing Lab'' in a higher semester, (2) respective program
   238 development with as little advice from the {\sisac}-team and (3) records
   239 and comments for the main steps of development in an Isabelle theory;
   240 this theory should provide guidelines for future programmers. An
   241 excerpt from this theory is the main part of this paper.
   242 \par
   243 The paper will use the problem in Fig.\ref{fig-interactive} as a
   244 running example:
   245 \begin{figure} [htb]
   246 \begin{center}
   247 \includegraphics[width=140mm]{fig/isac-Ztrans-math-3}
   248 %\includegraphics[width=140mm]{fig/isac-Ztrans-math}
   249 \caption{Step-wise problem solving guided by the TP-based program}
   250 \label{fig-interactive}
   251 \end{center}
   252 \end{figure}
   253 
   254 \paragraph{The problem is} from the domain of Signal Processing and requests to
   255 determine the inverse Z-transform for a given term. Fig.\ref{fig-interactive}
   256 also shows the beginning of the interactive construction of a solution
   257 for the problem. This construction is done in the right window named
   258 ``Worksheet''.
   259 \par
   260 User-interaction on the Worksheet is {\em checked} and {\em guided} by
   261 TP services:
   262 \begin{enumerate}
   263 \item Formulas input by the user are {\em checked} by TP: such a
   264 formula establishes a proof situation --- the prover has to derive the
   265 formula from the logical context. The context is built up from the
   266 formal specification of the problem (here hidden from the user) by the
   267 Lucas-Interpreter.
   268 \item If the user gets stuck, the program developed below in this
   269 paper ``knows the next step'' from behind the scenes. How the latter
   270 TP-service is exploited by dialogue authoring is out of scope of this
   271 paper and can be studied in~\cite{gdaroczy-EP-13}.
   272 \end{enumerate} It should be noted that the programmer using the
   273 TP-based language is not concerned with interaction at all; we will
   274 see that the program contains neither input-statements nor
   275 output-statements. Rather, interaction is handled by services
   276 generated automatically.
   277 \par
   278 So there is a clear separation of concerns: Dialogues are
   279 adapted by dialogue authors (in Java-based tools), using automatically
   280 generated TP services, while the TP-based program is written by
   281 mathematics experts (in Isabelle/ML). The latter is concern of this
   282 paper.
   283 
   284 \paragraph{The paper is structed} as follows: The introduction
   285 \S\ref{intro} is followed by a brief re-introduction of the TP-based
   286 programming language in \S\ref{PL}, which extends the executable
   287 fragment of Isabelle's language (\S\ref{PL-isab}) by tactics which
   288 play a specific role in Lucas-Interpretation and in providing the TP
   289 services (\S\ref{PL-tacs}). The main part in \S\ref{trial} describes
   290 the main steps in developing the program for the running example:
   291 prepare domain knowledge, implement the formal specification of the
   292 problem, prepare the environment for the program, implement the
   293 program. The workflow of programming, debugging and testing is
   294 described in \S\ref{workflow}. The conclusion \S\ref{conclusion} will
   295 give directions identified for future development. 
   296 
   297 
   298 \section{\isac's Prototype for a Programming Language}\label{PL} 
   299 The prototype's language extends the executable fragment in the
   300 language of the theorem prover
   301 Isabelle~\cite{Nipkow-Paulson-Wenzel:2002}\footnote{http://isabelle.in.tum.de/}
   302 by tactics which have a specific role in Lucas-Interpretation.
   303 
   304 \subsection{The Executable Fragment of Isabelle's Language}\label{PL-isab}
   305 The executable fragment consists of data-type and function
   306 definitions.  It's usability even suggests that fragment for
   307 introductory courses \cite{nipkow-prog-prove}. HOL is a typed logic
   308 whose type system resembles that of functional programming
   309 languages. Thus there are
   310 \begin{description}
   311 \item[base types,] in particular \textit{bool}, the type of truth
   312 values, \textit{nat}, \textit{int}, \textit{complex}, and the types of
   313 natural, integer and complex numbers respectively in mathematics.
   314 \item[type constructors] allow to define arbitrary types, from
   315 \textit{set}, \textit{list} to advanced data-structures like
   316 \textit{trees}, red-black-trees etc.
   317 \item[function types,] denoted by $\Rightarrow$.
   318 \item[type variables,] denoted by $^\prime a, ^\prime b$ etc, provide
   319 type polymorphism. Isabelle automatically computes the type of each
   320 variable in a term by use of Hindley-Milner type inference
   321 \cite{pl:hind97,Milner-78}.
   322 \end{description}
   323 
   324 \textbf{Terms} are formed as in functional programming by applying
   325 functions to arguments. If $f$ is a function of type
   326 $\tau_1\Rightarrow \tau_2$ and $t$ is a term of type $\tau_1$ then
   327 $f\;t$ is a term of type~$\tau_2$. $t\;::\;\tau$ means that term $t$
   328 has type $\tau$. There are many predefined infix symbols like $+$ and
   329 $\leq$ most of which are overloaded for various types.
   330 
   331 HOL also supports some basic constructs from functional programming:
   332 {\it\label{isabelle-stmts}
   333 \begin{tabbing} 123\=\kill
   334 \>$( \; {\tt if} \; b \; {\tt then} \; t_1 \; {\tt else} \; t_2 \;)$\\
   335 \>$( \; {\tt let} \; x=t \; {\tt in} \; u \; )$\\
   336 \>$( \; {\tt case} \; t \; {\tt of} \; {\it pat}_1
   337   \Rightarrow t_1 \; |\dots| \; {\it pat}_n\Rightarrow t_n \; )$
   338 \end{tabbing} }
   339 \noindent \textit{The running example's program uses some of these elements
   340 (marked by {\tt tt-font} on p.\pageref{expl-program}): ${\tt
   341 let}\dots{\tt in}$ in lines $02 \dots 11$, as well as {\tt last} for
   342 lists and {\tt o} for functional (forward) composition in line
   343 $10$. In fact, the whole program is an Isabelle term with specific
   344 function constants like {\sc program}, {\sc Substitute} and {\sc
   345 Rewrite\_Set\_Inst} in lines $01$ and $10$ respectively.}
   346 
   347 % Terms may also contain $\lambda$-abstractions. For example, $\lambda
   348 % x. \; x$ is the identity function.
   349 
   350 %JR warum auskommentiert? WN2...
   351 %WN2 weil ein Punkt wie dieser in weiteren Zusammenh"angen innerhalb
   352 %WN2 des Papers auftauchen m"usste; nachdem ich einen solchen
   353 %WN2 Zusammenhang _noch_ nicht sehe, habe ich den Punkt _noch_ nicht
   354 %WN2 gel"oscht.
   355 %WN2 Wenn der Punkt nicht weiter gebraucht wird, nimmt er nur wertvollen
   356 %WN2 Platz f"ur Anderes weg.
   357 
   358 \textbf{Formulae} are terms of type \textit{bool}. There are the basic
   359 constants \textit{True} and \textit{False} and the usual logical
   360 connectives (in decreasing order of precedence): $\neg, \land, \lor,
   361 \rightarrow$.
   362 
   363 \textbf{Equality} is available in the form of the infix function $=$
   364 of type $a \Rightarrow a \Rightarrow {\it bool}$. It also works for
   365 formulas, where it means ``if and only if''.
   366 
   367 \textbf{Quantifiers} are written $\forall x. \; P$ and $\exists x. \;
   368 P$.  Quantifiers lead to non-executable functions, so functions do not
   369 always correspond to programs, for instance, if comprising \\$(
   370 \;{\it if} \; \exists x.\;P \; {\it then} \; e_1 \; {\it else} \; e_2
   371 \;)$.
   372 
   373 \subsection{\isac's Tactics for Lucas-Interpretation}\label{PL-tacs}
   374 The prototype extends Isabelle's language by specific statements
   375 called tactics~\footnote{{\sisac}'s tactics are different from
   376 Isabelle's tactics: the former concern steps in a calculation, the
   377 latter concern proof steps.}  and tacticals. For the programmer these
   378 statements are functions with the following signatures:
   379 
   380 \begin{description}
   381 \item[Rewrite:] ${\it theorem}\Rightarrow{\it term}\Rightarrow{\it
   382 term} * {\it term}\;{\it list}$:
   383 this tactic appplies {\it theorem} to a {\it term} yielding a {\it
   384 term} and a {\it term list}, the list are assumptions generated by
   385 conditional rewriting. For instance, the {\it theorem}
   386 $b\not=0\land c\not=0\Rightarrow\frac{a\cdot c}{b\cdot c}=\frac{a}{b}$
   387 applied to the {\it term} $\frac{2\cdot x}{3\cdot x}$ yields
   388 $(\frac{2}{3}, [x\not=0])$.
   389 
   390 \item[Rewrite\_Set:] ${\it ruleset}\Rightarrow{\it
   391 term}\Rightarrow{\it term} * {\it term}\;{\it list}$:
   392 this tactic appplies {\it ruleset} to a {\it term}; {\it ruleset} is
   393 a confluent and terminating term rewrite system, in general. If
   394 none of the rules ({\it theorem}s) is applicable on interpretation
   395 of this tactic, an exception is thrown.
   396 
   397 % \item[Rewrite\_Inst:] ${\it substitution}\Rightarrow{\it
   398 % theorem}\Rightarrow{\it term}\Rightarrow{\it term} * {\it term}\;{\it
   399 % list}$:
   400 % 
   401 % \item[Rewrite\_Set\_Inst:] ${\it substitution}\Rightarrow{\it
   402 % ruleset}\Rightarrow{\it term}\Rightarrow{\it term} * {\it term}\;{\it
   403 % list}$:
   404 
   405 \item[Substitute:] ${\it substitution}\Rightarrow{\it
   406 term}\Rightarrow{\it term}$:
   407 
   408 \item[Take:] ${\it term}\Rightarrow{\it term}$:
   409 this tactic has no effect in the program; but it creates a side-effect
   410 by Lucas-Interpretation (see below) and writes {\it term} to the
   411 Worksheet.
   412 
   413 \item[Subproblem:] ${\it theory} * {\it specification} * {\it
   414 method}\Rightarrow{\it argument}\;{\it list}\Rightarrow{\it term}$:
   415 this tactic allows to enter a phase of interactive specification
   416 of a theory ($\Re$, $\cal C$, etc), a formal specification (for instance,
   417 a specific type of equation) and a method (for instance, solving an
   418 equation symbolically or numerically).
   419 
   420 \end{description}
   421 The tactics play a specific role in
   422 Lucas-Interpretation~\cite{wn:lucas-interp-12}: they are treated as
   423 break-points and control is handed over to the user. The user is free
   424 to investigate underlying knowledge, applicable theorems, etc.  And
   425 the user can proceed constructing a solution by input of a tactic to
   426 be applied or by input of a formula; in the latter case the
   427 Lucas-Interpreter has built up a logical context (initialised with the
   428 precondition of the formal specification) such that Isabelle can
   429 derive the formula from this context --- or give feedback, that no
   430 derivation can be found.
   431 
   432 \subsection{Tacticals for Control of Interpretation}
   433 The flow of control in a program can be determined by {\tt if then else}
   434 and {\tt case of} as mentioned on p.\pageref{isabelle-stmts} and also
   435 by additional tacticals:
   436 \begin{description}
   437 \item[Repeat:] ${\it tactic}\Rightarrow{\it term}\Rightarrow{\it
   438 term}$: iterates over tactics which take a {\it term} as argument as
   439 long as a tactic is applicable (for instance, {\it Rewrite\_Set} might
   440 not be applicable).
   441 
   442 \item[Try:] ${\it tactic}\Rightarrow{\it term}\Rightarrow{\it term}$:
   443 if {\it tactic} is applicable, then it is applied to {\it term},
   444 otherwise {\it term} is passed on unchanged.
   445 
   446 \item[Or:] ${\it tactic}\Rightarrow{\it tactic}\Rightarrow{\it
   447 term}\Rightarrow{\it term}$:
   448 
   449 
   450 \item[@@:] ${\it tactic}\Rightarrow{\it tactic}\Rightarrow{\it
   451 term}\Rightarrow{\it term}$:
   452 
   453 \item[While:] ${\it term::bool}\Rightarrow{\it tactic}\Rightarrow{\it
   454 term}\Rightarrow{\it term}$:
   455 
   456 \end{description}
   457 
   458 no input / output --- Lucas-Interpretation
   459 
   460 .\\.\\.\\TODO\\.\\.\\
   461 
   462 \section{Development of a Program on Trial}\label{trial} 
   463 As mentioned above, {\sisac} is an experimental system for a proof of
   464 concept for Lucas-Interpretation~\cite{wn:lucas-interp-12}.  The
   465 latter interprets a specific kind of TP-based programming language,
   466 which is as experimental as the whole prototype --- so programming in
   467 this language can be only ``on trial'', presently.  However, as a
   468 prototype, the language addresses essentials described below.
   469 
   470 \subsection{Mechanization of Math --- Domain Engineering}\label{isabisac}
   471 
   472 %WN was Fachleute unter obigem Titel interessiert findet sich
   473 %WN unterhalb des auskommentierten Textes.
   474 
   475 %WN der Text unten spricht Benutzer-Aspekte anund ist nicht speziell
   476 %WN auf Computer-Mathematiker fokussiert.
   477 % \paragraph{As mentioned in the introduction,} a prototype of an
   478 % educational math assistant called
   479 % {{\sisac}}\footnote{{{\sisac}}=\textbf{Isa}belle for
   480 % \textbf{C}alculations, see http://www.ist.tugraz.at/isac/.} bridges
   481 % the gap between (1) introducation and (2) application of mathematics:
   482 % {{\sisac}} is based on Computer Theorem Proving (TP), a technology which
   483 % requires each fact and each action justified by formal logic, so
   484 % {{{\sisac}{}}} makes justifications transparent to students in
   485 % interactive step-wise problem solving. By that way {{\sisac}} already
   486 % can serve both:
   487 % \begin{enumerate}
   488 %   \item Introduction of math stuff (in e.g. partial fraction
   489 % decomposition) by stepwise explaining and exercising respective
   490 % symbolic calculations with ``next step guidance (NSG)'' and rigorously
   491 % checking steps freely input by students --- this also in context with
   492 % advanced applications (where the stuff to be taught in higher
   493 % semesters can be skimmed through by NSG), and
   494 %   \item Application of math stuff in advanced engineering courses
   495 % (e.g. problems to be solved by inverse Z-transform in a Signal
   496 % Processing Lab) and now without much ado about basic math techniques
   497 % (like partial fraction decomposition): ``next step guidance'' supports
   498 % students in independently (re-)adopting such techniques.
   499 % \end{enumerate} 
   500 % Before the question is answers, how {{\sisac}}
   501 % accomplishes this task from a technical point of view, some remarks on
   502 % the state-of-the-art is given, therefor follow up Section~\ref{emas}.
   503 % 
   504 % \subsection{Educational Mathematics Assistants (EMAs)}\label{emas}
   505 % 
   506 % \paragraph{Educational software in mathematics} is, if at all, based
   507 % on Computer Algebra Systems (CAS, for instance), Dynamic Geometry
   508 % Systems (DGS, for instance \footnote{GeoGebra http://www.geogebra.org}
   509 % \footnote{Cinderella http://www.cinderella.de/}\footnote{GCLC
   510 % http://poincare.matf.bg.ac.rs/~janicic/gclc/}) or spread-sheets. These
   511 % base technologies are used to program math lessons and sometimes even
   512 % exercises. The latter are cumbersome: the steps towards a solution of
   513 % such an interactive exercise need to be provided with feedback, where
   514 % at each step a wide variety of possible input has to be foreseen by
   515 % the programmer - so such interactive exercises either require high
   516 % development efforts or the exercises constrain possible inputs.
   517 % 
   518 % \subparagraph{A new generation} of educational math assistants (EMAs)
   519 % is emerging presently, which is based on Theorem Proving (TP). TP, for
   520 % instance Isabelle and Coq, is a technology which requires each fact
   521 % and each action justified by formal logic. Pushed by demands for
   522 % \textit{proven} correctness of safety-critical software TP advances
   523 % into software engineering; from these advancements computer
   524 % mathematics benefits in general, and math education in particular. Two
   525 % features of TP are immediately beneficial for learning:
   526 % 
   527 % \paragraph{TP have knowledge in human readable format,} that is in
   528 % standard predicate calculus. TP following the LCF-tradition have that
   529 % knowledge down to the basic definitions of set, equality,
   530 % etc~\footnote{http://isabelle.in.tum.de/dist/library/HOL/HOL.html};
   531 % following the typical deductive development of math, natural numbers
   532 % are defined and their properties
   533 % proven~\footnote{http://isabelle.in.tum.de/dist/library/HOL/Number\_Theory/Primes.html},
   534 % etc. Present knowledge mechanized in TP exceeds high-school
   535 % mathematics by far, however by knowledge required in software
   536 % technology, and not in other engineering sciences.
   537 % 
   538 % \paragraph{TP can model the whole problem solving process} in
   539 % mathematical problem solving {\em within} a coherent logical
   540 % framework. This is already being done by three projects, by
   541 % Ralph-Johan Back, by ActiveMath and by Carnegie Mellon Tutor.
   542 % \par
   543 % Having the whole problem solving process within a logical coherent
   544 % system, such a design guarantees correctness of intermediate steps and
   545 % of the result (which seems essential for math software); and the
   546 % second advantage is that TP provides a wealth of theories which can be
   547 % exploited for mechanizing other features essential for educational
   548 % software.
   549 % 
   550 % \subsubsection{Generation of User Guidance in EMAs}\label{user-guid}
   551 % 
   552 % One essential feature for educational software is feedback to user
   553 % input and assistance in coming to a solution.
   554 % 
   555 % \paragraph{Checking user input} by ATP during stepwise problem solving
   556 % is being accomplished by the three projects mentioned above
   557 % exclusively. They model the whole problem solving process as mentioned
   558 % above, so all what happens between formalized assumptions (or formal
   559 % specification) and goal (or fulfilled postcondition) can be
   560 % mechanized. Such mechanization promises to greatly extend the scope of
   561 % educational software in stepwise problem solving.
   562 % 
   563 % \paragraph{NSG (Next step guidance)} comprises the system's ability to
   564 % propose a next step; this is a challenge for TP: either a radical
   565 % restriction of the search space by restriction to very specific
   566 % problem classes is required, or much care and effort is required in
   567 % designing possible variants in the process of problem solving
   568 % \cite{proof-strategies-11}.
   569 % \par
   570 % Another approach is restricted to problem solving in engineering
   571 % domains, where a problem is specified by input, precondition, output
   572 % and postcondition, and where the postcondition is proven by ATP behind
   573 % the scenes: Here the possible variants in the process of problem
   574 % solving are provided with feedback {\em automatically}, if the problem
   575 % is described in a TP-based programing language: \cite{plmms10} the
   576 % programmer only describes the math algorithm without caring about
   577 % interaction (the respective program is functional and even has no
   578 % input or output statements!); interaction is generated as a
   579 % side-effect by the interpreter --- an efficient separation of concern
   580 % between math programmers and dialog designers promising application
   581 % all over engineering disciplines.
   582 % 
   583 % 
   584 % \subsubsection{Math Authoring in Isabelle/ISAC\label{math-auth}}
   585 % Authoring new mathematics knowledge in {{\sisac}} can be compared with
   586 % ``application programing'' of engineering problems; most of such
   587 % programing uses CAS-based programing languages (CAS = Computer Algebra
   588 % Systems; e.g. Mathematica's or Maple's programing language).
   589 % 
   590 % \paragraph{A novel type of TP-based language} is used by {{\sisac}{}}
   591 % \cite{plmms10} for describing how to construct a solution to an
   592 % engineering problem and for calling equation solvers, integration,
   593 % etc~\footnote{Implementation of CAS-like functionality in TP is not
   594 % primarily concerned with efficiency, but with a didactic question:
   595 % What to decide for: for high-brow algorithms at the state-of-the-art
   596 % or for elementary algorithms comprehensible for students?} within TP;
   597 % TP can ensure ``systems that never make a mistake'' \cite{casproto} -
   598 % are impossible for CAS which have no logics underlying.
   599 % 
   600 % \subparagraph{Authoring is perfect} by writing such TP based programs;
   601 % the application programmer is not concerned with interaction or with
   602 % user guidance: this is concern of a novel kind of program interpreter
   603 % called Lucas-Interpreter. This interpreter hands over control to a
   604 % dialog component at each step of calculation (like a debugger at
   605 % breakpoints) and calls automated TP to check user input following
   606 % personalized strategies according to a feedback module.
   607 % \par
   608 % However ``application programing with TP'' is not done with writing a
   609 % program: according to the principles of TP, each step must be
   610 % justified. Such justifications are given by theorems. So all steps
   611 % must be related to some theorem, if there is no such theorem it must
   612 % be added to the existing knowledge, which is organized in so-called
   613 % \textbf{theories} in Isabelle. A theorem must be proven; fortunately
   614 % Isabelle comprises a mechanism (called ``axiomatization''), which
   615 % allows to omit proofs. Such a theorem is shown in
   616 % Example~\ref{eg:neuper1}.
   617 
   618 The running example, introduced by Fig.\ref{fig-interactive} on
   619 p.\pageref{fig-interactive}, requires to determine the inverse $\cal
   620 Z$-transform for a class of functions. The domain of Signal Processing
   621 is accustomed to specific notation for the resulting functions, which
   622 are absolutely summable and are called TODO: $u[n]$, where $u$ is the
   623 function, $n$ is the argument and the brackets indicate that the
   624 arguments are TODO. Surprisingly, Isabelle accepts the rules for
   625 ${\cal Z}^{-1}$ in this traditional notation~\footnote{Isabelle
   626 experts might be particularly surprised, that the brackets do not
   627 cause errors in typing (as lists).}:
   628 %\vbox{
   629 % \begin{example}
   630   \label{eg:neuper1}
   631   {\small\begin{tabbing}
   632   123\=123\=123\=123\=\kill
   633   \hfill \\
   634   \>axiomatization where \\
   635   \>\>  rule1: ``${\cal Z}^{-1}\;1 = \delta [n]$'' and\\
   636   \>\>  rule2: ``$\vert\vert z \vert\vert > 1 \Rightarrow {\cal Z}^{-1}\;z / (z - 1) = u [n]$'' and\\
   637   \>\>  rule3: ``$\vert\vert$ z $\vert\vert$ < 1 ==> z / (z - 1) = -u [-n - 1]'' and \\
   638 %TODO
   639   \>\>  rule4: ``$\vert\vert$ z $\vert\vert$ > $\vert\vert$ $\alpha$ $\vert\vert$ ==> z / (z - $\alpha$) = $\alpha^n$ $\cdot$ u [n]'' and\\
   640 %TODO
   641   \>\>  rule5: ``$\vert\vert$ z $\vert\vert$ < $\vert\vert$ $\alpha$ $\vert\vert$ ==> z / (z - $\alpha$) = -($\alpha^n$) $\cdot$ u [-n - 1]'' and\\
   642 %TODO
   643   \>\>  rule6: ``$\vert\vert$ z $\vert\vert$ > 1 ==> z/(z - 1)$^2$ = n $\cdot$ u [n]''\\
   644 %TODO
   645   \end{tabbing}
   646   }
   647 % \end{example}
   648 %}
   649 These 6 rules can be used as conditional rewrite rules, depending on
   650 the respective convergence radius. Satisfaction from accordance with traditional notation
   651 contrasts with the above word {\em axiomatization}: As TP-based, the
   652 programming language expects these rules as {\em proved} theorems, and
   653 not as axioms implemented in the above brute force manner; otherwise
   654 all the verification efforts envisaged (like proof of the
   655 post-condition, see below) would be meaningless.
   656 
   657 Isabelle provides a large body of knowledge, rigorously proven from
   658 the basic axioms of mathematics~\footnote{This way of rigorously
   659 deriving all knowledge from first principles is called the
   660 LCF-paradigm in TP.}. In the case of the Z-Transform the most advanced
   661 knowledge can be found in the theoris on Multivariate
   662 Analysis~\footnote{http://isabelle.in.tum.de/dist/library/HOL/HOL-Multivariate\_Analysis}. However,
   663 building up knowledge such that a proof for the above rules would be
   664 reasonably short and easily comprehensible, still requires lots of
   665 work (and is definitely out of scope of our case study).
   666 
   667 \paragraph{At the state-of-the-art in mechanization of knowledge} in
   668 engineering sciences, the process does not stop with the mechanization of
   669 mathematics traditionally used in these sciences. Rather, ``Formal Methods''~\cite{TODO-formal-methods}
   670 are expected to proceed to formal and explicit description of physical items.  Signal Processing,
   671 for instance is concerned with physical devices for signal acquisition
   672 and reconstruction, which involve measuring a physical signal, storing
   673 it, and possibly later rebuilding the original signal or an
   674 approximation thereof. For digital systems, this typically includes
   675 sampling and quantization; devices for signal compression, including
   676 audio compression, image compression, and video compression, etc.
   677 ``Domain engineering''\cite{db-domain-engineering} is concerned with
   678 {\em specification} of these devices' components and features; this
   679 part in the process of mechanization is only at the beginning in domains
   680 like Signal Processing.
   681 
   682 \subparagraph{TP-based programming, concern of this paper,} is determined to
   683 add ``algorithmic knowledge'' in Fig.\ref{fig:mathuni} on
   684 p.\pageref{fig:mathuni}.  As we shall see below, TP-based programming
   685 starts with a formal {\em specification} of the problem to be solved.
   686 
   687 
   688 \subsection{Specification of the Problem}\label{spec}
   689 %WN <--> \chapter 7 der Thesis
   690 %WN die Argumentation unten sollte sich NUR auf Verifikation beziehen..
   691 
   692 The problem of the running example is textually described in
   693 Fig.\ref{fig-interactive} on p.\pageref{fig-interactive}. The {\em
   694 formal} specification of this problem, in traditional mathematical
   695 notation, could look lik is this:
   696 
   697 %WN Hier brauchen wir die Spezifikation des 'running example' ...
   698 
   699 %JR Habe input, output und precond vom Beispiel eingefügt brauche aber Hilfe bei
   700 %JR der post condition - die existiert für uns ja eigentlich nicht aka
   701 %JR haben sie bis jetzt nicht beachtet WN...
   702 %WN2 Mein Vorschlag ist, das TODO zu lassen und deutlich zu kommentieren.
   703 
   704   \label{eg:neuper2}
   705   {\small\begin{tabbing}
   706   123,\=postcond \=: \= $\forall \,A^\prime\, u^\prime \,v^\prime.\,$\=\kill
   707   \hfill \\
   708   Specification:\\
   709     \>input    \>: filterExpression $X=\frac{3}{(z-\frac{1}{4}+-\frac{1}{8}*\frac{1}{z}}$\\
   710   \>precond  \>: filterExpression continius on $\mathbb{R}$ \\
   711   \>output   \>: stepResponse $x[n]$ \\
   712   \>postcond \>:{\small  $\;A=2uv-u^2 \;\land\; (\frac{u}{2})^2+(\frac{v}{2})^2=r^2 \;\land$}\\
   713   \>     \>\>{\small $\;\forall \;A^\prime\; u^\prime \;v^\prime.\;(A^\prime=2u^\prime v^\prime-(u^\prime)^2 \land
   714   (\frac{u^\prime}{2})^2+(\frac{v^\prime}{2})^2=r^2) \Longrightarrow A^\prime \leq A$} \\
   715   \end{tabbing}
   716   }
   717 
   718 And this is the implementation of the formal specification in the present
   719 prototype, still bar-bones without support for authoring:
   720 %WN Kopie von Inverse_Z_Transform.thy, leicht versch"onert:
   721 {\footnotesize
   722 \begin{verbatim}
   723    01  store_specification
   724    02    (prepare_specification
   725    03      ["Jan Rocnik"]
   726    04      "pbl_SP_Ztrans_inv"
   727    05      thy
   728    06      ( ["Inverse", "Z_Transform", "SignalProcessing"],
   729    07        [ ("#Given", ["filterExpression X_eq"]),
   730    08          ("#Pre"  , ["X_eq is_continuous"]),
   731    19          ("#Find" , ["stepResponse n_eq"]),
   732    10          ("#Post" , [" TODO "])],
   733    11        append_rls Erls [Calc ("Atools.is_continuous", eval_is_continuous)], 
   734    12        NONE, 
   735    13        [["SignalProcessing","Z_Transform","Inverse"]]));
   736 \end{verbatim}}
   737 Although the above details are partly very technical, we explain them
   738 in order to document some intricacies of TP-based programming in the
   739 present state of the {\sisac} prototype:
   740 \begin{description}
   741 \item[01..02]\textit{store\_specification:} stores the result of the
   742 function \textit{prep\_specification} in a global reference
   743 \textit{Unsynchronized.ref}, which causes principal conflicts with
   744 Isabelle's asyncronous document model~\cite{Wenzel-11:doc-orient} and
   745 parallel execution~\cite{Makarius-09:parall-proof} and is under
   746 reconstruction already.
   747 
   748 \textit{prep\_pbt:} translates the specification to an internal format
   749 which allows efficient processing; see for instance line {\rm 07}
   750 below.
   751 \item[03..04] are the ``mathematics author'' holding the copy-rights
   752 and a unique identifier for the specification within {\sisac},
   753 complare line {\rm 06}.
   754 \item[05] is the Isabelle \textit{theory} required to parse the
   755 specification in lines {\rm 07..10}.
   756 \item[06] is a key into the tree of all specifications as presented to
   757 the user (where some branches might be hidden by the dialog
   758 component).
   759 \item[07..10] are the specification with input, pre-condition, output
   760 and post-condition respectively; the post-condition is not handled in
   761 the prototype presently.
   762 \item[11] creates a term rewriting system (abbreviated \textit{rls} in
   763 {\sisac}) which evaluates the pre-condition for the actual input data.
   764 Only if the evaluation yields \textit{True}, a program con be started.
   765 \item[12]\textit{NONE:} could be \textit{SOME ``solve ...''} for a
   766 problem associated to a function from Computer Algebra (like an
   767 equation solver) which is not the case here.
   768 \item[13] is the specific key into the tree of programs addressing a
   769 method which is able to find a solution which satiesfies the
   770 post-condition of the specification.
   771 \end{description}
   772 
   773 
   774 %WN die folgenden Erkl"arungen finden sich durch "grep -r 'datatype pbt' *"
   775 %WN ...
   776 %  type pbt = 
   777 %     {guh  : guh,         (*unique within this isac-knowledge*)
   778 %      mathauthors: string list, (*copyright*)
   779 %      init  : pblID,      (*to start refinement with*)
   780 %      thy   : theory,     (* which allows to compile that pbt
   781 %			  TODO: search generalized for subthy (ref.p.69*)
   782 %      (*^^^ WN050912 NOT used during application of the problem,
   783 %       because applied terms may be from 'subthy' as well as from super;
   784 %       thus we take 'maxthy'; see match_ags !*)
   785 %      cas   : term option,(*'CAS-command'*)
   786 %      prls  : rls,        (* for preds in where_*)
   787 %      where_: term list,  (* where - predicates*)
   788 %      ppc   : pat list,
   789 %      (*this is the model-pattern; 
   790 %       it contains "#Given","#Where","#Find","#Relate"-patterns
   791 %       for constraints on identifiers see "fun cpy_nam"*)
   792 %      met   : metID list}; (* methods solving the pbt*)
   793 %
   794 %WN weil dieser Code sehr unaufger"aumt ist, habe ich die Erkl"arungen
   795 %WN oben selbst geschrieben.
   796 
   797 
   798 
   799 
   800 %WN das w"urde ich in \sec\label{progr} verschieben und
   801 %WN das SubProblem partial fractions zum Erkl"aren verwenden.
   802 % Such a specification is checked before the execution of a program is
   803 % started, the same applies for sub-programs. In the following example
   804 % (Example~\ref{eg:subprob}) shows the call of such a subproblem:
   805 % 
   806 % \vbox{
   807 %   \begin{example}
   808 %   \label{eg:subprob}
   809 %   \hfill \\
   810 %   {\ttfamily \begin{tabbing}
   811 %   ``(L\_L::bool list) = (\=SubProblem (\=Test','' \\
   812 %   ``\>\>[linear,univariate,equation,test],'' \\
   813 %   ``\>\>[Test,solve\_linear])'' \\
   814 %   ``\>[BOOL equ, REAL z])'' \\
   815 %   \end{tabbing}
   816 %   }
   817 %   {\small\textit{
   818 %     \noindent If a program requires a result which has to be
   819 % calculated first we can use a subproblem to do so. In our specific
   820 % case we wanted to calculate the zeros of a fraction and used a
   821 % subproblem to calculate the zeros of the denominator polynom.
   822 %     }}
   823 %   \end{example}
   824 % }
   825 
   826 \subsection{Implementation of the Method}\label{meth}
   827 %WN <--> \chapter 7 der Thesis
   828 TODO
   829 \subsection{Preparation of Simplifiers for the Program}\label{simp}
   830 TODO
   831 \subsection{Preparation of ML-Functions}\label{funs}
   832 %WN <--> Thesis 6.1 -- 6.3: jene ausw"ahlen, die Du f"ur \label{progr}
   833 %WN brauchst
   834 TODO
   835 \subsection{Implementation of the TP-based Program}\label{progr}
   836 %WN <--> \chapter 8 der Thesis
   837 .\\.\\.\\
   838 
   839 {\small\it\begin{tabbing}
   840 123l\=123\=123\=123\=123\=123\=123\=123\=123\=(x \=123\=123\=\kill
   841 \>{\rm 01}\>  {\tt Program} InverseZTransform (X\_eq::bool) =   \\
   842 \>{\rm 02}\>\>  {\tt LET}                                       \\
   843 \>{\rm 03}\>\>\>  X\_eq = {\tt Take} X\_eq ;   \\
   844 \>{\rm 04}\>\>\>  X\_eq = {\tt Rewrite} ruleZY X\_eq ; \\
   845 \>{\rm 05}\>\>\>  (X\_z::real) = lhs X\_eq ;       \\ %no inside-out evaluation
   846 \>{\rm 06}\>\>\>  (z::real) = argument\_in X\_z; \\
   847 \>{\rm 07}\>\>\>  (part\_frac::real) = {\tt SubProblem} \\
   848 \>{\rm 08}\>\>\>\>\>\>\>\>  ( \> Isac, \\
   849 \>{\rm 08}\>\>\>\>\>\>\>\>\>  [partial\_fraction, rational, simplification]\\
   850 \>{\rm 09}\>\>\>\>\>\>\>\>\>  [simplification,of\_rationals,to\_partial\_fraction] ) \\
   851 \>{\rm 10}\>\>\>\>\>\>\>\>  [ (rhs X\_eq)::real, z::real ]; \\
   852 \>{\rm 12}\>\>\>  (X'\_eq::bool) = {\tt Take} ((X'::real =$>$ bool) z = ZZ\_1 part\_frac) ; \\
   853 
   854 \>{\rm 12}\>\>\>  X'\_eq = {\tt Rewrite\_Set} ruleYZ X'\_eq ;   \\
   855 \>{\rm 15}\>\>\>  X'\_eq = {\tt Rewrite\_Set} inverse\_z X'\_eq \\
   856 \>{\rm 16}\>\>  {\tt IN } \\
   857 \>{\rm 15}\>\>\>  X'\_eq
   858 \end{tabbing}}
   859 % ORIGINAL FROM Inverse_Z_Transform.thy
   860 % "Script InverseZTransform (X_eq::bool) =            "^(*([], Frm), Problem (Isac, [Inverse, Z_Transform, SignalProcessing])*)
   861 % "(let X = Take X_eq;                                "^(*([1], Frm), X z = 3 / (z - 1 / 4 + -1 / 8 * (1 / z))*)
   862 % "  X' = Rewrite ruleZY False X;                     "^(*([1], Res), ?X' z = 3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))*)
   863 % "  (X'_z::real) = lhs X';                           "^(*            ?X' z*)
   864 % "  (zzz::real) = argument_in X'_z;                  "^(*            z *)
   865 % "  (funterm::real) = rhs X';                        "^(*            3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))*)
   866 %
   867 % "  (pbz::real) = (SubProblem (Isac',                "^(**)
   868 % "    [partial_fraction,rational,simplification],    "^
   869 % "    [simplification,of_rationals,to_partial_fraction]) "^
   870 % "    [REAL funterm, REAL zzz]);                     "^(*([2], Res), 4 / (z - 1 / 2) + -4 / (z - -1 / 4)*)
   871 %
   872 % "  (pbz_eq::bool) = Take (X'_z = pbz);              "^(*([3], Frm), ?X' z = 4 / (z - 1 / 2) + -4 / (z - -1 / 4)*)
   873 % "  pbz_eq = Rewrite ruleYZ False pbz_eq;            "^(*([3], Res), ?X' z = 4 * (?z / (z - 1 / 2)) + -4 * (?z / (z - -1 / 4))*)
   874 % "  pbz_eq = drop_questionmarks pbz_eq;              "^(*               4 * (z / (z - 1 / 2)) + -4 * (z / (z - -1 / 4))*)
   875 % "  (X_zeq::bool) = Take (X_z = rhs pbz_eq);         "^(*([4], Frm), X_z = 4 * (z / (z - 1 / 2)) + -4 * (z / (z - -1 / 4))*)
   876 % "  n_eq = (Rewrite_Set inverse_z False) X_zeq;      "^(*([4], Res), X_z = 4 * (1 / 2) ^^^ ?n * ?u [?n] + -4 * (-1 / 4) ^^^ ?n * ?u [?n]*)
   877 % "  n_eq = drop_questionmarks n_eq                   "^(*            X_z = 4 * (1 / 2) ^^^ n * u [n] + -4 * (-1 / 4) ^^^ n * u [n]*)
   878 % "in n_eq)"                                            (*([], Res), X_z = 4 * (1 / 2) ^^^ n * u [n] + -4 * (-1 / 4) ^^^ n * u [n]*)
   879 
   880 
   881 .\\.\\.\\
   882 
   883 \section{Workflow of Programming in the Prototype}\label{workflow}
   884 
   885 \cite{makar-jedit-12}
   886 
   887 \subsection{Preparations and Trials}\label{flow-prep}
   888 TODO Build\_Inverse\_Z\_Transform.thy ... ``imports PolyEq DiffApp Partial\_Fractions''
   889 .\\.\\.\\
   890 
   891 \subsection{Implementation in Isabelle/{\isac}}\label{flow-impl}
   892 TODO Build\_Inverse\_Z\_Transform.thy ... ``imports Isac''
   893 .\\.\\.\\
   894 
   895 below: calculation on the left; on the right are the tactics in the
   896 program which created the respective formula on the left.
   897 {\small\it
   898 \begin{tabbing}
   899 123l\=123\=123\=123\=123\=123\=123\=123\=123\=123\=123\=123\=\kill
   900 \>{\rm 01}\> $\bullet$  \> {\tt Problem } (Inverse\_Z\_Transform, [Inverse, Z\_Transform, SignalProcessing])       \`\\
   901 \>{\rm 02}\>\> $\vdash\;\;X z = \frac{3}{z - \frac{1}{4} - \frac{1}{8} \cdot z^{-1}}$       \`{\footnotesize {\tt Take} X\_eq}\\
   902 \>{\rm 03}\>\> $X z = \frac{3}{z + \frac{-1}{4} + \frac{-1}{8} \cdot \frac{1}{z}}$          \`{\footnotesize {\tt Rewrite} ruleZY X\_eq}\\
   903 \>{\rm 04}\>\> $\bullet$\> {\tt Problem } [partial\_fraction,rational,simplification]        \`{\footnotesize {\tt SubProblem} \dots}\\
   904 \>{\rm 05}\>\>\>  $\vdash\;\;\frac{3}{z + \frac{-1}{4} + \frac{-1}{8} \cdot \frac{1}{z}}=$    \`- - -\\
   905 \>{\rm 06}\>\>\>  $\frac{24}{-1 + -2 \cdot z + 8 \cdot z^2}$                                   \`- - -\\
   906 \>{\rm 07}\>\>\>  $\bullet$\> solve ($-1 + -2 \cdot z + 8 \cdot z^2,\;z$ )                      \`- - -\\
   907 \>{\rm 08}\>\>\>\>   $\vdash$ \> $\frac{3}{z + \frac{-1}{4} + \frac{-1}{8} \cdot \frac{1}{z}}=0$ \`- - -\\
   908 \>{\rm 09}\>\>\>\>   $z = \frac{2+\sqrt{-4+8}}{16}\;\lor\;z = \frac{2-\sqrt{-4+8}}{16}$           \`- - -\\
   909 \>{\rm 10}\>\>\>\>   $z = \frac{1}{2}\;\lor\;z =$ \_\_\_                                           \`- - -\\
   910 \>        \>\>\>\>   \_\_\_                                                                        \`- - -\\
   911 \>{\rm 11}\>\> \dots\> $\frac{4}{z - \frac{1}{2}} + \frac{-4}{z - \frac{-1}{4}}$                   \`\\
   912 \>{\rm 12}\>\> $X^\prime z = {\cal Z}^{-1} (\frac{4}{z - \frac{1}{2}} + \frac{-4}{z - \frac{-1}{4}})$ \`{\footnotesize {\tt Take} ((X'::real =$>$ bool) z = ZZ\_1 part\_frac)}\\
   913 \>{\rm 13}\>\> $X^\prime z = {\cal Z}^{-1} (4\cdot\frac{z}{z - \frac{1}{2}} + -4\cdot\frac{z}{z - \frac{-1}{4}})$ \`{\footnotesize{\tt Rewrite\_Set} ruleYZ X'\_eq }\\
   914 \>{\rm 14}\>\> $X^\prime z = 4\cdot(\frac{1}{2})^n \cdot u [n] + -4\cdot(\frac{-1}{4})^n \cdot u [n]$  \`{\footnotesize {\tt Rewrite\_Set} inverse\_z X'\_eq}\\
   915 \>{\rm 15}\> \dots\> $X^\prime z = 4\cdot(\frac{1}{2})^n \cdot u [n] + -4\cdot(\frac{-1}{4})^n \cdot u [n]$ \`{\footnotesize {\tt Check\_Postcond}}
   916 \end{tabbing}}
   917 % ORIGINAL FROM Inverse_Z_Transform.thy
   918 %    "Script InverseZTransform (X_eq::bool) =            "^(*([], Frm), Problem (Isac, [Inverse, Z_Transform, SignalProcessing])*)
   919 %    "(let X = Take X_eq;                                "^(*([1], Frm), X z = 3 / (z - 1 / 4 + -1 / 8 * (1 / z))*)
   920 %    "  X' = Rewrite ruleZY False X;                     "^(*([1], Res), ?X' z = 3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))*)
   921 %    "  (X'_z::real) = lhs X';                           "^(*            ?X' z*)
   922 %    "  (zzz::real) = argument_in X'_z;                  "^(*            z *)
   923 %    "  (funterm::real) = rhs X';                        "^(*            3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))*)
   924 % 
   925 %    "  (pbz::real) = (SubProblem (Isac',                "^(**)
   926 %    "    [partial_fraction,rational,simplification],    "^
   927 %    "    [simplification,of_rationals,to_partial_fraction]) "^
   928 %    "    [REAL funterm, REAL zzz]);                     "^(*([2], Res), 4 / (z - 1 / 2) + -4 / (z - -1 / 4)*)
   929 % 
   930 %    "  (pbz_eq::bool) = Take (X'_z = pbz);              "^(*([3], Frm), ?X' z = 4 / (z - 1 / 2) + -4 / (z - -1 / 4)*)
   931 %    "  pbz_eq = Rewrite ruleYZ False pbz_eq;            "^(*([3], Res), ?X' z = 4 * (?z / (z - 1 / 2)) + -4 * (?z / (z - -1 / 4))*)
   932 %    "  pbz_eq = drop_questionmarks pbz_eq;              "^(*               4 * (z / (z - 1 / 2)) + -4 * (z / (z - -1 / 4))*)
   933 %    "  (X_zeq::bool) = Take (X_z = rhs pbz_eq);         "^(*([4], Frm), X_z = 4 * (z / (z - 1 / 2)) + -4 * (z / (z - -1 / 4))*)
   934 %    "  n_eq = (Rewrite_Set inverse_z False) X_zeq;      "^(*([4], Res), X_z = 4 * (1 / 2) ^^^ ?n * ?u [?n] + -4 * (-1 / 4) ^^^ ?n * ?u [?n]*)
   935 %    "  n_eq = drop_questionmarks n_eq                   "^(*            X_z = 4 * (1 / 2) ^^^ n * u [n] + -4 * (-1 / 4) ^^^ n * u [n]*)
   936 %    "in n_eq)"                                            (*([], Res), X_z = 4 * (1 / 2) ^^^ n * u [n] + -4 * (-1 / 4) ^^^ n * u [n]*)
   937 
   938 .\\.\\.\\
   939 
   940 \subsection{Transfer into the Isabelle/{\isac} Knowledge}\label{flow-trans}
   941 TODO http://www.ist.tugraz.at/isac/index.php/Extend\_ISAC\_Knowledge\#Add\_an\_example ?
   942 
   943 
   944 
   945 
   946 \newpage
   947 -------------------------------------------------------------------
   948 
   949 Material, falls noch Platz bleibt ...
   950 
   951 -------------------------------------------------------------------
   952 
   953 
   954 \subsubsection{Trials on Notation and Termination}
   955 
   956 \paragraph{Technical notations} are a big problem for our piece of software,
   957 but the reason for that isn't a fault of the software itself, one of the
   958 troubles comes out of the fact that different technical subtopics use different
   959 symbols and notations for a different purpose. The most famous example for such
   960 a symbol is the complex number $i$ (in cassique math) or $j$ (in technical
   961 math). In the specific part of signal processing one of this notation issues is
   962 the use of brackets --- we use round brackets for analoge signals and squared
   963 brackets for digital samples. Also if there is no problem for us to handle this
   964 fact, we have to tell the machine what notation leads to wich meaning and that
   965 this purpose seperation is only valid for this special topic - signal
   966 processing.
   967 \subparagraph{In the programming language} itself it is not possible to declare
   968 fractions, exponents, absolutes and other operators or remarks in a way to make
   969 them pretty to read; our only posssiblilty were ASCII characters and a handfull
   970 greek symbols like: $\alpha, \beta, \gamma, \phi,\ldots$.
   971 \par
   972 With the upper collected knowledge it is possible to check if we were able to
   973 donate all required terms and expressions.
   974 
   975 \subsubsection{Definition and Usage of Rules}
   976 
   977 \paragraph{The core} of our implemented problem is the Z-Transformation, due
   978 the fact that the transformation itself would require higher math which isn't
   979 yet avaible in our system we decided to choose the way like it is applied in
   980 labratory and problem classes at our university - by applying transformation
   981 rules (collected in transformation tables).
   982 \paragraph{Rules,} in {\sisac{}}'s programming language can be designed by the
   983 use of axiomatizations like shown in Example~\ref{eg:ruledef}
   984 
   985 \begin{example}
   986   \label{eg:ruledef}
   987   \hfill\\
   988   \begin{verbatim}
   989   axiomatization where
   990     rule1: ``1 = $\delta$[n]'' and
   991     rule2: ``|| z || > 1 ==> z / (z - 1) = u [n]'' and
   992     rule3: ``|| z || < 1 ==> z / (z - 1) = -u [-n - 1]''
   993   \end{verbatim}
   994 \end{example}
   995 
   996 This rules can be collected in a ruleset and applied to a given expression as
   997 follows in Example~\ref{eg:ruleapp}.
   998 
   999 \begin{example}
  1000   \hfill\\
  1001   \label{eg:ruleapp}
  1002   \begin{enumerate}
  1003   \item Store rules in ruleset:
  1004   \begin{verbatim}
  1005   val inverse_Z = append_rls "inverse_Z" e_rls
  1006     [ Thm ("rule1",num_str @{thm rule1}),
  1007       Thm ("rule2",num_str @{thm rule2}),
  1008       Thm ("rule3",num_str @{thm rule3})
  1009     ];\end{verbatim}
  1010   \item Define exression:
  1011   \begin{verbatim}
  1012   val sample_term = str2term "z/(z-1)+z/(z-</delta>)+1";\end{verbatim}
  1013   \item Apply ruleset:
  1014   \begin{verbatim}
  1015   val SOME (sample_term', asm) = 
  1016     rewrite_set_ thy true inverse_Z sample_term;\end{verbatim}
  1017   \end{enumerate}
  1018 \end{example}
  1019 
  1020 The use of rulesets makes it much easier to develop our designated applications,
  1021 but the programmer has to be careful and patient. When applying rulesets
  1022 two important issues have to be mentionend:
  1023 \subparagraph{How often} the rules have to be applied? In case of
  1024 transformations it is quite clear that we use them once but other fields
  1025 reuqire to apply rules until a special condition is reached (e.g.
  1026 a simplification is finished when there is nothing to be done left).
  1027 \subparagraph{The order} in which rules are applied often takes a big effect
  1028 and has to be evaluated for each purpose once again.
  1029 \par
  1030 In our special case of Signal Processing and the rules defined in
  1031 Example~\ref{eg:ruledef} we have to apply rule~1 first of all to transform all
  1032 constants. After this step has been done it no mather which rule fit's next.
  1033 
  1034 \subsubsection{Helping Functions}
  1035 %get denom, argument in
  1036 \subsubsection{Trials on equation solving}
  1037 %simple eq and problem with double fractions/negative exponents
  1038 
  1039 
  1040 
  1041 \subsection{Formalization of missing knowledge in Isabelle}
  1042 
  1043 \paragraph{A problem} behind is the mechanization of mathematic
  1044 theories in TP-bases languages. There is still a huge gap between
  1045 these algorithms and this what we want as a solution - in Example
  1046 Signal Processing. 
  1047 
  1048 \vbox{
  1049   \begin{example}
  1050     \label{eg:gap}
  1051     \[
  1052       X\cdot(a+b)+Y\cdot(c+d)=aX+bX+cY+dY
  1053     \]
  1054     {\small\textit{
  1055       \noindent A very simple example on this what we call gap is the
  1056 simplification above. It is needles to say that it is correct and also
  1057 Isabelle for fills it correct - \emph{always}. But sometimes we don't
  1058 want expand such terms, sometimes we want another structure of
  1059 them. Think of a problem were we now would need only the coefficients
  1060 of $X$ and $Y$. This is what we call the gap between mechanical
  1061 simplification and the solution.
  1062     }}
  1063   \end{example}
  1064 }
  1065 
  1066 \paragraph{We are not able to fill this gap,} until we have to live
  1067 with it but first have a look on the meaning of this statement:
  1068 Mechanized math starts from mathematical models and \emph{hopefully}
  1069 proceeds to match physics. Academic engineering starts from physics
  1070 (experimentation, measurement) and then proceeds to mathematical
  1071 modeling and formalization. The process from a physical observance to
  1072 a mathematical theory is unavoidable bound of setting up a big
  1073 collection of standards, rules, definition but also exceptions. These
  1074 are the things making mechanization that difficult.
  1075 
  1076 \vbox{
  1077   \begin{example}
  1078     \label{eg:units}
  1079     \[
  1080       m,\ kg,\ s,\ldots
  1081     \]
  1082     {\small\textit{
  1083       \noindent Think about some units like that one's above. Behind
  1084 each unit there is a discerning and very accurate definition: One
  1085 Meter is the distance the light travels, in a vacuum, through the time
  1086 of 1 / 299.792.458 second; one kilogram is the weight of a
  1087 platinum-iridium cylinder in paris; and so on. But are these
  1088 definitions usable in a computer mechanized world?!
  1089     }}
  1090   \end{example}
  1091 }
  1092 
  1093 \paragraph{A computer} or a TP-System builds on programs with
  1094 predefined logical rules and does not know any mathematical trick
  1095 (follow up example \ref{eg:trick}) or recipe to walk around difficult
  1096 expressions. 
  1097 
  1098 \vbox{
  1099   \begin{example}
  1100     \label{eg:trick}
  1101   \[ \frac{1}{j\omega}\cdot\left(e^{-j\omega}-e^{j3\omega}\right)= \]
  1102   \[ \frac{1}{j\omega}\cdot e^{-j2\omega}\cdot\left(e^{j\omega}-e^{-j\omega}\right)=
  1103      \frac{1}{\omega}\, e^{-j2\omega}\cdot\colorbox{lgray}{$\frac{1}{j}\,\left(e^{j\omega}-e^{-j\omega}\right)$}= \]
  1104   \[ \frac{1}{\omega}\, e^{-j2\omega}\cdot\colorbox{lgray}{$2\, sin(\omega)$} \]
  1105     {\small\textit{
  1106       \noindent Sometimes it is also useful to be able to apply some
  1107 \emph{tricks} to get a beautiful and particularly meaningful result,
  1108 which we are able to interpret. But as seen in this example it can be
  1109 hard to find out what operations have to be done to transform a result
  1110 into a meaningful one.
  1111     }}
  1112   \end{example}
  1113 }
  1114 
  1115 \paragraph{The only possibility,} for such a system, is to work
  1116 through its known definitions and stops if none of these
  1117 fits. Specified on Signal Processing or any other application it is
  1118 often possible to walk through by doing simple creases. This creases
  1119 are in general based on simple math operational but the challenge is
  1120 to teach the machine \emph{all}\footnote{Its pride to call it
  1121 \emph{all}.} of them. Unfortunately the goal of TP Isabelle is to
  1122 reach a high level of \emph{all} but it in real it will still be a
  1123 survey of knowledge which links to other knowledge and {{\sisac}{}} a
  1124 trainer and helper but no human compensating calculator. 
  1125 \par
  1126 {{{\sisac}{}}} itself aims to adds \emph{Algorithmic Knowledge} (formal
  1127 specifications of problems out of topics from Signal Processing, etc.)
  1128 and \emph{Application-oriented Knowledge} to the \emph{deductive} axis of
  1129 physical knowledge. The result is a three-dimensional universe of
  1130 mathematics seen in Figure~\ref{fig:mathuni}.
  1131 
  1132 \begin{figure}
  1133   \begin{center}
  1134     \includegraphics{fig/universe}
  1135     \caption{Didactic ``Math-Universe'': Algorithmic Knowledge (Programs) is
  1136              combined with Application-oriented Knowledge (Specifications) and Deductive Knowledge (Axioms, Definitions, Theorems). The Result
  1137              leads to a three dimensional math universe.\label{fig:mathuni}}
  1138   \end{center}
  1139 \end{figure}
  1140 
  1141 %WN Deine aktuelle Benennung oben wird Dir kein Fachmann abnehmen;
  1142 %WN bitte folgende Bezeichnungen nehmen:
  1143 %WN 
  1144 %WN axis 1: Algorithmic Knowledge (Programs)
  1145 %WN axis 2: Application-oriented Knowledge (Specifications)
  1146 %WN axis 3: Deductive Knowledge (Axioms, Definitions, Theorems)
  1147 %WN 
  1148 %WN und bitte die R"ander von der Grafik wegschneiden (was ich f"ur *.pdf
  1149 %WN nicht hinkriege --- weshalb ich auch die eJMT-Forderung nicht ganz
  1150 %WN verstehe, separierte PDFs zu schicken; ich w"urde *.png schicken)
  1151 
  1152 %JR Ränder und beschriftung geändert. Keine Ahnung warum eJMT sich pdf's
  1153 %JR wünschen, würde ebenfalls png oder ähnliches verwenden, aber wenn pdf's
  1154 %JR gefordert werden WN2...
  1155 %WN2 meiner Meinung nach hat sich eJMT unklar ausgedr"uckt (z.B. kann
  1156 %WN2 man meines Wissens pdf-figures nicht auf eine bestimmte Gr"osse
  1157 %WN2 zusammenschneiden um die R"ander weg zu bekommen)
  1158 %WN2 Mein Vorschlag ist, in umserem tex-file bei *.png zu bleiben und
  1159 %WN2 png + pdf figures mitzuschicken.
  1160 
  1161 \subsection{Notes on Problems with Traditional Notation}
  1162 
  1163 \paragraph{During research} on these topic severely problems on
  1164 traditional notations have been discovered. Some of them have been
  1165 known in computer science for many years now and are still unsolved,
  1166 one of them aggregates with the so called \emph{Lambda Calculus},
  1167 Example~\ref{eg:lamda} provides a look on the problem that embarrassed
  1168 us.
  1169 
  1170 \vbox{
  1171   \begin{example}
  1172     \label{eg:lamda}
  1173 
  1174   \[ f(x)=\ldots\;  \quad R \rightarrow \quad R \]
  1175 
  1176 
  1177   \[ f(p)=\ldots\;  p \in \quad R \]
  1178 
  1179     {\small\textit{
  1180       \noindent Above we see two equations. The first equation aims to
  1181 be a mapping of an function from the reel range to the reel one, but
  1182 when we change only one letter we get the second equation which
  1183 usually aims to insert a reel point $p$ into the reel function. In
  1184 computer science now we have the problem to tell the machine (TP) the
  1185 difference between this two notations. This Problem is called
  1186 \emph{Lambda Calculus}.
  1187     }}
  1188   \end{example}
  1189 }
  1190 
  1191 \paragraph{An other problem} is that terms are not full simplified in
  1192 traditional notations, in {{\sisac}} we have to simplify them complete
  1193 to check weather results are compatible or not. in e.g. the solutions
  1194 of an second order linear equation is an rational in {{\sisac}} but in
  1195 tradition we keep fractions as long as possible and as long as they
  1196 aim to be \textit{beautiful} (1/8, 5/16,...).
  1197 \subparagraph{The math} which should be mechanized in Computer Theorem
  1198 Provers (\emph{TP}) has (almost) a problem with traditional notations
  1199 (predicate calculus) for axioms, definitions, lemmas, theorems as a
  1200 computer program or script is not able to interpret every Greek or
  1201 Latin letter and every Greek, Latin or whatever calculations
  1202 symbol. Also if we would be able to handle these symbols we still have
  1203 a problem to interpret them at all. (Follow up \hbox{Example
  1204 \ref{eg:symbint1}})
  1205 
  1206 \vbox{
  1207   \begin{example}
  1208     \label{eg:symbint1}
  1209     \[
  1210       u\left[n\right] \ \ldots \ unitstep
  1211     \]
  1212     {\small\textit{
  1213       \noindent The unitstep is something we need to solve Signal
  1214 Processing problem classes. But in {{{\sisac}{}}} the rectangular
  1215 brackets have a different meaning. So we abuse them for our
  1216 requirements. We get something which is not defined, but usable. The
  1217 Result is syntax only without semantic.
  1218     }}
  1219   \end{example}
  1220 }
  1221 
  1222 In different problems, symbols and letters have different meanings and
  1223 ask for different ways to get through. (Follow up \hbox{Example
  1224 \ref{eg:symbint2}}) 
  1225 
  1226 \vbox{
  1227   \begin{example}
  1228     \label{eg:symbint2}
  1229     \[
  1230       \widehat{\ }\ \widehat{\ }\ \widehat{\ } \  \ldots \  exponent
  1231     \]
  1232     {\small\textit{
  1233     \noindent For using exponents the three \texttt{widehat} symbols
  1234 are required. The reason for that is due the development of
  1235 {{{\sisac}{}}} the single \texttt{widehat} and also the double were
  1236 already in use for different operations.
  1237     }}
  1238   \end{example}
  1239 }
  1240 
  1241 \paragraph{Also the output} can be a problem. We are familiar with a
  1242 specified notations and style taught in university but a computer
  1243 program has no knowledge of the form proved by a professor and the
  1244 machines themselves also have not yet the possibilities to print every
  1245 symbol (correct) Recent developments provide proofs in a human
  1246 readable format but according to the fact that there is no money for
  1247 good working formal editors yet, the style is one thing we have to
  1248 live with.
  1249 
  1250 \section{Problems rising out of the Development Environment}
  1251 
  1252 fehlermeldungen! TODO
  1253 
  1254 \section{Conclusion}\label{conclusion}
  1255 
  1256 TODO
  1257 
  1258 \bibliographystyle{alpha}
  1259 \bibliography{references}
  1260 
  1261 \end{document}