doc-src/isac/jrocnik/eJMT-paper/jrocnik_eJMT.tex
author Jan Rocnik <jan.rocnik@student.tugraz.at>
Thu, 13 Sep 2012 22:44:56 +0200
changeset 42512 2dd662758ae2
parent 42511 8c892624d349
child 42513 f7aa38509a95
child 42514 5e8f68f7510c
permissions -rwxr-xr-x
jrocnik: paper: nearly finished formulars
jan@42463
     1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
     2
% Electronic Journal of Mathematics and Technology (eJMT) %
jan@42463
     3
% style sheet for LaTeX.  Please do not modify sections   %
jan@42463
     4
% or commands marked 'eJMT'.                              %
jan@42463
     5
%                                                         %
jan@42463
     6
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
     7
%                                                         %
jan@42463
     8
% eJMT commands                                           %
jan@42463
     9
%                                                         %
jan@42463
    10
\documentclass[12pt,a4paper]{article}%                    %
jan@42463
    11
\usepackage{times}                                        %
jan@42463
    12
\usepackage{amsfonts,amsmath,amssymb}                     %
jan@42463
    13
\usepackage[a4paper]{geometry}                            %
jan@42463
    14
\usepackage{fancyhdr}                                     %
jan@42463
    15
\usepackage{color}                                        %
jan@42463
    16
\usepackage[pdftex]{hyperref} % see note below            %
jan@42463
    17
\usepackage{graphicx}%                                    %
jan@42463
    18
\hypersetup{                                              %
jan@42463
    19
    a4paper,                                              %
jan@42463
    20
    breaklinks                                            %
jan@42463
    21
}                                                         %
jan@42463
    22
%                                                         %
jan@42463
    23
\newtheorem{theorem}{Theorem}                             %
jan@42463
    24
\newtheorem{acknowledgement}[theorem]{Acknowledgement}    %
jan@42463
    25
\newtheorem{algorithm}[theorem]{Algorithm}                %
jan@42463
    26
\newtheorem{axiom}[theorem]{Axiom}                        %
jan@42463
    27
\newtheorem{case}[theorem]{Case}                          %
jan@42463
    28
\newtheorem{claim}[theorem]{Claim}                        %
jan@42463
    29
\newtheorem{conclusion}[theorem]{Conclusion}              %
jan@42463
    30
\newtheorem{condition}[theorem]{Condition}                %
jan@42463
    31
\newtheorem{conjecture}[theorem]{Conjecture}              %
jan@42463
    32
\newtheorem{corollary}[theorem]{Corollary}                %
jan@42463
    33
\newtheorem{criterion}[theorem]{Criterion}                %
jan@42463
    34
\newtheorem{definition}[theorem]{Definition}              %
jan@42463
    35
\newtheorem{example}[theorem]{Example}                    %
jan@42463
    36
\newtheorem{exercise}[theorem]{Exercise}                  %
jan@42463
    37
\newtheorem{lemma}[theorem]{Lemma}                        %
jan@42463
    38
\newtheorem{notation}[theorem]{Notation}                  %
jan@42463
    39
\newtheorem{problem}[theorem]{Problem}                    %
jan@42463
    40
\newtheorem{proposition}[theorem]{Proposition}            %
jan@42463
    41
\newtheorem{remark}[theorem]{Remark}                      %
jan@42463
    42
\newtheorem{solution}[theorem]{Solution}                  %
jan@42463
    43
\newtheorem{summary}[theorem]{Summary}                    %
jan@42463
    44
