Louis . Granboulan on Sun, 30 Nov 1997 20:10:14 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Patch of gen3.c (it is working, now...)


The previous patch did not handle correctly the bug.
If you apply this patch instead of the other, you
can type :
  intformal([x^3,x^2],x)
and the answer will be
  [1/4*x^4, 1/3*x^3]

-- Ptitboul

*** ./src/basemath/gen3.c.orig  Fri Nov 14 04:53:28 1997
--- ./src/basemath/gen3.c       Sun Nov 30 17:28:46 1997
***************
*** 1407,1413 ****
        return gerepileupto(av,y);
  
      case t_VEC: case t_COL: case t_MAT:
!       return triv_integ(x,v,tx,lx);
    }
    err(typeer,"integ");
    return NULL; /* not reached */
--- 1407,1413 ----
        return gerepileupto(av,y);
  
      case t_VEC: case t_COL: case t_MAT:
!       lx=lg(x);
!       y=cgetg(lx,tx);
!       for (i=1; i<lg(x); i++) y[i]=linteg((GEN)x[i],v);
!       return y;
    }
    err(typeer,"integ");
    return NULL; /* not reached */