#282 EZGCD: ratsimp((%i+2)/(%i+1)),gcd:ez inf loop

open
nobody
5
2012-12-15
2003-03-29
No

ratsimp((%i+2)/(%i+1)),gcd:ez; gets into an infinite loop
inside ezgcd2.

Maxima 5.9.0 GCL 2.5.0 Windows 2000

Discussion

  • Stavros Macrakis

    Logged In: YES
    user_id=588346

    ratsimp((%i+2)/(%i+1)),gcd:mod;

    also gets into an infinite loop.

     
  • Stavros Macrakis

    Logged In: YES
    user_id=588346

    but ratsimp((%i+2)/(%i+1)),gcd:mod,algebraic does not

     
  • Robert Dodier

    Robert Dodier - 2006-07-06

    Logged In: YES
    user_id=501686

    Still present in 5.9.3cvs.

     
  • Robert Dodier

    Robert Dodier - 2006-07-06
    • labels: --> 460522
     
  • Robert Dodier

    Robert Dodier - 2006-09-09
    • labels: 460522 --> Lisp Core - Polynomials
     
  • Robert Dodier

    Robert Dodier - 2012-12-15

    Still present in 5.29.1. Gets stuck in PGCDM ("Brown's algorithm M"). Test at Step 12 always fails, so it loops back to Step 6. PGDM calls LAGRANGE3 with increasingly enormous integer arguments.

    (%i8) ratsimp((%i+2)/(%i+1)),gcd:ez;
    1. Trace: (PGCDA '(#:%I16214 1 1 0 2) '(#:%I16214 1 1 0 1) 'T)
    2. Trace: (EZGCD2 '(#:%I16214 1 1 0 2) '(#:%I16214 1 1 0 1))
    3. Trace: (NEWGCD '(#:%I16214 1 1 0 2) '(#:%I16214 1 1 0 1) 'NIL)
    4. Trace: (PGCDM '(#:%I16214 1 1 0 2) '(#:%I16214 1 1 0 1))
    5. Trace: (NEWGCD '(#:%I16214 1 1 0 2) '(#:%I16214 1 1 0 1) '8388593)
    6. Trace: (PGCDA '(#:%I16214 1 -1 0 1) '2 'T)
    6. Trace: PGCDA ==> (1)
    6. Trace: (PGCDA '1 '2 'T)
    6. Trace: PGCDA ==> (1 1 2)
    5. Trace: NEWGCD ==> ((#:%I16214 1 1 0 1) (#:%I16214 1 4194296 0 -4194295) 1)
    5. Trace: (NEWGCD '(#:%I16214 1 1 0 2) '(#:%I16214 1 1 0 1) '8388587)
    6. Trace: (PGCDA '(#:%I16214 1 -1 0 1) '2 'T)
    6. Trace: PGCDA ==> (1)
    6. Trace: (PGCDA '1 '2 'T)
    6. Trace: PGCDA ==> (1 1 2)
    5. Trace: NEWGCD ==> ((#:%I16214 1 1 0 1) (#:%I16214 1 4194293 0 -4194292) 1)
    5. Trace: (LAGRANGE3 '(#:%I16214 1 1 0 1) '(#:%I16214 1 1 0 1) '8388587 '8388593)
    5. Trace: LAGRANGE3 ==> (#:%I16214 1 1 0 1)
    5. Trace: (LAGRANGE3 '(#:%I16214 1 4194296 0 -4194295) '(#:%I16214 1 4194293 0 -4194292) '8388587 '8388593)
    5. Trace: LAGRANGE3 ==> (#:%I16214 1 35184221094045 0 -35184221094044)
    5. Trace: (LAGRANGE3 '1 '1 '8388587 '8388593)
    5. Trace: LAGRANGE3 ==> 1
    5. Trace: (NEWGCD '(#:%I16214 1 1 0 2) '(#:%I16214 1 1 0 1) '8388581)
    6. Trace: (PGCDA '(#:%I16214 1 -1 0 1) '2 'T)
    6. Trace: PGCDA ==> (1)
    6. Trace: (PGCDA '1 '2 'T)
    6. Trace: PGCDA ==> (1 1 2)
    5. Trace: NEWGCD ==> ((#:%I16214 1 1 0 1) (#:%I16214 1 4194290 0 -4194289) 1)
    5. Trace: (LAGRANGE3 '(#:%I16214 1 1 0 1) '(#:%I16214 1 1 0 1) '8388581 '70368442188091)
    5. Trace: LAGRANGE3 ==> (#:%I16214 1 1 0 1)
    5. Trace: (LAGRANGE3 '(#:%I16214 1 35184221094045 0 -35184221094044) '(#:%I16214 1 4194290 0 -4194289) '8388581 '70368442188091)
    5. Trace: LAGRANGE3 ==> (#:%I16214 1 295145688569309294435 0 -295145688569309294434)
    5. Trace: (LAGRANGE3 '1 '1 '8388581 '70368442188091)
    5. Trace: LAGRANGE3 ==> 1
    5. Trace: (NEWGCD '(#:%I16214 1 1 0 2) '(#:%I16214 1 1 0 1) '8388571)
    6. Trace: (PGCDA '(#:%I16214 1 -1 0 1) '2 'T)
    6. Trace: PGCDA ==> (1)
    6. Trace: (PGCDA '1 '2 'T)
    6. Trace: PGCDA ==> (1 1 2)
    5. Trace: NEWGCD ==> ((#:%I16214 1 1 0 1) (#:%I16214 1 4194285 0 -4194284) 1)
    5. Trace: (LAGRANGE3 '(#:%I16214 1 1 0 1) '(#:%I16214 1 1 0 1) '8388571 '590291377138618588871)
    5. Trace: LAGRANGE3 ==> (#:%I16214 1 1 0 1)
    5. Trace: (LAGRANGE3 '(#:%I16214 1 295145688569309294435 0 -295145688569309294434) '(#:%I16214 1 4194285 0 -4194284) '8388571 '590291377138618588871)
    5. Trace: LAGRANGE3 ==> (#:%I16214 1 2475850563907539437332096670 0 -2475850563907539437332096669)
    5. Trace: (LAGRANGE3 '1 '1 '8388571 '590291377138618588871)
    5. Trace: LAGRANGE3 ==> 1

     

Log in to post a comment.

Get latest updates about Open Source Projects, Conferences and News.

Sign up for the SourceForge newsletter:





No, thanks