Denis Simon on Wed, 06 Jun 2018 15:04:12 +0200


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

Re: Maximum of a rational function of many variables in the unit cube


Dear Jacques,

For your question 1: you can try forvec which handles variable numbers of loops.

Denis SIMON

----- Mail original -----
> De: "Jacques Gélinas" <jacquesg00@hotmail.com>
> À: "pari-users" <pari-users@pari.math.u-bordeaux.fr>
> Envoyé: Mercredi 6 Juin 2018 02:50:34
> Objet: Maximum of a rational function of many variables in the unit cube

> Now I have certain rational functions of the coefficients (a,b,c,...) of a real
> polynomial of degree n>1
> (expressions which look like the convergents of continued fractions) such as
> 
> f1(a,b) = 1/2 * (n-2)/(n-1) * (5-a*b)  /  (3-a);
> 
> f2(a,b,c) = b  *  (3-a)  /  (5-a*b)  * (1 - 1/3*(n-3)/(n-1)*(1-a*b*c/7)/(1-a/3))
> \
>          / (1 - 1/2*(n-2)/(n-1)*(1-a*b/5)/(1-a/3) );
> 
> The coefficients are in the unit cube, 0<=a<=1, 0<=b<=1, ..., and I would like
> to
> use Pari/GP to test the conjecture that the maxima over the cube is in (0,1] for
> all n.
> 
> 1. Could I determine the maxima over the vertices, say (0,0),(0,1),(1,0),(1,1)
> of f1(a,b),
> without using a variable number of "for loops" or constructors like "vector" ?
> 
> 2. What functions are available in GP for a Monte-Carlo search inside the cube
> for the maxima ?
> 
> Jacques Gelinas