aboutsummaryrefslogtreecommitdiff
path: root/src/interp/nag-c02.boot
diff options
context:
space:
mode:
Diffstat (limited to 'src/interp/nag-c02.boot')
-rw-r--r--src/interp/nag-c02.boot299
1 files changed, 0 insertions, 299 deletions
diff --git a/src/interp/nag-c02.boot b/src/interp/nag-c02.boot
deleted file mode 100644
index a8661a63..00000000
--- a/src/interp/nag-c02.boot
+++ /dev/null
@@ -1,299 +0,0 @@
--- Copyright (c) 1991-2002, The Numerical ALgorithms Group Ltd.
--- All rights reserved.
--- Copyright (C) 2007, Gabriel Dos Reis.
--- 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.
-
-
-import '"macros"
-)package "BOOT"
-
-c02aff() ==
- htInitPage('"C02AFF - All Zeros of a Complex Polynomial",nil)
- htMakePage '(
- (domainConditions
- (isDomain PI (PositiveInteger)))
- (text . "\windowlink{Manual Page}{manpageXXc02aff} for this routine ")
- (text . "\newline ")
- (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|c02aff| '|NagPolynomialRootsPackage|)} for this routine")
- (text . "\newline \horizontalline ")
- (text . "Finds all the roots of the complex polynomial equation ")
- (text . "\htbitmap{c02aff}, using a variant of Laguerre's method. ")
- (text . "\blankline ")
- (text . "\newline \menuitemstyle{} ")
- (text . "\tab{2} Enter the degree {\em n} of the polynomial:")
- (text . "\newline\tab{2} ")
- (bcStrings (5 5 n PI))
- (text . "\blankline")
- (text . "\newline")
- (text . "\newline \menuitemstyle{} \tab{2} Scale value:")
- (radioButtons scale
- ("" " True" true)
- ("" " False" false))
- (text . "\blankline ")
- (text . "\newline \menuitemstyle{} \tab{2} Ifail value:")
- (radioButtons ifail
- ("" " -1, Print error messages" minusOne)
- ("" " 1, Suppress error messages" one)))
- htMakeDoneButton('"Continue", 'c02affSolve)
- htShowPage()
-
-c02affSolve htPage ==
- n :=
- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
- objValUnwrap htpLabelSpadValue(htPage, 'n)
- logical := htpButtonValue(htPage,'scale)
- scale :=
- logical = 'true => '"true"
- '"false"
- error := htpButtonValue(htPage,'ifail)
- ifail :=
- error = 'one => '1
- '-1
- n = '5 => c02affDefaultSolve(htPage,scale,ifail)
- labelList :=
- "append"/[f(i,n) for i in 1..(n+1)] where f(i,n) ==
- prefix := ('"\newline z**")
- prefix := STRCONC(prefix,STRINGIMAGE (n-i+1),'"\space{1}")
- post := ('"\tab{30} ")
- post := STRCONC(post,'"\space{1}")
- rnam := INTERN STRCONC ('"r",STRINGIMAGE i)
- inam := INTERN STRCONC ('"i",STRINGIMAGE i)
- [['text,:prefix],['bcStrings,[10, 0.0, rnam, 'P]],
- ['text,:post],['bcStrings,[10, 0.0, inam, 'P]]]
- equationPart := [
- '(domainConditions
- (isDomain P (Polynomial $EmptyMode))
- (isDomain S (String))
- (isDomain PI (PositiveInteger))),
- :labelList]
- page := htInitPage("C02AFF - All Zeros of a Complex Polynomial", htpPropertyList htPage)
- htSay '"\menuitemstyle{} \tab{2} Enter the coefficients of the polynomial: "
- htSay '"\blankline "
- htSay '"Real parts \tab{30} Imaginary parts "
- htSay '"\newline "
- htMakePage equationPart
- htMakeDoneButton('"Continue",'c02affGen)
- htpSetProperty(page,'n,n)
- htpSetProperty(page,'scale,scale)
- htpSetProperty(page,'ifail,ifail)
- htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
- htShowPage()
-
-c02affDefaultSolve (htPage, scale, ifail) ==
- n := '5
- page := htInitPage('"C02AFF - All Zeros of a Complex Polynomial",nil)
- htMakePage '(
- (domainConditions
- (isDomain F (Float)))
- (text . "\newline \menuitemstyle{} \tab{2} ")
- (text . "Enter the coefficients of the polynomial: ")
- (text . "\blankline ")
- (text . "Real parts \tab{30} Imaginary parts ")
- (text . "\newline z**5 \space{1} ")
- (bcStrings (10 "5.0" r1 F))
- (text . "\tab{30} ")
- (text . "\space{1} ")
- (bcStrings (10 "6.0" i1 F))
- (text . "\newline ")
- (text . "z**4 \space{1} ")
- (bcStrings (10 "30.0" r2 F))
- (text . "\tab{30} ")
- (text . "\space{1} ")
- (bcStrings (10 "20.0" i2 F))
- (text . "\newline ")
- (text . "z**3 \space{1} ")
- (bcStrings (10 "-0.2" r3 F))
- (text . "\tab{30} ")
- (text . "\space{1} ")
- (bcStrings (10 "-6.0" i3 F))
- (text . "\newline ")
- (text . "z**2 \space{1} ")
- (bcStrings (10 "50.0" r4 F))
- (text . "\tab{30} ")
- (text . "\space{1} ")
- (bcStrings (10 "100000.0" i4 F))
- (text . "\newline ")
- (text . "z**1 \space{1} ")
- (bcStrings (10 "-2.0" r5 F))
- (text . "\tab{30} ")
- (text . "\space{1} ")
- (bcStrings (10 "40.0" i5 F))
- (text . "\newline ")
- (text . "z**0 \space{1} ")
- (bcStrings (10 "10.0" r6 F))
- (text . "\tab{30} ")
- (text . "\space{1} ")
- (bcStrings (10 "1.0" i6 F))
- (text . "\newline ")
- (text . "\blankline"))
- htMakeDoneButton('"Continue",'c02affGen)
- htpSetProperty(page,'n,n)
- htpSetProperty(page,'scale,scale)
- htpSetProperty(page,'ifail,ifail)
- htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
- htShowPage()
-
-
-c02affGen htPage ==
- n := htpProperty(htPage,'n)
- scale := htpProperty(htPage,'scale)
- ifail := htpProperty(htPage,'ifail)
- alist := htpInputAreaAlist htPage
- y := alist
- while y repeat
- right := STRCONC((first y).1," ")
- y := rest y
- left := STRCONC((first y).1," ")
- y := rest y
- reallist := [left,:reallist]
- imaglist := [right,:imaglist]
- realstring := bcwords2liststring reallist
- imagstring := bcwords2liststring imaglist
- linkGen STRCONC ('"c02aff([",realstring,",",imagstring,"],",STRINGIMAGE n,",",scale,",",STRINGIMAGE ifail,")")
-
-c02agf() ==
- htInitPage('"C02AGF - All Zeros of a Real Polynomial",nil)
- htMakePage '(
- (domainConditions
- (isDomain PI (PositiveInteger)))
- (text . "\windowlink{Manual Page}{manpageXXc02agf} for this routine ")
- (text . "\newline ")
- (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|c02agf| '|NagPolynomialRootsPackage|)} for this routine")
- (text . "\newline \horizontalline ")
- (text . "Finds all the roots of the real polynomial equation ")
- (text . "\htbitmap{c02aff}, using a variant of Laguerre's method. ")
- (text . "\blankline ")
- (text . "\newline \menuitemstyle{} ")
- (text . "\tab{2} Enter the degree {\em n} of the polynomial:")
- (text . "\newline\tab{2} ")
- (bcStrings (5 5 n PI))
- (text . "\blankline")
- (text . "\newline")
- (text . "\newline \menuitemstyle{} \tab{2} Scale value:")
- (radioButtons scale
- ("" " True" true)
- ("" " False" false))
- (text . "\blankline ")
- (text . "\newline \menuitemstyle{} \tab{2} Ifail value:")
- (radioButtons ifail
- ("" " -1, Print error messages" minusOne)
- ("" " 1, Suppress error messages" one)))
- htMakeDoneButton('"Continue", 'c02agfSolve)
- htShowPage()
-
-c02agfSolve htPage ==
- n :=
- $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
- objValUnwrap htpLabelSpadValue(htPage, 'n)
- logical := htpButtonValue(htPage,'scale)
- scale :=
- logical = 'true => '"true"
- '"false"
- error := htpButtonValue(htPage,'ifail)
- ifail :=
- error = 'one => '1
- '-1
- n = '5 => c02agfDefaultSolve(htPage,scale,ifail)
- labelList :=
- "append"/[f(i,n) for i in 1..(n+1)] where f(i,n) ==
- prefix := ('"\newline z**")
- prefix := STRCONC(prefix,STRINGIMAGE (n-i+1),'"\space{1}")
- rnam := INTERN STRCONC ('"r",STRINGIMAGE i)
- [['text,:prefix],['bcStrings,[10, 0.0, rnam, 'P]]]
- equationPart := [
- '(domainConditions
- (isDomain P (Polynomial $EmptyMode))
- (isDomain S (String))
- (isDomain PI (PositiveInteger))),
- :labelList]
- page := htInitPage("C02AGF - All Zeros of a Real Polynomial", htpPropertyList htPage)
- htSay '"\menuitemstyle{} \tab{2} "
- htSay '"Enter the coefficients of the polynomial: "
- htSay '"\newline "
- htMakePage equationPart
- htMakeDoneButton('"Continue",'c02agfGen)
- htpSetProperty(page,'n,n)
- htpSetProperty(page,'scale,scale)
- htpSetProperty(page,'ifail,ifail)
- htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
- htShowPage()
-
-
-c02agfDefaultSolve (htPage, scale, ifail) ==
- n := '5
- page := htInitPage('"C02AGF - All Zeros of a Real Polynomial",nil)
- htMakePage '(
- (domainConditions
- (isDomain F (Float)))
- (text . "\newline \menuitemstyle{} \tab{2} ")
- (text . "Enter the coefficients of the polynomial: ")
- (text . "\newline ")
- (text . "z**5 \space{1} ")
- (bcStrings (10 "1.0" r1 F))
- (text . "\newline ")
- (text . "z**4 \space{1} ")
- (bcStrings (10 "2.0" r2 F))
- (text . "\newline ")
- (text . "z**3 \space{1} ")
- (bcStrings (10 "3.0" r3 F))
- (text . "\newline ")
- (text . "z**2 \space{1} ")
- (bcStrings (10 "4.0" r4 F))
- (text . "\newline ")
- (text . "z**1 \space{1} ")
- (bcStrings (10 "5.0" r5 F))
- (text . "\newline ")
- (text . "z**0 \space{1} ")
- (bcStrings (10 "6.0" r6 F))
- (text . "\newline ")
- (text . "\blankline"))
- htMakeDoneButton('"Continue",'c02agfGen)
- htpSetProperty(page,'n,n)
- htpSetProperty(page,'scale,scale)
- htpSetProperty(page,'ifail,ifail)
- htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
- htShowPage()
-
-
-c02agfGen htPage ==
- n := htpProperty(htPage,'n)
- scale := htpProperty(htPage,'scale)
- ifail := htpProperty(htPage,'ifail)
- alist := htpInputAreaAlist htPage
- y := alist
- while y repeat
- left := STRCONC((first y).1," ")
- y := rest y
- reallist := [left,:reallist]
- realstring := bcwords2liststring reallist
- linkGen STRCONC ('"c02agf([",realstring,"],",STRINGIMAGE n,",",scale,",",STRINGIMAGE ifail,")")
-
-