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