Karim Belabas on Fri, 02 Nov 2007 13:10:34 +0100


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

Re: Feature request: forvec with an initial point


* Olivier Ramare [2007-11-02 12:30]:
>  When computations are aborted in the middle
> of a forvec, it would be extremely interesting to
> be able to start them again from where they
> stopped.
>  I guess it is not so difficult to implement,
> though it will force the third argument of
> forvec to be present.

It is easy to implement, but the user interface would be cumbersome.
Maybe the following is enough ?

  forvec(X = v,
    if (lex(X, init) < 0, next);
    ...
  )

If also update v so that v[1] = init[1], the extra work should be
negligible.

Cheers,

    K.B.
--
Karim Belabas                  Tel: (+33) (0)5 40 00 26 17
Universite Bordeaux 1          Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation    http://www.math.u-bordeaux.fr/~belabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux.fr/  [PARI/GP]
`