\newenvironment{proof}[1][Proof]{\noindent\textbf{#1.} }  %
jan@42463
    45
{\ \rule{0.5em}{0.5em}}                                   %
jan@42463
    46
%                                                         %
jan@42463
    47
% eJMT page dimensions                                    %
jan@42463
    48
%                                                         %
jan@42463
    49
\geometry{left=2cm,right=2cm,top=3.2cm,bottom=4cm}        %
jan@42463
    50
%                                                         %
jan@42463
    51
% eJMT header & footer                                    %
jan@42463
    52
%                                                         %
jan@42463
    53
\newcounter{ejmtFirstpage}                                %
jan@42463
    54
\setcounter{ejmtFirstpage}{1}                             %
jan@42463
    55
\pagestyle{empty}                                         %
jan@42463
    56
\setlength{\headheight}{14pt}                             %
jan@42463
    57
\geometry{left=2cm,right=2cm,top=3.2cm,bottom=4cm}        %
jan@42463
    58
\pagestyle{fancyplain}                                    %
jan@42463
    59
\fancyhf{}                                                %
jan@42463
    60
\fancyhead[c]{\small The Electronic Journal of Mathematics%
jan@42463
    61
\ and Technology, Volume 1, Number 1, ISSN 1933-2823}     %
jan@42463
    62
\cfoot{%                                                  %
jan@42463
    63
  \ifnum\value{ejmtFirstpage}=0%                          %
jan@42463
    64
    {\vtop to\hsize{\hrule\vskip .2cm\thepage}}%          %
jan@42463
    65
  \else\setcounter{ejmtFirstpage}{0}\fi%                  %
jan@42463
    66
}                                                         %
jan@42463
    67
%                                                         %
jan@42463
    68
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
    69
%
jan@42463
    70
% Please place your own definitions here
jan@42463
    71
%
jan@42463
    72
\def\isac{${\cal I}\mkern-2mu{\cal S}\mkern-5mu{\cal AC}$}
jan@42463
    73
\def\sisac{\footnotesize${\cal I}\mkern-2mu{\cal S}\mkern-5mu{\cal AC}$}
jan@42463
    74
jan@42463
    75
\usepackage{color}
jan@42463
    76
\definecolor{lgray}{RGB}{238,238,238}
jan@42463
    77
jan@42463
    78
%
jan@42463
    79
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
    80
%                                                         %
jan@42463
    81
% How to use hyperref                                     %
jan@42463
    82
% -------------------                                     %
jan@42463
    83
%                                                         %
jan@42463
    84
% Probably the only way you will need to use the hyperref %
jan@42463
    85
% package is as follows.  To make some text, say          %
jan@42463
    86
% "My Text Link", into a link to the URL                  %
jan@42463
    87
% http://something.somewhere.com/mystuff, use             %
jan@42463
    88
%                                                         %
jan@42463
    89
% \href{http://something.somewhere.com/mystuff}{My Text Link}
jan@42463
    90
%                                                         %
jan@42463
    91
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
    92
%
jan@42463
    93
\begin{document}
jan@42463
    94
%
jan@42463
    95
% document title
jan@42463
    96
%
neuper@42464
    97
\title{Trials with TP-based Programming
neuper@42464
    98
\\
neuper@42464
    99
for Interactive Course Material}%
jan@42463
   100
%
jan@42463
   101
% Single author.  Please supply at least your name,
jan@42463
   102
% email address, and affiliation here.
jan@42463
   103
%
jan@42463
   104
\author{\begin{tabular}{c}
jan@42463
   105
\textit{Jan Ro\v{c}nik} \\
jan@42463
   106
jan.rocnik@student.tugraz.at \\
jan@42463
   107
IST, SPSC\\
neuper@42464
   108
Graz University of Technologie\\
jan@42463
   109
Austria\end{tabular}
jan@42463
   110
}%
jan@42463
   111
%
jan@42463
   112
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
   113
%                                                         %
jan@42463
   114
% eJMT commands - do not change these                     %
jan@42463
   115
%                                                         %
jan@42463
   116
\date{}                                                   %
jan@42463
   117
\maketitle                                                %
jan@42463
   118
%                                                         %
jan@42463
   119
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
   120
%
jan@42463
   121
% abstract
jan@42463
   122
%
jan@42463
   123
\begin{abstract}
jan@42463
   124
jan@42463
   125
Traditional course material in engineering disciplines lacks an
jan@42463
   126
important component, interactive support for step-wise problem
neuper@42464
   127
solving. Theorem-Proving (TP) technology is appropriate for one part
jan@42463
   128
of such support, in checking user-input. For the other part of such
jan@42463
   129
support, guiding the learner towards a solution, another kind of
jan@42463
   130
technology is required. %TODO ... connect to prototype ...
jan@42463
   131
jan@42511
   132
Both kinds of support can be achieved by so-called
neuper@42504
   133
Lucas-Interpretation which combines deduction and computation and, for
neuper@42504
   134
the latter, uses a novel kind of programming language. This language
neuper@42504
   135
is based on (Computer) Theorem Proving (TP), thus called a ``TP-based
neuper@42504
   136
programming language''.
jan@42463
   137
neuper@42504
   138
This paper is the experience report of the first ``application
neuper@42507
   139
programmer'' using this language for creating exercises in step-wise
neuper@42507
   140
problem solving for an advanced lab in Signal Processing. The tasks
neuper@42507
   141
involved in TP-based programming are described together with the
neuper@42507
   142
experience gained from a prototype of the programming language and of
neuper@42507
   143
it's interpreter.
neuper@42504
   144
neuper@42504
   145
The report concludes with a positive proof of concept, states
jan@42512
   146
insufficiency usability of the prototype and captures the requirements
neuper@42504
   147
for further development of both, the programming language and the
neuper@42504
   148
interpreter.
jan@42463
   149
%
jan@42463
   150
\end{abstract}%
jan@42463
   151
%
jan@42463
   152
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
   153
%                                                         %
jan@42463
   154
% eJMT command                                            %
jan@42463
   155
%                                                         %
jan@42463
   156
\thispagestyle{fancy}                                     %
jan@42463
   157
%                                                         %
jan@42463
   158
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
jan@42463
   159
%
jan@42463
   160
% Please use the following to indicate sections, subsections,
jan@42463
   161
% etc.  Please also use \subsubsection{...}, \paragraph{...}
jan@42463
   162
% and \subparagraph{...} as necessary.
jan@42463
   163
%
jan@42463
   164
neuper@42464
   165
\section{Introduction}\label{intro}
jan@42463
   166
jan@42466
   167
% \paragraph{Didactics of mathematics} 
jan@42466
   168
%WN: wenn man in einem high-quality paper von 'didactics' spricht, 
jan@42466
   169
%WN muss man am state-of-the-art ankn"upfen -- siehe
jan@42466
   170
%WN W.Neuper, On the Emergence of TP-based Educational Math Assistants
neuper@42464
   171
% faces a specific issue, a gap
neuper@42464
   172
% between (1) introduction of math concepts and skills and (2)
neuper@42464
   173
% application of these concepts and skills, which usually are separated
neuper@42464
   174
% into different units in curricula (for good reasons). For instance,
neuper@42464
   175
% (1) teaching partial fraction decomposition is separated from (2)
neuper@42464
   176
% application for inverse Z-transform in signal processing.
neuper@42464
   177
% 
neuper@42464
   178
% \par This gap is an obstacle for applying math as an fundamental
neuper@42464
   179
% thinking technology in engineering: In (1) motivation is lacking
neuper@42464
   180
% because the question ``What is this stuff good for?'' cannot be
neuper@42464
   181
% treated sufficiently, and in (2) the ``stuff'' is not available to
neuper@42464
   182
% students in higher semesters as widespread experience shows.
neuper@42464
   183
% 
neuper@42464
   184
% \paragraph{Motivation} taken by this didactic issue on the one hand,
neuper@42464
   185
% and ongoing research and development on a novel kind of educational
neuper@42464
   186
% mathematics assistant at Graz University of
neuper@42464
   187
% Technology~\footnote{http://www.ist.tugraz.at/isac/} promising to
neuper@42464
   188
% scope with this issue on the other hand, several institutes are
neuper@42464
   189
% planning to join their expertise: the Institute for Information
neuper@42464
   190
% Systems and Computer Media (IICM), the Institute for Software
neuper@42464
   191
% Technology (IST), the Institutes for Mathematics, the Institute for
neuper@42464
   192
% Signal Processing and Speech Communication (SPSC), the Institute for
neuper@42464
   193
% Structural Analysis and the Institute of Electrical Measurement and
neuper@42464
   194
% Measurement Signal Processing.
jan@42466
   195
%WN diese Information ist f"ur das Paper zu spezielle, zu aktuell 
jan@42466
   196
%WN und damit zu verg"anglich.
neuper@42464
   197
% \par This thesis is the first attempt to tackle the above mentioned
neuper@42464
   198
% issue, it focuses on Telematics, because these specific studies focus
neuper@42464
   199
% on mathematics in \emph{STEOP}, the introductory orientation phase in
neuper@42464
   200
% Austria. \emph{STEOP} is considered an opportunity to investigate the
neuper@42464
   201
% impact of {\sisac}'s prototype on the issue and others.
neuper@42464
   202
% 
jan@42466
   203
jan@42502
   204
Traditional course material in engineering disciplines lacks an
neuper@42464
   205
important component, interactive support for step-wise problem
neuper@42464
   206
solving. Theorem-Proving (TP) technology can provide such support by
neuper@42464
   207
specific services. An important part of such services is called
neuper@42464
   208
``next-step-guidance'', generated by a specific kind of ``TP-based
neuper@42464
   209
programming language''. In the
neuper@42464
   210
{\sisac}-project~\footnote{http://www.ist.tugraz.at/projects/isac/} such
neuper@42464
   211
a language is prototyped in line with~\cite{plmms10} and built upon
jan@42511
   212
the theorem prover.
neuper@42464
   213
Isabelle~\cite{Nipkow-Paulson-Wenzel:2002}\footnote{http://isabelle.in.tum.de/}.
neuper@42464
   214
The TP services are coordinated by a specific interpreter for the
neuper@42464
   215
programming language, called
neuper@42464
   216
Lucas-Interpreter~\cite{wn:lucas-interp-12}. The language and the
neuper@42464
   217
interpreter will be briefly re-introduced in order to make the paper
neuper@42464
   218
self-contained.
jan@42463
   219
neuper@42504
   220
The main part of the paper is an account of first experiences
neuper@42464
   221
with programming in this TP-based language. The experience was gained
neuper@42464
   222
in a case study by the author. The author was considered an ideal
neuper@42464
   223
candidate for this study for the following reasons: as a student in
neuper@42464
   224
Telematics (computer science with focus on Signal Processing) he had
neuper@42464
   225
general knowledge in programming as well as specific domain knowledge
neuper@42504
   226
in Signal Processing; and he was {\em not} involved in the development of
jan@42511
   227
{\sisac}'s programming language and interpreter, thus a novice to the
neuper@42464
   228
language.
jan@42463
   229
neuper@42504
   230
The goal of the case study was (1) some TP-based programs for
jan@42511
   231
interactive course material for a specific ``Advanced Signal
neuper@42464
   232
Processing Lab'' in a higher semester, (2) respective program
neuper@42464
   233
development with as little advice from the {\sisac}-team and (3) records
neuper@42464
   234
and comments for the main steps of development in an Isabelle theory;
neuper@42464
   235
this theory should provide guidelines for future programmers. An
neuper@42464
   236
excerpt from this theory is the main part of this paper.
jan@42466
   237
\par
jan@42466
   238
The paper will use the problem in Fig.\ref{fig-interactive} as a
jan@42463
   239
running example:
jan@42463
   240
\begin{figure} [htb]
jan@42463
   241
\begin{center}
neuper@42468
   242
\includegraphics[width=140mm]{fig/isac-Ztrans-math-3}
neuper@42468
   243
%\includegraphics[width=140mm]{fig/isac-Ztrans-math}
jan@42512
   244
\caption{Step-wise problem solving guided by the TP-based program
jan@42512
   245
\label{fig-interactive}}
jan@42463
   246
\end{center}
jan@42463
   247
\end{figure}
jan@42466
   248
jan@42502
   249
The problem is from the domain of Signal Processing and requests to
jan@42512
   250
determine the inverse ${\cal z}$-transform for a given term.
jan@42511
   251
Fig.\ref{fig-interactive}
neuper@42464
   252
also shows the beginning of the interactive construction of a solution
neuper@42464
   253
for the problem. This construction is done in the right window named
neuper@42464
   254
``Worksheet''.
jan@42466
   255
\par
neuper@42464
   256
User-interaction on the Worksheet is {\em checked} and {\em guided} by
neuper@42464
   257
TP services:
neuper@42464
   258
\begin{enumerate}
neuper@42464
   259
\item Formulas input by the user are {\em checked} by TP: such a
neuper@42464
   260
formula establishes a proof situation --- the prover has to derive the
neuper@42464
   261
formula from the logical context. The context is built up from the
neuper@42464
   262
formal specification of the problem (here hidden from the user) by the
neuper@42464
   263
Lucas-Interpreter.
neuper@42464
   264
\item If the user gets stuck, the program developed below in this
neuper@42504
   265
paper ``knows the next step'' and Lucas-Interpretation provides services
neuper@42504
   266
featuring so-called ``next-step-guidance''; this is out of scope of this
neuper@42464
   267
paper and can be studied in~\cite{gdaroczy-EP-13}.
neuper@42464
   268
\end{enumerate} It should be noted that the programmer using the
neuper@42464
   269
TP-based language is not concerned with interaction at all; we will
neuper@42464
   270
see that the program contains neither input-statements nor
neuper@42504
   271
output-statements. Rather, interaction is handled by the interpreter
jan@42511
   272
of the language.
neuper@42504
   273
neuper@42504
   274
So there is a clear separation of concerns: Dialogues are adapted by
neuper@42504
   275
dialogue authors (in Java-based tools), using TP services provided by
neuper@42504
   276
Lucas-Interpretation. The latter acts on programs developed by
neuper@42504
   277
mathematics-authors (in Isabelle/ML); their task is concern of this
neuper@42464
   278
paper.
jan@42466
   279
jan@42511
   280
\paragraph{The paper is structured} as follows: The introduction
neuper@42464
   281
\S\ref{intro} is followed by a brief re-introduction of the TP-based
neuper@42464
   282
programming language in \S\ref{PL}, which extends the executable
neuper@42464
   283
fragment of Isabelle's language (\S\ref{PL-isab}) by tactics which
neuper@42464
   284
play a specific role in Lucas-Interpretation and in providing the TP
neuper@42504
   285
services (\S\ref{PL-tacs}). The main part \S\ref{trial} describes
neuper@42464
   286
the main steps in developing the program for the running example:
neuper@42464
   287
prepare domain knowledge, implement the formal specification of the
neuper@42504
   288
problem, prepare the environment for the interpreter, implement the
neuper@42504
   289
program in \S\ref{isabisac} to \S\ref{progr} respectively. 
jan@42511
   290
The work-flow of programming, debugging and testing is
neuper@42464
   291
described in \S\ref{workflow}. The conclusion \S\ref{conclusion} will
neuper@42464
   292
give directions identified for future development. 
neuper@42464
   293
jan@42463
   294
jan@42463
   295
\section{\isac's Prototype for a Programming Language}\label{PL} 
neuper@42504
   296
The prototype of the language and of the Lucas-Interpreter are briefly
neuper@42504
   297
described from the point of view of a programmer. The language extends
neuper@42504
   298
the executable fragment in the language of the theorem prover
neuper@42504
   299
Isabelle~\cite{Nipkow-Paulson-Wenzel:2002}\footnote{http://isabelle.in.tum.de/}.
jan@42463
   300
jan@42463
   301
\subsection{The Executable Fragment of Isabelle's Language}\label{PL-isab}
jan@42463
   302
The executable fragment consists of data-type and function
jan@42463
   303
definitions.  It's usability even suggests that fragment for
jan@42463
   304
introductory courses \cite{nipkow-prog-prove}. HOL is a typed logic
jan@42463
   305
whose type system resembles that of functional programming
jan@42463
   306
languages. Thus there are
jan@42463
   307
\begin{description}
jan@42463
   308
\item[base types,] in particular \textit{bool}, the type of truth
jan@42463
   309
values, \textit{nat}, \textit{int}, \textit{complex}, and the types of
jan@42463
   310
natural, integer and complex numbers respectively in mathematics.
jan@42463
   311
\item[type constructors] allow to define arbitrary types, from
jan@42463
   312
\textit{set}, \textit{list} to advanced data-structures like
jan@42463
   313
\textit{trees}, red-black-trees etc.
jan@42463
   314
\item[function types,] denoted by $\Rightarrow$.
jan@42463
   315
\item[type variables,] denoted by $^\prime a, ^\prime b$ etc, provide
jan@42463
   316
type polymorphism. Isabelle automatically computes the type of each
jan@42463
   317
variable in a term by use of Hindley-Milner type inference
jan@42463
   318
\cite{pl:hind97,Milner-78}.
jan@42463
   319
\end{description}
jan@42463
   320
jan@42463
   321
\textbf{Terms} are formed as in functional programming by applying
jan@42463
   322
functions to arguments. If $f$ is a function of type
jan@42463
   323
$\tau_1\Rightarrow \tau_2$ and $t$ is a term of type $\tau_1$ then
jan@42463
   324
$f\;t$ is a term of type~$\tau_2$. $t\;::\;\tau$ means that term $t$
jan@42463
   325
has type $\tau$. There are many predefined infix symbols like $+$ and
jan@42463
   326
$\leq$ most of which are overloaded for various types.
jan@42463
   327
jan@42463
   328
HOL also supports some basic constructs from functional programming:
jan@42512
   329
{\footnotesize\it\label{isabelle-stmts}
jan@42463
   330
\begin{tabbing} 123\=\kill
jan@42463
   331
\>$( \; {\tt if} \; b \; {\tt then} \; t_1 \; {\tt else} \; t_2 \;)$\\
jan@42463
   332
\>$( \; {\tt let} \; x=t \; {\tt in} \; u \; )$\\
jan@42463
   333
\>$( \; {\tt case} \; t \; {\tt of} \; {\it pat}_1
jan@42463
   334
  \Rightarrow t_1 \; |\dots| \; {\it pat}_n\Rightarrow t_n \; )$
jan@42512
   335
\end{tabbing}}
neuper@42482
   336
\noindent The running example's program uses some of these elements
neuper@42482
   337
(marked by {\tt tt-font} on p.\pageref{s:impl}): for instance {\tt
neuper@42482
   338
let}\dots{\tt in} in lines {\rm 02} \dots {\rm 13}. In fact, the whole program
neuper@42482
   339
is an Isabelle term with specific function constants like {\tt
neuper@42482
   340
program}, {\tt Take}, {\tt Rewrite}, {\tt Subproblem} and {\tt
neuper@42482
   341
Rewrite\_Set} in lines {\rm 01, 03. 04, 07, 10} and {\rm 11, 12}
neuper@42482
   342
respectively.
jan@42463
   343
jan@42463
   344
% Terms may also contain $\lambda$-abstractions. For example, $\lambda
jan@42463
   345
% x. \; x$ is the identity function.
jan@42463
   346
neuper@42467
   347
%JR warum auskommentiert? WN2...
neuper@42467
   348
%WN2 weil ein Punkt wie dieser in weiteren Zusammenh"angen innerhalb
neuper@42467
   349
%WN2 des Papers auftauchen m"usste; nachdem ich einen solchen
neuper@42467
   350
%WN2 Zusammenhang _noch_ nicht sehe, habe ich den Punkt _noch_ nicht
neuper@42467
   351
%WN2 gel"oscht.
neuper@42467
   352
%WN2 Wenn der Punkt nicht weiter gebraucht wird, nimmt er nur wertvollen
neuper@42467
   353
%WN2 Platz f"ur Anderes weg.
jan@42466
   354
neuper@42464
   355
\textbf{Formulae} are terms of type \textit{bool}. There are the basic
jan@42463
   356
constants \textit{True} and \textit{False} and the usual logical
jan@42463
   357
connectives (in decreasing order of precedence): $\neg, \land, \lor,
jan@42463
   358
\rightarrow$.
jan@42463
   359
neuper@42464
   360
\textbf{Equality} is available in the form of the infix function $=$
neuper@42464
   361
of type $a \Rightarrow a \Rightarrow {\it bool}$. It also works for
neuper@42464
   362
formulas, where it means ``if and only if''.
jan@42463
   363
jan@42463
   364
\textbf{Quantifiers} are written $\forall x. \; P$ and $\exists x. \;
jan@42463
   365
P$.  Quantifiers lead to non-executable functions, so functions do not
jan@42463
   366
always correspond to programs, for instance, if comprising \\$(
jan@42463
   367
\;{\it if} \; \exists x.\;P \; {\it then} \; e_1 \; {\it else} \; e_2
jan@42463
   368
\;)$.
jan@42463
   369
jan@42463
   370
\subsection{\isac's Tactics for Lucas-Interpretation}\label{PL-tacs}
jan@42463
   371
The prototype extends Isabelle's language by specific statements
neuper@42464
   372
called tactics~\footnote{{\sisac}'s tactics are different from
jan@42463
   373
Isabelle's tactics: the former concern steps in a calculation, the
jan@42511
   374
latter concern proofs.}  and tactics. For the programmer these
jan@42463
   375
statements are functions with the following signatures:
jan@42463
   376
jan@42463
   377
\begin{description}
jan@42463
   378
\item[Rewrite:] ${\it theorem}\Rightarrow{\it term}\Rightarrow{\it
jan@42463
   379
term} * {\it term}\;{\it list}$:
jan@42511
   380
this tactic applies {\it theorem} to a {\it term} yielding a {\it
jan@42463
   381
term} and a {\it term list}, the list are assumptions generated by
jan@42463
   382
conditional rewriting. For instance, the {\it theorem}
jan@42463
   383
$b\not=0\land c\not=0\Rightarrow\frac{a\cdot c}{b\cdot c}=\frac{a}{b}$
jan@42463
   384
applied to the {\it term} $\frac{2\cdot x}{3\cdot x}$ yields
jan@42463
   385
$(\frac{2}{3}, [x\not=0])$.
jan@42463
   386
jan@42463
   387
\item[Rewrite\_Set:] ${\it ruleset}\Rightarrow{\it
jan@42463
   388
term}\Rightarrow{\it term} * {\it term}\;{\it list}$:
jan@42511
   389
this tactic applies {\it ruleset} to a {\it term}; {\it ruleset} is
jan@42463
   390
a confluent and terminating term rewrite system, in general. If
jan@42463
   391
none of the rules ({\it theorem}s) is applicable on interpretation
jan@42463
   392
of this tactic, an exception is thrown.
jan@42463
   393
jan@42463
   394
% \item[Rewrite\_Inst:] ${\it substitution}\Rightarrow{\it
jan@42463
   395
% theorem}\Rightarrow{\it term}\Rightarrow{\it term} * {\it term}\;{\it
jan@42463
   396
% list}$:
jan@42463
   397
% 
jan@42463
   398
% \item[Rewrite\_Set\_Inst:] ${\it substitution}\Rightarrow{\it
jan@42463
   399
% ruleset}\Rightarrow{\it term}\Rightarrow{\it term} * {\it term}\;{\it
jan@42463
   400
% list}$:
jan@42463
   401
neuper@42504
   402
%SPACEvvv
jan@42463
   403
\item[Substitute:] ${\it substitution}\Rightarrow{\it
neuper@42482
   404
term}\Rightarrow{\it term}$: allows to access sub-terms.
neuper@42504
   405
%SPACE^^^
jan@42463
   406
jan@42463
   407
\item[Take:] ${\it term}\Rightarrow{\it term}$:
jan@42463
   408
this tactic has no effect in the program; but it creates a side-effect
jan@42463
   409
by Lucas-Interpretation (see below) and writes {\it term} to the
jan@42463
   410
Worksheet.
jan@42463
   411
jan@42463
   412
\item[Subproblem:] ${\it theory} * {\it specification} * {\it
jan@42463
   413
method}\Rightarrow{\it argument}\;{\it list}\Rightarrow{\it term}$:
neuper@42482
   414
this tactic is a generalisation of a function call: it takes an
neuper@42482
   415
\textit{argument list} as usual, and additionally a triple consisting
neuper@42482
   416
of an Isabelle \textit{theory}, an implicit \textit{specification} of the
neuper@42482
   417
program and a \textit{method} containing data for Lucas-Interpretation,
neuper@42482
   418
last not least a program (as an explicit specification)~\footnote{In
neuper@42482
   419
interactive tutoring these three items can be determined explicitly
neuper@42482
   420
by the user.}.
jan@42463
   421
\end{description}
jan@42463
   422
The tactics play a specific role in
jan@42463
   423
Lucas-Interpretation~\cite{wn:lucas-interp-12}: they are treated as
neuper@42482
   424
break-points where, as a side-effect, a line is added to a calculation
neuper@42483
   425
as a protocol for proceeding towards a solution in step-wise problem
neuper@42483
   426
solving. At the same points Lucas-Interpretation serves interactive
neuper@42504
   427
tutoring and hands over control to the user. The user is free to
neuper@42483
   428
investigate underlying knowledge, applicable theorems, etc.  And the
neuper@42483
   429
user can proceed constructing a solution by input of a tactic to be
neuper@42483
   430
applied or by input of a formula; in the latter case the
jan@42463
   431
Lucas-Interpreter has built up a logical context (initialised with the
jan@42463
   432
precondition of the formal specification) such that Isabelle can
jan@42463
   433
derive the formula from this context --- or give feedback, that no
jan@42463
   434
derivation can be found.
jan@42463
   435
jan@42511
   436
\subsection{Tactics as Control Flow Statements}
jan@42463
   437
The flow of control in a program can be determined by {\tt if then else}
jan@42463
   438
and {\tt case of} as mentioned on p.\pageref{isabelle-stmts} and also
jan@42511
   439
by additional tactics:
jan@42463
   440
\begin{description}
jan@42463
   441
\item[Repeat:] ${\it tactic}\Rightarrow{\it term}\Rightarrow{\it
jan@42463
   442
term}$: iterates over tactics which take a {\it term} as argument as
neuper@42482
   443
long as a tactic is applicable (for instance, {\tt Rewrite\_Set} might
jan@42463
   444
not be applicable).
jan@42463
   445
jan@42463
   446
\item[Try:] ${\it tactic}\Rightarrow{\it term}\Rightarrow{\it term}$:
jan@42463
   447
if {\it tactic} is applicable, then it is applied to {\it term},
neuper@42483
   448
otherwise {\it term} is passed on without changes.
jan@42463
   449
jan@42463
   450
\item[Or:] ${\it tactic}\Rightarrow{\it tactic}\Rightarrow{\it
neuper@42483
   451
term}\Rightarrow{\it term}$: If the first {\it tactic} is applicable,
neuper@42483
   452
it is applied to the first {\it term} yielding another {\it term},
neuper@42483
   453
otherwise the second {\it tactic} is applied; if none is applicable an
neuper@42483
   454
exception is raised.
jan@42463
   455
jan@42463
   456
\item[@@:] ${\it tactic}\Rightarrow{\it tactic}\Rightarrow{\it
neuper@42483
   457
term}\Rightarrow{\it term}$: applies the first {\it tactic} to the
neuper@42483
   458
first {\it term} yielding an intermediate term (not appearing in the
neuper@42483
   459
signature) to which the second {\it tactic} is applied.
jan@42463
   460
jan@42463
   461
\item[While:] ${\it term::bool}\Rightarrow{\it tactic}\Rightarrow{\it
neuper@42483
   462
term}\Rightarrow{\it term}$: if the first {\it term} is true, then the
neuper@42483
   463
{\it tactic} is applied to the first {\it term} yielding an
neuper@42483
   464
intermediate term (not appearing in the signature); the intermediate
neuper@42483
   465
term is added to the environment the first {\it term} is evaluated in
jan@42511
   466
etc. as long as the first {\it term} is true.
jan@42463
   467
\end{description}
jan@42511
   468
The tactics are not treated as break-points by Lucas-Interpretation
neuper@42504
   469
and thus do neither contribute to the calculation nor to interaction.
jan@42463
   470
neuper@42498
   471
\section{Concepts and Tasks in TP-based Programming}\label{trial}
neuper@42498
   472
%\section{Development of a Program on Trial}
neuper@42498
   473
neuper@42498
   474
This section presents all the concepts involved in TP-based
neuper@42498
   475
programming and all the tasks to be accomplished by programmers. The
neuper@42504
   476
presentation uses the running example from
neuper@42498
   477
Fig.\ref{fig-interactive} on p.\pageref{fig-interactive}.
jan@42466
   478
jan@42466
   479
\subsection{Mechanization of Math --- Domain Engineering}\label{isabisac}
jan@42466
   480
neuper@42467
   481
%WN was Fachleute unter obigem Titel interessiert findet sich
jan@42466
   482
%WN unterhalb des auskommentierten Textes.
jan@42466
   483
jan@42466
   484
%WN der Text unten spricht Benutzer-Aspekte anund ist nicht speziell
jan@42466
   485
%WN auf Computer-Mathematiker fokussiert.
neuper@42464
   486
% \paragraph{As mentioned in the introduction,} a prototype of an
neuper@42464
   487
% educational math assistant called
neuper@42464
   488
% {{\sisac}}\footnote{{{\sisac}}=\textbf{Isa}belle for
neuper@42464
   489
% \textbf{C}alculations, see http://www.ist.tugraz.at/isac/.} bridges
neuper@42464
   490
% the gap between (1) introducation and (2) application of mathematics:
neuper@42464
   491
% {{\sisac}} is based on Computer Theorem Proving (TP), a technology which
neuper@42464
   492
% requires each fact and each action justified by formal logic, so
neuper@42464
   493
% {{{\sisac}{}}} makes justifications transparent to students in
neuper@42464
   494
% interactive step-wise problem solving. By that way {{\sisac}} already
neuper@42464
   495
% can serve both:
neuper@42464
   496
% \begin{enumerate}
neuper@42464
   497
%   \item Introduction of math stuff (in e.g. partial fraction
neuper@42464
   498
% decomposition) by stepwise explaining and exercising respective
neuper@42464
   499
% symbolic calculations with ``next step guidance (NSG)'' and rigorously
neuper@42464
   500
% checking steps freely input by students --- this also in context with
neuper@42464
   501
% advanced applications (where the stuff to be taught in higher
neuper@42464
   502
% semesters can be skimmed through by NSG), and
neuper@42464
   503
%   \item Application of math stuff in advanced engineering courses
neuper@42464
   504
% (e.g. problems to be solved by inverse Z-transform in a Signal
neuper@42464
   505
% Processing Lab) and now without much ado about basic math techniques
neuper@42464
   506
% (like partial fraction decomposition): ``next step guidance'' supports
neuper@42464
   507
% students in independently (re-)adopting such techniques.
neuper@42464
   508
% \end{enumerate} 
neuper@42464
   509
% Before the question is answers, how {{\sisac}}
neuper@42464
   510
% accomplishes this task from a technical point of view, some remarks on
neuper@42464
   511
% the state-of-the-art is given, therefor follow up Section~\ref{emas}.
neuper@42464
   512
% 
neuper@42464
   513
% \subsection{Educational Mathematics Assistants (EMAs)}\label{emas}
neuper@42464
   514
% 
jan@42466
   515
% \paragraph{Educational software in mathematics} is, if at all, based
jan@42466
   516
% on Computer Algebra Systems (CAS, for instance), Dynamic Geometry
jan@42466
   517
% Systems (DGS, for instance \footnote{GeoGebra http://www.geogebra.org}
jan@42466
   518
% \footnote{Cinderella http://www.cinderella.de/}\footnote{GCLC
jan@42466
   519
% http://poincare.matf.bg.ac.rs/~janicic/gclc/}) or spread-sheets. These
jan@42466
   520
% base technologies are used to program math lessons and sometimes even
jan@42466
   521
% exercises. The latter are cumbersome: the steps towards a solution of
jan@42466
   522
% such an interactive exercise need to be provided with feedback, where
jan@42466
   523
% at each step a wide variety of possible input has to be foreseen by
jan@42466
   524
% the programmer - so such interactive exercises either require high
neuper@42464
   525
% development efforts or the exercises constrain possible inputs.
neuper@42464
   526
% 
jan@42466
   527
% \subparagraph{A new generation} of educational math assistants (EMAs)
jan@42466
   528
% is emerging presently, which is based on Theorem Proving (TP). TP, for
jan@42466
   529
% instance Isabelle and Coq, is a technology which requires each fact
jan@42466
   530
% and each action justified by formal logic. Pushed by demands for
jan@42466
   531
% \textit{proven} correctness of safety-critical software TP advances
jan@42466
   532
% into software engineering; from these advancements computer
jan@42466
   533
% mathematics benefits in general, and math education in particular. Two
neuper@42464
   534
% features of TP are immediately beneficial for learning:
neuper@42464
   535
% 
jan@42466
   536
% \paragraph{TP have knowledge in human readable format,} that is in
jan@42466
   537
% standard predicate calculus. TP following the LCF-tradition have that
jan@42466
   538
% knowledge down to the basic definitions of set, equality,
jan@42466
   539
% etc~\footnote{http://isabelle.in.tum.de/dist/library/HOL/HOL.html};
jan@42466
   540
% following the typical deductive development of math, natural numbers
jan@42466
   541
% are defined and their properties
jan@42466
   542
% proven~\footnote{http://isabelle.in.tum.de/dist/library/HOL/Number\_Theory/Primes.html},
jan@42466
   543
% etc. Present knowledge mechanized in TP exceeds high-school
jan@42466
   544
% mathematics by far, however by knowledge required in software
neuper@42464
   545
% technology, and not in other engineering sciences.
neuper@42464
   546
% 
jan@42466
   547
% \paragraph{TP can model the whole problem solving process} in
jan@42466
   548
% mathematical problem solving {\em within} a coherent logical
jan@42466
   549
% framework. This is already being done by three projects, by
neuper@42464
   550
% Ralph-Johan Back, by ActiveMath and by Carnegie Mellon Tutor.
neuper@42464
   551
% \par
jan@42466
   552
% Having the whole problem solving process within a logical coherent
jan@42466
   553
% system, such a design guarantees correctness of intermediate steps and
jan@42466
   554
% of the result (which seems essential for math software); and the
jan@42466
   555
% second advantage is that TP provides a wealth of theories which can be
jan@42466
   556
% exploited for mechanizing other features essential for educational
neuper@42464
   557
% software.
neuper@42464
   558
% 
neuper@42464
   559
% \subsubsection{Generation of User Guidance in EMAs}\label{user-guid}
neuper@42464
   560
% 
jan@42466
   561
% One essential feature for educational software is feedback to user
neuper@42464
   562
% input and assistance in coming to a solution.
neuper@42464
   563
% 
jan@42466
   564
% \paragraph{Checking user input} by ATP during stepwise problem solving
jan@42466
   565
% is being accomplished by the three projects mentioned above
jan@42466
   566
% exclusively. They model the whole problem solving process as mentioned
jan@42466
   567
% above, so all what happens between formalized assumptions (or formal
jan@42466
   568
% specification) and goal (or fulfilled postcondition) can be
jan@42466
   569
% mechanized. Such mechanization promises to greatly extend the scope of
neuper@42464
   570
% educational software in stepwise problem solving.
neuper@42464
   571
% 
jan@42466
   572
% \paragraph{NSG (Next step guidance)} comprises the system's ability to
jan@42466
   573
% propose a next step; this is a challenge for TP: either a radical
jan@42466
   574
% restriction of the search space by restriction to very specific
jan@42466
   575
% problem classes is required, or much care and effort is required in
jan@42466
   576
% designing possible variants in the process of problem solving
neuper@42464
   577
% \cite{proof-strategies-11}.
neuper@42464
   578
% \par
jan@42466
   579
% Another approach is restricted to problem solving in engineering
jan@42466
   580
% domains, where a problem is specified by input, precondition, output
jan@42466
   581
% and postcondition, and where the postcondition is proven by ATP behind
jan@42466
   582
% the scenes: Here the possible variants in the process of problem
jan@42466
   583
% solving are provided with feedback {\em automatically}, if the problem
jan@42466
   584
% is described in a TP-based programing language: \cite{plmms10} the
jan@42466
   585
% programmer only describes the math algorithm without caring about
jan@42466
   586
% interaction (the respective program is functional and even has no
jan@42466
   587
% input or output statements!); interaction is generated as a
jan@42466
   588
% side-effect by the interpreter --- an efficient separation of concern
jan@42466
   589
% between math programmers and dialog designers promising application
neuper@42464
   590
% all over engineering disciplines.
neuper@42464
   591
% 
neuper@42464
   592
% 
neuper@42464
   593
% \subsubsection{Math Authoring in Isabelle/ISAC\label{math-auth}}
jan@42466
   594
% Authoring new mathematics knowledge in {{\sisac}} can be compared with
jan@42466
   595
% ``application programing'' of engineering problems; most of such
jan@42466
   596
% programing uses CAS-based programing languages (CAS = Computer Algebra
neuper@42464
   597
% Systems; e.g. Mathematica's or Maple's programing language).
neuper@42464
   598
% 
jan@42466
   599
% \paragraph{A novel type of TP-based language} is used by {{\sisac}{}}
jan@42466
   600
% \cite{plmms10} for describing how to construct a solution to an
jan@42466
   601
% engineering problem and for calling equation solvers, integration,
jan@42466
   602
% etc~\footnote{Implementation of CAS-like functionality in TP is not
jan@42466
   603
% primarily concerned with efficiency, but with a didactic question:
jan@42466
   604
% What to decide for: for high-brow algorithms at the state-of-the-art
jan@42466
   605
% or for elementary algorithms comprehensible for students?} within TP;
jan@42466
   606
% TP can ensure ``systems that never make a mistake'' \cite{casproto} -
neuper@42464
   607
% are impossible for CAS which have no logics underlying.
neuper@42464
   608
% 
jan@42466
   609
% \subparagraph{Authoring is perfect} by writing such TP based programs;
jan@42466
   610
% the application programmer is not concerned with interaction or with
jan@42466
   611
% user guidance: this is concern of a novel kind of program interpreter
jan@42466
   612
% called Lucas-Interpreter. This interpreter hands over control to a
jan@42466
   613
% dialog component at each step of calculation (like a debugger at
jan@42466
   614
% breakpoints) and calls automated TP to check user input following
neuper@42464
   615
% personalized strategies according to a feedback module.
neuper@42464
   616
% \par
jan@42466
   617
% However ``application programing with TP'' is not done with writing a
jan@42466
   618
% program: according to the principles of TP, each step must be
jan@42466
   619
% justified. Such justifications are given by theorems. So all steps
jan@42466
   620
% must be related to some theorem, if there is no such theorem it must
jan@42466
   621
% be added to the existing knowledge, which is organized in so-called
jan@42466
   622
% \textbf{theories} in Isabelle. A theorem must be proven; fortunately
jan@42466
   623
% Isabelle comprises a mechanism (called ``axiomatization''), which
jan@42466
   624
% allows to omit proofs. Such a theorem is shown in
neuper@42464
   625
% Example~\ref{eg:neuper1}.
jan@42466
   626
neuper@42498
   627
The running example requires to determine the inverse $\cal
jan@42466
   628
Z$-transform for a class of functions. The domain of Signal Processing
jan@42466
   629
is accustomed to specific notation for the resulting functions, which
jan@42511
   630
are absolutely capable of being totalled and are called step-response: $u[n]$, where $u$ is the
jan@42466
   631
function, $n$ is the argument and the brackets indicate that the
neuper@42504
   632
arguments are discrete. Surprisingly, Isabelle accepts the rules for
jan@42512
   633
${\cal z}^{-1}$ in this traditional notation~\footnote{Isabelle
jan@42466
   634
experts might be particularly surprised, that the brackets do not
jan@42466
   635
cause errors in typing (as lists).}:
neuper@42464
   636
%\vbox{
neuper@42464
   637
% \begin{example}
jan@42463
   638
  \label{eg:neuper1}
jan@42509
   639
  {\footnotesize\begin{tabbing}
jan@42463
   640
  123\=123\=123\=123\=\kill
jan@42509
   641
jan@42463
   642
  \>axiomatization where \\
jan@42512
   643
  \>\>  rule1: ``${\cal z}^{-1}\;1 = \delta [n]$'' and\\
jan@42512
   644
  \>\>  rule2: ``$\vert\vert z \vert\vert > 1 \Rightarrow {\cal z}^{-1}\;z / (z - 1) = u [n]$'' and\\
jan@42466
   645
  \>\>  rule3: ``$\vert\vert$ z $\vert\vert$ < 1 ==> z / (z - 1) = -u [-n - 1]'' and \\
jan@42466
   646
  \>\>  rule4: ``$\vert\vert$ z $\vert\vert$ > $\vert\vert$ $\alpha$ $\vert\vert$ ==> z / (z - $\alpha$) = $\alpha^n$ $\cdot$ u [n]'' and\\
jan@42466
   647
  \>\>  rule5: ``$\vert\vert$ z $\vert\vert$ < $\vert\vert$ $\alpha$ $\vert\vert$ ==> z / (z - $\alpha$) = -($\alpha^n$) $\cdot$ u [-n - 1]'' and\\
jan@42509
   648
  \>\>  rule6: ``$\vert\vert$ z $\vert\vert$ > 1 ==> z/(z - 1)$^2$ = n $\cdot$ u [n]''
jan@42509
   649
  \end{tabbing}}
neuper@42464
   650
% \end{example}
jan@42466
   651
%}
jan@42466
   652
These 6 rules can be used as conditional rewrite rules, depending on
jan@42466
   653
the respective convergence radius. Satisfaction from accordance with traditional notation
jan@42466
   654
contrasts with the above word {\em axiomatization}: As TP-based, the
jan@42466
   655
programming language expects these rules as {\em proved} theorems, and
jan@42466
   656
not as axioms implemented in the above brute force manner; otherwise
jan@42466
   657
all the verification efforts envisaged (like proof of the
jan@42466
   658
post-condition, see below) would be meaningless.
jan@42466
   659
jan@42466
   660
Isabelle provides a large body of knowledge, rigorously proven from
jan@42466
   661
the basic axioms of mathematics~\footnote{This way of rigorously
jan@42466
   662
deriving all knowledge from first principles is called the
jan@42512
   663
LCF-paradigm in TP.}. In the case of the ${\cal z}$-Transform the most advanced
jan@42511
   664
knowledge can be found in the theories on Multivariate
jan@42466
   665
Analysis~\footnote{http://isabelle.in.tum.de/dist/library/HOL/HOL-Multivariate\_Analysis}. However,
jan@42466
   666
building up knowledge such that a proof for the above rules would be
jan@42466
   667
reasonably short and easily comprehensible, still requires lots of
jan@42466
   668
work (and is definitely out of scope of our case study).
jan@42466
   669
neuper@42508
   670
%REMOVED DUE TO SPACE CONSTRAINTS
neuper@42508
   671
%At the state-of-the-art in mechanization of knowledge in engineering
neuper@42508
   672
%sciences, the process does not stop with the mechanization of
neuper@42508
   673
%mathematics traditionally used in these sciences. Rather, ``Formal
neuper@42508
   674
%Methods''~\cite{ fm-03} are expected to proceed to formal and explicit
neuper@42508
   675
%description of physical items.  Signal Processing, for instance is
neuper@42508
   676
%concerned with physical devices for signal acquisition and
neuper@42508
   677
%reconstruction, which involve measuring a physical signal, storing it,
neuper@42508
   678
%and possibly later rebuilding the original signal or an approximation
neuper@42508
   679
%thereof. For digital systems, this typically includes sampling and
neuper@42508
   680
%quantization; devices for signal compression, including audio
neuper@42508
   681
%compression, image compression, and video compression, etc.  ``Domain
neuper@42508
   682
%engineering''\cite{db:dom-eng} is concerned with {\em specification}
neuper@42508
   683
%of these devices' components and features; this part in the process of
neuper@42508
   684
%mechanization is only at the beginning in domains like Signal
neuper@42508
   685
%Processing.
neuper@42508
   686
%
neuper@42508
   687
%TP-based programming, concern of this paper, is determined to
neuper@42508
   688
%add ``algorithmic knowledge'' to the mechanised body of knowledge.
neuper@42508
   689
%% in Fig.\ref{fig:mathuni} on
neuper@42508
   690
%% p.\pageref{fig:mathuni}.  As we shall see below, TP-based programming
neuper@42508
   691
%% starts with a formal {\em specification} of the problem to be solved.
neuper@42508
   692
%% \begin{figure}
neuper@42508
   693
%%   \begin{center}
neuper@42508
   694
%%     \includegraphics[width=110mm]{../../fig/jrocnik/math-universe-small}
neuper@42508
   695
%%     \caption{The three-dimensional universe of mathematics knowledge}
neuper@42508
   696
%%     \label{fig:mathuni}
neuper@42508
   697
%%   \end{center}
neuper@42508
   698
%% \end{figure}
neuper@42508
   699
%% The language for both axes is defined in the axis at the bottom, deductive
neuper@42508
   700
%% knowledge, in {\sisac} represented by Isabelle's theories.
jan@42466
   701
jan@42466
   702
\subsection{Preparation of Simplifiers for the Program}\label{simp}
jan@42469
   703
jan@42511
   704
All evaluation in the prototype's Lucas-Interpreter is done by term rewriting on
neuper@42507
   705
Isabelle's terms, see \S\ref{meth} below; in this section some of respective
jan@42505
   706
preparations are described. In order to work reliably with term rewriting, the
jan@42505
   707
respective rule-sets must be confluent and terminating~\cite{nipk:rew-all-that},
jan@42505
   708
then they are called (canonical) simplifiers. These properties do not go without
jan@42505
   709
saying, their establishment is a difficult task for the programmer; this task is
neuper@42508
   710
not yet supported in the prototype.
jan@42505
   711
jan@42505
   712
The prototype rewrites using theorems only. Axioms which are theorems as well 
jan@42505
   713
have been already shown in \S\ref{eg:neuper1} on p.\pageref{eg:neuper1} , we
jan@42512
   714
assemble them in a rule-set and apply them in ML as follows:
jan@42505
   715
neuper@42508
   716
{\footnotesize
neuper@42508
   717
\begin{verbatim}
jan@42512
   718
   01  val inverse_z = Rls 
jan@42512
   719
   02      {id       = "inverse_z",
jan@42512
   720
   03       rew_ord  = dummy_ord,
jan@42512
   721
   04       erls     = Erls,
jan@42512
   722
   05       rules    = [Thm ("rule1", @{thm rule1}), Thm ("rule2", @{thm rule1}), 
jan@42512
   723
   06                   Thm ("rule3", @{thm rule3}), Thm ("rule4", @{thm rule4}), 
jan@42512
   724
   07                   Thm ("rule5", @{thm rule5}), Thm ("rule6", @{thm rule6})],
jan@42512
   725
   08       errpatts = [],
jan@42512
   726
   09       scr      = ""}
neuper@42508
   727
\end{verbatim}}
jan@42505
   728
neuper@42508
   729
\noindent The items, line by line, in the above record have the following purpose:
neuper@42508
   730
\begin{description}
jan@42512
   731
\item[01..02] the ML-value \textit{inverse\_z} stores it's identifier
neuper@42508
   732
as a string for ``reflection'' when switching between the language
neuper@42508
   733
layers of Isabelle/ML (like in the Lucas-Interpreter) and
neuper@42508
   734
Isabelle/Isar (like in the example program on p.\pageref{s:impl} on
neuper@42508
   735
line {\rm 12}).
jan@42475
   736
jan@42512
   737
\item[03..04] both, (a) the rewrite-order~\cite{nipk:rew-all-that}
neuper@42508
   738
\textit{rew\_ord} and (b) the rule-set \textit{erls} are trivial here:
neuper@42508
   739
(a) the \textit{rules} in {\rm 07..12} don't need ordered rewriting
neuper@42508
   740
and (b) the assumptions of the \textit{rules} need not be evaluated
neuper@42508
   741
(they just go into the context during rewriting).
jan@42505
   742
jan@42512
   743
\item[05..07] the \textit{rules} are the axioms from p.\pageref{eg:neuper1};
neuper@42508
   744
also ML-functions (\S\ref{funs}) can come into this list as shown in
neuper@42508
   745
\S\ref{flow-prep}; so they are distinguished by type-constructors \textit{Thm}
neuper@42508
   746
and \textit{Calc} respectively; for the purpose of reflection both
neuper@42508
   747
contain their identifiers.
jan@42502
   748
jan@42512
   749
\item[08..09] are error-patterns not discussed here and \textit{scr}
neuper@42508
   750
is prepared to get a program, automatically generated by {\sisac} for
neuper@42508
   751
producing intermediate rewrites when requested by the user.
jan@42502
   752
neuper@42508
   753
\end{description}
jan@42505
   754
neuper@42508
   755
\noindent It is advisable to immediately test rule-sets; for that
neuper@42508
   756
purpose an appropriate term has to be created; \textit{parse} takes a
neuper@42508
   757
context \textit{ctxt} and a string (with \textit{ZZ\_1} denoting ${\cal
neuper@42508
   758
Z}^{-1}$) and creates a term:
jan@42505
   759
neuper@42508
   760
{\footnotesize
neuper@42508
   761
\begin{verbatim}
neuper@42508
   762
   01 ML {*
neuper@42508
   763
   02   val t = parse ctxt "ZZ_1 (z / (z - 1) + z / (z - </alpha>) + 1)";
neuper@42508
   764
   03 *}
neuper@42508
   765
   04 val t = Const ("Build_Inverse_Z_Transform.ZZ_1", 
neuper@42508
   766
   05   "RealDef.real => RealDef.real => RealDef.real") $
neuper@42508
   767
   06     (Const (...) $ (Const (...) $ Free (...) $ (Const (...) $ Free (...) 
neuper@42508
   768
\end{verbatim}}
jan@42505
   769
neuper@42508
   770
\noindent The internal representation of the term, as required for
neuper@42508
   771
rewriting, consists of \textit{Const}ants, a pair of a string
neuper@42508
   772
\textit{"Groups.plus\_class.plus"} for $+$ and a type, variables
neuper@42508
   773
\textit{Free} and the respective constructor \textit{\$}. Now the
neuper@42508
   774
term can be rewritten by the rule-set \textit{inverse\_z}:
jan@42505
   775
neuper@42508
   776
{\footnotesize
neuper@42508
   777
\begin{verbatim}
neuper@42508
   778
   01 ML {*
neuper@42508
   779
   02   val SOME (t', asm) = rewrite_set_ @{theory} inverse\_z t;
neuper@42508
   780
   03   term2str t';
neuper@42508
   781
   04   terms2str asm;
neuper@42508
   782
   05 *}
neuper@42508
   783
   06 val it = "u[n] + </alpha> ^ n * u[n] + </delta>[n]" : string
neuper@42508
   784
   07 val it = "|| z || > 1 & || z || > </alpha>" : string
neuper@42508
   785
\end{verbatim}}
jan@42505
   786
neuper@42508
   787
\noindent The resulting term \textit{t} and the assumptions
jan@42511
   788
\textit{asm} are converted to readable strings by \textit{term2str}
neuper@42508
   789
and \textit{terms2str}.
jan@42505
   790
jan@42466
   791
\subsection{Preparation of ML-Functions}\label{funs}
neuper@42504
   792
Some functionality required in programming, cannot be accomplished by
neuper@42504
   793
rewriting. So the prototype has a mechanism to call functions within
jan@42511
   794
the rewrite-engine: certain regexes in Isabelle terms call these
neuper@42504
   795
functions written in SML~\cite{pl:milner97}, the implementation {\em
neuper@42504
   796
and} meta-language of Isabelle. The programmer has to use this
neuper@42504
   797
mechanism.
jan@42469
   798
neuper@42498
   799
In the running example's program on p.\pageref{s:impl} the lines {\rm
neuper@42498
   800
05} and {\rm 06} contain such functions; we go into the details with
neuper@42498
   801
\textit{argument\_in X\_z;}. This function fetches the argument from a
neuper@42498
   802
function application: Line {\rm 03} in the example calculation on
neuper@42498
   803
p.\pageref{exp-calc} is created by line {\rm 06} of the example
neuper@42498
   804
program on p.\pageref{s:impl} where the program's environment assigns
neuper@42498
   805
the value \textit{X z} to the variable \textit{X\_z}; so the function
neuper@42498
   806
shall extract the argument \textit{z}.
jan@42469
   807
neuper@42498
   808
\medskip In order to be recognised as a function constant in the
neuper@42499
   809
program source the constant needs to be declared in a theory, here in
neuper@42498
   810
\textit{Build\_Inverse\_Z\_Transform.thy}; then it can be parsed in
neuper@42498
   811
the context \textit{ctxt} of that theory:
neuper@42504
   812
neuper@42498
   813
{\footnotesize
neuper@42498
   814
\begin{verbatim}
neuper@42498
   815
   consts
neuper@42504
   816
     argument'_in :: "real => real" ("argument'_in _" 10)
neuper@42507
   817
\end{verbatim}}
neuper@42498
   818
   
neuper@42507
   819
%^3.2^    ML {* val SOME t = parse ctxt "argument_in (X z)"; *}
neuper@42507
   820
%^3.2^    val t = Const ("Build_Inverse_Z_Transform.argument'_in", "RealDef.real ⇒ RealDef.real") 
neuper@42507
   821
%^3.2^              $ (Free ("X", "RealDef.real ⇒ RealDef.real") $ Free ("z", "RealDef.real")): term
neuper@42507
   822
%^3.2^ \end{verbatim}}
neuper@42507
   823
%^3.2^ 
neuper@42507
   824
%^3.2^ \noindent Parsing produces a term \texttt{t} in internal
neuper@42507
   825
%^3.2^ representation~\footnote{The attentive reader realizes the 
neuper@42507
   826
%^3.2^ differences between interal and extermal representation even in the
neuper@42507
   827
%^3.2^ strings, i.e \texttt{'\_}}, consisting of \texttt{Const
neuper@42507
   828
%^3.2^ ("argument'\_in", type)} and the two variables \texttt{Free ("X",
neuper@42507
   829
%^3.2^ type)} and \texttt{Free ("z", type)}, \texttt{\$} is the term
neuper@42507
   830
%^3.2^ constructor. 
neuper@42507
   831
The function body below is implemented directly in SML,
neuper@42499
   832
i.e in an \texttt{ML \{* *\}} block; the function definition provides
neuper@42499
   833
a unique prefix \texttt{eval\_} to the function name:
jan@42473
   834
neuper@42498
   835
{\footnotesize
jan@42470
   836
\begin{verbatim}
neuper@42498
   837
   ML {*
neuper@42498
   838
     fun eval_argument_in _ 
neuper@42498
   839
       "Build_Inverse_Z_Transform.argument'_in" 
jan@42512
   840
       (t as (Const ("Build_Inverse_Z_Transform.argument'_in", _) $(f $arg))) _ =
neuper@42498
   841
         if is_Free arg (*could be something to be simplified before*)
jan@42512
   842
         then SOME (term2str t ^"="^ term2str arg, Trueprop $(mk_equality (t, arg)))
neuper@42498
   843
         else NONE
neuper@42498
   844
     | eval_argument_in _ _ _ _ = NONE;
neuper@42498
   845
   *}
neuper@42498
   846
\end{verbatim}}
jan@42469
   847
neuper@42498
   848
\noindent The function body creates either creates \texttt{NONE}
jan@42511
   849
telling the rewrite-engine to search for the next regex, or creates an
neuper@42498
   850
ad-hoc theorem for rewriting, thus the programmer needs to adopt many
neuper@42498
   851
technicalities of Isabelle, for instance, the \textit{Trueprop}
neuper@42498
   852
constant.
jan@42469
   853
neuper@42498
   854
\bigskip This sub-task particularly sheds light on basic issues in the
jan@42511
   855
design of a programming language, the integration of differential language
neuper@42498
   856
layers, the layer of Isabelle/Isar and Isabelle/ML.
jan@42469
   857
neuper@42498
   858
Another point of improvement for the prototype is the rewrite-engine: The
neuper@42498
   859
program on p.\pageref{s:impl} would not allow to contract the two lines {\rm 05}
neuper@42498
   860
and {\rm 06} to
jan@42469
   861
neuper@42498
   862
{\small\it\label{s:impl}
neuper@42498
   863
\begin{tabbing}
neuper@42498
   864
123l\=123\=123\=123\=123\=123\=123\=((x\=123\=(x \=123\=123\=\kill
jan@42512
   865
\>{\rm 05/06}\>\>\>  (z::real) = argument\_in (lhs X\_eq) ;
neuper@42498
   866
\end{tabbing}}
jan@42469
   867
neuper@42498
   868
\noindent because nested function calls would require creating redexes
neuper@42498
   869
inside-out; however, the prototype's rewrite-engine only works top down
neuper@42498
   870
from the root of a term down to the leaves.
jan@42469
   871
neuper@42504
   872
How all these technicalities are to be checked in the prototype is 
neuper@42498
   873
shown in \S\ref{flow-prep} below.
jan@42473
   874
neuper@42498
   875
% \paragraph{Explicit Problems} require explicit methods to solve them, and within
neuper@42498
   876
% this methods we have some explicit steps to do. This steps can be unique for
neuper@42498
   877
% a special problem or refindable in other problems. No mather what case, such
neuper@42498
   878
% steps often require some technical functions behind. For the solving process
neuper@42498
   879
% of the Inverse Z Transformation and the corresponding partial fraction it was
neuper@42498
   880
% neccessary to build helping functions like \texttt{get\_denominator},
neuper@42498
   881
% \texttt{get\_numerator} or \texttt{argument\_in}. First two functions help us
neuper@42498
   882
% to filter the denominator or numerator out of a fraction, last one helps us to
neuper@42498
   883
% get to know the bound variable in a equation.
neuper@42498
   884
% \par
neuper@42498
   885
% By taking \texttt{get\_denominator} as an example, we want to explain how to 
neuper@42498
   886
% implement new functions into the existing system and how we can later use them
neuper@42498
   887
% in our program.
neuper@42498
   888
% 
neuper@42498
   889
% \subsubsection{Find a place to Store the Function}
neuper@42498
   890
% 
neuper@42498
   891
% The whole system builds up on a well defined structure of Knowledge. This
neuper@42498
   892
% Knowledge sets up at the Path:
neuper@42498
   893
% \begin{center}\ttfamily src/Tools/isac/Knowledge\normalfont\end{center}
neuper@42498
   894
% For implementing the Function \texttt{get\_denominator} (which let us extract
neuper@42498
   895
% the denominator out of a fraction) we have choosen the Theory (file)
neuper@42498
   896
% \texttt{Rational.thy}.
neuper@42498
   897
% 
neuper@42498
   898
% \subsubsection{Write down the new Function}
neuper@42498
   899
% 
neuper@42498
   900
% In upper Theory we now define the new function and its purpose:
neuper@42498
   901
% \begin{verbatim}
neuper@42498
   902
%   get_denominator :: "real => real"
neuper@42498
   903
% \end{verbatim}
neuper@42498
   904
% This command tells the machine that a function with the name
neuper@42498
   905
% \texttt{get\_denominator} exists which gets a real expression as argument and
neuper@42498
   906
% returns once again a real expression. Now we are able to implement the function
neuper@42498
   907
% itself, upcoming example now shows the implementation of
neuper@42498
   908
% \texttt{get\_denominator}.
neuper@42498
   909
% 
neuper@42498
   910
% %\begin{example}
neuper@42498
   911
%   \label{eg:getdenom}
neuper@42498
   912
%   \begin{verbatim}
neuper@42498
   913
% 
neuper@42498
   914
% 01  (*
neuper@42498
   915
% 02   *("get_denominator",
neuper@42498
   916
% 03   *  ("Rational.get_denominator", eval_get_denominator ""))
neuper@42498
   917
% 04   *)
neuper@42498
   918
% 05  fun eval_get_denominator (thmid:string) _ 
neuper@42498
   919
% 06            (t as Const ("Rational.get_denominator", _) $
neuper@42498
   920
% 07                (Const ("Rings.inverse_class.divide", _) $num 
neuper@42498
   921
% 08                  $denom)) thy = 
neuper@42498
   922
% 09          SOME (mk_thmid thmid "" 
neuper@42498
   923
% 10              (Print_Mode.setmp [] 
neuper@42498
   924
% 11                (Syntax.string_of_term (thy2ctxt thy)) denom) "", 
neuper@42498
   925
% 12              Trueprop $ (mk_equality (t, denom)))
neuper@42498
   926
% 13    | eval_get_denominator _ _ _ _ = NONE;\end{verbatim}
neuper@42498
   927
% %\end{example}
neuper@42498
   928
% 
neuper@42498
   929
% Line \texttt{07} and \texttt{08} are describing the mode of operation the best -
neuper@42498
   930
% there is a fraction\\ (\ttfamily Rings.inverse\_class.divide\normalfont) 
neuper@42498
   931
% splittet
neuper@42498
   932
% into its two parts (\texttt{\$num \$denom}). The lines before are additionals
neuper@42498
   933
% commands for declaring the function and the lines after are modeling and 
neuper@42498
   934
% returning a real variable out of \texttt{\$denom}.
neuper@42498
   935
% 
neuper@42498
   936
% \subsubsection{Add a test for the new Function}
neuper@42498
   937
% 
neuper@42498
   938
% \paragraph{Everytime when adding} a new function it is essential also to add
neuper@42498
   939
% a test for it. Tests for all functions are sorted in the same structure as the
neuper@42498
   940
% knowledge it self and can be found up from the path:
neuper@42498
   941
% \begin{center}\ttfamily test/Tools/isac/Knowledge\normalfont\end{center}
neuper@42498
   942
% This tests are nothing very special, as a first prototype the functionallity
neuper@42498
   943
% of a function can be checked by evaluating the result of a simple expression
neuper@42498
   944
% passed to the function. Example~\ref{eg:getdenomtest} shows the test for our
neuper@42498
   945
% \textit{just} created function \texttt{get\_denominator}.
neuper@42498
   946
% 
neuper@42498
   947
% %\begin{example}
neuper@42498
   948
% \label{eg:getdenomtest}
neuper@42498
   949
% \begin{verbatim}
neuper@42498
   950
% 
neuper@42498
   951
% 01 val thy = @{theory Isac};
neuper@42498
   952
% 02 val t = term_of (the (parse thy "get_denominator ((a +x)/b)"));
neuper@42498
   953
% 03 val SOME (_, t') = eval_get_denominator "" 0 t thy;
neuper@42498
   954
% 04 if term2str t' = "get_denominator ((a + x) / b) = b" then ()
neuper@42498
   955
% 05 else error "get_denominator ((a + x) / b) = b" \end{verbatim}
neuper@42498
   956
% %\end{example}
neuper@42498
   957
% 
neuper@42498
   958
% \begin{description}
neuper@42498
   959
% \item[01] checks if the proofer set up on our {\sisac{}} System.
neuper@42498
   960
% \item[02] passes a simple expression (fraction) to our suddenly created
neuper@42498
   961
%           function.
neuper@42498
   962
% \item[04] checks if the resulting variable is the correct one (in this case
neuper@42498
   963
%           ``b'' the denominator) and returns.
neuper@42498
   964
% \item[05] handels the error case and reports that the function is not able to
neuper@42498
   965
%           solve the given problem.
neuper@42498
   966
% \end{description}
jan@42469
   967
jan@42491
   968
\subsection{Specification of the Problem}\label{spec}
jan@42491
   969
%WN <--> \chapter 7 der Thesis
jan@42491
   970
%WN die Argumentation unten sollte sich NUR auf Verifikation beziehen..
jan@42491
   971
neuper@42504
   972
Mechanical treatment requires to translate a textual problem
neuper@42504
   973
description like in Fig.\ref{fig-interactive} on
neuper@42504
   974
p.\pageref{fig-interactive} into a {\em formal} specification. The
neuper@42504
   975
formal specification of the running example could look like is this:
jan@42491
   976
jan@42491
   977
%WN Hier brauchen wir die Spezifikation des 'running example' ...
jan@42491
   978
%JR Habe input, output und precond vom Beispiel eingefügt brauche aber Hilfe bei
jan@42491
   979
%JR der post condition - die existiert für uns ja eigentlich nicht aka
jan@42491
   980
%JR haben sie bis jetzt nicht beachtet WN...
jan@42491
   981
%WN2 Mein Vorschlag ist, das TODO zu lassen und deutlich zu kommentieren.
jan@42491
   982
%JR2 done
jan@42491
   983
neuper@42504
   984
\label{eg:neuper2}
neuper@42504
   985
{\small\begin{tabbing}
neuper@42504
   986
  123\=123\=postcond \=: \= $\forall \,A^\prime\, u^\prime \,v^\prime.\,$\=\kill
neuper@42504
   987
  %\hfill \\
neuper@42504
   988
  \>Specification:\\
neuper@42507
   989
  \>  \>input    \>: ${\it filterExpression} \;\;X\;z=\frac{3}{z-\frac{1}{4}+-\frac{1}{8}*\frac{1}{z}}, \;{\it domain}\;\mathbb{R}-\{\frac{1}{2}, \frac{-1}{4}\}$\\
neuper@42504
   990
  \>\>precond  \>: $\frac{3}{z-\frac{1}{4}+-\frac{1}{8}*\frac{1}{z}}\;\; {\it continuous\_on}\; \mathbb{R}-\{\frac{1}{2}, \frac{-1}{4}\}$ \\
neuper@42504
   991
  \>\>output   \>: stepResponse $x[n]$ \\
neuper@42504
   992
  \>\>postcond \>: TODO
neuper@42504
   993
\end{tabbing}}
jan@42491
   994
jan@42500
   995
%JR wie besprochen, kein remark, keine begründung, nur simples "nicht behandelt"
jan@42500
   996
jan@42500
   997
% \begin{remark}
jan@42500
   998
%    Defining the postcondition requires a high amount mathematical 
jan@42500
   999
%    knowledge, the difficult part in our case is not to set up this condition 
jan@42500
  1000
%    nor it is more to define it in a way the interpreter is able to handle it. 
jan@42500
  1001
%    Due the fact that implementing that mechanisms is quite the same amount as 
jan@42500
  1002
%    creating the programm itself, it is not avaible in our prototype.
jan@42500
  1003
%    \label{rm:postcond}
jan@42500
  1004
% \end{remark}
jan@42491
  1005
neuper@42504
  1006
The implementation of the formal specification in the present
neuper@42504
  1007
prototype, still bar-bones without support for authoring, is done
neuper@42504
  1008
like that:
jan@42491
  1009
%WN Kopie von Inverse_Z_Transform.thy, leicht versch"onert:
neuper@42504
  1010
jan@42491
  1011
{\footnotesize\label{exp-spec}
jan@42491
  1012
\begin{verbatim}
neuper@42504
  1013
   00 ML {*
jan@42491
  1014
   01  store_specification
jan@42491
  1015
   02    (prepare_specification
neuper@42504
  1016
   03      "pbl_SP_Ztrans_inv"
neuper@42504
  1017
   04      ["Jan Rocnik"]
jan@42491
  1018
   05      thy
jan@42491
  1019
   06      ( ["Inverse", "Z_Transform", "SignalProcessing"],
neuper@42507
  1020
   07        [ ("#Given", ["filterExpression X_eq", "domain D"]),
neuper@42507
  1021
   08          ("#Pre"  , ["(rhs X_eq) is_continuous_in D"]),
jan@42494
  1022
   09          ("#Find" , ["stepResponse n_eq"]),
neuper@42507
  1023
   10          ("#Post" , [" TODO "])])
neuper@42507
  1024
   11        prls
neuper@42507
  1025
   12        NONE
neuper@42507
  1026
   13        [["SignalProcessing","Z_Transform","Inverse"]]);
neuper@42504
  1027
   14 *}
jan@42491
  1028
\end{verbatim}}
neuper@42504
  1029
jan@42491
  1030
Although the above details are partly very technical, we explain them
jan@42491
  1031
in order to document some intricacies of TP-based programming in the
jan@42491
  1032
present state of the {\sisac} prototype:
jan@42491
  1033
\begin{description}
jan@42491
  1034
\item[01..02]\textit{store\_specification:} stores the result of the
jan@42491
  1035
function \textit{prep\_specification} in a global reference
jan@42491
  1036
\textit{Unsynchronized.ref}, which causes principal conflicts with
jan@42511
  1037
Isabelle's asynchronous document model~\cite{Wenzel-11:doc-orient} and
jan@42491
  1038
parallel execution~\cite{Makarius-09:parall-proof} and is under
jan@42491
  1039
reconstruction already.
jan@42491
  1040
neuper@42504
  1041
\textit{prep\_specification:} translates the specification to an internal format
jan@42491
  1042
which allows efficient processing; see for instance line {\rm 07}
jan@42491
  1043
below.
neuper@42504
  1044
\item[03..04] are a unique identifier for the specification within {\sisac}
neuper@42504
  1045
and the ``mathematics author'' holding the copy-rights.
jan@42491
  1046
\item[05] is the Isabelle \textit{theory} required to parse the
jan@42491
  1047
specification in lines {\rm 07..10}.
jan@42491
  1048
\item[06] is a key into the tree of all specifications as presented to
jan@42511
  1049
the user (where some branches might be hidden by the dialogue
jan@42491
  1050
component).
jan@42491
  1051
\item[07..10] are the specification with input, pre-condition, output
neuper@42507
  1052
and post-condition respectively; note that the specification contains
neuper@42507
  1053
variables to be instantiated with concrete values for a concrete problem ---
neuper@42507
  1054
thus the specification actually captures a class of problems. The post-condition is not handled in
neuper@42504
  1055
the prototype presently.
neuper@42507
  1056
\item[11] is a rule-set (defined elsewhere) for evaluation of the pre-condition: \textit{(rhs X\_eq) is\_continuous\_in D}, instantiated with the values of a concrete problem, evaluates to true or false --- and all evaluation is done by
neuper@42507
  1057
rewriting determined by rule-sets.
jan@42491
  1058
\item[12]\textit{NONE:} could be \textit{SOME ``solve ...''} for a
jan@42491
  1059
problem associated to a function from Computer Algebra (like an
jan@42491
  1060
equation solver) which is not the case here.
neuper@42504
  1061
\item[13] is a list of methods solving the specified problem (here
neuper@42504
  1062
only one list item) represented analogously to {\rm 06}.
jan@42491
  1063
\end{description}
jan@42491
  1064
jan@42491
  1065
jan@42491
  1066
%WN die folgenden Erkl"arungen finden sich durch "grep -r 'datatype pbt' *"
jan@42491
  1067
%WN ...
jan@42491
  1068
%  type pbt = 
jan@42491
  1069
%     {guh  : guh,         (*unique within this isac-knowledge*)
jan@42491
  1070
%      mathauthors: string list, (*copyright*)
jan@42491
  1071
%      init  : pblID,      (*to start refinement with*)
jan@42491
  1072
%      thy   : theory,     (* which allows to compile that pbt
jan@42491
  1073
%			  TODO: search generalized for subthy (ref.p.69*)
jan@42491
  1074
%      (*^^^ WN050912 NOT used during application of the problem,
jan@42491
  1075
%       because applied terms may be from 'subthy' as well as from super;
jan@42491
  1076
%       thus we take 'maxthy'; see match_ags !*)
jan@42491
  1077
%      cas   : term option,(*'CAS-command'*)
jan@42491
  1078
%      prls  : rls,        (* for preds in where_*)
jan@42491
  1079
%      where_: term list,  (* where - predicates*)
jan@42491
  1080
%      ppc   : pat list,
jan@42491
  1081
%      (*this is the model-pattern; 
jan@42491
  1082
%       it contains "#Given","#Where","#Find","#Relate"-patterns
jan@42491
  1083
%       for constraints on identifiers see "fun cpy_nam"*)
jan@42491
  1084
%      met   : metID list}; (* methods solving the pbt*)
jan@42491
  1085
%
jan@42491
  1086
%WN weil dieser Code sehr unaufger"aumt ist, habe ich die Erkl"arungen
jan@42491
  1087
%WN oben selbst geschrieben.
jan@42491
  1088
jan@42491
  1089
jan@42491
  1090
jan@42491
  1091
jan@42491
  1092
%WN das w"urde ich in \sec\label{progr} verschieben und
jan@42491
  1093
%WN das SubProblem partial fractions zum Erkl"aren verwenden.
jan@42491
  1094
% Such a specification is checked before the execution of a program is
jan@42491
  1095
% started, the same applies for sub-programs. In the following example
jan@42491
  1096
% (Example~\ref{eg:subprob}) shows the call of such a subproblem:
jan@42491
  1097
% 
jan@42491
  1098
% \vbox{
jan@42491
  1099
%   \begin{example}
jan@42491
  1100
%   \label{eg:subprob}
jan@42491
  1101
%   \hfill \\
jan@42491
  1102
%   {\ttfamily \begin{tabbing}
jan@42491
  1103
%   ``(L\_L::bool list) = (\=SubProblem (\=Test','' \\
jan@42491
  1104
%   ``\>\>[linear,univariate,equation,test],'' \\
jan@42491
  1105
%   ``\>\>[Test,solve\_linear])'' \\
jan@42491
  1106
%   ``\>[BOOL equ, REAL z])'' \\
jan@42491
  1107
%   \end{tabbing}
jan@42491
  1108
%   }
jan@42491
  1109
%   {\small\textit{
jan@42491
  1110
%     \noindent If a program requires a result which has to be
jan@42491
  1111
% calculated first we can use a subproblem to do so. In our specific
jan@42491
  1112
% case we wanted to calculate the zeros of a fraction and used a
jan@42491
  1113
% subproblem to calculate the zeros of the denominator polynom.
jan@42491
  1114
%     }}
jan@42491
  1115
%   \end{example}
jan@42491
  1116
% }
jan@42491
  1117
jan@42491
  1118
\subsection{Implementation of the Method}\label{meth}
neuper@42504
  1119
A method collects all data required to interpret a certain program by
neuper@42504
  1120
Lucas-Interpretation. The \texttt{program} from p.\pageref{s:impl} of
neuper@42507
  1121
the running example is embedded on the last line in the following method:
neuper@42504
  1122
%The methods represent the different ways a problem can be solved. This can
neuper@42504
  1123
%include mathematical tactics as well as tactics taught in different courses.
neuper@42504
  1124
%Declaring the Method itself gives us the possibilities to describe the way of 
neuper@42504
  1125
%calculation in deep, as well we get the oppertunities to build in different
neuper@42504
  1126
%rulesets.
jan@42491
  1127
jan@42502
  1128
{\footnotesize
jan@42491
  1129
\begin{verbatim}
neuper@42504
  1130
   00 ML {*
neuper@42504
  1131
   01  store_method
neuper@42504
  1132
   02    (prep_method
neuper@42504
  1133
   03      "SP_InverseZTransformation_classic" 
neuper@42504
  1134
   04      ["Jan Rocnik"]
neuper@42504
  1135
   05      thy 
neuper@42507
  1136
   06      ( ["SignalProcessing", "Z_Transform", "Inverse"], 
neuper@42507
  1137
   07        [ ("#Given", ["filterExpression X_eq", "domain D"]),
neuper@42507
  1138
   08          ("#Pre"  , ["(rhs X_eq) is_continuous_in D"]),
neuper@42507
  1139
   09          ("#Find" , ["stepResponse n_eq"]),
neuper@42507
  1140
   10        rew_ord  erls
neuper@42507
  1141
   11        srls  prls  nrls
neuper@42507
  1142
   12        errpats 
neuper@42507
  1143
   13        program);
neuper@42507
  1144
   14 *}
neuper@42504
  1145
\end{verbatim}}
jan@42494
  1146
neuper@42504
  1147
\noindent The above code stores the whole structure analogously to a
neuper@42507
  1148
specification as described above:
neuper@42504
  1149
\begin{description}
neuper@42504
  1150
\item[01..06] are identical to those for the example specification on
neuper@42504
  1151
p.\pageref{exp-spec}.
jan@42494
  1152
neuper@42504
  1153
\item[07..09] show something looking like the specification; this is a
neuper@42507
  1154
{\em guard}: as long as not all \textit{Given} items are present and
neuper@42507
  1155
the \textit{Pre}-conditions is not true, interpretation of the program
neuper@42504
  1156
is not started.
neuper@42504
  1157
neuper@42507
  1158
\item[10..11] all concern rewriting (the respective data are defined elsewhere): \textit{rew\_ord} is the rewrite order~\cite{nipk:rew-all-that} in case
neuper@42507
  1159
\textit{program} contains a \textit{Rewrite} tactic; and in case the respective rule is a conditional rewrite-rule, \textit{erls} features evaluating the conditions. The rule-sets 
neuper@42507
  1160
\textit{srls, prls, nrls} feature evaluating (a) the ML-functions in the program (e.g.
jan@42511
  1161
\textit{lhs, argument\_in, rhs} in the program on p.\pageref{s:impl}, (b) the pre-condition analogous to the specification in line 11 on p.\pageref{exp-spec}
neuper@42507
  1162
and (c) is required for the derivation-machinery checking user-input formulas.
neuper@42504
  1163
neuper@42507
  1164
\item[12..13] \textit{errpats} are error-patterns~\cite{gdaroczy-EP-13} for this method and \textit{program} is the variable holding the example from p.\pageref {s:impl}.
jan@42494
  1165
\end{description}
neuper@42507
  1166
The many rule-sets above cause considerable efforts for the
neuper@42507
  1167
programmers, in particular, because there are no tools for checking
neuper@42507
  1168
essential features of rule-sets.
neuper@42504
  1169
neuper@42504
  1170
% is again very technical and goes hard in detail. Unfortunataly
neuper@42504
  1171
% most declerations are not essential for a basic programm but leads us to a huge
neuper@42504
  1172
% range of powerful possibilities.
neuper@42504
  1173
% 
neuper@42504
  1174
% \begin{description}
neuper@42504
  1175
% \item[01..02] stores the method with the given name into the system under a global
neuper@42504
  1176
% reference.
neuper@42504
  1177
% \item[03] specifies the topic within which context the method can be found.
neuper@42504
  1178
% \item[04..05] as the requirements for different methods can be deviant we 
neuper@42504
  1179
% declare what is \emph{given} and and what to \emph{find} for this specific method.
neuper@42504
  1180
% The code again helds on the topic of the case studie, where the inverse 
neuper@42504
  1181
% z-transformation does a switch between a term describing a electrical filter into
neuper@42504
  1182
% its step response. Also the datatype has to be declared (bool - due the fact that 
neuper@42504
  1183
% we handle equations).
neuper@42504
  1184
% \item[06] \emph{rewrite order} is the order of this rls (ruleset), where one 
neuper@42504
  1185
% theorem of it is used for rewriting one single step.
neuper@42504
  1186
% \item[07] \texttt{rls} is the currently used ruleset for this method. This set
neuper@42504
  1187
% has already been defined before.
neuper@42504
  1188
% \item[08] we would have the possiblitiy to add this method to a predefined tree of
neuper@42504
  1189
% calculations, i.eg. if it would be a sub of a bigger problem, here we leave it
neuper@42504
  1190
% independend.
neuper@42504
  1191
% \item[09] The \emph{source ruleset}, can be used to evaluate list expressions in 
neuper@42504
  1192
% the source.
neuper@42504
  1193
% \item[10] \emph{predicates ruleset} can be used to indicates predicates within 
neuper@42504
  1194
% model patterns.
neuper@42504
  1195
% \item[11] The \emph{check ruleset} summarizes rules for checking formulas 
neuper@42504
  1196
% elementwise.
neuper@42504
  1197
% \item[12] \emph{error patterns} which are expected in this kind of method can be
neuper@42504
  1198
% pre-specified to recognize them during the method.
neuper@42504
  1199
% \item[13] finally the \emph{canonical ruleset}, declares the canonical simplifier 
neuper@42504
  1200
% of the specific method.
neuper@42504
  1201
% \item[14] for this code snipset we don't specify the programm itself and keep it 
neuper@42504
  1202
% empty. Follow up \S\ref{progr} for informations on how to implement this
neuper@42504
  1203
% \textit{main} part.
neuper@42504
  1204
% \end{description}
neuper@42504
  1205
neuper@42478
  1206
\subsection{Implementation of the TP-based Program}\label{progr} 
neuper@42507
  1207
So finally all the prerequisites are described and the final task can
neuper@42480
  1208
be addressed. The program below comes back to the running example: it
neuper@42480
  1209
computes a solution for the problem from Fig.\ref{fig-interactive} on
neuper@42480
  1210
p.\pageref{fig-interactive}. The reader is reminded of
neuper@42480
  1211
\S\ref{PL-isab}, the introduction of the programming language:
jan@42502
  1212
jan@42502
  1213
{\footnotesize\it\label{s:impl}
neuper@42482
  1214
\begin{tabbing}
neuper@42478
  1215
123l\=123\=123\=123\=123\=123\=123\=((x\=123\=(x \=123\=123\=\kill
neuper@42507
  1216
\>{\rm 00}\>ML \{*\\
neuper@42480
  1217
\>{\rm 00}\>val program =\\
neuper@42480
  1218
\>{\rm 01}\>  "{\tt Program} InverseZTransform (X\_eq::bool) =   \\
neuper@42482
  1219
\>{\rm 02}\>\>  {\tt let}                                       \\
neuper@42468
  1220
\>{\rm 03}\>\>\>  X\_eq = {\tt Take} X\_eq ;   \\
neuper@42507
  1221
\>{\rm 04}\>\>\>  X\_eq = {\tt Rewrite} prep\_for\_part\_frac X\_eq ; \\
neuper@42468
  1222
\>{\rm 05}\>\>\>  (X\_z::real) = lhs X\_eq ;       \\ %no inside-out evaluation
neuper@42468
  1223
\>{\rm 06}\>\>\>  (z::real) = argument\_in X\_z; \\
neuper@42468
  1224
\>{\rm 07}\>\>\>  (part\_frac::real) = {\tt SubProblem} \\
neuper@42478
  1225
\>{\rm 08}\>\>\>\>\>\>\>\>  ( Isac, [partial\_fraction, rational, simplification], [] )\\
neuper@42478
  1226
%\>{\rm 10}\>\>\>\>\>\>\>\>\>  [simplification, of\_rationals, to\_partial\_fraction] ) \\
neuper@42478
  1227
\>{\rm 09}\>\>\>\>\>\>\>\>  [ (rhs X\_eq)::real, z::real ]; \\
neuper@42478
  1228
\>{\rm 10}\>\>\>  (X'\_eq::bool) = {\tt Take} ((X'::real =$>$ bool) z = ZZ\_1 part\_frac) ; \\
neuper@42507
  1229
\>{\rm 11}\>\>\>  X'\_eq = (({\tt Rewrite\_Set} prep\_for\_inverse\_z) @@   \\
neuper@42478
  1230
\>{\rm 12}\>\>\>\>\>  $\;\;$ ({\tt Rewrite\_Set} inverse\_z)) X'\_eq \\
neuper@42482
  1231
\>{\rm 13}\>\>  {\tt in } \\
neuper@42504
  1232
\>{\rm 14}\>\>\>  X'\_eq"\\
neuper@42507
  1233
\>{\rm 15}\>*\}
neuper@42478
  1234
\end{tabbing}}
neuper@42468
  1235
% ORIGINAL FROM Inverse_Z_Transform.thy
neuper@42468
  1236
% "Script InverseZTransform (X_eq::bool) =            "^(*([], Frm), Problem (Isac, [Inverse, Z_Transform, SignalProcessing])*)
neuper@42468
  1237
% "(let X = Take X_eq;                                "^(*([1], Frm), X z = 3 / (z - 1 / 4 + -1 / 8 * (1 / z))*)
neuper@42468
  1238
% "  X' = Rewrite ruleZY False X;                     "^(*([1], Res), ?X' z = 3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))*)
neuper@42468
  1239
% "  (X'_z::real) = lhs X';                           "^(*            ?X' z*)
neuper@42468
  1240
% "  (zzz::real) = argument_in X'_z;                  "^(*            z *)
neuper@42468
  1241
% "  (funterm::real) = rhs X';                        "^(*            3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))*)
neuper@42468
  1242
%
neuper@42468
  1243
% "  (pbz::real) = (SubProblem (Isac',                "^(**)
neuper@42468
  1244
% "    [partial_fraction,rational,simplification],    "^
neuper@42468
  1245
% "    [simplification,of_rationals,to_partial_fraction]) "^
neuper@42468
  1246
% "    [REAL funterm, REAL zzz]);                     "^(*([2], Res), 4 / (z - 1 / 2) + -4 / (z - -1 / 4)*)
neuper@42468
  1247
%
neuper@42468
  1248
% "  (pbz_eq::bool) = Take (X'_z = pbz);              "^(*([3], Frm), ?X' z = 4 / (z - 1 / 2) + -4 / (z - -1 / 4)*)
neuper@42468
  1249
% "  pbz_eq = Rewrite ruleYZ False pbz_eq;            "^(*([3], Res), ?X' z = 4 * (?z / (z - 1 / 2)) + -4 * (?z / (z - -1 / 4))*)
neuper@42468
  1250
% "  pbz_eq = drop_questionmarks pbz_eq;              "^(*               4 * (z / (z - 1 / 2)) + -4 * (z / (z - -1 / 4))*)
neuper@42468
  1251
% "  (X_zeq::bool) = Take (X_z = rhs pbz_eq);         "^(*([4], Frm), X_z = 4 * (z / (z - 1 / 2)) + -4 * (z / (z - -1 / 4))*)
neuper@42468
  1252
% "  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]*)
neuper@42468
  1253
% "  n_eq = drop_questionmarks n_eq                   "^(*            X_z = 4 * (1 / 2) ^^^ n * u [n] + -4 * (-1 / 4) ^^^ n * u [n]*)
neuper@42468
  1254
% "in n_eq)"                                            (*([], Res), X_z = 4 * (1 / 2) ^^^ n * u [n] + -4 * (-1 / 4) ^^^ n * u [n]*)
neuper@42480
  1255
The program is represented as a string and part of the method in
neuper@42480
  1256
\S\ref{meth}.  As mentioned in \S\ref{PL} the program is purely
neuper@42480
  1257
functional and lacks any input statements and output statements. So
neuper@42480
  1258
the steps of calculation towards a solution (and interactive tutoring
neuper@42480
  1259
in step-wise problem solving) are created as a side-effect by
neuper@42480
  1260
Lucas-Interpretation.  The side-effects are triggered by the tactics
neuper@42482
  1261
\texttt{Take}, \texttt{Rewrite}, \texttt{SubProblem} and
neuper@42482
  1262
\texttt{Rewrite\_Set} in the above lines {\rm 03, 04, 07, 10, 11} and
neuper@42507
  1263
{\rm 12} respectively. These tactics produce the respective lines in the
neuper@42480
  1264
calculation on p.\pageref{flow-impl}.
neuper@42478
  1265
neuper@42480
  1266
The above lines {\rm 05, 06} do not contain a tactics, so they do not
neuper@42480
  1267
immediately contribute to the calculation on p.\pageref{flow-impl};
neuper@42482
  1268
rather, they compute actual arguments for the \texttt{SubProblem} in
neuper@42480
  1269
line {\rm 09}~\footnote{The tactics also are break-points for the
neuper@42480
  1270
interpreter, where control is handed over to the user in interactive
neuper@42482
  1271
tutoring.}. Line {\rm 11} contains tactical \textit{@@}.
neuper@42480
  1272
neuper@42480
  1273
\medskip The above program also indicates the dominant role of interactive
neuper@42478
  1274
selection of knowledge in the three-dimensional universe of
neuper@42478
  1275
mathematics as depicted in Fig.\ref{fig:mathuni} on
neuper@42482
  1276
p.\pageref{fig:mathuni}, The \texttt{SubProblem} in the above lines
neuper@42478
  1277
{\rm 07..09} is more than a function call with the actual arguments
neuper@42478
  1278
\textit{[ (rhs X\_eq)::real, z::real ]}. The programmer has to determine
neuper@42478
  1279
three items:
neuper@42480
  1280
neuper@42478
  1281
\begin{enumerate}
neuper@42478
  1282
\item the theory, in the example \textit{Isac} because different
neuper@42478
  1283
methods can be selected in Pt.3 below, which are defined in different
neuper@42478
  1284
theories with \textit{Isac} collecting them.
neuper@42480
  1285
\item the specification identified by \textit{[partial\_fraction,
neuper@42480
  1286
rational, simplification]} in the tree of specifications; this
neuper@42480
  1287
specification is analogous to the specification of the main program
neuper@42480
  1288
described in \S\ref{spec}; the problem is to find a ``partial fraction
neuper@42480
  1289
decomposition'' for a univariate rational polynomial.
neuper@42480
  1290
\item the method in the above example is \textit{[ ]}, i.e. empty,
neuper@42480
  1291
which supposes the interpreter to select one of the methods predefined
neuper@42480
  1292
in the specification, for instance in line {\rm 13} in the running
neuper@42480
  1293
example's specification on p.\pageref{exp-spec}~\footnote{The freedom
neuper@42480
  1294
(or obligation) for selection carries over to the student in
neuper@42480
  1295
interactive tutoring.}.
neuper@42478
  1296
\end{enumerate}
neuper@42478
  1297
neuper@42480
  1298
The program code, above presented as a string, is parsed by Isabelle's
neuper@42480
  1299
parser --- the program is an Isabelle term. This fact is expected to
neuper@42480
  1300
simplify verification tasks in the future; on the other hand, this
jan@42511
  1301
fact causes troubles in error detection which are discussed as part
neuper@42480
  1302
of the workflow in the subsequent section.
neuper@42467
  1303
jan@42463
  1304
\section{Workflow of Programming in the Prototype}\label{workflow}
neuper@42498
  1305
The new prover IDE Isabelle/jEdit~\cite{makar-jedit-12} is a great
neuper@42498
  1306
step forward for interactive theory and proof development. The
neuper@42498
  1307
{\sisac}-prototype re-uses this IDE as a programming environment.  The
neuper@42498
  1308
experiences from this re-use show, that the essential components are
neuper@42498
  1309
available from Isabelle/jEdit. However, additional tools and features
jan@42511
  1310
are required to achieve acceptable usability.
neuper@42498
  1311
neuper@42498
  1312
So notable experiences are reported here, also as a requirement
neuper@42498
  1313
capture for further development of TP-based languages and respective
neuper@42498
  1314
IDEs.
neuper@42468
  1315
jan@42466
  1316
\subsection{Preparations and Trials}\label{flow-prep}
neuper@42499
  1317
The many sub-tasks to be accomplished {\em before} the first line of
neuper@42499
  1318
program code can be written and tested suggest an approach which
neuper@42499
  1319
step-wise establishes the prerequisites. The case study underlying
neuper@42499
  1320
this paper~\cite{jrocnik-bakk} documents the approach in a separate
neuper@42499
  1321
Isabelle theory,
neuper@42499
  1322
\textit{Build\_Inverse\_Z\_Transform.thy}~\footnote{http://www.ist.tugraz.at/projects/isac/publ/Build\_Inverse\_Z\_Transform.thy}. Part
neuper@42499
  1323
II in the study comprises this theory, \LaTeX ed from the theory by
neuper@42499
  1324
use of Isabelle's document preparation system. This paper resembles
neuper@42499
  1325
the approach in \S\ref{isabisac} to \S\ref{meth}, which in actual
neuper@42499
  1326
implementation work involves several iterations.
neuper@42498
  1327
neuper@42499
  1328
\bigskip For instance, only the last step, implementing the program
neuper@42499
  1329
described in \S\ref{meth}, reveals details required. Let us assume,
neuper@42499
  1330
this is the ML-function \textit{argument\_in} required in line {\rm 06}
neuper@42499
  1331
of the example program on p.\pageref{s:impl}; how this function needs
neuper@42499
  1332
to be implemented in the prototype has been discussed in \S\ref{funs}
neuper@42499
  1333
already.
neuper@42498
  1334
neuper@42499
  1335
Now let us assume, that calling this function from the program code
neuper@42499
  1336
does not work; so testing this function is required in order to find out
neuper@42499
  1337
the reason: type errors, a missing entry of the function somewhere or
neuper@42499
  1338
even more nasty technicalities \dots
neuper@42498
  1339
neuper@42499
  1340
{\footnotesize
neuper@42482
  1341
\begin{verbatim}
neuper@42482
  1342
   ML {*
neuper@42499
  1343
     val SOME t = parseNEW ctxt "argument_in (X (z::real))";
neuper@42499
  1344
     val SOME (str, t') = eval_argument_in "" 
neuper@42499
  1345
       "Build_Inverse_Z_Transform.argument'_in" t 0;
neuper@42499
  1346
     term2str t';
neuper@42499
  1347
   *}
neuper@42499
  1348
   val it = "(argument_in X z) = z": string
neuper@42482
  1349
\end{verbatim}}
neuper@42499
  1350
neuper@42499
  1351
\noindent So, this works: we get an ad-hoc theorem, which used in
neuper@42499
  1352
rewriting would reduce \texttt{argument\_in X z} to \texttt{z}. Now we check this
neuper@42499
  1353
reduction and create a rule-set \texttt{rls} for that purpose:
neuper@42499
  1354
neuper@42499
  1355
{\footnotesize
neuper@42482
  1356
\begin{verbatim}
neuper@42482
  1357
   ML {*
neuper@42499
  1358
     val rls = append_rls "test" e_rls 
neuper@42499
  1359
       [Calc ("Build_Inverse_Z_Transform.argument'_in", eval_argument_in "")]
neuper@42499
  1360
     val SOME (t', asm) = rewrite_set_ @{theory} rls t;
neuper@42499
  1361
   *}
neuper@42499
  1362
   val t' = Free ("z", "RealDef.real"): term
neuper@42499
  1363
   val asm = []: term list
neuper@42482
  1364
\end{verbatim}}
neuper@42499
  1365
neuper@42499
  1366
\noindent The resulting term \texttt{t'} is \texttt{Free ("z",
neuper@42499
  1367
"RealDef.real")}, i.e the variable \texttt{z}, so all is
neuper@42499
  1368
perfect. Probably we have forgotten to store this function correctly~?
neuper@42499
  1369
We review the respective \texttt{calclist} (again an
neuper@42499
  1370
\textit{Unsynchronized.ref} to be removed in order to adjust to
jan@42511
  1371
IsabelleIsar's asynchronous document model):
neuper@42499
  1372
neuper@42499
  1373
{\footnotesize
neuper@42499
  1374
\begin{verbatim}
neuper@42499
  1375
   calclist:= overwritel (! calclist, 
neuper@42499
  1376
    [("argument_in",("Build_Inverse_Z_Transform.argument'_in", eval_argument_in "")),
neuper@42499
  1377
     ...
neuper@42499
  1378
     ]);
neuper@42499
  1379
\end{verbatim}}
neuper@42499
  1380
neuper@42499
  1381
\noindent The entry is perfect. So what is the reason~? Ah, probably there
neuper@42499
  1382
is something messed up with the many rule-sets in the method, see \S\ref{meth} ---
neuper@42499
  1383
right, the function \texttt{argument\_in} is not contained in the respective
neuper@42499
  1384
rule-set \textit{srls} \dots this just as an example of the intricacies in
neuper@42499
  1385
debugging a program in the present state of the prototype.
neuper@42499
  1386
neuper@42499
  1387
\subsection{Implementation in Isabelle/{\isac}}\label{flow-impl}
neuper@42499
  1388
Given all the prerequisites from \S\ref{isabisac} to \S\ref{meth},
neuper@42499
  1389
usually developed within several iterations, the program can be
neuper@42499
  1390
assembled; on p.\pageref{s:impl} there is the complete program of the
neuper@42499
  1391
running example.
neuper@42499
  1392
neuper@42499
  1393
The completion of this program required efforts for several weeks
neuper@42499
  1394
(after some months of familiarisation with {\sisac}), caused by the
neuper@42499
  1395
abundance of intricacies indicated above. Also writing the program is
neuper@42499
  1396
not pleasant, given Isabelle/Isar/ without add-ons for
neuper@42499
  1397
programming. Already writing and parsing a few lines of program code
neuper@42499
  1398
is a challenge: the program is an Isabelle term; Isabelle's parser,
neuper@42499
  1399
however, is not meant for huge terms like the program of the running
neuper@42499
  1400
example. So reading out the specific error (usually type errors) from
neuper@42499
  1401
Isabelle's message is difficult.
neuper@42499
  1402
neuper@42499
  1403
\medskip Testing the evaluation of the program has to rely on very
neuper@42499
  1404
simple tools. Step-wise execution is modelled by a function
neuper@42499
  1405
\texttt{me}, short for mathematics-engine~\footnote{The interface used
neuper@42499
  1406
by the fron-end which created the calculation on
neuper@42499
  1407
p.\pageref{fig-interactive} is different from this function}:
neuper@42499
  1408
%the following is a simplification of the actual function 
neuper@42499
  1409
neuper@42499
  1410
{\footnotesize
neuper@42499
  1411
\begin{verbatim}
neuper@42499
  1412
   ML {* me; *}
neuper@42499
  1413
   val it = tac -> ctree * pos -> mout * tac * ctree * pos
neuper@42499
  1414
\end{verbatim}} 
neuper@42499
  1415
neuper@42499
  1416
\noindent This function takes as arguments a tactic \texttt{tac} which
neuper@42499
  1417
determines the next step, the step applied to the interpreter-state
neuper@42499
  1418
\texttt{ctree * pos} as last argument taken. The interpreter-state is
neuper@42499
  1419
a pair of a tree \texttt{ctree} representing the calculation created
neuper@42499
  1420
(see the example below) and a position \texttt{pos} in the
jan@42511
  1421
calculation. The function delivers a quadruple, beginning with the new
neuper@42499
  1422
formula \texttt{mout} and the next tactic followed by the new
neuper@42499
  1423
interpreter-state.
neuper@42499
  1424
neuper@42499
  1425
This function allows to stepwise check the program:
neuper@42499
  1426
neuper@42499
  1427
{\footnotesize
neuper@42482
  1428
\begin{verbatim}
neuper@42482
  1429
   ML {*
neuper@42499
  1430
     val fmz =
neuper@42499
  1431
       ["filterExpression (X z = 3 / ((z::real) + 1/10 - 1/50*(1/z)))",
neuper@42499
  1432
        "stepResponse (x[n::real]::bool)"];     
neuper@42499
  1433
     val (dI,pI,mI) =
neuper@42499
  1434
       ("Isac", 
neuper@42499
  1435
        ["Inverse", "Z_Transform", "SignalProcessing"], 
neuper@42499
  1436
        ["SignalProcessing","Z_Transform","Inverse"]);
neuper@42499
  1437
     val (mout, tac, ctree, pos)  = CalcTreeTEST [(fmz, (dI, pI, mI))];
neuper@42499
  1438
     val (mout, tac, ctree, pos)  = me tac (ctree, pos);
neuper@42499
  1439
     val (mout, tac, ctree, pos)  = me tac (ctree, pos);
neuper@42499
  1440
     val (mout, tac, ctree, pos)  = me tac (ctree, pos);
jan@42512
  1441
     ...\end{verbatim}} 
neuper@42481
  1442
jan@42511
  1443
\noindent Several dozens of calls for \texttt{me} are required to
neuper@42499
  1444
create the lines in the calculation below (including the sub-problems
neuper@42499
  1445
not shown). When an error occurs, the reason might be located
neuper@42499
  1446
many steps before: if evaluation by rewriting, as done by the prototype,
neuper@42499
  1447
fails, then first nothing happens --- the effects come later and
neuper@42499
  1448
cause unpleasant checks.
neuper@42481
  1449
neuper@42499
  1450
The checks comprise watching the rewrite-engine for many different
neuper@42499
  1451
kinds of rule-sets (see \S\ref{meth}), the interpreter-state, in
neuper@42499
  1452
particular the environment and the context at the states position ---
neuper@42499
  1453
all checks have to rely on simple functions accessing the
neuper@42499
  1454
\texttt{ctree}. So getting the calculation below (which resembles the
neuper@42499
  1455
calculation in Fig.\ref{fig-interactive} on p.\pageref{fig-interactive})
neuper@42507
  1456
is the result of several weeks of development:
jan@42469
  1457
neuper@42498
  1458
{\small\it\label{exp-calc}
neuper@42468
  1459
\begin{tabbing}
neuper@42468
  1460
123l\=123\=123\=123\=123\=123\=123\=123\=123\=123\=123\=123\=\kill
neuper@42468
  1461
\>{\rm 01}\> $\bullet$  \> {\tt Problem } (Inverse\_Z\_Transform, [Inverse, Z\_Transform, SignalProcessing])       \`\\
neuper@42468
  1462
\>{\rm 02}\>\> $\vdash\;\;X z = \frac{3}{z - \frac{1}{4} - \frac{1}{8} \cdot z^{-1}}$       \`{\footnotesize {\tt Take} X\_eq}\\
neuper@42507
  1463
\>{\rm 03}\>\> $X z = \frac{3}{z + \frac{-1}{4} + \frac{-1}{8} \cdot \frac{1}{z}}$          \`{\footnotesize {\tt Rewrite} prep\_for\_part\_frac X\_eq}\\
neuper@42468
  1464
\>{\rm 04}\>\> $\bullet$\> {\tt Problem } [partial\_fraction,rational,simplification]        \`{\footnotesize {\tt SubProblem} \dots}\\
neuper@42468
  1465
\>{\rm 05}\>\>\>  $\vdash\;\;\frac{3}{z + \frac{-1}{4} + \frac{-1}{8} \cdot \frac{1}{z}}=$    \`- - -\\
neuper@42468
  1466
\>{\rm 06}\>\>\>  $\frac{24}{-1 + -2 \cdot z + 8 \cdot z^2}$                                   \`- - -\\
neuper@42468
  1467
\>{\rm 07}\>\>\>  $\bullet$\> solve ($-1 + -2 \cdot z + 8 \cdot z^2,\;z$ )                      \`- - -\\
neuper@42468
  1468
\>{\rm 08}\>\>\>\>   $\vdash$ \> $\frac{3}{z + \frac{-1}{4} + \frac{-1}{8} \cdot \frac{1}{z}}=0$ \`- - -\\
neuper@42468
  1469
\>{\rm 09}\>\>\>\>   $z = \frac{2+\sqrt{-4+8}}{16}\;\lor\;z = \frac{2-\sqrt{-4+8}}{16}$           \`- - -\\
neuper@42468
  1470
\>{\rm 10}\>\>\>\>   $z = \frac{1}{2}\;\lor\;z =$ \_\_\_                                           \`- - -\\
neuper@42468
  1471
\>        \>\>\>\>   \_\_\_                                                                        \`- - -\\
neuper@42468
  1472
\>{\rm 11}\>\> \dots\> $\frac{4}{z - \frac{1}{2}} + \frac{-4}{z - \frac{-1}{4}}$                   \`\\
jan@42512
  1473
\>{\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)}\\
jan@42512
  1474
\>{\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} prep\_for\_inverse\_z X'\_eq }\\
neuper@42468
  1475
\>{\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}\\
neuper@42468
  1476
\>{\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}}
neuper@42468
  1477
\end{tabbing}}
neuper@42507
  1478
The tactics on the right margin of the above calculation are those in
neuper@42507
  1479
the program on p.\pageref{s:impl} which create the respective formulas
neuper@42507
  1480
on the left.
neuper@42468
  1481
% ORIGINAL FROM Inverse_Z_Transform.thy
neuper@42468
  1482
%    "Script InverseZTransform (X_eq::bool) =            "^(*([], Frm), Problem (Isac, [Inverse, Z_Transform, SignalProcessing])*)
neuper@42468
  1483
%    "(let X = Take X_eq;                                "^(*([1], Frm), X z = 3 / (z - 1 / 4 + -1 / 8 * (1 / z))*)
neuper@42468
  1484
%    "  X' = Rewrite ruleZY False X;                     "^(*([1], Res), ?X' z = 3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))*)
neuper@42468
  1485
%    "  (X'_z::real) = lhs X';                           "^(*            ?X' z*)
neuper@42468
  1486
%    "  (zzz::real) = argument_in X'_z;                  "^(*            z *)
neuper@42468
  1487
%    "  (funterm::real) = rhs X';                        "^(*            3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))*)
neuper@42468
  1488
% 
neuper@42468
  1489
%    "  (pbz::real) = (SubProblem (Isac',                "^(**)
neuper@42468
  1490
%    "    [partial_fraction,rational,simplification],    "^
neuper@42468
  1491
%    "    [simplification,of_rationals,to_partial_fraction]) "^
neuper@42468
  1492
%    "    [REAL funterm, REAL zzz]);                     "^(*([2], Res), 4 / (z - 1 / 2) + -4 / (z - -1 / 4)*)
neuper@42468
  1493
% 
neuper@42468
  1494
%    "  (pbz_eq::bool) = Take (X'_z = pbz);              "^(*([3], Frm), ?X' z = 4 / (z - 1 / 2) + -4 / (z - -1 / 4)*)
neuper@42468
  1495
%    "  pbz_eq = Rewrite ruleYZ False pbz_eq;            "^(*([3], Res), ?X' z = 4 * (?z / (z - 1 / 2)) + -4 * (?z / (z - -1 / 4))*)
neuper@42468
  1496
%    "  pbz_eq = drop_questionmarks pbz_eq;              "^(*               4 * (z / (z - 1 / 2)) + -4 * (z / (z - -1 / 4))*)
neuper@42468
  1497
%    "  (X_zeq::bool) = Take (X_z = rhs pbz_eq);         "^(*([4], Frm), X_z = 4 * (z / (z - 1 / 2)) + -4 * (z / (z - -1 / 4))*)
neuper@42468
  1498
%    "  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]*)
neuper@42468
  1499
%    "  n_eq = drop_questionmarks n_eq                   "^(*            X_z = 4 * (1 / 2) ^^^ n * u [n] + -4 * (-1 / 4) ^^^ n * u [n]*)
neuper@42468
  1500
%    "in n_eq)"                                            (*([], Res), X_z = 4 * (1 / 2) ^^^ n * u [n] + -4 * (-1 / 4) ^^^ n * u [n]*)
neuper@42468
  1501
neuper@42499
  1502
\subsection{Transfer into the Isabelle/{\isac} Knowledge}\label{flow-trans}
neuper@42499
  1503
Finally \textit{Build\_Inverse\_Z\_Transform.thy} has got the job done
neuper@42499
  1504
and the knowledge accumulated in it can be distributed to appropriate
neuper@42499
  1505
theories: the program to \textit{Inverse\_Z\_Transform.thy}, the
neuper@42499
  1506
sub-problem accomplishing the partial fraction decomposition to
neuper@42499
  1507
\textit{Partial\_Fractions.thy}. Since there are hacks into Isabelle's
neuper@42499
  1508
internals, this kind of distribution is not trivial. For instance, the
neuper@42499
  1509
function \texttt{argument\_in} in \S\ref{funs} explicitly contains a
neuper@42499
  1510
string with the theory it has been defined in, so this string needs to
neuper@42499
  1511
be updated from \texttt{Build\_Inverse\_Z\_Transform} to
neuper@42499
  1512
\texttt{Atools} if that function is transferred to theory
neuper@42499
  1513
\textit{Atools.thy}.
neuper@42468
  1514
neuper@42499
  1515
In order to obtain the functionality presented in Fig.\ref{fig-interactive} on p.\pageref{fig-interactive} data must be exported from SML-structures to XML.
neuper@42499
  1516
This process is also rather bare-bones without authoring tools and is
neuper@42499
  1517
described in detail in the {\sisac} wiki~\footnote{http://www.ist.tugraz.at/isac/index.php/Generate\_representations\_for\_ISAC\_Knowledge}.
neuper@42468
  1518
neuper@42478
  1519
% \newpage
neuper@42478
  1520
% -------------------------------------------------------------------
neuper@42478
  1521
% 
neuper@42478
  1522
% Material, falls noch Platz bleibt ...
neuper@42478
  1523
% 
neuper@42478
  1524
% -------------------------------------------------------------------
neuper@42478
  1525
% 
neuper@42478
  1526
% 
neuper@42478
  1527
% \subsubsection{Trials on Notation and Termination}
neuper@42478
  1528
% 
neuper@42478
  1529
% \paragraph{Technical notations} are a big problem for our piece of software,
neuper@42478
  1530
% but the reason for that isn't a fault of the software itself, one of the
neuper@42478
  1531
% troubles comes out of the fact that different technical subtopics use different
neuper@42478
  1532
% symbols and notations for a different purpose. The most famous example for such
neuper@42478
  1533
% a symbol is the complex number $i$ (in cassique math) or $j$ (in technical
neuper@42478
  1534
% math). In the specific part of signal processing one of this notation issues is
neuper@42478
  1535
% the use of brackets --- we use round brackets for analoge signals and squared
neuper@42478
  1536
% brackets for digital samples. Also if there is no problem for us to handle this
neuper@42478
  1537
% fact, we have to tell the machine what notation leads to wich meaning and that
neuper@42478
  1538
% this purpose seperation is only valid for this special topic - signal
neuper@42478
  1539
% processing.
neuper@42478
  1540
% \subparagraph{In the programming language} itself it is not possible to declare
neuper@42478
  1541
% fractions, exponents, absolutes and other operators or remarks in a way to make
neuper@42478
  1542
% them pretty to read; our only posssiblilty were ASCII characters and a handfull
neuper@42478
  1543
% greek symbols like: $\alpha, \beta, \gamma, \phi,\ldots$.
neuper@42478
  1544
% \par
neuper@42478
  1545
% With the upper collected knowledge it is possible to check if we were able to
neuper@42478
  1546
% donate all required terms and expressions.
neuper@42478
  1547
% 
neuper@42478
  1548
% \subsubsection{Definition and Usage of Rules}
neuper@42478
  1549
% 
neuper@42478
  1550
% \paragraph{The core} of our implemented problem is the Z-Transformation, due
neuper@42478
  1551
% the fact that the transformation itself would require higher math which isn't
neuper@42478
  1552
% yet avaible in our system we decided to choose the way like it is applied in
neuper@42478
  1553
% labratory and problem classes at our university - by applying transformation
neuper@42478
  1554
% rules (collected in transformation tables).
neuper@42478
  1555
% \paragraph{Rules,} in {\sisac{}}'s programming language can be designed by the
neuper@42478
  1556
% use of axiomatizations like shown in Example~\ref{eg:ruledef}
neuper@42478
  1557
% 
neuper@42478
  1558
% \begin{example}
neuper@42478
  1559
%   \label{eg:ruledef}
neuper@42478
  1560
%   \hfill\\
neuper@42478
  1561
%   \begin{verbatim}
neuper@42478
  1562
%   axiomatization where
neuper@42478
  1563
%     rule1: ``1 = $\delta$[n]'' and
neuper@42478
  1564
%     rule2: ``|| z || > 1 ==> z / (z - 1) = u [n]'' and
neuper@42478
  1565
%     rule3: ``|| z || < 1 ==> z / (z - 1) = -u [-n - 1]''
neuper@42478
  1566
%   \end{verbatim}
neuper@42478
  1567
% \end{example}
neuper@42478
  1568
% 
neuper@42478
  1569
% This rules can be collected in a ruleset and applied to a given expression as
neuper@42478
  1570
% follows in Example~\ref{eg:ruleapp}.
neuper@42478
  1571
% 
neuper@42478
  1572
% \begin{example}
neuper@42478
  1573
%   \hfill\\
neuper@42478
  1574
%   \label{eg:ruleapp}
neuper@42478
  1575
%   \begin{enumerate}
neuper@42478
  1576
%   \item Store rules in ruleset:
neuper@42478
  1577
%   \begin{verbatim}
neuper@42478
  1578
%   val inverse_Z = append_rls "inverse_Z" e_rls
neuper@42478
  1579
%     [ Thm ("rule1",num_str @{thm rule1}),
neuper@42478
  1580
%       Thm ("rule2",num_str @{thm rule2}),
neuper@42478
  1581
%       Thm ("rule3",num_str @{thm rule3})
neuper@42478
  1582
%     ];\end{verbatim}
neuper@42478
  1583
%   \item Define exression:
neuper@42478
  1584
%   \begin{verbatim}
neuper@42478
  1585
%   val sample_term = str2term "z/(z-1)+z/(z-</delta>)+1";\end{verbatim}
neuper@42478
  1586
%   \item Apply ruleset:
neuper@42478
  1587
%   \begin{verbatim}
neuper@42478
  1588
%   val SOME (sample_term', asm) = 
neuper@42478
  1589
%     rewrite_set_ thy true inverse_Z sample_term;\end{verbatim}
neuper@42478
  1590
%   \end{enumerate}
neuper@42478
  1591
% \end{example}
neuper@42478
  1592
% 
neuper@42478
  1593
% The use of rulesets makes it much easier to develop our designated applications,
neuper@42478
  1594
% but the programmer has to be careful and patient. When applying rulesets
neuper@42478
  1595
% two important issues have to be mentionend:
neuper@42478
  1596
% \subparagraph{How often} the rules have to be applied? In case of
neuper@42478
  1597
% transformations it is quite clear that we use them once but other fields
neuper@42478
  1598
% reuqire to apply rules until a special condition is reached (e.g.
neuper@42478
  1599
% a simplification is finished when there is nothing to be done left).
neuper@42478
  1600
% \subparagraph{The order} in which rules are applied often takes a big effect
neuper@42478
  1601
% and has to be evaluated for each purpose once again.
neuper@42478
  1602
% \par
neuper@42478
  1603
% In our special case of Signal Processing and the rules defined in
neuper@42478
  1604
% Example~\ref{eg:ruledef} we have to apply rule~1 first of all to transform all
neuper@42478
  1605
% constants. After this step has been done it no mather which rule fit's next.
neuper@42478
  1606
% 
neuper@42478
  1607
% \subsubsection{Helping Functions}
neuper@42478
  1608
% 
neuper@42478
  1609
% \paragraph{New Programms require,} often new ways to get through. This new ways
neuper@42478
  1610
% means that we handle functions that have not been in use yet, they can be 
neuper@42478
  1611
% something special and unique for a programm or something famous but unneeded in
neuper@42478
  1612
% the system yet. In our dedicated example it was for example neccessary to split
neuper@42478
  1613
% a fraction into numerator and denominator; the creation of such function and
neuper@42478
  1614
% even others is described in upper Sections~\ref{simp} and \ref{funs}.
neuper@42478
  1615
% 
neuper@42478
  1616
% \subsubsection{Trials on equation solving}
neuper@42478
  1617
% %simple eq and problem with double fractions/negative exponents
neuper@42478
  1618
% \paragraph{The Inverse Z-Transformation} makes it neccessary to solve
neuper@42478
  1619
% equations degree one and two. Solving equations in the first degree is no 
neuper@42478
  1620
% problem, wether for a student nor for our machine; but even second degree
neuper@42478
  1621
% equations can lead to big troubles. The origin of this troubles leads from
neuper@42478
  1622
% the build up process of our equation solving functions; they have been
neuper@42478
  1623
% implemented some time ago and of course they are not as good as we want them to
neuper@42478
  1624
% be. Wether or not following we only want to show how cruel it is to build up new
neuper@42478
  1625
% work on not well fundamentials.
neuper@42478
  1626
% \subparagraph{A simple equation solving,} can be set up as shown in the next
neuper@42478
  1627
% example:
neuper@42478
  1628
% 
neuper@42478
  1629
% \begin{example}
neuper@42478
  1630
% \begin{verbatim}
neuper@42478
  1631
%   
neuper@42478
  1632
%   val fmz =
neuper@42478
  1633
%     ["equality (-1 + -2 * z + 8 * z ^^^ 2 = (0::real))",
neuper@42478
  1634
%      "solveFor z",
neuper@42478
  1635
%      "solutions L"];                                    
neuper@42478
  1636
% 
neuper@42478
  1637
%   val (dI',pI',mI') =
neuper@42478
  1638
%     ("Isac", 
neuper@42478
  1639
%       ["abcFormula","degree_2","polynomial","univariate","equation"],
neuper@42478
  1640
%       ["no_met"]);\end{verbatim}
neuper@42478
  1641
% \end{example}
neuper@42478
  1642
% 
neuper@42478
  1643
% Here we want to solve the equation: $-1+-2\cdot z+8\cdot z^{2}=0$. (To give
neuper@42478
  1644
% a short overview on the commands; at first we set up the equation and tell the
neuper@42478
  1645
% machine what's the bound variable and where to store the solution. Second step 
neuper@42478
  1646
% is to define the equation type and determine if we want to use a special method
neuper@42478
  1647
% to solve this type.) Simple checks tell us that the we will get two results for
neuper@42478
  1648
% this equation and this results will be real.
neuper@42478
  1649
% So far it is easy for us and for our machine to solve, but
neuper@42478
  1650
% mentioned that a unvariate equation second order can have three different types
neuper@42478
  1651
% of solutions it is getting worth.
neuper@42478
  1652
% \subparagraph{The solving of} all this types of solutions is not yet supported.
neuper@42478
  1653
% Luckily it was needed for us; but something which has been needed in this 
neuper@42478
  1654
% context, would have been the solving of an euation looking like:
neuper@42478
  1655
% $-z^{-2}+-2\cdot z^{-1}+8=0$ which is basically the same equation as mentioned
neuper@42478
  1656
% before (remember that befor it was no problem to handle for the machine) but
neuper@42478
  1657
% now, after a simple equivalent transformation, we are not able to solve
neuper@42478
  1658
% it anymore.
neuper@42478
  1659
% \subparagraph{Error messages} we get when we try to solve something like upside
neuper@42478
  1660
% were very confusing and also leads us to no special hint about a problem.
neuper@42478
  1661
% \par The fault behind is, that we have no well error handling on one side and
neuper@42478
  1662
% no sufficient formed equation solving on the other side. This two facts are
neuper@42478
  1663
% making the implemention of new material very difficult.
neuper@42478
  1664
% 
neuper@42478
  1665
% \subsection{Formalization of missing knowledge in Isabelle}
neuper@42478
  1666
% 
neuper@42478
  1667
% \paragraph{A problem} behind is the mechanization of mathematic
neuper@42478
  1668
% theories in TP-bases languages. There is still a huge gap between
neuper@42478
  1669
% these algorithms and this what we want as a solution - in Example
neuper@42478
  1670
% Signal Processing. 
neuper@42478
  1671
% 
neuper@42478
  1672
% \vbox{
neuper@42478
  1673
%   \begin{example}
neuper@42478
  1674
%     \label{eg:gap}
neuper@42478
  1675
%     \[
neuper@42478
  1676
%       X\cdot(a+b)+Y\cdot(c+d)=aX+bX+cY+dY
neuper@42478
  1677
%     \]
neuper@42478
  1678
%     {\small\textit{
neuper@42478
  1679
%       \noindent A very simple example on this what we call gap is the
neuper@42478
  1680
% simplification above. It is needles to say that it is correct and also
neuper@42478
  1681
% Isabelle for fills it correct - \emph{always}. But sometimes we don't
neuper@42478
  1682
% want expand such terms, sometimes we want another structure of
neuper@42478
  1683
% them. Think of a problem were we now would need only the coefficients
neuper@42478
  1684
% of $X$ and $Y$. This is what we call the gap between mechanical
neuper@42478
  1685
% simplification and the solution.
neuper@42478
  1686
%     }}
neuper@42478
  1687
%   \end{example}
neuper@42478
  1688
% }
neuper@42478
  1689
% 
neuper@42478
  1690
% \paragraph{We are not able to fill this gap,} until we have to live
neuper@42478
  1691
% with it but first have a look on the meaning of this statement:
neuper@42478
  1692
% Mechanized math starts from mathematical models and \emph{hopefully}
neuper@42478
  1693
% proceeds to match physics. Academic engineering starts from physics
neuper@42478
  1694
% (experimentation, measurement) and then proceeds to mathematical
neuper@42478
  1695
% modeling and formalization. The process from a physical observance to
neuper@42478
  1696
% a mathematical theory is unavoidable bound of setting up a big
neuper@42478
  1697
% collection of standards, rules, definition but also exceptions. These
neuper@42478
  1698
% are the things making mechanization that difficult.
neuper@42478
  1699
% 
neuper@42478
  1700
% \vbox{
neuper@42478
  1701
%   \begin{example}
neuper@42478
  1702
%     \label{eg:units}
neuper@42478
  1703
%     \[
neuper@42478
  1704
%       m,\ kg,\ s,\ldots
neuper@42478
  1705
%     \]
neuper@42478
  1706
%     {\small\textit{
neuper@42478
  1707
%       \noindent Think about some units like that one's above. Behind
neuper@42478
  1708
% each unit there is a discerning and very accurate definition: One
neuper@42478
  1709
% Meter is the distance the light travels, in a vacuum, through the time
neuper@42478
  1710
% of 1 / 299.792.458 second; one kilogram is the weight of a
neuper@42478
  1711
% platinum-iridium cylinder in paris; and so on. But are these
neuper@42478
  1712
% definitions usable in a computer mechanized world?!
neuper@42478
  1713
%     }}
neuper@42478
  1714
%   \end{example}
neuper@42478
  1715
% }
neuper@42478
  1716
% 
neuper@42478
  1717
% \paragraph{A computer} or a TP-System builds on programs with
neuper@42478
  1718
% predefined logical rules and does not know any mathematical trick
neuper@42478
  1719
% (follow up example \ref{eg:trick}) or recipe to walk around difficult
neuper@42478
  1720
% expressions. 
neuper@42478
  1721
% 
neuper@42478
  1722
% \vbox{
neuper@42478
  1723
%   \begin{example}
neuper@42478
  1724
%     \label{eg:trick}
neuper@42478
  1725
%   \[ \frac{1}{j\omega}\cdot\left(e^{-j\omega}-e^{j3\omega}\right)= \]
neuper@42478
  1726
%   \[ \frac{1}{j\omega}\cdot e^{-j2\omega}\cdot\left(e^{j\omega}-e^{-j\omega}\right)=
neuper@42478
  1727
%      \frac{1}{\omega}\, e^{-j2\omega}\cdot\colorbox{lgray}{$\frac{1}{j}\,\left(e^{j\omega}-e^{-j\omega}\right)$}= \]
neuper@42478
  1728
%   \[ \frac{1}{\omega}\, e^{-j2\omega}\cdot\colorbox{lgray}{$2\, sin(\omega)$} \]
neuper@42478
  1729
%     {\small\textit{
neuper@42478
  1730
%       \noindent Sometimes it is also useful to be able to apply some
neuper@42478
  1731
% \emph{tricks} to get a beautiful and particularly meaningful result,
neuper@42478
  1732
% which we are able to interpret. But as seen in this example it can be
neuper@42478
  1733
% hard to find out what operations have to be done to transform a result
neuper@42478
  1734
% into a meaningful one.
neuper@42478
  1735
%     }}
neuper@42478
  1736
%   \end{example}
neuper@42478
  1737
% }
neuper@42478
  1738
% 
neuper@42478
  1739
% \paragraph{The only possibility,} for such a system, is to work
neuper@42478
  1740
% through its known definitions and stops if none of these
neuper@42478
  1741
% fits. Specified on Signal Processing or any other application it is
neuper@42478
  1742
% often possible to walk through by doing simple creases. This creases
neuper@42478
  1743
% are in general based on simple math operational but the challenge is
neuper@42478
  1744
% to teach the machine \emph{all}\footnote{Its pride to call it
neuper@42478
  1745
% \emph{all}.} of them. Unfortunately the goal of TP Isabelle is to
neuper@42478
  1746
% reach a high level of \emph{all} but it in real it will still be a
neuper@42478
  1747
% survey of knowledge which links to other knowledge and {{\sisac}{}} a
neuper@42478
  1748
% trainer and helper but no human compensating calculator. 
neuper@42478
  1749
% \par
neuper@42478
  1750
% {{{\sisac}{}}} itself aims to adds \emph{Algorithmic Knowledge} (formal
neuper@42478
  1751
% specifications of problems out of topics from Signal Processing, etc.)
neuper@42478
  1752
% and \emph{Application-oriented Knowledge} to the \emph{deductive} axis of
neuper@42478
  1753
% physical knowledge. The result is a three-dimensional universe of
neuper@42478
  1754
% mathematics seen in Figure~\ref{fig:mathuni}.
neuper@42478
  1755
% 
neuper@42478
  1756
% \begin{figure}
neuper@42478
  1757
%   \begin{center}
neuper@42478
  1758
%     \includegraphics{fig/universe}
neuper@42478
  1759
%     \caption{Didactic ``Math-Universe'': Algorithmic Knowledge (Programs) is
neuper@42478
  1760
%              combined with Application-oriented Knowledge (Specifications) and Deductive Knowledge (Axioms, Definitions, Theorems). The Result
neuper@42478
  1761
%              leads to a three dimensional math universe.\label{fig:mathuni}}
neuper@42478
  1762
%   \end{center}
neuper@42478
  1763
% \end{figure}
neuper@42478
  1764
% 
neuper@42478
  1765
% %WN Deine aktuelle Benennung oben wird Dir kein Fachmann abnehmen;
neuper@42478
  1766
% %WN bitte folgende Bezeichnungen nehmen:
neuper@42478
  1767
% %WN 
neuper@42478
  1768
% %WN axis 1: Algorithmic Knowledge (Programs)
neuper@42478
  1769
% %WN axis 2: Application-oriented Knowledge (Specifications)
neuper@42478
  1770
% %WN axis 3: Deductive Knowledge (Axioms, Definitions, Theorems)
neuper@42478
  1771
% %WN 
neuper@42478
  1772
% %WN und bitte die R"ander von der Grafik wegschneiden (was ich f"ur *.pdf
neuper@42478
  1773
% %WN nicht hinkriege --- weshalb ich auch die eJMT-Forderung nicht ganz
neuper@42478
  1774
% %WN verstehe, separierte PDFs zu schicken; ich w"urde *.png schicken)
neuper@42478
  1775
% 
neuper@42478
  1776
% %JR Ränder und beschriftung geändert. Keine Ahnung warum eJMT sich pdf's
neuper@42478
  1777
% %JR wünschen, würde ebenfalls png oder ähnliches verwenden, aber wenn pdf's
neuper@42478
  1778
% %JR gefordert werden WN2...
neuper@42478
  1779
% %WN2 meiner Meinung nach hat sich eJMT unklar ausgedr"uckt (z.B. kann
neuper@42478
  1780
% %WN2 man meines Wissens pdf-figures nicht auf eine bestimmte Gr"osse
neuper@42478
  1781
% %WN2 zusammenschneiden um die R"ander weg zu bekommen)
neuper@42478
  1782
% %WN2 Mein Vorschlag ist, in umserem tex-file bei *.png zu bleiben und
neuper@42478
  1783
% %WN2 png + pdf figures mitzuschicken.
neuper@42478
  1784
% 
neuper@42478
  1785
% \subsection{Notes on Problems with Traditional Notation}
neuper@42478
  1786
% 
neuper@42478
  1787
% \paragraph{During research} on these topic severely problems on
neuper@42478
  1788
% traditional notations have been discovered. Some of them have been
neuper@42478
  1789
% known in computer science for many years now and are still unsolved,
neuper@42478
  1790
% one of them aggregates with the so called \emph{Lambda Calculus},
neuper@42478
  1791
% Example~\ref{eg:lamda} provides a look on the problem that embarrassed
neuper@42478
  1792
% us.
neuper@42478
  1793
% 
neuper@42478
  1794
% \vbox{
neuper@42478
  1795
%   \begin{example}
neuper@42478
  1796
%     \label{eg:lamda}
neuper@42478
  1797
% 
neuper@42478
  1798
%   \[ f(x)=\ldots\;  \quad R \rightarrow \quad R \]
neuper@42478
  1799
% 
neuper@42478
  1800
% 
neuper@42478
  1801
%   \[ f(p)=\ldots\;  p \in \quad R \]
neuper@42478
  1802
% 
neuper@42478
  1803
%     {\small\textit{
neuper@42478
  1804
%       \noindent Above we see two equations. The first equation aims to
neuper@42478
  1805
% be a mapping of an function from the reel range to the reel one, but
neuper@42478
  1806
% when we change only one letter we get the second equation which
neuper@42478
  1807
% usually aims to insert a reel point $p$ into the reel function. In
neuper@42478
  1808
% computer science now we have the problem to tell the machine (TP) the
neuper@42478
  1809
% difference between this two notations. This Problem is called
neuper@42478
  1810
% \emph{Lambda Calculus}.
neuper@42478
  1811
%     }}
neuper@42478
  1812
%   \end{example}
neuper@42478
  1813
% }
neuper@42478
  1814
% 
neuper@42478
  1815
% \paragraph{An other problem} is that terms are not full simplified in
neuper@42478
  1816
% traditional notations, in {{\sisac}} we have to simplify them complete
neuper@42478
  1817
% to check weather results are compatible or not. in e.g. the solutions
neuper@42478
  1818
% of an second order linear equation is an rational in {{\sisac}} but in
neuper@42478
  1819
% tradition we keep fractions as long as possible and as long as they
neuper@42478
  1820
% aim to be \textit{beautiful} (1/8, 5/16,...).
neuper@42478
  1821
% \subparagraph{The math} which should be mechanized in Computer Theorem
neuper@42478
  1822
% Provers (\emph{TP}) has (almost) a problem with traditional notations
neuper@42478
  1823
% (predicate calculus) for axioms, definitions, lemmas, theorems as a
neuper@42478
  1824
% computer program or script is not able to interpret every Greek or
neuper@42478
  1825
% Latin letter and every Greek, Latin or whatever calculations
neuper@42478
  1826
% symbol. Also if we would be able to handle these symbols we still have
neuper@42478
  1827
% a problem to interpret them at all. (Follow up \hbox{Example
neuper@42478
  1828
% \ref{eg:symbint1}})
neuper@42478
  1829
% 
neuper@42478
  1830
% \vbox{
neuper@42478
  1831
%   \begin{example}
neuper@42478
  1832
%     \label{eg:symbint1}
neuper@42478
  1833
%     \[
neuper@42478
  1834
%       u\left[n\right] \ \ldots \ unitstep
neuper@42478
  1835
%     \]
neuper@42478
  1836
%     {\small\textit{
neuper@42478
  1837
%       \noindent The unitstep is something we need to solve Signal
neuper@42478
  1838
% Processing problem classes. But in {{{\sisac}{}}} the rectangular
neuper@42478
  1839
% brackets have a different meaning. So we abuse them for our
neuper@42478
  1840
% requirements. We get something which is not defined, but usable. The
neuper@42478
  1841
% Result is syntax only without semantic.
neuper@42478
  1842
%     }}
neuper@42478
  1843
%   \end{example}
neuper@42478
  1844
% }
neuper@42478
  1845
% 
neuper@42478
  1846
% In different problems, symbols and letters have different meanings and
neuper@42478
  1847
% ask for different ways to get through. (Follow up \hbox{Example
neuper@42478
  1848
% \ref{eg:symbint2}}) 
neuper@42478
  1849
% 
neuper@42478
  1850
% \vbox{
neuper@42478
  1851
%   \begin{example}
neuper@42478
  1852
%     \label{eg:symbint2}
neuper@42478
  1853
%     \[
neuper@42478
  1854
%       \widehat{\ }\ \widehat{\ }\ \widehat{\ } \  \ldots \  exponent
neuper@42478
  1855
%     \]
neuper@42478
  1856
%     {\small\textit{
neuper@42478
  1857
%     \noindent For using exponents the three \texttt{widehat} symbols
neuper@42478
  1858
% are required. The reason for that is due the development of
neuper@42478
  1859
% {{{\sisac}{}}} the single \texttt{widehat} and also the double were
neuper@42478
  1860
% already in use for different operations.
neuper@42478
  1861
%     }}
neuper@42478
  1862
%   \end{example}
neuper@42478
  1863
% }
neuper@42478
  1864
% 
neuper@42478
  1865
% \paragraph{Also the output} can be a problem. We are familiar with a
neuper@42478
  1866
% specified notations and style taught in university but a computer
neuper@42478
  1867
% program has no knowledge of the form proved by a professor and the
neuper@42478
  1868
% machines themselves also have not yet the possibilities to print every
neuper@42478
  1869
% symbol (correct) Recent developments provide proofs in a human
neuper@42478
  1870
% readable format but according to the fact that there is no money for
neuper@42478
  1871
% good working formal editors yet, the style is one thing we have to
neuper@42478
  1872
% live with.
neuper@42478
  1873
% 
neuper@42478
  1874
% \section{Problems rising out of the Development Environment}
neuper@42478
  1875
% 
neuper@42478
  1876
% fehlermeldungen! TODO
jan@42463
  1877
neuper@42492
  1878
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\end{verbatim}
neuper@42492
  1879
neuper@42464
  1880
\section{Conclusion}\label{conclusion}
jan@42463
  1881
jan@42512
  1882
%JR obvious
jan@42512
  1883
jan@42512
  1884
%This paper gives a first experience report about programming with a
jan@42512
  1885
%TP-based programming language.
jan@42512
  1886
jan@42512
  1887
A brief re-introduction of the novel kind of programming
neuper@42492
  1888
language by example of the {\sisac}-prototype makes the paper
neuper@42492
  1889
self-contained. The main section describes all the main concepts
neuper@42492
  1890
involved in TP-based programming and all the sub-tasks concerning
neuper@42492
  1891
respective implementation: mechanisation of mathematics and domain
neuper@42492
  1892
modelling, implementation of term rewriting systems for the
neuper@42492
  1893
rewriting-engine, formal (implicit) specification of the problem to be
neuper@42507
  1894
(explicitly) described by the program, implementation of the many components
neuper@42492
  1895
required for Lucas-Interpretation and finally implementation of the
neuper@42492
  1896
program itself.
neuper@42492
  1897
neuper@42492
  1898
The many concepts and sub-tasks involved in programming require a
neuper@42492
  1899
comprehensive workflow; first experiences with the workflow as
neuper@42492
  1900
supported by the present prototype are described as well: Isabelle +
neuper@42492
  1901
Isar + jEdit provide appropriate components for establishing an
neuper@42492
  1902
efficient development environment integrating computation and
neuper@42492
  1903
deduction. However, the present state of the prototype is far off a
neuper@42492
  1904
state appropriate for wide-spread use: the prototype of the program
neuper@42492
  1905
language lacks expressiveness and elegance, the prototype of the
neuper@42492
  1906
development environment is hardly usable: error messages still address
neuper@42492
  1907
the developer of the prototype's interpreter rather than the
neuper@42492
  1908
application programmer, implementation of the many settings for the
neuper@42492
  1909
Lucas-Interpreter is cumbersome.
neuper@42492
  1910
neuper@42492
  1911
From these experiences a successful proof of concept can be concluded:
neuper@42492
  1912
programming arbitrary problems from engineering sciences is possible,
neuper@42492
  1913
in principle even in the prototype. Furthermore the experiences allow
neuper@42492
  1914
to conclude detailed requirements for further development:
neuper@42492
  1915
\begin{itemize}
neuper@42492
  1916
\item Clarify underlying logics such that programming is smoothly
neuper@42492
  1917
integrated with verification of the program; the post-condition should
neuper@42492
  1918
be proved more or less automatically, otherwise working engineers
neuper@42492
  1919
would not encounter such programming.
neuper@42492
  1920
\item Combine the prototype's programming language with Isabelle's
neuper@42492
  1921
powerful function package and probably with more of SML's
neuper@42492
  1922
pattern-matching features; include parallel execution on multi-core
jan@42511
  1923
machines into the language design.
neuper@42492
  1924
\item Extend the prototype's Lucas-Interpreter such that it also
neuper@42492
  1925
handles functions defined by use of Isabelle's functions package; and
neuper@42492
  1926
generalize Isabelle's code generator such that efficient code for the
neuper@42507
  1927
whole definition of the programming language can be generated (for
neuper@42492
  1928
multi-core machines).
neuper@42492
  1929
\item Develop an efficient development environment with
neuper@42492
  1930
integration of programming and proving, with management not only of
neuper@42492
  1931
Isabelle theories, but also of large collections of specifications and
neuper@42492
  1932
of programs.
neuper@42492
  1933
\end{itemize} 
neuper@42492
  1934
Provided successful accomplishment, these points provide distinguished
jan@42511
  1935
components for virtual workbenches appealing to practitioner of
neuper@42492
  1936
engineering in the near future.
neuper@42492
  1937
jan@42511
  1938
\medskip Interactive course material, as addressed by the title, then
neuper@42507
  1939
can comprise step-wise problem solving created as a side-effect of a
neuper@42507
  1940
TP-based program: Lucas-Interpretation not only provides an
neuper@42507
  1941
interactive programming environment, Lucas-Interpretation also can
jan@42511
  1942
provide TP-based services for a flexible dialogue component with
neuper@42507
  1943
adaptive user guidance for independent and inquiry-based learning.
neuper@42492
  1944
jan@42463
  1945
jan@42463
  1946
\bibliographystyle{alpha}
neuper@42507
  1947
{\small\bibliography{references}}
jan@42463
  1948
jan@42463
  1949
\end{document}