• Solving system of non-linear equations

    From David Dooner@david.dooner@upr.edu to comp.lang.apl on Tue Mar 30 05:21:35 2021
    From Newsgroup: comp.lang.apl

    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
    ZraE(FN ZERO_NEWTON)X;rAoIO;X0;FX;TOL
    riY
    riY Newton's method to find roots to system of eqs. (see ZERO_SECANT for scalar version)
    riY
    riY Form of function FN is the following (sample: FCT_CUTTER_SETUP_I)
    riY FN[1;] -- 1st equation
    riY FN[2;] -- 2nd equation
    riY FN[3;] -- 3rd equation
    riY FN[n;] -- nth equation
    riY
    rAoIOraE1
    TOLraE0.00001
    ZraE(0,2|uri|,X)ri|0
    L1:X0raEX
    FXraEFN X0
    ZraEZ,[1]FX,X0
    XraEX0-FXri|FN JACOBI X
    raA(TOL<+/(X-X0)*2)/L1
    ZraEX
    ZraE(FN JACOBI)X;rAoIO;rea;N;reaX
    riY
    riY Jacobian of function FN evaluated at 'X' (see ZERO_NEWTON, OPTI_NEWTON, HESSIAN)
    riY
    rAoIOraE1
    rearaE0.0001
    NraEri|ri|XraEreeX
    reaXraEree[1]rea|uNrey.=N
    ZraE-rierea((reeFN X)-FN-?reaX+reeX)|+rea
    --- Synchronet 3.21d-Linux NewsLink 1.2
  • From eki...@gmail.com@ekimlpa@gmail.com to comp.lang.apl on Tue Mar 30 11:27:43 2021
    From Newsgroup: comp.lang.apl

    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,
    David

    ZraE(FN ZERO_NEWTON)X;rAoIO;X0;FX;TOL
    riY
    riY Newton's method to find roots to system of eqs. (see ZERO_SECANT for scalar version)
    riY
    riY Form of function FN is the following (sample: FCT_CUTTER_SETUP_I)
    riY FN[1;] -- 1st equation
    riY FN[2;] -- 2nd equation
    riY FN[3;] -- 3rd equation
    riY FN[n;] -- nth equation
    riY

    rAoIOraE1
    TOLraE0.00001
    ZraE(0,2|uri|,X)ri|0
    L1:X0raEX
    FXraEFN X0
    ZraEZ,[1]FX,X0
    XraEX0-FXri|FN JACOBI X
    raA(TOL<+/(X-X0)*2)/L1
    ZraEX


    ZraE(FN JACOBI)X;rAoIO;rea;N;reaX
    riY
    riY Jacobian of function FN evaluated at 'X' (see ZERO_NEWTON, OPTI_NEWTON, HESSIAN)
    riY
    rAoIOraE1
    rearaE0.0001
    NraEri|ri|XraEreeX
    reaXraEree[1]rea|uNrey.=N
    ZraE-rierea((reeFN X)-FN-?reaX+reeX)|+rea
    --- Synchronet 3.21d-Linux NewsLink 1.2
  • From eki...@gmail.com@ekimlpa@gmail.com to comp.lang.apl on Tue Mar 30 12:00:21 2021
    From Newsgroup: comp.lang.apl

    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,
    David
    [[ snipped ]]
    The dl.acm.org link in your message points to the Table of Contents for the APL93 Conference Proceedings.
    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 student members.
    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 end up needing as few as 2 articles.
    --- Synchronet 3.21d-Linux NewsLink 1.2
  • From David Dooner@david.dooner@upr.edu to comp.lang.apl on Tue Mar 30 12:56:36 2021
    From Newsgroup: comp.lang.apl

    On Tuesday, March 30, 2021 at 3:00:22 PM UTC-4, eki...@gmail.com wrote:
    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,
    David
    [[ snipped ]]


    The dl.acm.org link in your message points to the Table of Contents for the APL93 Conference Proceedings.

    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 student members.

    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 end up needing as few as 2 articles.
    Thank you, I was able to purchase a copy.
    --- Synchronet 3.21d-Linux NewsLink 1.2