All,
I have a system of non-linear equations. Below is my current operator. Hmm, I have doubts. I look online for related code and find below by Stephan Mansour:
https://www.deepdyve.com/lp/association-for-computing-machinery/using-defined-operators-and-function-arrays-to-solve-non-linear-0CRfFEIce6
It's dated Sept 1993. I look in SIGAPL for the article and Sept. is missing.
https://dl.acm.org/toc/sigapl/1993/24/1
Please share if you have APL code for solving non-linear equations. FYI, I'm using APL2 (ex to be).
Thanks,
David
All,
I have a system of non-linear equations. Below is my current operator. Hmm, I have doubts. I look online for related code and find below by Stephan Mansour:
https://www.deepdyve.com/lp/association-for-computing-machinery/using-defined-operators-and-function-arrays-to-solve-non-linear-0CRfFEIce6
It's dated Sept 1993. I look in SIGAPL for the article and Sept. is missing.
https://dl.acm.org/toc/sigapl/1993/24/1
Please share if you have APL code for solving non-linear equations. FYI, I'm using APL2 (ex to be).
Thanks,
David
Z←(FN ZERO_NEWTON)X;⎕IO;X0;FX;TOL
⍝
⍝ Newton's method to find roots to system of eqs. (see ZERO_SECANT for scalar version)
⍝
⍝ Form of function FN is the following (sample: FCT_CUTTER_SETUP_I)
⍝ FN[1;] -- 1st equation
⍝ FN[2;] -- 2nd equation
⍝ FN[3;] -- 3rd equation
⍝ FN[n;] -- nth equation
⍝
⎕IO←1
TOL←0.00001
Z←(0,2×⍴,X)⍴0
L1:X0←X
FX←FN X0
Z←Z,[1]FX,X0
X←X0-FX⌹FN JACOBI X
→(TOL<+/(X-X0)*2)/L1
Z←X
Z←(FN JACOBI)X;⎕IO;∆;N;∆X
⍝
⍝ Jacobian of function FN evaluated at 'X' (see ZERO_NEWTON, OPTI_NEWTON, HESSIAN)
⍝
⎕IO←1
∆←0.0001
N←⍳⍴X←∊X
∆X←⊂[1]∆×N∘.=N
Z←-⍉⊃((⊂FN X)-FN¨∆X+⊂X)÷∆
On Tuesday, March 30, 2021 at 8:21:36 AM UTC-4, David Dooner wrote:
All,
I have a system of non-linear equations. Below is my current operator. Hmm, I have doubts. I look online for related code and find below by Stephan Mansour:
https://www.deepdyve.com/lp/association-for-computing-machinery/using-defined-operators-and-function-arrays-to-solve-non-linear-0CRfFEIce6
It's dated Sept 1993. I look in SIGAPL for the article and Sept. is missing.
https://dl.acm.org/toc/sigapl/1993/24/1
Please share if you have APL code for solving non-linear equations. FYI, I'm using APL2 (ex to be).
Thanks,[[ snipped ]]
David
The dl.acm.org link in your message points to the Table of Contents for the APL93 Conference Proceedings.student members.
If you scroll down through the list of articles, the one you want is the 18th. Clicking on the article title links to the abstract and some metadata. The full text is behind ACM's paywall and costs $15 for anybody, $10 for ACM members, or $5 for ACM
ACM SIGPLAN members have free access to all SIGPLAN proceedings, which now includes all SIGAPL conference proceedings. SIGPLAN membership ( https://www.sigplan.org/Membership ) for a year is $25 ($15 for students), and is thus less expensive if you endup needing as few as 2 articles.
Sysop: | Keyop |
---|---|
Location: | Huddersfield, West Yorkshire, UK |
Users: | 113 |
Nodes: | 8 (1 / 7) |
Uptime: | 123:53:16 |
Calls: | 2,501 |
Files: | 8,692 |
Messages: | 1,924,232 |