aboutsummaryrefslogtreecommitdiff
path: root/src/interp/nag-f01.boot.pamphlet
diff options
context:
space:
mode:
authordos-reis <gdr@axiomatics.org>2007-08-14 05:14:52 +0000
committerdos-reis <gdr@axiomatics.org>2007-08-14 05:14:52 +0000
commitab8cc85adde879fb963c94d15675783f2cf4b183 (patch)
treec202482327f474583b750b2c45dedfc4e4312b1d /src/interp/nag-f01.boot.pamphlet
downloadopen-axiom-ab8cc85adde879fb963c94d15675783f2cf4b183.tar.gz
Initial population.
Diffstat (limited to 'src/interp/nag-f01.boot.pamphlet')
-rw-r--r--src/interp/nag-f01.boot.pamphlet2252
1 files changed, 2252 insertions, 0 deletions
diff --git a/src/interp/nag-f01.boot.pamphlet b/src/interp/nag-f01.boot.pamphlet
new file mode 100644
index 00000000..0fa5c044
--- /dev/null
+++ b/src/interp/nag-f01.boot.pamphlet
@@ -0,0 +1,2252 @@
+\documentclass{article}
+\usepackage{axiom}
+\begin{document}
+\title{\$SPAD/src/interp nag-f01.boot}
+\author{The Axiom Team}
+\maketitle
+\begin{abstract}
+\end{abstract}
+\eject
+\tableofcontents
+\eject
+\section{License}
+<<license>>=
+-- Copyright (c) 1991-2002, The Numerical ALgorithms Group Ltd.
+-- All rights reserved.
+--
+-- Redistribution and use in source and binary forms, with or without
+-- modification, are permitted provided that the following conditions are
+-- met:
+--
+-- - Redistributions of source code must retain the above copyright
+-- notice, this list of conditions and the following disclaimer.
+--
+-- - Redistributions in binary form must reproduce the above copyright
+-- notice, this list of conditions and the following disclaimer in
+-- the documentation and/or other materials provided with the
+-- distribution.
+--
+-- - Neither the name of The Numerical ALgorithms Group Ltd. nor the
+-- names of its contributors may be used to endorse or promote products
+-- derived from this software without specific prior written permission.
+--
+-- THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS
+-- IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+-- TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
+-- PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
+-- OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
+-- EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
+-- PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
+-- PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
+-- LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
+-- NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
+-- SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+@
+<<*>>=
+<<license>>
+
+f01brf() ==
+ htInitPage("F01BRF - LU factorization of real sparse matrix",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float)))
+ (text . "\windowlink{Manual Page}{manpageXXf01brf} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01brf| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "Factorizes a real sparse matrix A of order n. The routine forms ")
+ (text . "the {\it LU} factorization of the entire matrix, or ,")
+ (text . "optionally, first permutes the matrix to block lower ")
+ (text . "triangular form and then only factorizes the diagonal block. ")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Enter the order {\em n} of the matrix A ")
+ (text . "\htbitmap{great=} 1:")
+ (text . "\newline\tab{2} ")
+ (bcStrings (8 6 n PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Number of non-zero elements {\it nz}:")
+ (text . "\tab{32} \menuitemstyle{}\tab{34}")
+ (text . "{\it pivot}:")
+ (text . "\newline \tab{2} ")
+ (bcStrings (8 15 nz PI))
+ (text . "\tab{34} ")
+ (bcStrings (8 "0.1" pivot PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Dimension of A & ICN {\it licn}:")
+ (text . "\tab{32} \menuitemstyle{}\tab{34}")
+ (text . "Dimension of IRN {\it lirn}:")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 150 licn PI))
+ (text . "\tab{34} ")
+ (bcStrings (6 75 lirn PI))
+ (text . "\blankline")
+ (text . "\menuitemstyle{}\tab{2} Grow value:")
+ (radioButtons grow
+ ("" " True" gr_true)
+ ("" " False" gr_false))
+ (text . "\blankline")
+ (text . "\menuitemstyle{}\tab{2} Lblock value:")
+ (radioButtons lblock
+ ("" " True" lb_true)
+ ("" " False" lb_false))
+ (text . "\blankline ")
+ (text . "\newline \tab{2} ")
+ (text . "Ifail is input in three components: ")
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "{\it a} ")
+ (radioButtons afail
+ ("" " 0, hard failure" azero)
+ ("" " 1, soft failure" aone))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "{\it b} ")
+ (radioButtons bfail
+ ("" " 1, print error messages" bone)
+ ("" " 0, suppress error messages" bzero))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "{\it c} ")
+ (radioButtons cfail
+ ("" " 1, print warning messages" cone)
+ ("" " 0, suppress warning messages" czero)))
+ htMakeDoneButton('"Continue", 'f01brfSolve)
+ htShowPage()
+
+f01brfSolve htPage ==
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ nz :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nz)
+ objValUnwrap htpLabelSpadValue(htPage, 'nz)
+ licn :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'licn)
+ objValUnwrap htpLabelSpadValue(htPage, 'licn)
+ lirn :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lirn)
+ objValUnwrap htpLabelSpadValue(htPage, 'lirn)
+ pivot := htpLabelInputString(htPage, 'pivot)
+ gr := htpButtonValue(htPage,'grow)
+ grow :=
+ gr = 'gr_true => '"true"
+ '"false"
+ lb := htpButtonValue(htPage,'lblock)
+ lblock :=
+ lb = 'lb_true => '"true"
+ '"false"
+ aerror := htpButtonValue(htPage,'afail)
+ afail :=
+ aerror = 'azero => '0
+ '1
+ berror := htpButtonValue(htPage,'bfail)
+ bfail :=
+ berror = 'bone => '1
+ '0
+ cerror := htpButtonValue(htPage,'cfail)
+ cfail :=
+ cerror = 'cone => '1
+ '0
+ ifail := 100*cfail + 10*bfail + afail
+ ((n = '6 and nz = '15) and (licn = '150 and lirn = '75))
+ => f01brfDefaultSolve(htPage,pivot,grow,lblock,ifail)
+ labelList :=
+ "append"/[f(i) for i in 1..nz] where f(i) ==
+ prefix := ('"\newline \tab{2} ")
+ anam := INTERN STRCONC ('"a",STRINGIMAGE i)
+ mid := ('"\tab{32} ")
+ rnam := INTERN STRCONC ('"irn",STRINGIMAGE i)
+ end := ('"\tab{42} ")
+ cnam := INTERN STRCONC ('"icn",STRINGIMAGE i)
+ [['text,:prefix],['bcStrings,[8, 0.0, anam, 'F]],
+ ['text,:mid],['bcStrings,[4, 0, rnam, 'PI]],
+ ['text,:end],['bcStrings,[4, 0, cnam, 'PI]]]
+ abortList :=
+ [['bcStrings,[6, '"true", 'abortone, 'EM]],
+ ['bcStrings,[6, '"true", 'aborttwo, 'EM]],
+ ['bcStrings,[6, '"false", 'abortthree, 'EM]],
+ ['bcStrings,[6, '"true", 'abortfour, 'EM]]]
+ prefix := ('"\blankline \menuitemstyle{}\tab{2} Abort: ")
+ abortList := [['text,:prefix],:abortList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain EM (EmptyMode))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :labelList,:abortList]
+ page := htInitPage("F01BRF - LU factorization of real sparse matrix",nil)
+ htSay '"\menuitemstyle{}\tab{2} Non-zero elements of A: "
+ htSay '"\tab{30} \menuitemstyle{}\tab{32} Row: "
+ htSay '"\tab{40} \menuitemstyle{}\tab{42} Column: "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01brfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'nz,nz)
+ htpSetProperty(page,'licn,licn)
+ htpSetProperty(page,'lirn,lirn)
+ htpSetProperty(page,'pivot,pivot)
+ htpSetProperty(page,'grow,grow)
+ htpSetProperty(page,'lblock,lblock)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01brfDefaultSolve(htPage,pivot,grow,lblock,ifail) ==
+ n := '6
+ nz := '15
+ licn := '150
+ lirn := '75
+ page := htInitPage("F01BRF - LU factorization of real sparse matrix",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain PI (Positive Integer))
+ (isDomain EM $EmptyMode)
+ (isDomain F (Float)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Non-zero elements of A: ")
+ (text . "\tab{30} \menuitemstyle{}\tab{32} Row: ")
+ (text . "\tab{40} \menuitemstyle{}\tab{42} Column: ")
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "5.0" a1 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 1 irn1 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 1 icn1 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "2.0" a2 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 2 irn2 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 2 icn2 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-1.0" a3 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 2 irn3 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 3 icn3 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "2.0" a4 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 2 irn4 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 4 icn4 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "3.0" a5 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 3 irn5 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 3 icn5 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-2.0" a6 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 4 irn6 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 1 icn6 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a7 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 4 irn7 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 4 icn7 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a8 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 4 irn8 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 5 icn8 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-1.0" a9 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 5 irn9 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 1 icn9 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-1.0" a10 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 5 irn10 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 4 icn10 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "2.0" a11 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 5 irn11 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 5 icn11 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-3.0" a12 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 5 irn12 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 6 icn12 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-1.0" a13 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 6 irn13 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 1 icn13 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-1.0" a14 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 6 irn14 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 2 icn14 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "6.0" a15 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 6 irn15 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 6 icn15 PI))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{} \tab{2} Abort :")
+ (bcStrings (8 "true" abort_one EM))
+ (bcStrings (8 "true" abort_two EM))
+ (bcStrings (8 "false" abort_three EM))
+ (bcStrings (8 "true" abort_four EM)))
+ htMakeDoneButton('"Continue",'f01brfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'nz,nz)
+ htpSetProperty(page,'licn,licn)
+ htpSetProperty(page,'lirn,lirn)
+ htpSetProperty(page,'pivot,pivot)
+ htpSetProperty(page,'grow,grow)
+ htpSetProperty(page,'lblock,lblock)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01brfGen htPage ==
+ n := htpProperty(htPage,'n)
+ nz := htpProperty(htPage,'nz)
+ licn := htpProperty(htPage,'licn)
+ lirn := htpProperty(htPage,'lirn)
+ pivot := htpProperty(htPage,'pivot)
+ grow := htpProperty(htPage,'grow)
+ lblock := htpProperty(htPage,'lblock)
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ for i in 1..4 repeat
+ abort := STRCONC((first y).1," ")
+ y := rest y
+ abortList := [abort,:abortList]
+ astring := bcwords2liststring abortList
+ while y repeat
+ end := STRCONC ((first y).1," ")
+ y := rest y
+ mid := STRCONC ((first y).1," ")
+ y := rest y
+ top := STRCONC ((first y).1," ")
+ y := rest y
+ cList := [end,:cList]
+ rList := [mid,:rList]
+ matList := [top,:matList]
+ for i in 1..(licn-nz) repeat
+ cList := [:cList,'"0 "]
+ matList := [:matList,'"0 "]
+ for i in 1..(lirn-nz) repeat
+ rList := [:rList,'"0 "]
+ cstring := bcwords2liststring cList
+ rstring := bcwords2liststring rList
+ matstring := bcwords2liststring matList
+ prefix := STRCONC('"f01brf(",STRINGIMAGE n,", ",STRINGIMAGE nz,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE licn,", ",STRINGIMAGE lirn,", ",pivot)
+ prefix := STRCONC(prefix,", ",lblock,", ",grow,", ",astring,",[",matstring)
+ prefix := STRCONC(prefix,"],[",rstring,"],[",cstring,"], ")
+ linkGen STRCONC(prefix,STRINGIMAGE ifail,")")
+
+f01bsf() ==
+ htInitPage("F01BSF - LU factorization of real sparse matrix with known sparsity pattern",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float)))
+ (text . "\windowlink{Manual Page}{manpageXXf01bsf} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01bsf| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "Factorizes a real sparse matrix A of order n using the pivotal ")
+ (text . "sequence previously obtained by F01BRF when a matrix of the ")
+ (text . "same sparsity pattern was factorized. ")
+ (text . "\blankline ")
+ (text . "Read the input file to see the example program. ")
+ (text . "\spadpaste{)read f01bsf \bound{s0}} ")
+ (text . "\blankline")
+ (text . "\newline "))
+ htShowPage()
+
+f01maf() ==
+ htInitPage("F01MAF - \htbitmap{llt} factorization of real sparse symmetric positive-definite matrix",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float)))
+ (text . "\windowlink{Manual Page}{manpageXXf01maf} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01maf| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "Computes an incomplete Cholesky factorization of a real ")
+ (text . "sparse symmetric positive-definite matrix A of order n. ")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Enter the order {\em n} of the matrix A ")
+ (text . "\htbitmap{great=} 1:")
+ (text . "\newline\tab{2} ")
+ (bcStrings (8 16 n PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Number of non-zero elements {\it nz}:")
+ (text . "\newline \tab{2} ")
+ (bcStrings (8 40 nz PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Dimension of A & ICN {\it licn}:")
+ (text . "\tab{32} \menuitemstyle{}\tab{34}")
+ (text . "Dimension of IRN {\it lirn}:")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 90 licn PI))
+ (text . "\tab{34} ")
+ (bcStrings (6 50 lirn PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Tolerance {\it droptl}: ")
+ (text . "\tab{32} \menuitemstyle{}\tab{34}")
+ (text . "{\it densw}:")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "0.1" droptl F))
+ (text . "\tab{34} ")
+ (bcStrings (6 "0.8" densw F))
+ (text . "\blankline ")
+ (text . "\newline \tab{2} ")
+ (text . "Ifail is input in three components: ")
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "{\it a} ")
+ (radioButtons afail
+ ("" " 0, hard failure" azero)
+ ("" " 1, soft failure" aone))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "{\it b} ")
+ (radioButtons bfail
+ ("" " 1, print error messages" bone)
+ ("" " 0, suppress error messages" bzero))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "{\it c} ")
+ (radioButtons cfail
+ ("" " 1, print warning messages" cone)
+ ("" " 0, suppress warning messages" czero)))
+ htMakeDoneButton('"Continue", 'f01mafSolve)
+ htShowPage()
+
+f01mafSolve htPage ==
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ nz :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nz)
+ objValUnwrap htpLabelSpadValue(htPage, 'nz)
+ licn :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'licn)
+ objValUnwrap htpLabelSpadValue(htPage, 'licn)
+ lirn :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lirn)
+ objValUnwrap htpLabelSpadValue(htPage, 'lirn)
+ aerror := htpButtonValue(htPage,'afail)
+ afail :=
+ aerror = 'azero => '0
+ '1
+ berror := htpButtonValue(htPage,'bfail)
+ bfail :=
+ berror = 'bone => '1
+ '0
+ cerror := htpButtonValue(htPage,'cfail)
+ cfail :=
+ cerror = 'cone => '1
+ '0
+ ifail := 100*cfail + 10*bfail + afail
+ droptl := htpLabelInputString(htPage, 'droptl)
+ densw := htpLabelInputString(htPage, 'densw)
+ ((n = '16 and nz = '40) and (licn = '90 and lirn = '50))
+ => f01mafDefaultSolve(htPage,droptl,densw,ifail)
+ labelList :=
+ "append"/[f(i) for i in 1..nz] where f(i) ==
+ prefix := ('"\newline \tab{2} ")
+ anam := INTERN STRCONC ('"a",STRINGIMAGE i)
+ mid := ('"\tab{32} ")
+ rnam := INTERN STRCONC ('"irn",STRINGIMAGE i)
+ end := ('"\tab{42} ")
+ cnam := INTERN STRCONC ('"icn",STRINGIMAGE i)
+ [['text,:prefix],['bcStrings,[8, 0.0, anam, 'F]],
+ ['text,:mid],['bcStrings,[4, 0, rnam, 'PI]],
+ ['text,:end],['bcStrings,[4, 0, cnam, 'PI]]]
+ abortList :=
+ [['bcStrings,[6, '"true", 'abortone, 'EM]],
+ ['bcStrings,[6, '"true", 'aborttwo, 'EM]],
+ ['bcStrings,[6, '"true", 'abortthree, 'EM]]]
+ prefix := ('"\blankline \menuitemstyle{}\tab{2} Abort: ")
+ abortList := [['text,:prefix],:abortList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain EM (EmptyMode))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :labelList,:abortList]
+ page := htInitPage("F01MAF - \htbitmap{llt} factorization of real sparse symmetric positive-definite matrix",nil)
+ htSay '"\menuitemstyle{}\tab{2} Non-zero elements of A: "
+ htSay '"\tab{30} \menuitemstyle{}\tab{32} Row: "
+ htSay '"\tab{40} \menuitemstyle{}\tab{42} Column: "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01mafGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'nz,nz)
+ htpSetProperty(page,'licn,licn)
+ htpSetProperty(page,'lirn,lirn)
+ htpSetProperty(page,'droptl,droptl)
+ htpSetProperty(page,'densw,densw)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01mafDefaultSolve(htPage,droptl,densw,ifail) ==
+ n := '16
+ nz := '40
+ licn := '90
+ lirn := '50
+ page := htInitPage("F01MAF - \htbitmap{llt} factorization of real sparse symmetric positive-definite matrix",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain PI (Positive Integer))
+ (isDomain EM $EmptyMode)
+ (isDomain F (Float)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Non-zero elements of A: ")
+ (text . "\tab{30} \menuitemstyle{}\tab{32} Row: ")
+ (text . "\tab{40} \menuitemstyle{}\tab{42} Column: ")
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a1 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 1 irn1 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 1 icn1 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a2 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 2 irn2 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 2 icn2 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a3 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 3 irn3 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 3 icn3 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a4 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 4 irn4 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 4 icn4 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a5 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 5 irn5 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 5 icn5 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a6 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 6 irn6 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 6 icn6 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a7 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 7 irn7 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 7 icn7 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a8 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 8 irn8 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 8 icn8 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a9 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 9 irn9 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 9 icn9 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a10 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 10 irn10 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 10 icn10 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a11 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 11 irn11 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 11 icn11 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a12 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 12 irn12 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 12 icn12 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a13 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 13 irn13 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 13 icn13 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a14 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 14 irn14 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 14 icn14 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a15 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 15 irn15 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 15 icn15 PI))
+ (text . "\blankline ")
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "1.0" a16 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 16 irn16 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 16 icn16 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a17 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 1 irn17 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 2 icn17 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a18 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 2 irn18 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 3 icn18 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a19 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 3 irn19 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 4 icn19 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a20 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 5 irn20 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 6 icn20 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a21 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 6 irn21 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 7 icn21 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a22 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 7 irn22 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 8 icn22 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a23 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 9 irn23 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 10 icn23 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a24 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 10 irn24 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 11 icn24 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a25 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 11 irn25 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 12 icn25 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a26 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 13 irn26 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 14 icn26 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a27 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 14 irn27 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 15 icn27 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a28 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 15 irn28 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 16 icn28 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a29 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 1 irn29 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 5 icn29 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a30 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 2 irn30 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 6 icn30 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a31 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 3 irn31 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 7 icn31 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a32 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 4 irn32 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 8 icn32 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a33 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 5 irn33 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 9 icn33 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a34 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 6 irn34 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 10 icn34 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a35 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 7 irn35 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 11 icn35 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a36 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 8 irn36 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 12 icn36 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a37 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 9 irn37 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 13 icn37 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a38 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 10 irn38 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 14 icn38 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a39 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 11 irn39 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 15 icn39 PI))
+ (text . "\newline \tab{2}")
+ (bcStrings (8 "-0.25" a40 F))
+ (text . "\tab{32} ")
+ (bcStrings (4 12 irn40 PI))
+ (text . "\tab{42} ")
+ (bcStrings (4 16 icn40 PI))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{} \tab{2} Abort :")
+ (bcStrings (8 "true" abort_one EM))
+ (bcStrings (8 "true" abort_two EM))
+ (bcStrings (8 "true" abort_three EM)))
+ htMakeDoneButton('"Continue",'f01mafGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'nz,nz)
+ htpSetProperty(page,'licn,licn)
+ htpSetProperty(page,'lirn,lirn)
+ htpSetProperty(page,'droptl,droptl)
+ htpSetProperty(page,'densw,densw)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01mafGen htPage ==
+ n := htpProperty(htPage,'n)
+ nz := htpProperty(htPage,'nz)
+ licn := htpProperty(htPage,'licn)
+ lirn := htpProperty(htPage,'lirn)
+ droptl := htpProperty(htPage,'droptl)
+ densw := htpProperty(htPage,'densw)
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ for i in 1..3 repeat
+ abort := STRCONC((first y).1," ")
+ y := rest y
+ abortList := [abort,:abortList]
+ astring := bcwords2liststring abortList
+ while y repeat
+ end := STRCONC ((first y).1," ")
+ y := rest y
+ mid := STRCONC ((first y).1," ")
+ y := rest y
+ top := STRCONC ((first y).1," ")
+ y := rest y
+ cList := [end,:cList]
+ rList := [mid,:rList]
+ matList := [top,:matList]
+ for i in 1..(licn-nz) repeat
+ cList := [:cList,'"0 "]
+ matList := [:matList,'"0 "]
+ for i in 1..(lirn-nz) repeat
+ rList := [:rList,'"0 "]
+ cstring := bcwords2liststring cList
+ rstring := bcwords2liststring rList
+ matstring := bcwords2liststring matList
+ prefix := STRCONC('"f01maf(",STRINGIMAGE n,", ",STRINGIMAGE nz,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE licn,", ",STRINGIMAGE lirn,", ")
+ prefix := STRCONC(prefix,astring,",[",matstring)
+ prefix := STRCONC(prefix,"],[",rstring,"],[",cstring,"], ",droptl,", ",densw)
+ linkGen STRCONC(prefix,", ",STRINGIMAGE ifail,")")
+
+
+
+
+f01mcf() ==
+ htInitPage("F01MCF - \htbitmap{ldlt} factorization of real symmetric positive-definite variable-bandwidth matrix",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float)))
+ (text . "\windowlink{Manual Page}{manpageXXf01mcf} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01mcf| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "Computes the Cholesky factorization of a real symmetric positive")
+ (text . "-definite variable-bandwidth matrix {\it A} or order {\it n}. ")
+ (text . "That is, {\it A = }\htbitmap{ldlt}, where {\it L} is ")
+ (text . "a unit lower triangular matrix and {\it D} is a diagonal matrix.")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Enter the order {\em n} of the matrix A ")
+ (text . "\htbitmap{great=} 1:")
+ (text . "\newline\tab{2} ")
+ (bcStrings (9 6 n PI))
+ (text . "\blankline")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "\newline Enter the number of elements: ")
+ (text . "\newline\tab{2} ")
+ (bcStrings (9 14 lal PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2}")
+ (text . "Ifail value:")
+ (radioButtons ifail
+ ("" " -1, Print error messages" minusOne)
+ ("" " 1, Suppress error messages" one)))
+ htMakeDoneButton('"Continue", 'f01mcfSolve)
+ htShowPage()
+
+f01mcfSolve htPage ==
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ lal :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lal)
+ objValUnwrap htpLabelSpadValue(htPage, 'lal)
+ error := htpButtonValue(htPage,'ifail)
+ ifail :=
+ error = 'one => '1
+ '-1
+ (n = '6 and lal = '14) => f01mcfDefaultSolve(htPage,ifail)
+ labelList :=
+ "append"/[f(i) for i in 1..lal] where f(i) ==
+ xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
+ [['bcStrings,[6, 0.0, xnam, 'F]]]
+ nrowList :=
+ "append"/[g(j) for j in 1..n] where g(j) ==
+ nam := INTERN STRCONC ('"n",STRINGIMAGE j)
+ [['bcStrings,[6, 0, nam, 'PI]]]
+ prefix := ('"\blankline \menuitemstyle{}\tab{2} {\it NROW(n)} the width ")
+ prefix := STRCONC(prefix,"of the ith row of A: \newline \tab{2} ")
+ nrowList := [['text,:prefix],:nrowList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :labelList,:nrowList]
+ page := htInitPage("F01MCF - \htbitmap{ldlt} factorization of real symmetric positive-definite variable-bandwidth matrix",nil)
+ htSay '"\menuitemstyle{}\tab{2} Elements of matrix {\it A} in row by row "
+ htSay '"order: \newline \tab{2} "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01mcfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'lal,lal)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01mcfDefaultSolve (htPage,ifail) ==
+ n := '6
+ lal := '14
+ page := htInitPage("F01MCF - \htbitmap{ldlt} factorization of real symmetric positive-definite variable-bandwidth matrix",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain PI (Positive Integer))
+ (isDomain F (Float)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Elements of matrix {\it A} in row by ")
+ (text . "row order: ")
+ (text . "\newline ")
+ (bcStrings (6 "1.0" x1 F))
+ (bcStrings (6 "2.0" x2 F))
+ (bcStrings (6 "5.0" x3 F))
+ (bcStrings (6 "3.0" x4 F))
+ (bcStrings (6 "13.0" x5 F))
+ (bcStrings (6 "16.0" x6 F))
+ (bcStrings (6 "5.0" x7 F))
+ (bcStrings (6 "14.0" x8 F))
+ (bcStrings (6 "18.0" x9 F))
+ (bcStrings (6 "8.0" x10 F))
+ (bcStrings (6 "55.0" x11 F))
+ (bcStrings (6 "24.0" x12 F))
+ (bcStrings (6 "17.0" x13 F))
+ (bcStrings (6 "77.0" x14 F))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2} {\it NROW(n)} the width of the ith row ")
+ (text . "of A: ")
+ (text . "\newline ")
+ (bcStrings (6 1 n1 PI))
+ (bcStrings (6 2 n2 PI))
+ (bcStrings (6 2 n3 PI))
+ (bcStrings (6 1 n4 PI))
+ (bcStrings (6 5 n5 PI))
+ (bcStrings (6 3 n6 PI))
+ (text . "\blankline "))
+ htMakeDoneButton('"Continue",'f01mcfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'lal,lal)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01mcfGen htPage ==
+ n := htpProperty(htPage,'n)
+ lal := htpProperty(htPage,'lal)
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ for i in 1..n repeat
+ right := STRCONC ((first y).1," ")
+ y := rest y
+ nrowList := [right,:nrowList]
+ nrowstring := bcwords2liststring nrowList
+ while y repeat
+ right := STRCONC ((first y).1," ")
+ y := rest y
+ matList := [right,:matList]
+ matstring := bcwords2liststring matList
+ prefix := STRCONC('"f01mcf(",STRINGIMAGE n,", [",matstring,"], ")
+ prefix := STRCONC(prefix,STRINGIMAGE lal,", [",nrowstring,"], ")
+ linkGen STRCONC(prefix,STRINGIMAGE ifail,")")
+
+
+f01qcf() ==
+ htInitPage('"F01QCF - QR factorization or real m by n matrix (m \htbitmap{great=} n)",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger)))
+ (text . "\windowlink{Manual Page}{manpageXXf01qcf} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01qcf| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "\newline ")
+ (text . "Finds the QR factorization of a real {\it m} by {it n} ({\it m ")
+ (text . "\htbitmap{great=} n}) matrix {\it A}, which ")
+ (text . "is factorized as \htbitmap{f01qcf}, ")
+ (text . "where {\it m > n} and {\it A = QR } when {\it m = n }, where ")
+ (text . "{\it Q} is an {\it m} by {\it m } orthogonal matrix and {\it R} ")
+ (text . "is an {\it n} by {\it n} upper triangular matrix. The {\it k}th ")
+ (text . "transformation matrix,{\it Qk}, ")
+ (text . "which is used to introduce zeros into the {\it k}th column of ")
+ (text . "{\it A}, is given in the form ")
+ (text . "\htbitmap{f01qcf1}, ")
+ (text . "where \htbitmap{f01qcf2}, ")
+ (text . "\htbitmap{f01qcf3}, ")
+ (text . "\htbitmap{zetak} is a scalar and ")
+ (text . "\htbitmap{zk} is an (m-k) element vector. ")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline Rows of matrix A, {\it m}: ")
+ (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+ (text . "Columns of matrix A, {\it n}: \newline \tab{2} ")
+ (bcStrings (6 5 m PI))
+ (text . "\tab{34} ")
+ (bcStrings (6 3 n PI))
+-- (text . "\blankline ")
+-- (text . "\newline ")
+-- (text . "\menuitemstyle{} \tab{2} ")
+-- (text . "\newline First dimension of A, {\it lda} ")
+-- (text . "\htbitmap{great=} m: ")
+-- (text . "\newline \tab{2} ")
+-- (bcStrings (6 5 lda PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Ifail value: ")
+ (radioButtons ifail
+ ("" " -1, Print error messages" minusOne)
+ ("" " 1, Suppress error messages" one)))
+ htMakeDoneButton('"Continue", 'f01qcfSolve)
+ htShowPage()
+
+f01qcfSolve htPage ==
+ m :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
+ objValUnwrap htpLabelSpadValue(htPage, 'm)
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ lda := m
+-- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lda)
+-- objValUnwrap htpLabelSpadValue(htPage, 'lda)
+ error := htpButtonValue(htPage,'ifail)
+ ifail :=
+ error = 'one => '1
+ '-1
+ (m = '5 and n = '3) => f01qcfDefaultSolve(htPage,lda,ifail)
+ matList :=
+ "append"/[f(i,n) for i in 1..lda] where f(i,n) ==
+ labelList :=
+ "append"/[g(i,j) for j in 1..n] where g(i,j) ==
+ anam := INTERN STRCONC ('"a",STRINGIMAGE i, STRINGIMAGE j)
+ [['bcStrings,[6, "0.0", anam, 'F]]]
+ prefix := ('"\newline \tab{2} ")
+ labelList := [['text,:prefix],:labelList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :matList]
+ page := htInitPage("F01QCF - QR factorization or real m by n matrix (m \htbitmap{great=} n)",nil)
+ htSay '"\newline \menuitemstyle{}\tab{2} Enter values of {\it A}: "
+ htSay '"\newline \tab{2} "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01qcfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01qcfDefaultSolve (htPage,lda,ifail) ==
+ n := '3
+ m := '5
+ page := htInitPage('"F01QCF - QR factorization or real m by n matrix (m \htbitmap{great=} n)",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain I (Integer)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of {\it A}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "2.0" a11 F))
+ (bcStrings (6 "2.5" a12 F))
+ (bcStrings (6 "2.5" a13 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "2.0" a21 F))
+ (bcStrings (6 "2.5" a22 F))
+ (bcStrings (6 "2.5" a23 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "1.6" a31 F))
+ (bcStrings (6 "-0.4" a32 F))
+ (bcStrings (6 "2.8" a33 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "2.0" a41 F))
+ (bcStrings (6 "-0.5" a42 F))
+ (bcStrings (6 "0.5" a43 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "1.2" a51 F))
+ (bcStrings (6 "-0.3" a52 F))
+ (bcStrings (6 "-2.9" a53 F)))
+ htMakeDoneButton('"Continue",'f01qcfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01qcfGen htPage ==
+ n := htpProperty(htPage,'n)
+ m := htpProperty(htPage,'m)
+-- lda := htpProperty(htPage,'lda)
+ lda := m
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ y := REVERSE y
+ for i in 1..lda repeat
+ for j in 1..n repeat
+ elm := STRCONC((first y).1," ")
+ rowList := [:rowList,elm]
+ y := rest y
+ matform := [:matform,rowList]
+ rowList := []
+ matstring := bcwords2liststring [bcwords2liststring x for x in matform]
+ prefix := STRCONC('"f01qcf(",STRINGIMAGE m,", ",STRINGIMAGE n,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE lda,", ",matstring,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE ifail,")")
+ linkGen prefix
+
+f01qdf() ==
+ htInitPage('"F01QDF - Operations with orthogonal matrices, compute {\it QB} or \htbitmap{f01qdf} after factorization by F01QCF or F01QDF",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger)))
+ (text . "\windowlink{Manual Page}{manpageXXf01qdf} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01qdf| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "\newline ")
+ (text . "Performs one of the transformations {\it B = QB or B = }")
+ (text . "\htbitmap{f01qdf}, where {\it B} is a real {\it m} ")
+ (text . "by {\it ncolb} matrix and {\it Q} is an {\it m} by {\it m} ")
+ (text . "orthogonal matrix assumed to be given by {\it Q = }")
+ (text . "\htbitmap{f01qdf1}, \htbitmap{f01qdf2} ")
+ (text . "being given in the form ")
+ (text . "\htbitmap{f01qcf1}, ")
+ (text . "where \htbitmap{f01qcf2}, ")
+ (text . "\htbitmap{f01qcf3}, ")
+ (text . "\htbitmap{zetak} is a scalar and ")
+ (text . "\htbitmap{zk} is an (m-k) element vector. ")
+ (text . "The routine is intended for use following F01QCF or F01QFF. ")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline Rows of matrix A, {\it m}: ")
+ (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+ (text . "Columns of matrix A, {\it n}: \newline \tab{2} ")
+ (bcStrings (6 5 m PI))
+ (text . "\tab{34} ")
+ (bcStrings (6 3 n PI))
+-- (text . "\blankline ")
+-- (text . "\newline ")
+-- (text . "\menuitemstyle{} \tab{2} ")
+-- (text . "\newline First dimension of A, {\it lda} ")
+-- (text . "\htbitmap{great=} m: ")
+-- (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+-- (text . "First dimension of B, {\it ldb} ")
+-- (text . "\htbitmap{great=} m: ")
+-- (text . "\newline \tab{2} ")
+-- (bcStrings (6 5 lda PI))
+-- (text . "\tab{34} ")
+-- (bcStrings (6 5 ldb PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Number of columns of matrix B {\it ncolb}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 2 ncolb PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Transformation to be performed: ")
+ (radioButtons trans
+ (" " " {\it B = QB}" no_trans)
+ (" " " {\it B =} \htbitmap{f01qdf}" trans))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Where the elements can be found: ")
+ (radioButtons wheret
+ (" " " the elements of \zeta are in A" in_a)
+ (" " " the elements of \zeta are in ZETA, returned by F01QCF/F01QFF" seperate))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Ifail value: ")
+ (radioButtons ifail
+ ("" " -1, Print error messages" minusOne)
+ ("" " 1, Suppress error messages" one)))
+ htMakeDoneButton('"Continue", 'f01qdfSolve)
+ htShowPage()
+
+f01qdfSolve htPage ==
+ m :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
+ objValUnwrap htpLabelSpadValue(htPage, 'm)
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ lda := m
+-- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lda)
+-- objValUnwrap htpLabelSpadValue(htPage, 'lda)
+ ldb := m
+-- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ldb)
+-- objValUnwrap htpLabelSpadValue(htPage, 'ldb)
+ ncolb :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ncolb)
+ objValUnwrap htpLabelSpadValue(htPage, 'ncolb)
+ operation := htpButtonValue(htPage,'trans)
+ trans :=
+ operation = 'no_trans => '"n"
+ '"t"
+ elements := htpButtonValue(htPage,'wheret)
+ wheret :=
+ elements = 'in_a => '"i"
+ '"s"
+ error := htpButtonValue(htPage,'ifail)
+ ifail :=
+ error = 'one => '1
+ '-1
+ ((m = '5 and n = '3) and ncolb = '2) => f01qdfDefaultSolve(htPage,lda,ldb,trans,wheret,ifail)
+ matList :=
+ "append"/[fa(i,n) for i in 1..lda] where fa(i,n) ==
+ labelList :=
+ "append"/[ga(i,j) for j in 1..n] where ga(i,j) ==
+ anam := INTERN STRCONC ('"a",STRINGIMAGE i, STRINGIMAGE j)
+ [['bcStrings,[6, "0.0", anam, 'F]]]
+ prefix := ('"\newline \tab{2} ")
+ labelList := [['text,:prefix],:labelList]
+ bList :=
+ "append"/[fb(i,ncolb) for i in 1..ldb] where fb(i,ncolb) ==
+ labelList :=
+ "append"/[gb(i,j) for j in 1..ncolb] where gb(i,j) ==
+ bnam := INTERN STRCONC ('"b",STRINGIMAGE i, STRINGIMAGE j)
+ [['bcStrings,[6, "0.0", bnam, 'F]]]
+ prefix := ('"\newline \tab{2} ")
+ labelList := [['text,:prefix],:labelList]
+ prefix := ("\blankline \menuitemstyle{}\tab{2} Enter values of {\it B}: ")
+ bList := [['text,:prefix],:bList]
+ zList :=
+ "append"/[fz(i) for i in 1..n] where fz(i) ==
+ znam := INTERN STRCONC ('"z",STRINGIMAGE i)
+ [['bcStrings,[6, "0.0", znam, 'F]]]
+ prefix := ("\blankline \menuitemstyle{}\tab{2} Enter values of \zeta ")
+ prefix := STRCONC(prefix,"(if required): \newline \tab{2}")
+ zList := [['text,:prefix],:zList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :matList,:bList,:zList]
+ page := htInitPage('"F01QDF - Operations with orthogonal matrices, compute {\it QB} or \htbitmap{f01qdf} after factorization by F01QCF or F01QDF",nil)
+ htSay '"\newline \menuitemstyle{}\tab{2} Enter values of {\it A}: "
+ htSay '"\newline \tab{2} "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01qdfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+-- htpSetProperty(page,'ldb,ldb)
+ htpSetProperty(page,'ncolb,ncolb)
+ htpSetProperty(page,'trans,trans)
+ htpSetProperty(page,'wheret,wheret)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01qdfDefaultSolve (htPage,lda,ldb,trans,wheret,ifail) ==
+ n := '3
+ m := '5
+ ncolb := '2
+ page := htInitPage('"F01QDF - Operations with orthogonal matrices, compute {\it QB} or \htbitmap{f01qdf} after factorization by F01QCF or F01QDF",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain I (Integer)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of {\it A}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "2.0" a11 F))
+ (bcStrings (6 "2.5" a12 F))
+ (bcStrings (6 "2.5" a13 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "2.0" a21 F))
+ (bcStrings (6 "2.5" a22 F))
+ (bcStrings (6 "2.5" a23 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "1.6" a31 F))
+ (bcStrings (6 "-0.4" a32 F))
+ (bcStrings (6 "2.8" a33 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "2.0" a41 F))
+ (bcStrings (6 "-0.5" a42 F))
+ (bcStrings (6 "0.5" a43 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "1.2" a51 F))
+ (bcStrings (6 "-0.3" a52 F))
+ (bcStrings (6 "-2.9" a53 F))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of {\it B}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "1.1" b11 F))
+ (bcStrings (6 "0.0" b12 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "0.9" b21 F))
+ (bcStrings (6 "0.0" b22 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "0.6" b31 F))
+ (bcStrings (6 "1.32" b32 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "0.0" b41 F))
+ (bcStrings (6 "1.1" b42 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 "-0.8" b51 F))
+ (bcStrings (6 "-0.26" b52 F))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of \zeta (if required): ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (10 "0.0" z1 F))
+ (bcStrings (10 "0.0" z2 F))
+ (bcStrings (10 "0.0" z3 F))
+ (text . "\blankline "))
+ htMakeDoneButton('"Continue",'f01qdfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+-- htpSetProperty(page,'ldb,ldb)
+ htpSetProperty(page,'ncolb,ncolb)
+ htpSetProperty(page,'trans,trans)
+ htpSetProperty(page,'wheret,wheret)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01qdfGen htPage ==
+ n := htpProperty(htPage,'n)
+ m := htpProperty(htPage,'m)
+-- lda := htpProperty(htPage,'lda)
+-- ldb := htpProperty(htPage,'ldb)
+ lda := m
+ ldb := m
+ ncolb := htpProperty(htPage,'ncolb)
+ trans := htpProperty(htPage,'trans)
+ wheret := htpProperty(htPage,'wheret)
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ for i in 1..n repeat
+ left := STRCONC((first y).1," ")
+ y := rest y
+ zetalist := [left,:zetalist]
+ zetastring := bcwords2liststring zetalist
+ y := REVERSE y
+ for i in 1..lda repeat
+ for j in 1..n repeat
+ elm := STRCONC((first y).1," ")
+ rowList := [:rowList,elm]
+ y := rest y
+ matform := [:matform,rowList]
+ rowList := []
+ matstring := bcwords2liststring [bcwords2liststring x for x in matform]
+ for i in 1..ldb repeat
+ for j in 1..ncolb repeat
+ elm := STRCONC((first y).1," ")
+ rowList := [:rowList,elm]
+ y := rest y
+ bform := [:bform,rowList]
+ rowList := []
+ bstring := bcwords2liststring [bcwords2liststring x for x in bform]
+ prefix := STRCONC('"f01qdf(_"",trans,"_",_"",wheret,"_",",STRINGIMAGE m,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE n,", ",matstring,", ",STRINGIMAGE lda)
+ prefix := STRCONC(prefix,",[",zetastring,"],",STRINGIMAGE ncolb,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE ldb,", ",bstring,", ",STRINGIMAGE ifail,")")
+ linkGen prefix
+
+
+f01qef() ==
+ htInitPage('"F01QEF - Operations with orthogonal matrices, form columns of {\it Q} after factorization by F01QCF or F01QFF",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger)))
+ (text . "\windowlink{Manual Page}{manpageXXf01qef} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01qef| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "\newline ")
+ (text . "Returns the first {\it ncolq} columns of the real {\it m} by ")
+ (text . "{\it n} orthogonal matrix {\it Q}, where {\it Q} is assumed ")
+ (text . "to be given by {\it Q = }\htbitmap{f01qdf1}, ")
+ (text . "\htbitmap{f01qdf2} being given in the form ")
+ (text . "\htbitmap{f01qcf1}, ")
+ (text . "where \htbitmap{f01qcf2}, ")
+ (text . "\htbitmap{f01qcf3}, ")
+ (text . "\htbitmap{zetak} is a scalar and ")
+ (text . "\htbitmap{zk} is an (m-k) element vector. ")
+ (text . "The routine is intended for use following F01QCF or F01QFF. ")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline Rows of matrix A, {\it m}: ")
+ (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+ (text . "Columns of matrix A, {\it n}: \newline \tab{2} ")
+ (bcStrings (6 5 m PI))
+ (text . "\tab{34} ")
+ (bcStrings (6 3 n PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+-- (text . "\newline First dimension of A, {\it lda} ")
+-- (text . "\htbitmap{great=} m: ")
+-- (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+ (text . "Number columns of matrix Q {\it ncolq}: ")
+ (text . "\newline \tab{2} ")
+-- (bcStrings (6 5 lda PI))
+-- (text . "\tab{34} ")
+ (bcStrings (6 5 ncolq PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Where the elements can be found: ")
+ (radioButtons wheret
+ (" " " the elements of \zeta are in ZETA, returned by F01QCF/F01QFF" subsequent)
+ (" " " the elements of \zeta are in A" initial))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Ifail value: ")
+ (radioButtons ifail
+ ("" " -1, Print error messages" minusOne)
+ ("" " 1, Suppress error messages" one)))
+ htMakeDoneButton('"Continue", 'f01qefSolve)
+ htShowPage()
+
+f01qefSolve htPage ==
+ m :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
+ objValUnwrap htpLabelSpadValue(htPage, 'm)
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ lda := m
+-- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lda)
+-- objValUnwrap htpLabelSpadValue(htPage, 'lda)
+ ncolq :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ncolq)
+ objValUnwrap htpLabelSpadValue(htPage, 'ncolq)
+ elements := htpButtonValue(htPage,'wheret)
+ wheret :=
+ elements = 'initial => '"i"
+ '"s"
+ error := htpButtonValue(htPage,'ifail)
+ ifail :=
+ error = 'one => '1
+ '-1
+ ((m = '5 and n = '3) and ncolq = '5) => f01qefDefaultSolve(htPage,lda,wheret,ifail)
+ matList :=
+ "append"/[fa(i,ncolq) for i in 1..lda] where fa(i,ncolq) ==
+ labelList :=
+ "append"/[ga(i,j) for j in 1..ncolq] where ga(i,j) ==
+ anam := INTERN STRCONC ('"a",STRINGIMAGE i, STRINGIMAGE j)
+ [['bcStrings,[7, "0.0", anam, 'F]]]
+ prefix := ('"\newline ")
+ labelList := [['text,:prefix],:labelList]
+ zList :=
+ "append"/[fz(i) for i in 1..n] where fz(i) ==
+ znam := INTERN STRCONC ('"z",STRINGIMAGE i)
+ [['bcStrings,[7, "0.0", znam, 'F]]]
+ prefix := ("\blankline \menuitemstyle{}\tab{2} Enter values of \zeta ")
+ prefix := STRCONC(prefix,"(if required): \newline ")
+ zList := [['text,:prefix],:zList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :matList,:zList]
+ page := htInitPage('"F01QEF - Operations with orthogonal matrices, form columns of {\it Q} after factorization by F01QCF or F01QFF",nil)
+ htSay '"\newline \menuitemstyle{}\tab{2} Enter values of {\it Q}: "
+ htSay '"\newline "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01qefGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+ htpSetProperty(page,'ncolq,ncolq)
+ htpSetProperty(page,'wheret,wheret)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01qefDefaultSolve (htPage,lda,wheret,ifail) ==
+ n := '3
+ m := '5
+ ncolq := '5
+ page := htInitPage('"F01QEF - Operations with orthogonal matrices, form columns of {\it Q} after factorization by F01QCF or F01QFF",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain I (Integer)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of {\it Q}")
+ (text . "(in this case returned by the default entries of F01QCF) : ")
+ (text . "\newline ")
+ (bcStrings (7 "-4.0" a11 F))
+ (bcStrings (7 "-2.0" a12 F))
+ (bcStrings (7 "-3.0" a13 F))
+ (bcStrings (7 "0.0" a14 F))
+ (bcStrings (7 "0.0" a15 F))
+ (text . "\newline ")
+ (bcStrings (7 "0.4085" a21 F))
+ (bcStrings (7 "-3.0" a22 F))
+ (bcStrings (7 "-2.0" a23 F))
+ (bcStrings (7 "0.0" a24 F))
+ (bcStrings (7 "0.0" a25 F))
+ (text . "\newline ")
+ (bcStrings (7 "0.3266" a31 F))
+ (bcStrings (7 "-0.4619" a32 F))
+ (bcStrings (7 "-4.0" a33 F))
+ (bcStrings (7 "0.0" a34 F))
+ (bcStrings (7 "0.0" a35 F))
+ (text . "\newline ")
+ (bcStrings (7 "0.4082" a41 F))
+ (bcStrings (7 "-0.5774" a42 F))
+ (bcStrings (7 "0.0" a43 F))
+ (bcStrings (7 "0.0" a44 F))
+ (bcStrings (7 "0.0" a45 F))
+ (text . "\newline ")
+ (bcStrings (7 "0.2449" a51 F))
+ (bcStrings (7 "-0.3464" a52 F))
+ (bcStrings (7 "-0.6326" a53 F))
+ (bcStrings (7 "0.0" a54 F))
+ (bcStrings (7 "0.0" a55 F))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of \zeta: ")
+ (text . "\newline ")
+ (bcStrings (10 "1.2247" z1 F))
+ (bcStrings (10 "1.1547" z2 F))
+ (bcStrings (10 "1.2649" z3 F))
+ (text . "\blankline "))
+ htMakeDoneButton('"Continue",'f01qefGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+ htpSetProperty(page,'ncolq,ncolq)
+ htpSetProperty(page,'wheret,wheret)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01qefGen htPage ==
+ n := htpProperty(htPage,'n)
+ m := htpProperty(htPage,'m)
+-- lda := htpProperty(htPage,'lda)
+ lda := m
+ ncolq := htpProperty(htPage,'ncolq)
+ wheret := htpProperty(htPage,'wheret)
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ for i in 1..n repeat
+ left := STRCONC((first y).1," ")
+ y := rest y
+ zetalist := [left,:zetalist]
+ zetastring := bcwords2liststring zetalist
+ y := REVERSE y
+ for i in 1..lda repeat
+ for j in 1..ncolq repeat
+ elm := STRCONC((first y).1," ")
+ rowList := [:rowList,elm]
+ y := rest y
+ matform := [:matform,rowList]
+ rowList := []
+ matstring := bcwords2liststring [bcwords2liststring x for x in matform]
+ prefix := STRCONC('"f01qef(_"",wheret,"_",",STRINGIMAGE m,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE n,", ",STRINGIMAGE lda,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE ncolq,",[",zetastring,"],")
+ prefix := STRCONC(prefix,matstring,", ",STRINGIMAGE ifail,")")
+ linkGen prefix
+
+f01rcf() ==
+ htInitPage('"F01RCF - {\it QR} factorization of complex {\it m} by {\it n} matrix (m \htbitmap{great=} n)",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger)))
+ (text . "\windowlink{Manual Page}{manpageXXf01rcf} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01rcf| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "\newline ")
+ (text . "Finds the QR factorization of the complex m by n matrix {\it A},")
+ (text . " which is factorized as \htbitmap{f01qcf}, where m > n")
+ (text . " and A = QR when m = n , where Q is an m by m unitary matrix ")
+ (text . "and R is an n by n upper triangular matrix with real diagonal ")
+ (text . "elements. The {\it k}th transformation matrix,{\it Qk}, ")
+ (text . "which is used to introduce zeros into the {\it k}th column of ")
+ (text . "{\it A}, is given in the form ")
+ (text . "\htbitmap{f01qcf1}, ")
+ (text . "where \htbitmap{f01rdf2}, ")
+ (text . "\htbitmap{f01qcf3}, ")
+ (text . "\htbitmap{gammak} is a scalar for which Re ")
+ (text . "\htbitmap{gammak} = 1.0, \htbitmap{zetak} ")
+ (text . "is a real scalar and \htbitmap{zk} is an ")
+ (text . "(m-k) element vector. ")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline Rows of matrix A, {\it m}: ")
+ (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+ (text . "Columns of matrix A, {\it n}: \newline \tab{2} ")
+ (bcStrings (6 5 m PI))
+ (text . "\tab{34} ")
+ (bcStrings (6 3 n PI))
+-- (text . "\blankline ")
+-- (text . "\newline ")
+-- (text . "\menuitemstyle{} \tab{2} ")
+-- (text . "\newline First dimension of A, {\it lda} ")
+-- (text . "\htbitmap{great=} m: ")
+-- (text . "\newline \tab{2} ")
+-- (bcStrings (6 5 lda PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Ifail value: ")
+ (radioButtons ifail
+ ("" " -1, Print error messages" minusOne)
+ ("" " 1, Suppress error messages" one)))
+ htMakeDoneButton('"Continue", 'f01rcfSolve)
+ htShowPage()
+
+f01rcfSolve htPage ==
+ m :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
+ objValUnwrap htpLabelSpadValue(htPage, 'm)
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ lda := m
+-- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lda)
+-- objValUnwrap htpLabelSpadValue(htPage, 'lda)
+ error := htpButtonValue(htPage,'ifail)
+ ifail :=
+ error = 'one => '1
+ '-1
+ (m = '5 and n = '3) => f01rcfDefaultSolve(htPage,ifail)
+ matList :=
+ "append"/[fa(i,n) for i in 1..lda] where fa(i,n) ==
+ labelList :=
+ "append"/[ga(i,j) for j in 1..n] where ga(i,j) ==
+ anam := INTERN STRCONC ('"a",STRINGIMAGE i, STRINGIMAGE j)
+ [['bcStrings,[16, "0.0 + 0.0*%i", anam, 'F]]]
+ prefix := ('"\newline \tab{2} ")
+ labelList := [['text,:prefix],:labelList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :matList]
+ page := htInitPage('"F01RCF - {\it QR} factorization of complex {\it m} by {\it n} matrix (m \htbitmap{great=} n)",nil)
+ htSay '"\newline \menuitemstyle{}\tab{2} Enter values of {\it A}: "
+ htSay '"\newline \tab{2} "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01rcfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01rcfDefaultSolve (htPage,ifail) ==
+ n := '3
+ m := '5
+ lda := '5
+ page := htInitPage('"F01RCF - {\it QR} factorization of complex {\it m} by {\it n} matrix (m \htbitmap{great=} n)",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain I (Integer)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of {\it A}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.5*%i" a11 F))
+ (bcStrings (15 "-0.5 + 1.5*%i" a12 F))
+ (bcStrings (15 "-1.0 + 1.0*%i" a13 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.4 + 0.3*%i" a21 F))
+ (bcStrings (15 "0.9 + 1.3*%i" a22 F))
+ (bcStrings (15 "0.2 + 1.4*%i" a23 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.4" a31 F))
+ (bcStrings (15 "-0.4 + 0.4*%i" a32 F))
+ (bcStrings (15 "1.8" a33 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.3 - 0.4*%i" a41 F))
+ (bcStrings (15 "0.1 + 0.7*%i" a42 F))
+ (bcStrings (15 "0.0" a43 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "-0.3*%i" a51 F))
+ (bcStrings (15 "0.3 + 0.3*%i" a52 F))
+ (bcStrings (15 "2.4*%i" a53 F))
+ (text . "\blankline "))
+ htMakeDoneButton('"Continue",'f01rcfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01rcfGen htPage ==
+ n := htpProperty(htPage,'n)
+ m := htpProperty(htPage,'m)
+-- lda := htpProperty(htPage,'lda)
+ lda := m
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ y := REVERSE y
+ for i in 1..lda repeat
+ for j in 1..n repeat
+ elm := STRCONC((first y).1," ")
+ rowList := [:rowList,elm]
+ y := rest y
+ matform := [:matform,rowList]
+ rowList := []
+ matstring := bcwords2liststring [bcwords2liststring x for x in matform]
+ prefix := STRCONC('"f01rcf(",STRINGIMAGE m,", ",STRINGIMAGE n,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE lda,", ",matstring)
+ linkGen STRCONC(prefix,", ",STRINGIMAGE ifail,")")
+
+f01rdf() ==
+ htInitPage('"F01RDF - Operations with unitary matrices, compute {\it QB} or \htbitmap{f01rdf} after factorization by F01QCF",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger)))
+ (text . "\windowlink{Manual Page}{manpageXXf01rdf} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01rdf| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "\newline ")
+ (text . "Performs one of the transformations B = QB or B = ")
+ (text . "\htbitmap{f01rdf}, where B is an m ")
+ (text . "by ncolb matrix and Q is an m by m ")
+ (text . "unitary matrix assumed to be given by Q = ")
+ (text . "\htbitmap{f01rdf1}, \htbitmap{f01qdf2} ")
+ (text . "being given in the form \htbitmap{f01qcf1}, ")
+ (text . "where \htbitmap{f01rdf2}, \htbitmap{f01qcf3}")
+ (text . ", \htbitmap{gammak} is a scalar for which Re ")
+ (text . "\htbitmap{gammak} = 1.0, \htbitmap{zetak} ")
+ (text . "is a real scalar and \htbitmap{zk} is an ")
+ (text . "(m-k) element vector. ")
+ (text . "The routine is intended for use following F01QCF or F01QFF. ")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline Rows of matrix A, {\it m}: ")
+ (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+ (text . "Columns of matrix A, {\it n}: \newline \tab{2} ")
+ (bcStrings (6 5 m PI))
+ (text . "\tab{34} ")
+ (bcStrings (6 3 n PI))
+-- (text . "\blankline ")
+-- (text . "\newline ")
+-- (text . "\menuitemstyle{} \tab{2} ")
+-- (text . "\newline First dimension of A, {\it lda} ")
+-- (text . "\htbitmap{great=} m: ")
+-- (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+-- (text . "First dimension of B, {\it ldb} ")
+-- (text . "\htbitmap{great=} m: ")
+-- (text . "\newline \tab{2} ")
+-- (bcStrings (6 5 lda PI))
+-- (text . "\tab{34} ")
+-- (bcStrings (6 5 ldb PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Number of columns of matrix B {\it ncolb}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 2 ncolb PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Transformation to be performed: ")
+ (radioButtons trans
+ (" " " {\it B = QB}" no_trans)
+ (" " " {\it B =} \htbitmap{f01rdf} (Conjugate Transpose)" trans))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Where the elements can be found: ")
+ (radioButtons wheret
+ (" " " the elements of \theta are in A" in_a)
+ (" " " the elements of \theta are in THETA" seperate))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Ifail value: ")
+ (radioButtons ifail
+ ("" " -1, Print error messages" minusOne)
+ ("" " 1, Suppress error messages" one)))
+ htMakeDoneButton('"Continue", 'f01rdfSolve)
+ htShowPage()
+
+f01rdfSolve htPage ==
+ m :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
+ objValUnwrap htpLabelSpadValue(htPage, 'm)
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ lda := m
+-- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lda)
+-- objValUnwrap htpLabelSpadValue(htPage, 'lda)
+ ldb := m
+-- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ldb)
+-- objValUnwrap htpLabelSpadValue(htPage, 'ldb)
+ ncolb :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ncolb)
+ objValUnwrap htpLabelSpadValue(htPage, 'ncolb)
+ operation := htpButtonValue(htPage,'trans)
+ trans :=
+ operation = 'no_trans => '"n"
+ '"c"
+ elements := htpButtonValue(htPage,'wheret)
+ wheret :=
+ elements = 'in_a => '"i"
+ '"c"
+ error := htpButtonValue(htPage,'ifail)
+ ifail :=
+ error = 'one => '1
+ '-1
+ ((m = '5 and n = '3) and ncolb = '2) => f01rdfDefaultSolve(htPage,lda,ldb,trans,wheret,ifail)
+ matList :=
+ "append"/[fa(i,n) for i in 1..lda] where fa(i,n) ==
+ labelList :=
+ "append"/[ga(i,j) for j in 1..n] where ga(i,j) ==
+ anam := INTERN STRCONC ('"a",STRINGIMAGE i, STRINGIMAGE j)
+ [['bcStrings,[16, "0.0 + 0.0*%i", anam, 'F]]]
+ prefix := ('"\newline \tab{2} ")
+ labelList := [['text,:prefix],:labelList]
+ bList :=
+ "append"/[fb(i,ncolb) for i in 1..ldb] where fb(i,ncolb) ==
+ labelList :=
+ "append"/[gb(i,j) for j in 1..ncolb] where gb(i,j) ==
+ bnam := INTERN STRCONC ('"b",STRINGIMAGE i, STRINGIMAGE j)
+ [['bcStrings,[16, "0.0 + 0.0*%i", bnam, 'F]]]
+ prefix := ('"\newline \tab{2} ")
+ labelList := [['text,:prefix],:labelList]
+ prefix := ("\blankline \menuitemstyle{}\tab{2} Enter values of {\it B}: ")
+ bList := [['text,:prefix],:bList]
+ zList :=
+ "append"/[fz(i) for i in 1..n] where fz(i) ==
+ znam := INTERN STRCONC ('"z",STRINGIMAGE i)
+ [['bcStrings,[16, "0.0", znam, 'F]]]
+ prefix := ("\blankline \menuitemstyle{}\tab{2} Enter values of \theta ")
+ prefix := STRCONC(prefix,"(if required): \newline \tab{2}")
+ zList := [['text,:prefix],:zList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :matList,:bList,:zList]
+ page := htInitPage('"F01RDF - Operations with orthogonal matrices, compute {\it QB} or \htbitmap{f01rdf} after factorization by F01QCF or F01RDF",nil)
+ htSay '"\newline \menuitemstyle{}\tab{2} Enter values of {\it A}: "
+ htSay '"\newline \tab{2} "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01rdfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+-- htpSetProperty(page,'ldb,ldb)
+ htpSetProperty(page,'ncolb,ncolb)
+ htpSetProperty(page,'trans,trans)
+ htpSetProperty(page,'wheret,wheret)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01rdfDefaultSolve (htPage,lda,ldb,trans,wheret,ifail) ==
+ n := '3
+ m := '5
+ ncolb := '2
+ page := htInitPage('"F01RDF - Operations with orthogonal matrices, compute {\it QB} or \htbitmap{f01rdf} after factorization by F01QCF or F01RDF",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain I (Integer)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of {\it A}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.5*%i" a11 F))
+ (bcStrings (15 "-0.5 + 1.5*%i" a12 F))
+ (bcStrings (15 "-1.0 + 1.0*%i" a13 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.4 + 0.3*%i" a21 F))
+ (bcStrings (15 "0.9 + 1.3*%i" a22 F))
+ (bcStrings (15 "0.2 + 1.4*%i" a23 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.4" a31 F))
+ (bcStrings (15 "-0.4 + 0.4*%i" a32 F))
+ (bcStrings (15 "1.8" a33 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.3 - 0.4*%i" a41 F))
+ (bcStrings (15 "0.1 + 0.7*%i" a42 F))
+ (bcStrings (15 "0.0" a43 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "-0.3*%i" a51 F))
+ (bcStrings (15 "0.3 + 0.3*%i" a52 F))
+ (bcStrings (15 "2.4" a53 F))
+ (text . "\blankline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of {\it B}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "-0.55 + 1.05*%i" b11 F))
+ (bcStrings (15 "0.45 + 1.05*%i" b12 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.49 + 0.93*%i" b21 F))
+ (bcStrings (15 "1.09 + 0.13*%i" b22 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.56 - 0.16*%i" b31 F))
+ (bcStrings (15 "0.64 + 0.16*%i" b32 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.39 + 0.23*%i" b41 F))
+ (bcStrings (15 "-0.39 - 0.23*%i" b42 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "1.13 + 0.83*%i" b51 F))
+ (bcStrings (15 "-1.13 + 0.77*%i" b52 F))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of \theta (if required): ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (15 "0.0" z1 F))
+ (bcStrings (15 "0.0" z2 F))
+ (bcStrings (15 "0.0" z3 F))
+ (text . "\blankline "))
+ htMakeDoneButton('"Continue",'f01rdfGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+-- htpSetProperty(page,'ldb,ldb)
+ htpSetProperty(page,'ncolb,ncolb)
+ htpSetProperty(page,'trans,trans)
+ htpSetProperty(page,'wheret,wheret)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01rdfGen htPage ==
+ n := htpProperty(htPage,'n)
+ m := htpProperty(htPage,'m)
+-- lda := htpProperty(htPage,'lda)
+-- ldb := htpProperty(htPage,'ldb)
+ lda := m
+ ldb := m
+ ncolb := htpProperty(htPage,'ncolb)
+ trans := htpProperty(htPage,'trans)
+ wheret := htpProperty(htPage,'wheret)
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ for i in 1..n repeat
+ left := STRCONC((first y).1," ")
+ y := rest y
+ zetalist := [left,:zetalist]
+ zetastring := bcwords2liststring zetalist
+ y := REVERSE y
+ for i in 1..lda repeat
+ for j in 1..n repeat
+ elm := STRCONC((first y).1," ")
+ rowList := [:rowList,elm]
+ y := rest y
+ matform := [:matform,rowList]
+ rowList := []
+ matstring := bcwords2liststring [bcwords2liststring x for x in matform]
+ for i in 1..ldb repeat
+ for j in 1..ncolb repeat
+ elm := STRCONC((first y).1," ")
+ rowList := [:rowList,elm]
+ y := rest y
+ bform := [:bform,rowList]
+ rowList := []
+ bstring := bcwords2liststring [bcwords2liststring x for x in bform]
+ prefix := STRCONC('"f01rdf(_"",trans,"_",_"",wheret,"_",",STRINGIMAGE m,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE n,", ",matstring,", ",STRINGIMAGE lda)
+ prefix := STRCONC(prefix,",[",zetastring,"],",STRINGIMAGE ncolb,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE ldb,", ",bstring,", ",STRINGIMAGE ifail,")")
+ linkGen prefix
+
+
+f01ref() ==
+ htInitPage('"F01REF - Operations with unitary matrices, form columns of {\it Q} after factorization by F01RCF",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger)))
+ (text . "\windowlink{Manual Page}{manpageXXf01ref} for this routine ")
+ (text . "\newline ")
+ (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|f01ref| '|NagMatrixOperationsPackage|)} for this routine")
+ (text . "\newline \horizontalline ")
+ (text . "\newline ")
+ (text . "Returns the first {\it ncolq} columns of the real {\it m} by ")
+ (text . "{\it m} unitary matrix {\it Q}, where {\it Q} is assumed ")
+ (text . "to be given by {\it Q = }\htbitmap{f01rdf1}, ")
+ (text . "\htbitmap{f01qdf2} being given in the form ")
+ (text . "\htbitmap{f01qcf1}, ")
+ (text . "where \htbitmap{f01rdf2}, ")
+ (text . "\htbitmap{f01qcf3}, ")
+ (text . "\htbitmap{gammak} is a scalar for which Re ")
+ (text . "\htbitmap{gammak} = 1.0, \htbitmap{zetak} ")
+ (text . "is a real scalar and \htbitmap{zk} is an ")
+ (text . "(m-k) element vector. ")
+ (text . "The routine is intended for use following F01RCF or F01RFF. ")
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline Rows of matrix A, {\it m}: ")
+ (text . "\tab{32} \menuitemstyle{} \tab{34} ")
+ (text . "Columns of matrix A, {\it n}: \newline \tab{2} ")
+ (bcStrings (6 5 m PI))
+ (text . "\tab{34} ")
+ (bcStrings (6 3 n PI))
+-- (text . "\blankline ")
+-- (text . "\newline ")
+-- (text . "\menuitemstyle{} \tab{2} ")
+-- (text . "\newline First dimension of A, {\it lda} ")
+-- (text . "\htbitmap{great=} m: ")
+-- (text . "\newline \tab{2} ")
+-- (bcStrings (6 5 lda PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Required number of columns of matrix Q {\it ncolq}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (6 2 ncolq PI))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Where the elements can be found: ")
+ (radioButtons wheret
+ (" " " the elements of \theta are in THETA" seperate)
+ (" " " the elements of \theta are in A" in_a))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{} \tab{2} ")
+ (text . "\newline \tab{2} ")
+ (text . "Ifail value: ")
+ (radioButtons ifail
+ ("" " -1, Print error messages" minusOne)
+ ("" " 1, Suppress error messages" one)))
+ htMakeDoneButton('"Continue", 'f01refSolve)
+ htShowPage()
+
+f01refSolve htPage ==
+ m :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
+ objValUnwrap htpLabelSpadValue(htPage, 'm)
+ n :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
+ objValUnwrap htpLabelSpadValue(htPage, 'n)
+ lda := m
+-- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lda)
+-- objValUnwrap htpLabelSpadValue(htPage, 'lda)
+ ncolq :=
+ $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ncolq)
+ objValUnwrap htpLabelSpadValue(htPage, 'ncolq)
+ elements := htpButtonValue(htPage,'wheret)
+ wheret :=
+ elements = 'in_a => '"i"
+ '"s"
+ error := htpButtonValue(htPage,'ifail)
+ ifail :=
+ error = 'one => '1
+ '-1
+ ((m = '5 and n = '3) and ncolq = '2) => f01refDefaultSolve(htPage,lda,wheret,ifail)
+ matList :=
+ "append"/[fa(i,n) for i in 1..lda] where fa(i,n) ==
+ labelList :=
+ "append"/[ga(i,j) for j in 1..n] where ga(i,j) ==
+ anam := INTERN STRCONC ('"a",STRINGIMAGE i, STRINGIMAGE j)
+ [['bcStrings,[20, "0.0 + 0.0*%i", anam, 'F]]]
+ prefix := ('"\newline \tab{2} ")
+ labelList := [['text,:prefix],:labelList]
+ zList :=
+ "append"/[fz(i) for i in 1..n] where fz(i) ==
+ znam := INTERN STRCONC ('"z",STRINGIMAGE i)
+ [['bcStrings,[20, "0.0", znam, 'F]]]
+ prefix := ("\blankline \menuitemstyle{}\tab{2} Enter values of \theta ")
+ prefix := STRCONC(prefix,"(if required): \newline \tab{2}")
+ zList := [['text,:prefix],:zList]
+ equationPart := [
+ '(domainConditions
+ (isDomain P (Polynomial $EmptyMode))
+ (isDomain S (String))
+ (isDomain F (Float))
+ (isDomain PI (PositiveInteger))),
+ :matList,:zList]
+ page := htInitPage('"F01REF - Operations with unitary matrices, form columns of {\it Q} after factorization by F01RCF",nil)
+ htSay '"\newline \menuitemstyle{}\tab{2} Enter values of {\it A}: "
+ htSay '"\newline \tab{2} "
+ htMakePage equationPart
+ htSay '"\blankline "
+ htMakeDoneButton('"Continue",'f01refGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+ htpSetProperty(page,'ncolq,ncolq)
+ htpSetProperty(page,'wheret,wheret)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+
+f01refDefaultSolve (htPage,lda,wheret,ifail) ==
+ n := '3
+ m := '5
+ ncolq := '2
+ page := htInitPage('"F01REF - Operations with unitary matrices, form columns of {\it Q} after factorization by F01RCF",nil)
+ htMakePage '(
+ (domainConditions
+ (isDomain F (Float))
+ (isDomain I (Integer)))
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of {\it A}: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (16 "1" a11 F))
+ (bcStrings (16 "1 + %i" a12 F))
+ (bcStrings (16 "1 + %i" a13 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (16 "-0.2-0.4*%i" a21 F))
+ (bcStrings (16 "-2" a22 F))
+ (bcStrings (16 "-1 - %i" a23 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (16 "-0.32 - 0.16*%i" a31 F))
+ (bcStrings (16 "-0.3505+0.263*%i" a32 F))
+ (bcStrings (16 "-3" a33 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (16 "-0.4 + 0.2*%i" a41 F))
+ (bcStrings (16 "0.5477*%i" a42 F))
+ (bcStrings (16 "0.0" a43 F))
+ (text . "\newline \tab{2} ")
+ (bcStrings (16 "-0.12 + 0.24*%i" a51 F))
+ (bcStrings (16 "0.1972+0.2629*%i" a52 F))
+ (bcStrings (16 "0.6325" a53 F))
+ (text . "\blankline ")
+ (text . "\newline ")
+ (text . "\menuitemstyle{}\tab{2} Enter values of \theta: ")
+ (text . "\newline \tab{2} ")
+ (bcStrings (16 "1 + 0.5*%i" z1 F))
+ (bcStrings (16 "1.0954-0.3333*%i" z2 F))
+ (bcStrings (16 "1.2649-1.1565*%i" z3 F))
+ (text . "\blankline "))
+ htMakeDoneButton('"Continue",'f01refGen)
+ htpSetProperty(page,'n,n)
+ htpSetProperty(page,'m,m)
+-- htpSetProperty(page,'lda,lda)
+ htpSetProperty(page,'ncolq,ncolq)
+ htpSetProperty(page,'wheret,wheret)
+ htpSetProperty(page,'ifail,ifail)
+ htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
+ htShowPage()
+
+f01refGen htPage ==
+ n := htpProperty(htPage,'n)
+ m := htpProperty(htPage,'m)
+-- lda := htpProperty(htPage,'lda)
+ lda := m
+ ncolq := htpProperty(htPage,'ncolq)
+ wheret := htpProperty(htPage,'wheret)
+ ifail := htpProperty(htPage,'ifail)
+ alist := htpInputAreaAlist htPage
+ y := alist
+ for i in 1..n repeat
+ left := STRCONC((first y).1," ")
+ y := rest y
+ thetalist := [left,:thetalist]
+ thetastring := bcwords2liststring thetalist
+ y := REVERSE y
+ for i in 1..lda repeat
+ for j in 1..n repeat
+ elm := STRCONC((first y).1," ")
+ rowList := [:rowList,elm]
+ y := rest y
+ matform := [:matform,rowList]
+ rowList := []
+ matstring := bcwords2liststring [bcwords2liststring x for x in matform]
+ prefix := STRCONC('"f01ref(_"",wheret,"_",",STRINGIMAGE m,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE n,", ",STRINGIMAGE ncolq,", ")
+ prefix := STRCONC(prefix,STRINGIMAGE lda,",[",thetastring,"],")
+ prefix := STRCONC(prefix,matstring,", ",STRINGIMAGE ifail,")")
+ linkGen prefix
+
+@
+\eject
+\begin{thebibliography}{99}
+\bibitem{1} nothing
+\end{thebibliography}
+\end{document}