Initial population.

1 files changed, 69 insertions, 0 deletions

diff --git a/src/input/f04maf.input.pamphlet b/src/input/f04maf.input.pamphlet
new file mode 100644
index 00000000..b5e02a31
--- /dev/null
+++ b/src/input/f04maf.input.pamphlet
@@ -0,0 +1,69 @@
+\documentclass{article}
+\usepackage{axiom}
+\begin{document}
+\title{\$SPAD/src/input f04maf.input}
+\author{The Axiom Team}
+\maketitle
+\begin{abstract}
+\end{abstract}
+\eject
+\tableofcontents
+\eject
+\section{License}
+<<license>>=
+--Copyright The Numerical Algorithms Group Limited 1994.
+@
+<<*>>=
+<<license>>
+)cl all -- first run f01maf
+showArrayValues true
+showScalarValues true
+
+n := 4;
+m := 4;
+nn := n*m;
+licn := 90;
+lirn := 50;
+a := new(1,licn,0.0)$Matrix DoubleFloat;
+irn := new(1,lirn,0)$Matrix Integer;
+icn := new(1,licn,0)$Matrix Integer;
+k := 0;
+for i in 1..nn repeat
+  k := k + 1
+  a(1,k) := 1
+  irn(1,k) := i
+  icn(1,k) := i
+for i in 1..n repeat
+  for j in 1..(m-1) repeat
+    k := k + 1
+    jj := (i-1)*n + j
+    a(1,k) := -0.25
+    irn(1,k) := jj
+    icn(1,k) := jj + 1
+for i in (n+1)..nn repeat
+  k := k + 1
+  a(1,k) := -0.25
+  irn(1,k) := i-n
+  icn(1,k) := i
+nz := k;
+abort := [true,true,true];
+droptl := 0.1;
+densw := 0.8;
+ifail := 111;
+res1:=f01maf(nn,nz,licn,lirn,abort,a,irn,icn,droptl,densw,ifail)
+a := res1."avals"
+irn := res1."irn"
+icn := res1."icn"
+b := [[0.5,0.25,0.25,0.5,0.25,0,0,0.25,0.25,0,0,0.25,0.5,0.25,0.25,0.5]]
+acc := [[0.0001,0]]
+noits := [[50,0]]
+wkeep := res1."wkeep"
+ikeep := res1."ikeep"
+inform := res1."inform"
+result:=f04maf(nn,nz,a,licn,irn,lirn,icn,wkeep,ikeep,inform,b,acc,noits,ifail)
+@
+\eject
+\begin{thebibliography}{99}
+\bibitem{1} nothing
+\end{thebibliography}
+\end{document}