Code coverage tests

This page documents the degree to which the PARI/GP source code is tested by our public test suite, distributed with the source distribution in directory src/test/. This is measured by the gcov utility; we then process gcov output using the lcov frond-end.

We test a few variants depending on Configure flags on the pari.math.u-bordeaux1.fr machine (x86_64 architecture), and agregate them in the final report:

The target is 90% coverage for all mathematical modules (given that branches depending on DEBUGLEVEL or DEBUGMEM are not covered). This script is run to produce the results below.

LCOV - code coverage report
Current view: top level - language - eval.c (source / functions) Hit Total Coverage
Test: PARI/GP v2.8.0 lcov report (development 16962-5a32637) Lines: 801 1290 62.1 %
Date: 2014-10-29 Functions: 66 97 68.0 %
Legend: Lines: hit not hit | Branches: + taken - not taken # not executed Branches: 377 804 46.9 %

           Branch data     Line data    Source code
       1                 :            : /* Copyright (C) 2006  The PARI group.
       2                 :            : 
       3                 :            : This file is part of the PARI package.
       4                 :            : 
       5                 :            : PARI/GP is free software; you can redistribute it and/or modify it under the
       6                 :            : terms of the GNU General Public License as published by the Free Software
       7                 :            : Foundation. It is distributed in the hope that it will be useful, but WITHOUT
       8                 :            : ANY WARRANTY WHATSOEVER.
       9                 :            : 
      10                 :            : Check the License for details. You should have received a copy of it, along
      11                 :            : with the package; see the file 'COPYING'. If not, write to the Free Software
      12                 :            : Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. */
      13                 :            : 
      14                 :            : #include "pari.h"
      15                 :            : #include "paripriv.h"
      16                 :            : #include "anal.h"
      17                 :            : #include "opcode.h"
      18                 :            : 
      19                 :            : /********************************************************************/
      20                 :            : /*                                                                  */
      21                 :            : /*                   break/next/return handling                     */
      22                 :            : /*                                                                  */
      23                 :            : /********************************************************************/
      24                 :            : 
      25                 :            : static THREAD long br_status, br_count;
      26                 :            : static THREAD GEN br_res;
      27                 :            : 
      28                 :            : long
      29                 :   29246305 : loop_break(void)
      30                 :            : {
      31   [ -  +  +  +  :   29246305 :   switch(br_status)
                      + ]
      32                 :            :   {
      33                 :            :     case br_MULTINEXT :
      34         [ #  # ]:          0 :       if (! --br_count) br_status = br_NEXT;
      35                 :          0 :       return 1;
      36         [ +  + ]:      70140 :     case br_BREAK : if (! --br_count) br_status = br_NONE; /* fall through */
      37                 :      70154 :     case br_RETURN: return 1;
      38                 :      64694 :     case br_NEXT: br_status = br_NONE; /* fall through */
      39                 :            :   }
      40                 :   29246305 :   return 0;
      41                 :            : }
      42                 :            : 
      43                 :            : static void
      44                 :      21670 : reset_break(void)
      45                 :            : {
      46                 :      21670 :   br_status = br_NONE;
      47         [ +  + ]:      21670 :   if (br_res) { gunclone_deep(br_res); br_res = NULL; }
      48                 :      21670 : }
      49                 :            : 
      50                 :            : GEN
      51                 :        154 : return0(GEN x)
      52                 :            : {
      53                 :        154 :   GEN y = br_res;
      54 [ +  + ][ +  - ]:        154 :   br_res = (x && x != gnil)? gcloneref(x): NULL;
      55         [ -  + ]:        154 :   if (y) gunclone_deep(y);
      56                 :        154 :   br_status = br_RETURN; return NULL;
      57                 :            : }
      58                 :            : 
      59                 :            : GEN
      60                 :      65415 : next0(long n)
      61                 :            : {
      62         [ +  + ]:      65415 :   if (n < 1) pari_err_DOMAIN("next", "n", "<", gen_1, stoi(n));
      63         [ +  - ]:      65408 :   if (n == 1) br_status = br_NEXT;
      64                 :            :   else
      65                 :            :   {
      66                 :          0 :     br_count = n-1;
      67                 :          0 :     br_status = br_MULTINEXT;
      68                 :            :   }
      69                 :      65408 :   return NULL;
      70                 :            : }
      71                 :            : 
      72                 :            : GEN
      73                 :      70175 : break0(long n)
      74                 :            : {
      75         [ +  + ]:      70175 :   if (n < 1) pari_err_DOMAIN("break", "n", "<", gen_1, stoi(n));
      76                 :      70168 :   br_count = n;
      77                 :      70168 :   br_status = br_BREAK; return NULL;
      78                 :            : }
      79                 :            : 
      80                 :            : /*******************************************************************/
      81                 :            : /*                                                                 */
      82                 :            : /*                            VARIABLES                            */
      83                 :            : /*                                                                 */
      84                 :            : /*******************************************************************/
      85                 :            : 
      86                 :            : /* As a rule, ep->value is a clone (COPY). push_val and pop_val are private
      87                 :            :  * functions for use in sumiter: we want a temporary ep->value, which is NOT
      88                 :            :  * a clone (PUSH), to avoid unnecessary copies. */
      89                 :            : 
      90                 :            : enum {PUSH_VAL = 0, COPY_VAL = 1, DEFAULT_VAL = 2};
      91                 :            : 
      92                 :            : /* ep->args is the stack of old values (INITIAL if initial value, from
      93                 :            :  * installep) */
      94                 :            : typedef struct var_cell {
      95                 :            :   struct var_cell *prev; /* cell associated to previous value on stack */
      96                 :            :   GEN value; /* last value (not including current one, in ep->value) */
      97                 :            :   char flag; /* status of _current_ ep->value: PUSH or COPY ? */
      98                 :            :   long valence; /* valence of entree* associated to 'value', to be restored
      99                 :            :                     * by pop_val */
     100                 :            : } var_cell;
     101                 :            : #define INITIAL NULL
     102                 :            : 
     103                 :            : /* Push x on value stack associated to ep. */
     104                 :            : static void
     105                 :       6622 : new_val_cell(entree *ep, GEN x, char flag)
     106                 :            : {
     107                 :       6622 :   var_cell *v = (var_cell*) pari_malloc(sizeof(var_cell));
     108                 :       6622 :   v->value  = (GEN)ep->value;
     109                 :       6622 :   v->prev   = (var_cell*) ep->pvalue;
     110                 :       6622 :   v->flag   = flag;
     111                 :       6622 :   v->valence= ep->valence;
     112                 :            : 
     113                 :            :   /* beware: f(p) = Nv = 0
     114                 :            :    *         Nv = p; f(Nv) --> this call would destroy p [ isclone ] */
     115 [ +  - ][ #  # ]:       6622 :   ep->value = (flag == COPY_VAL)? gclone(x):
     116         [ #  # ]:          0 :                                   (x && isclone(x))? gcopy(x): x;
     117                 :            :   /* Do this last. In case the clone is <C-C>'ed before completion ! */
     118                 :       6622 :   ep->pvalue= (char*)v;
     119                 :       6622 :   ep->valence=EpVAR;
     120                 :       6622 : }
     121                 :            : 
     122                 :            : /* kill ep->value and replace by preceding one, poped from value stack */
     123                 :            : static void
     124                 :       6587 : pop_val(entree *ep)
     125                 :            : {
     126                 :       6587 :   var_cell *v = (var_cell*) ep->pvalue;
     127                 :            : 
     128         [ -  + ]:      13174 :   if (v == INITIAL) return;
     129         [ +  + ]:       6587 :   if (v->flag == COPY_VAL) gunclone_deep((GEN)ep->value);
     130                 :       6587 :   ep->value  = v->value;
     131                 :       6587 :   ep->pvalue = (char*) v->prev;
     132                 :       6587 :   ep->valence=v->valence;
     133                 :       6587 :   pari_free((void*)v);
     134                 :            : }
     135                 :            : 
     136                 :            : void
     137                 :      15081 : freeep(entree *ep)
     138                 :            : {
     139         [ -  + ]:      30162 :   if (EpSTATIC(ep)) return; /* gp function loaded at init time */
     140         [ +  + ]:      15081 :   if (ep->help) {pari_free((void*)ep->help); ep->help=NULL;}
     141         [ +  + ]:      15081 :   if (ep->code) {pari_free((void*)ep->code); ep->code=NULL;}
     142      [ +  +  + ]:      15081 :   switch(EpVALENCE(ep))
     143                 :            :   {
     144                 :            :     case EpVAR:
     145         [ +  + ]:      18537 :       while (ep->pvalue!=INITIAL) pop_val(ep);
     146                 :      11992 :       break;
     147                 :            :     case EpALIAS:
     148                 :         28 :       killblock((GEN)ep->value); ep->value=NULL; break;
     149                 :            :   }
     150                 :            : }
     151                 :            : 
     152                 :            : INLINE void
     153                 :         28 : pushvalue(entree *ep, GEN x) {
     154                 :         28 :   new_val_cell(ep, x, COPY_VAL);
     155                 :         28 : }
     156                 :            : 
     157                 :            : INLINE void
     158                 :         14 : zerovalue (entree *ep)
     159                 :            : {
     160                 :         14 :   var_cell *v = (var_cell*) pari_malloc(sizeof(var_cell));
     161                 :         14 :   v->value  = (GEN)ep->value;
     162                 :         14 :   v->prev   = (var_cell*) ep->pvalue;
     163                 :         14 :   v->flag   = PUSH_VAL;
     164                 :         14 :   v->valence= ep->valence;
     165                 :         14 :   ep->value = gen_0;
     166                 :         14 :   ep->pvalue= (char*)v;
     167                 :         14 :   ep->valence=EpVAR;
     168                 :         14 : }
     169                 :            : 
     170                 :            : 
     171                 :            : /* as above IF ep->value was PUSHed, or was created after block number 'loc'
     172                 :            :    return 0 if not deleted, 1 otherwise [for recover()] */
     173                 :            : int
     174                 :     562179 : pop_val_if_newer(entree *ep, long loc)
     175                 :            : {
     176                 :     562179 :   var_cell *v = (var_cell*) ep->pvalue;
     177                 :            : 
     178         [ +  + ]:     562179 :   if (v == INITIAL) return 0;
     179 [ +  - ][ +  + ]:      17990 :   if (v->flag == COPY_VAL && !pop_entree_block(ep, loc)) return 0;
     180                 :         49 :   ep->value = v->value;
     181                 :         49 :   ep->pvalue= (char*) v->prev;
     182                 :         49 :   ep->valence=v->valence;
     183                 :     562179 :   pari_free((void*)v); return 1;
     184                 :            : }
     185                 :            : 
     186                 :            : /* set new value of ep directly to val (COPY), do not save last value unless
     187                 :            :  * it's INITIAL. */
     188                 :            : void
     189                 :    5658289 : changevalue(entree *ep, GEN x)
     190                 :            : {
     191                 :    5658289 :   var_cell *v = (var_cell*) ep->pvalue;
     192         [ +  + ]:    5658289 :   if (v == INITIAL) new_val_cell(ep, x, COPY_VAL);
     193                 :            :   else
     194                 :            :   {
     195                 :    5651695 :     GEN old_val = (GEN) ep->value; /* beware: gunclone_deep may destroy old x */
     196                 :    5651695 :     ep->value = (void *) gclone(x);
     197         [ +  - ]:    5651695 :     if (v->flag == COPY_VAL) gunclone_deep(old_val); else v->flag = COPY_VAL;
     198                 :            :   }
     199                 :    5658289 : }
     200                 :            : 
     201                 :            : INLINE GEN
     202                 :     737520 : copyvalue(entree *ep)
     203                 :            : {
     204                 :     737520 :   var_cell *v = (var_cell*) ep->pvalue;
     205 [ +  + ][ -  + ]:     737520 :   if (v && v->flag != COPY_VAL)
     206                 :            :   {
     207                 :          0 :     ep->value = (void*) gclone((GEN)ep->value);
     208                 :          0 :     v->flag = COPY_VAL;
     209                 :            :   }
     210                 :     737520 :   return (GEN) ep->value;
     211                 :            : }
     212                 :            : 
     213                 :            : INLINE void
     214                 :          0 : err_var(GEN x) { pari_err_TYPE("evaluator [variable name expected]", x); }
     215                 :            : 
     216                 :            : INLINE void
     217                 :   37447585 : checkvalue(entree *ep)
     218                 :            : {
     219         [ -  + ]:    5349651 :   if (MT_IS_THREAD)
     220                 :          0 :     pari_err(e_MISC,"mt: global variable not supported: %s",ep->name);
     221         [ +  + ]:   37447585 :   if (ep->valence==EpNEW)
     222                 :            :   {
     223                 :      10780 :     pari_var_create(ep);
     224                 :      10780 :     ep->valence = EpVAR;
     225                 :      10780 :     ep->value = initial_value(ep);
     226                 :            :   }
     227         [ -  + ]:   37436805 :   else if (ep->valence!=EpVAR)
     228                 :          0 :     err_var(strtoGENstr(ep->name));
     229                 :   37447585 : }
     230                 :            : 
     231                 :            : /* make GP variables safe for avma = top */
     232                 :            : static void
     233                 :          0 : lvar_make_safe(void)
     234                 :            : {
     235                 :            :   long n;
     236                 :            :   entree *ep;
     237         [ #  # ]:          0 :   for (n = 0; n < functions_tblsz; n++)
     238         [ #  # ]:          0 :     for (ep = functions_hash[n]; ep; ep = ep->next)
     239         [ #  # ]:          0 :       if (EpVALENCE(ep) == EpVAR)
     240                 :            :       { /* make sure ep->value is a COPY */
     241                 :          0 :         var_cell *v = (var_cell*)ep->pvalue;
     242 [ #  # ][ #  # ]:          0 :         if (v && v->flag == PUSH_VAL) {
     243                 :          0 :           GEN x = (GEN)ep->value;
     244         [ #  # ]:          0 :           if (x) changevalue(ep, (GEN)ep->value); else pop_val(ep);
     245                 :            :         }
     246                 :            :       }
     247                 :          0 : }
     248                 :            : 
     249                 :            : static void
     250                 :   29029403 : check_array_index(long c, long l)
     251                 :            : {
     252         [ +  + ]:   29029403 :   if (c < 1) pari_err_COMPONENT("", "<", gen_1, stoi(c));
     253         [ +  + ]:   29029396 :   if (c >= l) pari_err_COMPONENT("", ">", stoi(l-1), stoi(c));
     254                 :   29029354 : }
     255                 :            : 
     256                 :            : GEN*
     257                 :          0 : safegel(GEN x, long l)
     258                 :            : {
     259         [ #  # ]:          0 :   if (!is_matvec_t(typ(x)))
     260                 :          0 :     pari_err_TYPE("safegel",x);
     261                 :          0 :   check_array_index(l, lg(x));
     262                 :          0 :   return &(gel(x,l));
     263                 :            : }
     264                 :            : 
     265                 :            : GEN*
     266                 :          0 : safelistel(GEN x, long l)
     267                 :            : {
     268                 :            :   GEN d;
     269         [ #  # ]:          0 :   if (typ(x)!=t_LIST)
     270                 :          0 :     pari_err_TYPE("safelistel",x);
     271                 :          0 :   d = list_data(x);
     272                 :          0 :   check_array_index(l, lg(d));
     273                 :          0 :   return &(gel(d,l));
     274                 :            : }
     275                 :            : 
     276                 :            : long*
     277                 :          0 : safeel(GEN x, long l)
     278                 :            : {
     279         [ #  # ]:          0 :   if (typ(x)!=t_VECSMALL)
     280                 :          0 :     pari_err_TYPE("safeel",x);
     281                 :          0 :   check_array_index(l, lg(x));
     282                 :          0 :   return &(x[l]);
     283                 :            : }
     284                 :            : 
     285                 :            : GEN*
     286                 :          0 : safegcoeff(GEN x, long a, long b)
     287                 :            : {
     288         [ #  # ]:          0 :   if (typ(x)!=t_MAT) pari_err_TYPE("safegcoeff", x);
     289                 :          0 :   check_array_index(b, lg(x));
     290                 :          0 :   check_array_index(a, lg(gel(x,b)));
     291                 :          0 :   return &(gcoeff(x,a,b));
     292                 :            : }
     293                 :            : 
     294                 :            : typedef struct matcomp
     295                 :            : {
     296                 :            :   GEN *ptcell;
     297                 :            :   GEN parent;
     298                 :            :   int full_col, full_row;
     299                 :            : } matcomp;
     300                 :            : 
     301                 :            : typedef struct gp_pointer
     302                 :            : {
     303                 :            :   matcomp c;
     304                 :            :   GEN x, ox;
     305                 :            :   entree *ep;
     306                 :            :   long vn;
     307                 :            :   long sp;
     308                 :            : } gp_pointer;
     309                 :            : 
     310                 :            : 
     311                 :            : /* assign res at *pt in "simple array object" p and return it, or a copy.*/
     312                 :            : static void
     313                 :    8193038 : change_compo(matcomp *c, GEN res)
     314                 :            : {
     315                 :    8193038 :   GEN p = c->parent, *pt = c->ptcell;
     316                 :            :   long i, t;
     317                 :            : 
     318         [ +  + ]:    8193038 :   if (typ(p) == t_VECSMALL)
     319                 :            :   {
     320 [ +  + ][ +  + ]:         21 :     if (typ(res) != t_INT || is_bigint(res))
     321                 :         14 :       pari_err_TYPE("t_VECSMALL assignment", res);
     322                 :          7 :     *pt = (GEN)itos(res); return;
     323                 :            :   }
     324                 :    8193017 :   t = typ(res);
     325         [ +  + ]:    8193017 :   if (c->full_row)
     326                 :            :   {
     327         [ +  + ]:         21 :     if (t != t_VEC) pari_err_TYPE("matrix row assignment", res);
     328         [ +  - ]:          7 :     if (lg(res) != lg(p)) pari_err_DIM("matrix row assignment");
     329         [ #  # ]:          0 :     for (i=1; i<lg(p); i++)
     330                 :            :     {
     331         [ #  # ]:          0 :       GEN p1 = gcoeff(p,c->full_row,i); if (isclone(p1)) gunclone_deep(p1);
     332                 :          0 :       gcoeff(p,c->full_row,i) = gclone(gel(res,i));
     333                 :            :     }
     334                 :          0 :     return;
     335                 :            :   }
     336         [ +  + ]:    8192996 :   if (c->full_col)
     337                 :            :   {
     338         [ +  + ]:         42 :     if (t != t_COL) pari_err_TYPE("matrix col assignment", res);
     339         [ +  + ]:         28 :     if (lg(res) != lg(*pt)) pari_err_DIM("matrix col assignment");
     340                 :            :   }
     341                 :            : 
     342                 :    8192975 :   res = gclone(res);
     343                 :    8192975 :   gunclone_deep(*pt);
     344                 :    8192982 :   *pt = res;
     345                 :            : }
     346                 :            : 
     347                 :            : /***************************************************************************
     348                 :            :  **                                                                       **
     349                 :            :  **                           Byte-code evaluator                         **
     350                 :            :  **                                                                       **
     351                 :            :  ***************************************************************************/
     352                 :            : 
     353                 :            : struct var_lex
     354                 :            : {
     355                 :            :   long flag;
     356                 :            :   GEN value;
     357                 :            : };
     358                 :            : 
     359                 :            : struct trace
     360                 :            : {
     361                 :            :   long pc;
     362                 :            :   GEN closure;
     363                 :            : };
     364                 :            : 
     365                 :            : static THREAD long sp, rp, dbg_level;
     366                 :            : static THREAD long *st, *precs;
     367                 :            : static THREAD gp_pointer *ptrs;
     368                 :            : static THREAD entree **lvars;
     369                 :            : static THREAD struct var_lex *var;
     370                 :            : static THREAD struct trace *trace;
     371                 :            : static THREAD pari_stack s_st, s_ptrs, s_var, s_lvars, s_trace, s_prec;
     372                 :            : 
     373                 :            : static void
     374                 :   47998909 : changelex(long vn, GEN x)
     375                 :            : {
     376                 :   47998909 :   struct var_lex *v=var+s_var.n+vn;
     377                 :   47998909 :   GEN old_val = v->value;
     378                 :   47998909 :   v->value = gclone(x);
     379         [ +  + ]:   47998909 :   if (v->flag == COPY_VAL) gunclone_deep(old_val); else v->flag = COPY_VAL;
     380                 :   47998909 : }
     381                 :            : 
     382                 :            : INLINE GEN
     383                 :    8190686 : copylex(long vn)
     384                 :            : {
     385                 :    8190686 :   struct var_lex *v = var+s_var.n+vn;
     386         [ +  + ]:    8190686 :   if (v->flag!=COPY_VAL)
     387                 :            :   {
     388                 :          7 :     v->value = gclone(v->value);
     389                 :          7 :     v->flag  = COPY_VAL;
     390                 :            :   }
     391                 :    8190686 :   return v->value;
     392                 :            : }
     393                 :            : 
     394                 :            : INLINE void
     395                 :   16797438 : pushlex(long vn, GEN x)
     396                 :            : {
     397                 :   16797438 :   struct var_lex *v=var+s_var.n+vn;
     398                 :   16797438 :   v->flag  = PUSH_VAL;
     399                 :   16797438 :   v->value = x;
     400                 :   16797438 : }
     401                 :            : 
     402                 :            : INLINE void
     403                 :   42393147 : freelex(long vn)
     404                 :            : {
     405                 :   42393147 :   struct var_lex *v=var+s_var.n+vn;
     406         [ +  + ]:   42393147 :   if (v->flag == COPY_VAL) gunclone_deep(v->value);
     407                 :   42393147 : }
     408                 :            : 
     409                 :            : INLINE void
     410                 :   76869364 : restore_vars(long nbmvar, long nblvar)
     411                 :            : {
     412                 :            :   long j;
     413         [ +  + ]:  116675583 :   for(j=1;j<=nbmvar;j++)
     414                 :   39806219 :     freelex(-j);
     415                 :   76869364 :   s_var.n-=nbmvar;
     416         [ +  + ]:   76869406 :   for(j=1;j<=nblvar;j++)
     417                 :         42 :     pop_val(lvars[s_lvars.n-j]);
     418                 :   76869364 :   s_lvars.n-=nblvar;
     419                 :   76869364 : }
     420                 :            : 
     421                 :            : INLINE void
     422                 :      20662 : restore_trace(long nbtrace)
     423                 :            : {
     424                 :            :   long j;
     425         [ +  + ]:      52648 :   for(j=1;j<=nbtrace;j++)
     426                 :            :   {
     427                 :      31986 :     GEN C = trace[s_trace.n-j].closure;
     428         [ +  + ]:      31986 :     if (isclone(C)) gunclone(C);
     429                 :            :   }
     430                 :      20662 :   s_trace.n-=nbtrace;
     431                 :      20662 : }
     432                 :            : 
     433                 :            : INLINE long
     434                 :   79467341 : trace_push(long pc, GEN C)
     435                 :            : {
     436                 :            :   long tr;
     437         [ +  - ]:   79467341 :   BLOCK_SIGINT_START
     438                 :   79467341 :   tr = pari_stack_new(&s_trace);
     439                 :   79467341 :   trace[tr].pc = pc;
     440                 :   79467341 :   trace[tr].closure = C;
     441 [ +  - ][ -  + ]:   79467341 :   BLOCK_SIGINT_END
         [ +  - ][ -  + ]
     442                 :   79467341 :   return tr;
     443                 :            : }
     444                 :            : 
     445                 :            : void
     446                 :    2587045 : push_lex(GEN a, GEN C)
     447                 :            : {
     448                 :    2587045 :   long vn=pari_stack_new(&s_var);
     449                 :    2587045 :   struct var_lex *v=var+vn;
     450                 :    2587045 :   v->flag  = PUSH_VAL;
     451                 :    2587045 :   v->value = a;
     452         [ +  + ]:    2587045 :   if (C) (void) trace_push(-1, C);
     453                 :    2587045 : }
     454                 :            : 
     455                 :            : GEN
     456                 :   22132198 : get_lex(long vn)
     457                 :            : {
     458                 :   22132198 :   struct var_lex *v=var+s_var.n+vn;
     459                 :   22132198 :   return v->value;
     460                 :            : }
     461                 :            : 
     462                 :            : void
     463                 :   23365784 : set_lex(long vn, GEN x)
     464                 :            : {
     465                 :   23365784 :   struct var_lex *v=var+s_var.n+vn;
     466         [ -  + ]:   23365784 :   if (v->flag == COPY_VAL) { gunclone_deep(v->value); v->flag = PUSH_VAL; }
     467                 :   23365784 :   v->value = x;
     468                 :   23365784 : }
     469                 :            : 
     470                 :            : void
     471                 :    2586627 : pop_lex(long n)
     472                 :            : {
     473                 :            :   long j;
     474         [ +  + ]:    5173555 :   for(j=1; j<=n; j++)
     475                 :    2586928 :     freelex(-j);
     476                 :    2586627 :   s_var.n-=n;
     477                 :    2586627 :   s_trace.n--;
     478                 :    2586627 : }
     479                 :            : 
     480                 :            : static THREAD pari_stack s_relocs;
     481                 :            : static THREAD entree **relocs;
     482                 :            : 
     483                 :            : void
     484                 :       1201 : pari_init_evaluator(void)
     485                 :            : {
     486                 :       1201 :   sp=0;
     487                 :       1201 :   pari_stack_init(&s_st,sizeof(*st),(void**)&st);
     488                 :       1201 :   pari_stack_alloc(&s_st,32);
     489                 :       1201 :   s_st.n=s_st.alloc;
     490                 :       1201 :   rp=0;
     491                 :       1201 :   pari_stack_init(&s_ptrs,sizeof(*ptrs),(void**)&ptrs);
     492                 :       1201 :   pari_stack_alloc(&s_ptrs,16);
     493                 :       1201 :   s_ptrs.n=s_ptrs.alloc;
     494                 :       1201 :   pari_stack_init(&s_var,sizeof(*var),(void**)&var);
     495                 :       1201 :   pari_stack_init(&s_lvars,sizeof(*lvars),(void**)&lvars);
     496                 :       1201 :   pari_stack_init(&s_trace,sizeof(*trace),(void**)&trace);
     497                 :       1201 :   br_res = NULL;
     498                 :       1201 :   pari_stack_init(&s_relocs,sizeof(*relocs),(void**)&relocs);
     499                 :       1201 :   pari_stack_init(&s_prec,sizeof(*precs),(void**)&precs);
     500                 :       1201 : }
     501                 :            : void
     502                 :       1189 : pari_close_evaluator(void)
     503                 :            : {
     504                 :       1189 :   pari_stack_delete(&s_st);
     505                 :       1189 :   pari_stack_delete(&s_ptrs);
     506                 :       1189 :   pari_stack_delete(&s_var);
     507                 :       1189 :   pari_stack_delete(&s_lvars);
     508                 :       1189 :   pari_stack_delete(&s_trace);
     509                 :       1189 :   pari_stack_delete(&s_relocs);
     510                 :       1189 :   pari_stack_delete(&s_prec);
     511                 :       1189 : }
     512                 :            : 
     513                 :            : static gp_pointer *
     514                 :   34052137 : new_ptr(void)
     515                 :            : {
     516         [ -  + ]:   34052137 :   if (rp==s_ptrs.n-1)
     517                 :            :   {
     518                 :            :     long i;
     519                 :          0 :     gp_pointer *old = ptrs;
     520                 :          0 :     (void)pari_stack_new(&s_ptrs);
     521         [ #  # ]:          0 :     if (old != ptrs)
     522         [ #  # ]:          0 :       for(i=0; i<rp; i++)
     523                 :            :       {
     524                 :          0 :         gp_pointer *g = &ptrs[i];
     525         [ #  # ]:          0 :         if(g->sp >= 0) gel(st,g->sp) = (GEN) &(g->x);
     526                 :            :       }
     527                 :            :   }
     528                 :   34052137 :   return &ptrs[rp++];
     529                 :            : }
     530                 :            : 
     531                 :            : void
     532                 :         28 : push_localprec(long p)
     533                 :            : {
     534                 :         28 :   long n = pari_stack_new(&s_prec);
     535                 :         28 :   precs[n] = p;
     536                 :         28 : }
     537                 :            : 
     538                 :            : void
     539                 :         21 : pop_localprec(void)
     540                 :            : {
     541                 :         21 :   s_prec.n--;
     542                 :         21 : }
     543                 :            : 
     544                 :            : long
     545                 :   10619174 : get_localprec(void)
     546                 :            : {
     547         [ +  + ]:   10619174 :   return s_prec.n? precs[s_prec.n-1]: precreal;
     548                 :            : }
     549                 :            : 
     550                 :            : void
     551                 :          7 : localprec(long p)
     552                 :            : {
     553         [ -  + ]:          7 :   if (p < 1) pari_err_DOMAIN("localprec", "p", "<", gen_1, stoi(p));
     554         [ -  + ]:          7 :   if (p > prec2ndec(LGBITS))
     555                 :          0 :     pari_err_DOMAIN("localprec", "p", "==", utoi(LONG_MAX), stoi(p));
     556                 :          7 :   push_localprec(ndec2prec(p));
     557                 :          7 : }
     558                 :            : 
     559                 :            : INLINE GEN
     560                 :   11492987 : copyupto(GEN z, GEN t)
     561                 :            : {
     562 [ +  + ][ +  + ]:   11492987 :   if (is_universal_constant(z) || (z>(GEN)pari_mainstack->bot && z<=t))
                 [ +  + ]
     563                 :   10121982 :     return z;
     564                 :            :   else
     565                 :   11492987 :     return gcopy(z);
     566                 :            : }
     567                 :            : 
     568                 :            : static void closure_eval(GEN C);
     569                 :            : 
     570                 :            : INLINE GEN
     571                 :   10961677 : closure_return(GEN C)
     572                 :            : {
     573                 :   10961677 :   pari_sp ltop=avma;
     574                 :   10961677 :   closure_eval(C);
     575         [ +  + ]:   10944468 :   if (br_status)
     576                 :            :   {
     577                 :            :     GEN z;
     578                 :        868 :     avma=ltop;
     579         [ +  + ]:        868 :     z=br_res?gcopy(br_res):gnil;
     580                 :        868 :     reset_break();
     581                 :        868 :     return z;
     582                 :            :   }
     583                 :   10944468 :   return gerepileupto(ltop,gel(st,--sp));
     584                 :            : }
     585                 :            : 
     586                 :            : /* for the break_loop debugger. Not memory clean */
     587                 :            : GEN
     588                 :        105 : closure_evalbrk(GEN C, long *status)
     589                 :            : {
     590                 :        105 :   closure_eval(C);
     591                 :         91 :   *status = br_status;
     592         [ +  + ]:         91 :   if (br_status)
     593                 :            :   {
     594         [ -  + ]:         35 :     GEN z = br_res? gcopy(br_res): gnil;
     595                 :         35 :     reset_break();
     596                 :         35 :     return z;
     597                 :            :   }
     598                 :         91 :   return gel(st,--sp);
     599                 :            : }
     600                 :            : 
     601                 :            : INLINE long
     602                 :      42245 : closure_varn(GEN x)
     603                 :            : {
     604         [ +  + ]:      42245 :   if (!x) return -1;
     605         [ -  + ]:      42098 :   if (!gequalX(x)) err_var(x);
     606                 :      42245 :   return varn(x);
     607                 :            : }
     608                 :            : 
     609                 :            : INLINE void
     610                 :   20795835 : closure_castgen(GEN z, long mode)
     611                 :            : {
     612   [ +  -  -  -  :   20795835 :   switch (mode)
                   -  - ]
     613                 :            :   {
     614                 :            :   case Ggen:
     615                 :   20795835 :     gel(st,sp++)=z;
     616                 :   20795835 :     break;
     617                 :            :   case Gsmall:
     618                 :          0 :     st[sp++]=gtos(z);
     619                 :          0 :     break;
     620                 :            :   case Gusmall:
     621                 :          0 :     st[sp++]=gtou(z);
     622                 :          0 :     break;
     623                 :            :   case Gvar:
     624                 :          0 :     st[sp++]=closure_varn(z);
     625                 :          0 :     break;
     626                 :            :   case Gvoid:
     627                 :          0 :     break;
     628                 :            :   default:
     629                 :          0 :     pari_err_BUG("closure_castgen, type unknown");
     630                 :            :   }
     631                 :   20795835 : }
     632                 :            : 
     633                 :            : INLINE void
     634                 :          0 : closure_castlong(long z, long mode)
     635                 :            : {
     636   [ #  #  #  # ]:          0 :   switch (mode)
     637                 :            :   {
     638                 :            :   case Gsmall:
     639                 :          0 :     st[sp++]=z;
     640                 :          0 :     break;
     641                 :            :   case Ggen:
     642                 :          0 :     gel(st,sp++)=stoi(z);
     643                 :          0 :     break;
     644                 :            :   case Gvar:
     645                 :          0 :     err_var(stoi(z));
     646                 :            :   default:
     647                 :          0 :     pari_err_BUG("closure_castlong, type unknown");
     648                 :            :   }
     649                 :          0 : }
     650                 :            : 
     651                 :            : const char *
     652                 :       6088 : closure_func_err(void)
     653                 :            : {
     654                 :       6088 :   long fun=s_trace.n-1, pc;
     655                 :            :   const char *code;
     656                 :            :   GEN C, oper;
     657 [ +  + ][ +  + ]:       6088 :   if (fun < 0 || trace[fun].pc < 0) return NULL;
     658                 :       5908 :   pc = trace[fun].pc; C  = trace[fun].closure;
     659                 :       5908 :   code = closure_codestr(C); oper = closure_get_oper(C);
     660 [ +  + ][ +  + ]:       5908 :   if (code[pc]==OCcallgen || code[pc]==OCcallgen2 ||
                 [ +  - ]
     661 [ +  + ][ +  + ]:       3379 :       code[pc]==OCcallint || code[pc]==OCcalllong || code[pc]==OCcallvoid)
     662                 :       2968 :     return ((entree*)oper[pc])->name;
     663                 :       6088 :   return NULL;
     664                 :            : }
     665                 :            : 
     666                 :            : /* return the next label for the call chain debugger closure_err(),
     667                 :            :  * incorporating the name of the user of member function. Return NULL for an
     668                 :            :  * anonymous (inline) closure. */
     669                 :            : static char *
     670                 :        126 : get_next_label(const char *s, int member, char **next_fun)
     671                 :            : {
     672                 :        126 :   const char *v, *t = s+1;
     673                 :            :   char *u, *next_label;
     674                 :            : 
     675         [ -  + ]:        126 :   if (!is_keyword_char(*s)) return NULL;
     676         [ +  + ]:        420 :   while (is_keyword_char(*t)) t++;
     677                 :            :   /* e.g. (x->1/x)(0) instead of (x)->1/x */
     678 [ +  + ][ +  - ]:        126 :   if (t[0] == '-' && t[1] == '>') return NULL;
     679                 :        119 :   next_label = (char*)pari_malloc(t - s + 32);
     680         [ +  + ]:        119 :   sprintf(next_label, "in %sfunction ", member? "member ": "");
     681                 :        119 :   u = *next_fun = next_label + strlen(next_label);
     682                 :        119 :   v = s;
     683         [ +  + ]:        532 :   while (v < t) *u++ = *v++;
     684                 :        126 :   *u++ = 0; return next_label;
     685                 :            : }
     686                 :            : 
     687                 :            : static const char *
     688                 :         14 : get_arg_name(GEN C, long i)
     689                 :            : {
     690                 :         14 :   GEN e = gmael(closure_get_dbg(C), 3, 1);
     691                 :         14 :   return ((entree*)e[i])->name;
     692                 :            : }
     693                 :            : 
     694                 :            : void
     695                 :       5913 : closure_err(long level)
     696                 :            : {
     697                 :            :   GEN base;
     698                 :       5913 :   const long lastfun = s_trace.n - 1 - level;
     699                 :            :   char *next_label, *next_fun;
     700                 :       5913 :   long i = maxss(0, lastfun - 19);
     701         [ +  - ]:      11826 :   if (lastfun < 0) return; /*e.g. when called by gp_main_loop's simplify */
     702 [ -  + ][ #  # ]:       5913 :   if (i > 0) while (lg(trace[i].closure)==6) i--;
     703                 :       5913 :   base = closure_get_text(trace[i].closure); /* gcc -Wall*/
     704         [ +  - ]:       5913 :   next_label = pari_strdup(i == 0? "at top-level": "[...] at");
     705                 :       5913 :   next_fun = next_label;
     706         [ +  - ]:      12139 :   for (; i <= lastfun; i++)
     707                 :            :   {
     708                 :       6226 :     GEN C = trace[i].closure;
     709         [ +  + ]:       6226 :     if (lg(C) >= 7) base=closure_get_text(C);
     710 [ +  + ][ +  + ]:       6226 :     if ((i==lastfun || lg(trace[i+1].closure)>=7))
     711                 :            :     {
     712                 :       6039 :       GEN dbg = gel(closure_get_dbg(C),1);
     713                 :            :       /* After a SIGINT, pc can be slightly off: ensure 0 <= pc < lg() */
     714         [ +  + ]:       6039 :       long pc = minss(lg(dbg)-1, trace[i].pc>=0 ? trace[i].pc: 1);
     715         [ +  + ]:       6039 :       long offset = pc? dbg[pc]: 0;
     716                 :            :       int member;
     717                 :            :       const char *s, *sbase;
     718         [ +  + ]:       6039 :       if (typ(base)!=t_VEC) sbase = GSTR(base);
     719         [ +  + ]:        112 :       else if (offset>=0)   sbase = GSTR(gel(base,2));
     720                 :         14 :       else { sbase = GSTR(gel(base,1)); offset += strlen(sbase); }
     721                 :       6039 :       s = sbase + offset;
     722 [ +  + ][ +  + ]:       6039 :       member = offset>0 && (s[-1] == '.');
     723                 :            :       /* avoid "in function foo: foo" */
     724 [ +  + ][ +  + ]:       6039 :       if (!next_fun || strcmp(next_fun, s)) {
     725                 :       6032 :         print_errcontext(pariErr, next_label, s, sbase);
     726                 :       6032 :         out_putc(pariErr, '\n');
     727                 :            :       }
     728                 :       6039 :       pari_free(next_label);
     729         [ +  + ]:       6039 :       if (i == lastfun) break;
     730                 :            : 
     731                 :        126 :       next_label = get_next_label(s, member, &next_fun);
     732         [ +  + ]:        126 :       if (!next_label) {
     733                 :          7 :         next_label = pari_strdup("in anonymous function");
     734                 :          7 :         next_fun = NULL;
     735                 :            :       }
     736                 :            :     }
     737                 :            :   }
     738                 :            : }
     739                 :            : 
     740                 :            : long
     741                 :         56 : closure_context(long start, long level)
     742                 :            : {
     743                 :         56 :   const long lastfun = s_trace.n - 1 - level;
     744                 :         56 :   long i, fun = lastfun;
     745         [ -  + ]:         56 :   if (fun<0) return lastfun;
     746 [ +  + ][ +  + ]:        147 :   while (fun>start && lg(trace[fun].closure)==6) fun--;
     747         [ +  + ]:        203 :   for (i=fun; i <= lastfun; i++)
     748                 :        147 :     push_frame(trace[i].closure, trace[i].pc,0);
     749         [ +  + ]:         84 :   for (  ; i < s_trace.n; i++)
     750                 :         28 :     push_frame(trace[i].closure, trace[i].pc,1);
     751                 :         56 :   return s_trace.n-level;
     752                 :            : }
     753                 :            : 
     754                 :            : INLINE void
     755                 :  894697922 : st_alloc(long n)
     756                 :            : {
     757         [ +  + ]:  894697922 :   if (sp+n>s_st.n)
     758                 :            :   {
     759                 :         21 :     pari_stack_alloc(&s_st,n+16);
     760                 :         21 :     s_st.n=s_st.alloc;
     761         [ -  + ]:         21 :     if (DEBUGMEM>=2) pari_warn(warner,"doubling evaluator stack");
     762                 :            :   }
     763                 :  894697922 : }
     764                 :            : 
     765                 :            : INLINE void
     766                 :    8193017 : ptr_proplock(gp_pointer *g, GEN C)
     767                 :            : {
     768                 :    8193017 :   g->x = C;
     769         [ +  + ]:    8193017 :   if (isclone(g->x))
     770                 :            :   {
     771                 :     297283 :     clone_unlock(g->ox);
     772                 :     297283 :     g->ox = g->x;
     773                 :     297283 :     ++bl_refc(g->ox);
     774                 :            :   }
     775                 :    8193017 : }
     776                 :            : 
     777                 :            : static void
     778                 :   76880597 : closure_eval(GEN C)
     779                 :            : {
     780                 :   76880597 :   const char *code=closure_codestr(C);
     781                 :   76880597 :   GEN oper=closure_get_oper(C);
     782                 :   76880597 :   GEN data=closure_get_data(C);
     783                 :   76880597 :   long loper=lg(oper);
     784                 :   76880597 :   long saved_sp=sp-closure_arity(C);
     785                 :   76880597 :   long saved_rp=rp, saved_prec=s_prec.n;
     786                 :   76880597 :   long j, nbmvar=0, nblvar=0;
     787                 :            :   long pc, t;
     788                 :   76880597 :   clone_lock(C);
     789                 :   76880597 :   t = trace_push(0, C);
     790         [ +  + ]:   76880597 :   if (lg(C)==8)
     791                 :            :   {
     792                 :       1022 :     GEN z=closure_get_frame(C);
     793                 :       1022 :     long l=lg(z)-1;
     794                 :       1022 :     pari_stack_alloc(&s_var,l);
     795                 :       1022 :     s_var.n+=l;
     796                 :       1022 :     nbmvar+=l;
     797         [ +  + ]:       2240 :     for(j=1;j<=l;j++)
     798                 :            :     {
     799                 :       1218 :       var[s_var.n-j].flag=PUSH_VAL;
     800                 :       1218 :       var[s_var.n-j].value=gel(z,j);
     801                 :            :     }
     802                 :            :   }
     803                 :            : 
     804         [ +  + ]:  953355871 :   for(pc=1;pc<loper;pc++)
     805                 :            :   {
     806                 :  876778727 :     op_code opcode=(op_code) code[pc];
     807                 :  876778727 :     long operand=oper[pc];
     808         [ -  + ]:  876778727 :     if (sp<0) pari_err_BUG("closure_eval, stack underflow");
     809                 :  876778727 :     st_alloc(16);
     810                 :  876778727 :     trace[t].pc = pc;
     811                 :            :     CHECK_CTRLC
     812   [ +  +  +  +  :  876778727 :     switch(opcode)
          +  +  +  +  +  
          +  +  +  +  +  
          +  +  +  +  +  
          +  +  +  +  +  
          -  +  -  +  +  
          +  +  +  +  +  
          +  +  +  +  +  
          +  +  +  +  +  
          +  +  +  +  +  
          +  +  +  +  +  
          +  +  -  -  +  
             +  +  +  +  
                      - ]
     813                 :            :     {
     814                 :            :     case OCpushlong:
     815                 :   64617841 :       st[sp++]=operand;
     816                 :   64617841 :       break;
     817                 :            :     case OCpushgnil:
     818                 :      22307 :       gel(st,sp++)=gnil;
     819                 :      22307 :       break;
     820                 :            :     case OCpushgen:
     821                 :   50100422 :       gel(st,sp++)=gel(data,operand);
     822                 :   50100422 :       break;
     823                 :            :     case OCpushreal:
     824                 :      93791 :       gel(st,sp++)=strtor(GSTR(data[operand]),get_localprec());
     825                 :      93791 :       break;
     826                 :            :     case OCpushstoi:
     827                 :   89952595 :       gel(st,sp++)=stoi(operand);
     828                 :   89952595 :       break;
     829                 :            :     case OCpushvar:
     830                 :            :       {
     831                 :       3955 :         entree *ep = (entree *)operand;
     832                 :       3955 :         pari_var_create(ep);
     833                 :       3955 :         gel(st,sp++)=(GEN)initial_value(ep);
     834                 :       3955 :         break;
     835                 :            :       }
     836                 :            :     case OCpushdyn:
     837                 :            :       {
     838                 :   31051685 :         entree *ep = (entree *)operand;
     839                 :   31051685 :         checkvalue(ep);
     840                 :   31051685 :         gel(st,sp++)=(GEN)ep->value;
     841                 :   31051685 :         break;
     842                 :            :       }
     843                 :            :     case OCpushlex:
     844                 :  155611448 :       gel(st,sp++)=var[s_var.n+operand].value;
     845                 :  155611448 :       break;
     846                 :            :     case OCsimpleptrdyn:
     847                 :            :       {
     848                 :     496223 :         gp_pointer *g = new_ptr();
     849                 :     496223 :         g->vn=0;
     850                 :     496223 :         g->ep = (entree*) operand;
     851                 :     496223 :         checkvalue(g->ep);
     852                 :     496223 :         g->x = (GEN) g->ep->value;
     853                 :     496223 :         g->ox = g->x; clone_lock(g->ox);
     854                 :     496223 :         g->sp = sp;
     855                 :     496223 :         gel(st,sp++) = (GEN)&(g->x);
     856                 :     496223 :         break;
     857                 :            :       }
     858                 :            :     case OCsimpleptrlex:
     859                 :            :       {
     860                 :   25362827 :         gp_pointer *g = new_ptr();
     861                 :   25362827 :         g->vn=operand;
     862                 :   25362827 :         g->ep=(entree *)0x1L;
     863                 :   25362827 :         g->x = (GEN) var[s_var.n+operand].value;
     864                 :   25362827 :         g->ox = g->x; clone_lock(g->ox);
     865                 :   25362827 :         g->sp = sp;
     866                 :   25362827 :         gel(st,sp++) = (GEN)&(g->x);
     867                 :   25362827 :         break;
     868                 :            :       }
     869                 :            :     case OCnewptrdyn:
     870                 :            :       {
     871                 :       2436 :         entree *ep = (entree *)operand;
     872                 :       2436 :         gp_pointer *g = new_ptr();
     873                 :            :         matcomp *C;
     874                 :       2436 :         checkvalue(ep);
     875                 :       2436 :         g->sp = -1;
     876                 :       2436 :         g->x = copyvalue(ep);
     877                 :       2436 :         g->ox = g->x; clone_lock(g->ox);
     878                 :       2436 :         g->vn=0;
     879                 :       2436 :         g->ep=NULL;
     880                 :       2436 :         C=&g->c;
     881                 :       2436 :         C->full_col = C->full_row = 0;
     882                 :       2436 :         C->parent   = (GEN)    g->x;
     883                 :       2436 :         C->ptcell   = (GEN *) &g->x;
     884                 :       2436 :         break;
     885                 :            :       }
     886                 :            :     case OCnewptrlex:
     887                 :            :       {
     888                 :    8190651 :         gp_pointer *g = new_ptr();
     889                 :            :         matcomp *C;
     890                 :    8190651 :         g->sp = -1;
     891                 :    8190651 :         g->x = copylex(operand);
     892                 :    8190651 :         g->ox = g->x; clone_lock(g->ox);
     893                 :    8190651 :         g->vn=0;
     894                 :    8190651 :         g->ep=NULL;
     895                 :    8190651 :         C=&g->c;
     896                 :    8190651 :         C->full_col = C->full_row = 0;
     897                 :    8190651 :         C->parent   = (GEN)     g->x;
     898                 :    8190651 :         C->ptcell   = (GEN *) &(g->x);
     899                 :    8190651 :         break;
     900                 :            :       }
     901                 :            :     case OCpushptr:
     902                 :            :       {
     903                 :         35 :         gp_pointer *g = &ptrs[rp-1];
     904                 :         35 :         g->sp = sp;
     905                 :         35 :         gel(st,sp++) = (GEN)&(g->x);
     906                 :            :       }
     907                 :         35 :       break;
     908                 :            :     case OCendptr:
     909         [ +  + ]:   51718072 :       for(j=0;j<operand;j++)
     910                 :            :       {
     911                 :   25859036 :         gp_pointer *g = &ptrs[--rp];
     912         [ +  + ]:   25859036 :         if (g->ep)
     913                 :            :         {
     914         [ +  + ]:   25859001 :           if (g->vn)
     915                 :   25362827 :             changelex(g->vn, g->x);
     916                 :            :           else
     917                 :     496174 :             changevalue(g->ep, g->x);
     918                 :            :         }
     919                 :         35 :         else change_compo(&(g->c), g->x);
     920                 :   25859036 :         clone_unlock(g->ox);
     921                 :            :       }
     922                 :   25859036 :       break;
     923                 :            :     case OCstoredyn:
     924                 :            :       {
     925                 :    5162115 :         entree *ep = (entree *)operand;
     926                 :    5162115 :         checkvalue(ep);
     927                 :    5162115 :         changevalue(ep, gel(st,--sp));
     928                 :    5162115 :         break;
     929                 :            :       }
     930                 :            :     case OCstorelex:
     931                 :   22636082 :       changelex(operand,gel(st,--sp));
     932                 :   22636082 :       break;
     933                 :            :     case OCstoreptr:
     934                 :            :       {
     935                 :    8193003 :         gp_pointer *g = &ptrs[--rp];
     936                 :    8193003 :         change_compo(&(g->c), gel(st,--sp));
     937                 :    8192947 :         clone_unlock(g->ox);
     938                 :    8192947 :         break;
     939                 :            :       }
     940                 :            :     case OCstackgen:
     941                 :            :       {
     942                 :   11478084 :         GEN z = gerepileupto(st[sp-2],gel(st,sp-1));
     943                 :   11478084 :         gmael(st,sp-3,operand) = copyupto(z,gel(st,sp-2));
     944                 :   11478084 :         st[sp-2] = avma;
     945                 :   11478084 :         sp--;
     946                 :   11478084 :         break;
     947                 :            :       }
     948                 :            :     case OCprecreal:
     949                 :   10525383 :       st[sp++]=get_localprec();
     950                 :   10525383 :       break;
     951                 :            :     case OCprecdl:
     952                 :        371 :       st[sp++]=precdl;
     953                 :        371 :       break;
     954                 :            :     case OCavma:
     955                 :       1477 :       st[sp++]=avma;
     956                 :       1477 :       break;
     957                 :            :     case OCcowvardyn:
     958                 :            :       {
     959                 :     735084 :         entree *ep = (entree *)operand;
     960                 :     735084 :         checkvalue(ep);
     961                 :     735084 :         (void)copyvalue(ep);
     962                 :     735084 :         break;
     963                 :            :       }
     964                 :            :     case OCcowvarlex:
     965                 :         35 :       (void)copylex(operand);
     966                 :         35 :       break;
     967                 :            :     case OCstoi:
     968                 :   10951081 :       gel(st,sp-1)=stoi(st[sp-1]);
     969                 :   10951081 :       break;
     970                 :            :     case OCutoi:
     971                 :          0 :       gel(st,sp-1)=utoi(st[sp-1]);
     972                 :          0 :       break;
     973                 :            :     case OCitos:
     974                 :   25765819 :       st[sp+operand]=gtos(gel(st,sp+operand));
     975                 :   25765798 :       break;
     976                 :            :     case OCitou:
     977                 :          0 :       st[sp+operand]=gtou(gel(st,sp+operand));
     978                 :          0 :       break;
     979                 :            :     case OCtostr:
     980                 :            :       {
     981                 :       6057 :         GEN z = gel(st,sp+operand);
     982                 :       6057 :         st[sp+operand] = (long)GENtostr_unquoted(z);
     983                 :       6057 :         break;
     984                 :            :       }
     985                 :            :     case OCvarn:
     986                 :      42245 :       st[sp+operand] = closure_varn(gel(st,sp+operand));
     987                 :      42245 :       break;
     988                 :            :     case OCcopy:
     989                 :      32400 :       gel(st,sp-1) = gcopy(gel(st,sp-1));
     990                 :      32400 :       break;
     991                 :            :     case OCgerepile:
     992                 :            :     {
     993                 :            :       pari_sp av;
     994                 :            :       GEN x;
     995                 :       1477 :       sp--;
     996                 :       1477 :       av = st[sp-1];
     997                 :       1477 :       x = gel(st,sp);
     998         [ +  - ]:       1477 :       if (isonstack(x))
     999                 :            :       {
    1000                 :       1477 :         pari_sp av2 = (pari_sp)(x + lg(x));
    1001         [ -  + ]:       1477 :         if ((long) (av - av2) > 1000000L)
    1002                 :            :         {
    1003         [ #  # ]:          0 :           if (DEBUGMEM>=2)
    1004                 :          0 :             pari_warn(warnmem,"eval: recovering %ld bytes", av - av2);
    1005                 :          0 :           x = gerepileupto(av, x);
    1006                 :            :         }
    1007                 :          0 :       } else avma = av;
    1008                 :       1477 :       gel(st,sp-1) = x;
    1009                 :       1477 :       break;
    1010                 :            :     }
    1011                 :            :     case OCcopyifclone:
    1012         [ +  + ]:     490319 :       if (isclone(gel(st,sp-1)))
    1013                 :      92845 :         gel(st,sp-1) = gcopy(gel(st,sp-1));
    1014                 :     490319 :       break;
    1015                 :            :     case OCcompo1:
    1016                 :            :       {
    1017                 :   20791663 :         GEN  p=gel(st,sp-2);
    1018                 :   20791663 :         long c=st[sp-1];
    1019                 :   20791663 :         sp-=2;
    1020   [ +  -  +  + ]:   20791663 :         switch(typ(p))
    1021                 :            :         {
    1022                 :            :         case t_VEC: case t_COL:
    1023                 :   20791614 :           check_array_index(c, lg(p));
    1024                 :   20791614 :           closure_castgen(gel(p,c),operand);
    1025                 :   20791614 :           break;
    1026                 :            :         case t_LIST:
    1027                 :            :           {
    1028                 :            :             long lx;
    1029         [ #  # ]:          0 :             p = list_data(p); lx = p? lg(p): 1;
    1030                 :          0 :             check_array_index(c, lx);
    1031                 :          0 :             closure_castgen(gel(p,c),operand);
    1032                 :          0 :             break;
    1033                 :            :           }
    1034                 :            :         case t_VECSMALL:
    1035                 :         14 :           check_array_index(c,lg(p));
    1036                 :          0 :           closure_castlong(p[c],operand);
    1037                 :          0 :           break;
    1038                 :            :         default:
    1039                 :         35 :           pari_err_TYPE("_[_] OCcompo1 [not a vector]", p);
    1040                 :          0 :           break;
    1041                 :            :         }
    1042                 :   20791614 :         break;
    1043                 :            :       }
    1044                 :            :     case OCcompo1ptr:
    1045                 :            :       {
    1046                 :    8192996 :         long c=st[sp-1];
    1047                 :            :         long lx;
    1048                 :    8192996 :         gp_pointer *g = &ptrs[rp-1];
    1049                 :    8192996 :         matcomp *C=&g->c;
    1050                 :    8192996 :         GEN p = g->x;
    1051                 :    8192996 :         sp--;
    1052   [ +  +  +  + ]:    8192996 :         switch(typ(p))
    1053                 :            :         {
    1054                 :            :         case t_VEC: case t_COL:
    1055                 :    8192947 :           check_array_index(c, lg(p));
    1056                 :    8192947 :           C->ptcell = (GEN *) p+c;
    1057                 :    8192947 :           ptr_proplock(g, *(C->ptcell));
    1058                 :    8192947 :           break;
    1059                 :            :         case t_VECSMALL:
    1060                 :         28 :           check_array_index(c, lg(p));
    1061                 :         21 :           C->ptcell = (GEN *) p+c;
    1062                 :         21 :           g->x = stoi(p[c]);
    1063                 :         21 :           break;
    1064                 :            :         case t_LIST:
    1065         [ +  - ]:         14 :           p = list_data(p); lx = p? lg(p): 1;
    1066                 :         14 :           check_array_index(c,lx);
    1067                 :         14 :           C->ptcell = (GEN *) p+c;
    1068                 :         14 :           ptr_proplock(g, *(C->ptcell));
    1069                 :         14 :           break;
    1070                 :            :         default:
    1071                 :          7 :           pari_err_TYPE("&_[_] OCcompo1ptr [not a vector]", p);
    1072                 :            :         }
    1073                 :    8192982 :         C->parent   = p;
    1074                 :    8192982 :         break;
    1075                 :            :       }
    1076                 :            :     case OCcompo2:
    1077                 :            :       {
    1078                 :       4228 :         GEN  p=gel(st,sp-3);
    1079                 :       4228 :         long c=st[sp-2];
    1080                 :       4228 :         long d=st[sp-1];
    1081         [ +  + ]:       4228 :         if (typ(p)!=t_MAT) pari_err_TYPE("_[_,_] OCcompo2 [not a matrix]", p);
    1082                 :       4221 :         check_array_index(d, lg(p));
    1083                 :       4221 :         check_array_index(c, lg(gel(p,d)));
    1084                 :       4221 :         sp-=3;
    1085                 :       4221 :         closure_castgen(gcoeff(p,c,d),operand);
    1086                 :       4221 :         break;
    1087                 :            :       }
    1088                 :            :     case OCcompo2ptr:
    1089                 :            :       {
    1090                 :         14 :         long c=st[sp-2];
    1091                 :         14 :         long d=st[sp-1];
    1092                 :         14 :         gp_pointer *g = &ptrs[rp-1];
    1093                 :         14 :         matcomp *C=&g->c;
    1094                 :         14 :         GEN p = g->x;
    1095                 :         14 :         sp-=2;
    1096         [ -  + ]:         14 :         if (typ(p)!=t_MAT)
    1097                 :          0 :           pari_err_TYPE("&_[_,_] OCcompo2ptr [not a matrix]", p);
    1098                 :         14 :         check_array_index(d, lg(p));
    1099                 :         14 :         check_array_index(c, lg(gel(p,d)));
    1100                 :         14 :         C->ptcell = (GEN *) gel(p,d)+c;
    1101                 :         14 :         C->parent   = p;
    1102                 :         14 :         ptr_proplock(g, *(C->ptcell));
    1103                 :         14 :         break;
    1104                 :            :       }
    1105                 :            :     case OCcompoC:
    1106                 :            :       {
    1107                 :      36239 :         GEN  p=gel(st,sp-2);
    1108                 :      36239 :         long c=st[sp-1];
    1109         [ +  + ]:      36239 :         if (typ(p)!=t_MAT)
    1110                 :          7 :           pari_err_TYPE("_[,_] OCcompoC [not a matrix]", p);
    1111                 :      36232 :         check_array_index(c, lg(p));
    1112                 :      36225 :         sp--;
    1113                 :      36225 :         gel(st,sp-1) = gel(p,c);
    1114                 :      36225 :         break;
    1115                 :            :       }
    1116                 :            :     case OCcompoCptr:
    1117                 :            :       {
    1118                 :         56 :         long c=st[sp-1];
    1119                 :         56 :         gp_pointer *g = &ptrs[rp-1];
    1120                 :         56 :         matcomp *C=&g->c;
    1121                 :         56 :         GEN p = g->x;
    1122                 :         56 :         sp--;
    1123         [ +  + ]:         56 :         if (typ(p)!=t_MAT)
    1124                 :          7 :           pari_err_TYPE("&_[,_] OCcompoCptr [not a matrix]", p);
    1125                 :         49 :         check_array_index(c, lg(p));
    1126                 :         42 :         C->ptcell = (GEN *) p+c;
    1127                 :         42 :         C->full_col = c;
    1128                 :         42 :         C->parent   = p;
    1129                 :         42 :         ptr_proplock(g, *(C->ptcell));
    1130                 :         42 :         break;
    1131                 :            :       }
    1132                 :            :     case OCcompoL:
    1133                 :            :       {
    1134                 :         14 :         GEN  p=gel(st,sp-2);
    1135                 :         14 :         long r=st[sp-1];
    1136                 :         14 :         sp--;
    1137         [ +  + ]:         14 :         if (typ(p)!=t_MAT)
    1138                 :          7 :           pari_err_TYPE("_[_,] OCcompoL [not a matrix]", p);
    1139         [ -  + ]:          7 :         check_array_index(r,lg(p) == 1? 1: lgcols(p));
    1140                 :          0 :         gel(st,sp-1) = row(p,r);
    1141                 :          0 :         break;
    1142                 :            :       }
    1143                 :            :     case OCcompoLptr:
    1144                 :            :       {
    1145                 :         35 :         long r=st[sp-1];
    1146                 :         35 :         gp_pointer *g = &ptrs[rp-1];
    1147                 :         35 :         matcomp *C=&g->c;
    1148                 :         35 :         GEN p = g->x, p2;
    1149                 :         35 :         sp--;
    1150         [ +  + ]:         35 :         if (typ(p)!=t_MAT)
    1151                 :          7 :           pari_err_TYPE("&_[_,] OCcompoLptr [not a matrix]", p);
    1152         [ +  - ]:         28 :         check_array_index(r,lg(p) == 1? 1: lgcols(p));
    1153                 :         21 :         p2 = rowcopy(p,r);
    1154                 :         21 :         C->full_row = r; /* record row number */
    1155                 :         21 :         C->ptcell = &p2;
    1156                 :         21 :         C->parent   = p;
    1157                 :         21 :         g->x = p2;
    1158                 :            :         break;
    1159                 :            :       }
    1160                 :            :     case OCdefaultarg:
    1161         [ +  + ]:       2170 :       if (var[s_var.n+operand].flag==DEFAULT_VAL)
    1162                 :            :       {
    1163                 :       1309 :         GEN z = gel(st,sp-1);
    1164         [ +  + ]:       1309 :         pushlex(operand,typ(z)== t_CLOSURE? closure_evalnobrk(z): z);
    1165                 :            :       }
    1166                 :       2170 :       sp--;
    1167                 :       2170 :       break;
    1168                 :            :     case OClocalvar:
    1169                 :            :       {
    1170                 :         28 :         long n = pari_stack_new(&s_lvars);
    1171                 :         28 :         entree *ep = (entree *)operand;
    1172                 :         28 :         checkvalue(ep);
    1173                 :         28 :         lvars[n] = ep;
    1174                 :         28 :         nblvar++;
    1175                 :         28 :         pushvalue(ep,gel(st,--sp));
    1176                 :         28 :         break;
    1177                 :            :       }
    1178                 :            :     case OClocalvar0:
    1179                 :            :       {
    1180                 :         14 :         long n = pari_stack_new(&s_lvars);
    1181                 :         14 :         entree *ep = (entree *)operand;
    1182                 :         14 :         checkvalue(ep);
    1183                 :         14 :         lvars[n] = ep;
    1184                 :         14 :         nblvar++;
    1185                 :         14 :         zerovalue(ep);
    1186                 :         14 :         break;
    1187                 :            :       }
    1188                 :            : 
    1189                 :            : #define EVAL_f(f) \
    1190                 :            :   switch (ep->arity) \
    1191                 :            :   { \
    1192                 :            :     case 0: f(); break; \
    1193                 :            :     case 1: sp--; f(st[sp]); break; \
    1194                 :            :     case 2: sp-=2; f(st[sp],st[sp+1]); break; \
    1195                 :            :     case 3: sp-=3; f(st[sp],st[sp+1],st[sp+2]); break; \
    1196                 :            :     case 4: sp-=4; f(st[sp],st[sp+1],st[sp+2],st[sp+3]); break; \
    1197                 :            :     case 5: sp-=5; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4]); break; \
    1198                 :            :     case 6: sp-=6; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5]); break; \
    1199                 :            :     case 7: sp-=7; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6]); break; \
    1200                 :            :     case 8: sp-=8; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7]); break; \
    1201                 :            :     case 9: sp-=9; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8]); break; \
    1202                 :            :     case 10: sp-=10; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9]); break; \
    1203                 :            :     case 11: sp-=11; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10]); break; \
    1204                 :            :     case 12: sp-=12; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11]); break; \
    1205                 :            :     case 13: sp-=13; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11],st[sp+12]); break; \
    1206                 :            :     case 14: sp-=14; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11],st[sp+12],st[sp+13]); break; \
    1207                 :            :     case 15: sp-=15; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11],st[sp+12],st[sp+13],st[sp+14]); break; \
    1208                 :            :     case 16: sp-=16; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11],st[sp+12],st[sp+13],st[sp+14],st[sp+15]); break; \
    1209                 :            :     case 17: sp-=17; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11],st[sp+12],st[sp+13],st[sp+14],st[sp+15],st[sp+16]); break; \
    1210                 :            :     case 18: sp-=18; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11],st[sp+12],st[sp+13],st[sp+14],st[sp+15],st[sp+16],st[sp+17]); break; \
    1211                 :            :     case 19: sp-=19; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11],st[sp+12],st[sp+13],st[sp+14],st[sp+15],st[sp+16],st[sp+17],st[sp+18]); break; \
    1212                 :            :     case 20: sp-=20; f(st[sp],st[sp+1],st[sp+2],st[sp+3],st[sp+4],st[sp+5],st[sp+6],st[sp+7],st[sp+8],st[sp+9],st[sp+10],st[sp+11],st[sp+12],st[sp+13],st[sp+14],st[sp+15],st[sp+16],st[sp+17],st[sp+18],st[sp+19]); break; \
    1213                 :            :     default: \
    1214                 :            :       pari_err_IMPL("functions with more than 20 parameters");\
    1215                 :            :       goto endeval; /*not reached*/ \
    1216                 :            :   }
    1217                 :            : 
    1218                 :            :     case OCcallgen:
    1219                 :            :       {
    1220                 :   34579278 :         entree *ep = (entree *)operand;
    1221                 :            :         GEN res;
    1222                 :            :         /* Macro Madness : evaluate function ep->value on arguments
    1223                 :            :          * st[sp-ep->arity .. sp]. Set res = result. */
    1224   [ +  +  -  +  :   34579278 :         EVAL_f(res = ((GEN (*)(ANYARG))ep->value));
          +  +  +  +  -  
          -  -  -  -  -  
          -  -  -  -  -  
                -  -  - ]
    1225         [ +  + ]:   34568295 :         if (br_status) goto endeval;
    1226                 :   34432544 :         gel(st,sp++)=res;
    1227                 :   34432544 :         break;
    1228                 :            :       }
    1229                 :            :     case OCcallgen2: /*same for ep->arity = 2. Is this optimization worth it ?*/
    1230                 :            :       {
    1231                 :  152745715 :         entree *ep = (entree *)operand;
    1232                 :            :         GEN res;
    1233                 :  152745715 :         sp-=2;
    1234                 :  152745715 :         res = ((GEN (*)(GEN,GEN))ep->value)(gel(st,sp),gel(st,sp+1));
    1235         [ +  + ]:  152733157 :         if (br_status) goto endeval;
    1236                 :  152733150 :         gel(st,sp++)=res;
    1237                 :  152733150 :         break;
    1238                 :            :       }
    1239                 :            :     case OCcalllong:
    1240                 :            :       {
    1241                 :   10779463 :         entree *ep = (entree *)operand;
    1242                 :            :         long res;
    1243   [ +  +  +  +  :   10779463 :         EVAL_f(res = ((long (*)(ANYARG))ep->value));
          +  -  -  -  -  
          -  -  -  -  -  
          -  -  -  -  -  
                -  -  - ]
    1244         [ -  + ]:   10778888 :         if (br_status) goto endeval;
    1245                 :   10778888 :         st[sp++] = res;
    1246                 :   10778888 :         break;
    1247                 :            :       }
    1248                 :            :     case OCcallint:
    1249                 :            :       {
    1250                 :     173404 :         entree *ep = (entree *)operand;
    1251                 :            :         long res;
    1252   [ -  +  +  -  :     173404 :         EVAL_f(res = ((int (*)(ANYARG))ep->value));
          -  -  -  -  -  
          -  -  -  -  -  
          -  -  -  -  -  
                -  -  - ]
    1253         [ -  + ]:     173397 :         if (br_status) goto endeval;
    1254                 :     173397 :         st[sp++] = res;
    1255                 :     173397 :         break;
    1256                 :            :       }
    1257                 :            :     case OCcallvoid:
    1258                 :            :       {
    1259                 :   19130374 :         entree *ep = (entree *)operand;
    1260   [ +  +  +  +  :   19130374 :         EVAL_f(((void (*)(ANYARG))ep->value));
          +  +  +  -  -  
          -  -  -  -  -  
          -  -  -  -  -  
                -  -  - ]
    1261         [ +  + ]:   19129833 :         if (br_status) goto endeval;
    1262                 :   18994138 :         break;
    1263                 :            :       }
    1264                 :            : #undef EVAL_f
    1265                 :            : 
    1266                 :            :     case OCcalluser:
    1267                 :            :       {
    1268                 :   10901772 :         long n=operand;
    1269                 :   10901772 :         GEN fun = gel(st,sp-1-n);
    1270                 :            :         long arity, isvar;
    1271                 :            :         GEN z;
    1272         [ +  + ]:   10901772 :         if (typ(fun)!=t_CLOSURE) pari_err(e_NOTFUNC, fun);
    1273                 :   10899063 :         isvar = closure_is_variadic(fun);
    1274                 :   10899063 :         arity = closure_arity(fun);
    1275 [ +  + ][ +  + ]:   10899063 :         if (!isvar || n < arity)
    1276                 :            :         {
    1277                 :   10899021 :           st_alloc(arity-n);
    1278         [ -  + ]:   10899021 :           if (n>arity)
    1279                 :          0 :             pari_err(e_MISC,"too many parameters in user-defined function call");
    1280         [ +  + ]:   10904292 :           for (j=n+1;j<=arity;j++)
    1281                 :       5271 :             gel(st,sp++)=0;
    1282         [ +  + ]:   10899021 :           if (isvar) gel(st,sp-1) = cgetg(1,t_VEC);
    1283                 :            :         }
    1284                 :            :         else
    1285                 :            :         {
    1286                 :            :           GEN v;
    1287                 :         42 :           long j, m = n-arity+1;
    1288                 :         42 :           v = cgetg(m+1,t_VEC);
    1289                 :         42 :           sp-=m;
    1290         [ +  + ]:        210 :           for (j=1; j<=m; j++)
    1291         [ +  + ]:        168 :             gel(v,j) = gel(st,sp+j-1)? gcopy(gel(st,sp+j-1)): gen_0;
    1292                 :         42 :           gel(st,sp++)=v;
    1293                 :            :         }
    1294                 :            : #ifdef STACK_CHECK
    1295 [ +  - ][ -  + ]:   10899063 :         if (PARI_stack_limit && (void*) &z <= PARI_stack_limit)
    1296                 :          0 :           pari_err(e_MISC, "deep recursion");
    1297                 :            : #endif
    1298                 :   10899063 :         z = closure_return(fun);
    1299         [ +  - ]:   10894653 :         if (br_status) goto endeval;
    1300                 :   10894653 :         gel(st, sp-1) = z;
    1301                 :            :         break;
    1302                 :            :       }
    1303                 :            :     case OCnewframe:
    1304         [ +  - ]:   15514464 :       if (operand>0) nbmvar+=operand;
    1305                 :          0 :       else operand=-operand;
    1306                 :   15514464 :       pari_stack_alloc(&s_var,operand);
    1307                 :   15514464 :       s_var.n+=operand;
    1308         [ +  + ]:   38521588 :       for(j=1;j<=operand;j++)
    1309                 :            :       {
    1310                 :   23007124 :         var[s_var.n-j].flag=PUSH_VAL;
    1311                 :   23007124 :         var[s_var.n-j].value=gen_0;
    1312                 :            :       }
    1313                 :   15514464 :       break;
    1314                 :            :     case OCsaveframe:
    1315                 :            :       {
    1316         [ +  + ]:        105 :         GEN cl = (operand?gcopy:shallowcopy)(gel(st,sp-1));
    1317                 :        105 :         long l = lg(gel(cl,7));
    1318                 :        105 :         GEN  v = cgetg(l, t_VEC);
    1319         [ +  + ]:        252 :         for(j=1; j<l; j++)
    1320                 :            :         {
    1321                 :        147 :           GEN val = var[s_var.n-j].value;
    1322         [ +  + ]:        147 :           gel(v,j) = operand?gcopy(val):val;
    1323                 :            :         }
    1324                 :        105 :         gel(cl,7) = v;
    1325                 :        105 :         gel(st,sp-1) = cl;
    1326                 :            :       }
    1327                 :        105 :       break;
    1328                 :            :     case OCgetargs:
    1329                 :   10839738 :       pari_stack_alloc(&s_var,operand);
    1330                 :   10839738 :       s_var.n+=operand;
    1331                 :   10839738 :       nbmvar+=operand;
    1332                 :   10839738 :       sp-=operand;
    1333         [ +  + ]:   27637498 :       for (j=0;j<operand;j++)
    1334                 :            :       {
    1335         [ +  + ]:   16797760 :         if (gel(st,sp+j))
    1336                 :   16796129 :           pushlex(j-operand,gel(st,sp+j));
    1337                 :            :         else
    1338                 :            :         {
    1339                 :       1631 :           var[s_var.n+j-operand].flag=DEFAULT_VAL;
    1340                 :       1631 :           var[s_var.n+j-operand].value=gen_0;
    1341                 :            :         }
    1342                 :            :       }
    1343                 :   10839738 :       break;
    1344                 :            :     case OCcheckuserargs:
    1345         [ +  + ]:         77 :       for (j=0; j<operand; j++)
    1346         [ +  + ]:         56 :         if (var[s_var.n-operand+j].flag==DEFAULT_VAL)
    1347                 :         14 :           pari_err(e_MISC,"missing mandatory argument"
    1348                 :            :                    " '%s' in user function",get_arg_name(C,j+1));
    1349                 :         21 :       break;
    1350                 :            :     case OCcheckargs:
    1351         [ +  + ]:     166950 :       for (j=sp-1;operand;operand>>=1UL,j--)
    1352 [ +  + ][ -  + ]:      97272 :         if ((operand&1L) && gel(st,j)==NULL)
    1353                 :          0 :           pari_err(e_MISC,"missing mandatory argument");
    1354                 :      69678 :       break;
    1355                 :            :     case OCcheckargs0:
    1356         [ +  + ]:        336 :       for (j=sp-1;operand;operand>>=1UL,j--)
    1357 [ +  - ][ -  + ]:        168 :         if ((operand&1L) && gel(st,j))
    1358                 :          0 :           pari_err(e_MISC,"argument type not implemented");
    1359                 :        168 :       break;
    1360                 :            :     case OCdefaultlong:
    1361                 :       2884 :       sp--;
    1362         [ +  + ]:       2884 :       if (st[sp+operand])
    1363                 :        910 :         st[sp+operand]=gtos(gel(st,sp+operand));
    1364                 :            :       else
    1365                 :       1974 :         st[sp+operand]=st[sp];
    1366                 :       2884 :       break;
    1367                 :            :     case OCdefaultulong:
    1368                 :          0 :       sp--;
    1369         [ #  # ]:          0 :       if (st[sp+operand])
    1370                 :          0 :         st[sp+operand]=gtou(gel(st,sp+operand));
    1371                 :            :       else
    1372                 :          0 :         st[sp+operand]=st[sp];
    1373                 :          0 :       break;
    1374                 :            :     case OCdefaultgen:
    1375                 :          0 :       sp--;
    1376         [ #  # ]:          0 :       if (!st[sp+operand])
    1377                 :          0 :         st[sp+operand]=st[sp];
    1378                 :          0 :       break;
    1379                 :            :     case OCvec:
    1380                 :    3764543 :       gel(st,sp++)=cgetg(operand,t_VEC);
    1381                 :    3764543 :       st[sp++]=avma;
    1382                 :    3764543 :       break;
    1383                 :            :     case OCcol:
    1384                 :        504 :       gel(st,sp++)=cgetg(operand,t_COL);
    1385                 :        504 :       st[sp++]=avma;
    1386                 :        504 :       break;
    1387                 :            :     case OCmat:
    1388                 :            :       {
    1389                 :            :         GEN z;
    1390                 :       5754 :         long l=st[sp-1];
    1391                 :       5754 :         z=cgetg(operand,t_MAT);
    1392         [ +  + ]:      33313 :         for(j=1;j<operand;j++)
    1393                 :      27559 :           gel(z,j) = cgetg(l,t_COL);
    1394                 :       5754 :         gel(st,sp-1) = z;
    1395                 :       5754 :         st[sp++]=avma;
    1396                 :            :       }
    1397                 :       5754 :       break;
    1398                 :            :     case OCpop:
    1399                 :   34844941 :       sp-=operand;
    1400                 :   34844941 :       break;
    1401                 :            :     case OCdup:
    1402                 :            :       {
    1403                 :    7012166 :         long i, s=st[sp-1];
    1404                 :    7012166 :         st_alloc(operand);
    1405         [ +  + ]:   14030583 :         for(i=1;i<=operand;i++)
    1406                 :    7018417 :           st[sp++]=s;
    1407                 :            :       }
    1408                 :    7012166 :       break;
    1409                 :            :     }
    1410                 :            :   }
    1411                 :            :   if (0)
    1412                 :            :   {
    1413                 :            : endeval:
    1414                 :     271453 :     sp = saved_sp;
    1415         [ -  + ]:     271453 :     for(  ; rp>saved_rp ;  )
    1416                 :            :     {
    1417                 :          0 :       gp_pointer *g = &ptrs[--rp];
    1418                 :          0 :       clone_unlock(g->ox);
    1419                 :            :     }
    1420                 :            :   }
    1421                 :   76848597 :   s_prec.n = saved_prec;
    1422                 :   76848597 :   s_trace.n--;
    1423                 :   76848597 :   restore_vars(nbmvar, nblvar);
    1424                 :   76848597 :   clone_unlock(C);
    1425                 :   76848597 : }
    1426                 :            : 
    1427                 :            : GEN
    1428                 :    9375428 : closure_evalgen(GEN C)
    1429                 :            : {
    1430                 :    9375428 :   pari_sp ltop=avma;
    1431                 :    9375428 :   closure_eval(C);
    1432         [ +  + ]:    9360714 :   if (br_status) { avma=ltop; return NULL; }
    1433                 :    9360714 :   return gerepileupto(ltop,gel(st,--sp));
    1434                 :            : }
    1435                 :            : 
    1436                 :            : void
    1437                 :    1278870 : evalstate_save(struct pari_evalstate *state)
    1438                 :            : {
    1439                 :    1278870 :   state->avma = avma;
    1440                 :    1278870 :   state->sp   = sp;
    1441                 :    1278870 :   state->rp   = rp;
    1442                 :    1278870 :   state->prec = s_prec.n;
    1443                 :    1278870 :   state->var  = s_var.n;
    1444                 :    1278870 :   state->lvars= s_lvars.n;
    1445                 :    1278870 :   state->trace= s_trace.n;
    1446                 :    1278870 :   compilestate_save(&state->comp);
    1447                 :    1278870 :   mtstate_save(&state->pending_threads);
    1448                 :    1278870 : }
    1449                 :            : 
    1450                 :            : void
    1451                 :      20662 : evalstate_restore(struct pari_evalstate *state)
    1452                 :            : {
    1453                 :      20662 :   avma = state->avma;
    1454                 :      20662 :   mtstate_restore(&state->pending_threads);
    1455                 :      20662 :   sp = state->sp;
    1456                 :      20662 :   rp = state->rp;
    1457                 :      20662 :   s_prec.n = state->prec;
    1458                 :      20662 :   restore_vars(s_var.n-state->var,s_lvars.n-state->lvars);
    1459                 :      20662 :   restore_trace(s_trace.n-state->trace);
    1460                 :      20662 :   reset_break();
    1461                 :      20662 :   compilestate_restore(&state->comp);
    1462                 :      20662 : }
    1463                 :            : 
    1464                 :            : GEN
    1465                 :      14700 : evalstate_restore_err(struct pari_evalstate *state)
    1466                 :            : {
    1467                 :      14700 :   GENbin* err = copy_bin(pari_err_last());
    1468                 :      14700 :   evalstate_restore(state);
    1469                 :      14700 :   return bin_copy(err);
    1470                 :            : }
    1471                 :            : 
    1472                 :            : void
    1473                 :        105 : evalstate_reset(void)
    1474                 :            : {
    1475                 :        105 :   mtstate_reset();
    1476                 :        105 :   sp = 0;
    1477                 :        105 :   rp = 0;
    1478                 :        105 :   dbg_level = 0;
    1479                 :        105 :   restore_vars(s_var.n, s_lvars.n);
    1480                 :        105 :   s_trace.n = 0;
    1481                 :        105 :   reset_break();
    1482                 :        105 :   compilestate_reset();
    1483                 :        105 :   parsestate_reset();
    1484                 :        105 :   avma = pari_mainstack->top;
    1485                 :        105 : }
    1486                 :            : 
    1487                 :            : void
    1488                 :          0 : evalstate_clone(void)
    1489                 :            : {
    1490                 :            :   long i;
    1491         [ #  # ]:          0 :   for (i = 1; i<=s_var.n; i++) copylex(-i);
    1492                 :          0 :   lvar_make_safe();
    1493         [ #  # ]:          0 :   for (i = 0; i< s_trace.n; i++)
    1494                 :            :   {
    1495                 :          0 :     GEN C = trace[i].closure;
    1496         [ #  # ]:          0 :     if (isonstack(C)) trace[i].closure = gclone(C);
    1497                 :            :   }
    1498                 :          0 : }
    1499                 :            : 
    1500                 :            : GEN
    1501                 :         21 : closure_trapgen(GEN C, long numerr)
    1502                 :            : {
    1503                 :            :   VOLATILE GEN x;
    1504                 :            :   struct pari_evalstate state;
    1505                 :         21 :   evalstate_save(&state);
    1506 [ +  + ][ +  - ]:         21 :   pari_CATCH(numerr) { x = (GEN)1L; }
                 [ +  + ]
    1507                 :         21 :   pari_TRY { x = closure_evalgen(C); } pari_ENDCATCH;
    1508         [ +  + ]:         14 :   if (x == (GEN)1L) evalstate_restore(&state);
    1509                 :         14 :   return x;
    1510                 :            : }
    1511                 :            : 
    1512                 :            : GEN
    1513                 :   21195798 : closure_evalnobrk(GEN C)
    1514                 :            : {
    1515                 :   21195798 :   pari_sp ltop=avma;
    1516                 :   21195798 :   closure_eval(C);
    1517         [ -  + ]:   21195791 :   if (br_status) pari_err(e_MISC, "break not allowed here");
    1518                 :   21195791 :   return gerepileupto(ltop,gel(st,--sp));
    1519                 :            : }
    1520                 :            : 
    1521                 :            : void
    1522                 :   35347589 : closure_evalvoid(GEN C)
    1523                 :            : {
    1524                 :   35347589 :   pari_sp ltop=avma;
    1525                 :   35347589 :   closure_eval(C);
    1526                 :   35347533 :   avma=ltop;
    1527                 :   35347533 : }
    1528                 :            : 
    1529                 :            : GEN
    1530                 :      52156 : closure_evalres(GEN C)
    1531                 :            : {
    1532                 :      52156 :   return closure_return(C);
    1533                 :            : }
    1534                 :            : 
    1535                 :            : INLINE GEN
    1536                 :      10458 : closure_returnupto(GEN C)
    1537                 :            : {
    1538                 :      10458 :   pari_sp av=avma;
    1539                 :      10458 :   return copyupto(closure_return(C),(GEN)av);
    1540                 :            : }
    1541                 :            : 
    1542                 :            : void
    1543                 :          0 : closure_callvoid1(GEN C, GEN x)
    1544                 :            : {
    1545                 :          0 :   long i, ar = closure_arity(C);
    1546                 :          0 :   gel(st,sp++) = x;
    1547         [ #  # ]:          0 :   for(i=2; i <= ar; i++) gel(st,sp++) = NULL;
    1548                 :          0 :   closure_evalvoid(C);
    1549                 :          0 : }
    1550                 :            : 
    1551                 :            : GEN
    1552                 :       2450 : closure_callgen1(GEN C, GEN x)
    1553                 :            : {
    1554                 :       2450 :   long i, ar = closure_arity(C);
    1555                 :       2450 :   gel(st,sp++) = x;
    1556         [ +  + ]:       2884 :   for(i=2; i<= ar; i++) gel(st,sp++) = NULL;
    1557                 :       2450 :   return closure_returnupto(C);
    1558                 :            : }
    1559                 :            : 
    1560                 :            : GEN
    1561                 :          0 : pareval_worker(GEN C)
    1562                 :            : {
    1563                 :          0 :   return closure_callgenall(C, 0);
    1564                 :            : }
    1565                 :            : 
    1566                 :            : GEN
    1567                 :          0 : pareval(GEN C)
    1568                 :            : {
    1569                 :          0 :   pari_sp av = avma;
    1570                 :          0 :   long l = lg(C), i, pending = 0, workid;
    1571                 :            :   struct pari_mt pt;
    1572                 :            :   GEN worker, V, done;
    1573         [ #  # ]:          0 :   if (!is_vec_t(typ(C))) pari_err_TYPE("pareval",C);
    1574         [ #  # ]:          0 :   for (i=1; i<l; i++)
    1575         [ #  # ]:          0 :     if (typ(gel(C,i))!=t_CLOSURE)
    1576                 :          0 :       pari_err_TYPE("pareval",gel(C,i));
    1577                 :          0 :   worker = snm_closure(is_entry("_pareval_worker"), NULL);
    1578                 :          0 :   V = cgetg(l, t_VEC);
    1579                 :          0 :   mt_queue_start(&pt, worker);
    1580 [ #  # ][ #  # ]:          0 :   for (i=1; i<l || pending; i++)
    1581                 :            :   {
    1582         [ #  # ]:          0 :     mt_queue_submit(&pt, i, i<l? mkvec(gel(C,i)): NULL);
    1583                 :          0 :     done = mt_queue_get(&pt, &workid, &pending);
    1584         [ #  # ]:          0 :     if (done) gel(V,workid) = done;
    1585                 :            :   }
    1586                 :          0 :   mt_queue_end(&pt);
    1587                 :          0 :   return gerepilecopy(av, V);
    1588                 :            : }
    1589                 :            : 
    1590                 :            : GEN
    1591                 :          0 : parvector_worker(GEN i, GEN C)
    1592                 :            : {
    1593                 :          0 :   return closure_callgen1(C, i);
    1594                 :            : }
    1595                 :            : 
    1596                 :            : GEN
    1597                 :          0 : parfor_worker(GEN i, GEN C)
    1598                 :            : {
    1599                 :          0 :   retmkvec2(gcopy(i), closure_callgen1(C, i));
    1600                 :            : }
    1601                 :            : 
    1602                 :            : GEN
    1603                 :          0 : parvector(long n, GEN code)
    1604                 :            : {
    1605                 :          0 :   long i, pending = 0, workid;
    1606                 :          0 :   GEN worker = snm_closure(is_entry("_parvector_worker"), mkvec(code));
    1607                 :            :   GEN a, V, done;
    1608                 :            :   struct pari_mt pt;
    1609                 :          0 :   mt_queue_start(&pt, worker);
    1610                 :          0 :   a = mkvec(cgetipos(3)); /* left on the stack */
    1611                 :          0 :   V = cgetg(n+1, t_VEC);
    1612 [ #  # ][ #  # ]:          0 :   for (i=1; i<=n || pending; i++)
    1613                 :            :   {
    1614                 :          0 :     mael(a,1,2) = i;
    1615         [ #  # ]:          0 :     mt_queue_submit(&pt, i, i<=n? a: NULL);
    1616                 :          0 :     done = mt_queue_get(&pt, &workid, &pending);
    1617         [ #  # ]:          0 :     if (done) gel(V,workid) = done;
    1618                 :            :   }
    1619                 :          0 :   mt_queue_end(&pt);
    1620                 :          0 :   return V;
    1621                 :            : }
    1622                 :            : 
    1623                 :            : GEN
    1624                 :          0 : parsum(GEN a, GEN b, GEN code, GEN x)
    1625                 :            : {
    1626                 :          0 :   pari_sp av = avma, av2;
    1627                 :          0 :   long pending = 0;
    1628                 :          0 :   GEN worker = snm_closure(is_entry("_parvector_worker"), mkvec(code));
    1629                 :            :   GEN done;
    1630                 :            :   struct pari_mt pt;
    1631         [ #  # ]:          0 :   if (typ(a) != t_INT) pari_err_TYPE("parsum",a);
    1632         [ #  # ]:          0 :   if (!x) x = gen_0;
    1633         [ #  # ]:          0 :   if (gcmp(b,a) < 0) return gcopy(x);
    1634                 :            : 
    1635                 :          0 :   mt_queue_start(&pt, worker);
    1636                 :          0 :   b = gfloor(b);
    1637                 :          0 :   a = mkvec(setloop(a));
    1638                 :          0 :   av2=avma;
    1639 [ #  # ][ #  # ]:          0 :   for (; cmpii(gel(a,1),b) <= 0 || pending; gel(a,1) = incloop(gel(a,1)))
    1640                 :            :   {
    1641         [ #  # ]:          0 :     mt_queue_submit(&pt, 0, cmpii(gel(a,1),b) <= 0? a: NULL);
    1642                 :          0 :     done = mt_queue_get(&pt, NULL, &pending);
    1643         [ #  # ]:          0 :     if (done)
    1644                 :            :     {
    1645                 :          0 :       x = gadd(x, done);
    1646         [ #  # ]:          0 :       if (gc_needed(av2,1))
    1647                 :            :       {
    1648         [ #  # ]:          0 :         if (DEBUGMEM>1) pari_warn(warnmem,"sum");
    1649                 :          0 :         x = gerepileupto(av2,x);
    1650                 :            :       }
    1651                 :            :     }
    1652                 :            :   }
    1653                 :          0 :   mt_queue_end(&pt);
    1654                 :          0 :   return gerepilecopy(av, x);
    1655                 :            : }
    1656                 :            : 
    1657                 :            : void
    1658                 :          0 : parfor(GEN a, GEN b, GEN code, GEN code2)
    1659                 :            : {
    1660                 :          0 :   pari_sp av = avma, av2;
    1661                 :          0 :   long running, pending = 0;
    1662                 :          0 :   long status = br_NONE;
    1663                 :          0 :   GEN worker = snm_closure(is_entry("_parfor_worker"), mkvec(code));
    1664                 :          0 :   GEN done, stop = NULL;
    1665                 :            :   struct pari_mt pt;
    1666         [ #  # ]:          0 :   if (typ(a) != t_INT) pari_err_TYPE("parfor",a);
    1667 [ #  # ][ #  # ]:          0 :   if (b && gcmp(b,a) < 0) return;
    1668                 :            : 
    1669                 :          0 :   mt_queue_start(&pt, worker);
    1670         [ #  # ]:          0 :   b = b ? gfloor(b): NULL;
    1671                 :          0 :   a = mkvec(setloop(a));
    1672                 :          0 :   av2 = avma;
    1673 [ #  # ][ #  # ]:          0 :   while ((running = (!stop && (!b || cmpii(gel(a,1),b) <= 0))) || pending)
         [ #  # ][ #  # ]
                 [ #  # ]
    1674                 :            :   {
    1675         [ #  # ]:          0 :     mt_queue_submit(&pt, 0, running ? a: NULL);
    1676                 :          0 :     done = mt_queue_get(&pt, NULL, &pending);
    1677 [ #  # ][ #  # ]:          0 :     if (code2 && done && (!stop || cmpii(gel(done,1),stop) < 0))
         [ #  # ][ #  # ]
    1678                 :            :     {
    1679                 :          0 :       push_lex(gel(done,1), code2);
    1680                 :          0 :       push_lex(gel(done,2), NULL);
    1681                 :          0 :       closure_evalvoid(code2);
    1682                 :          0 :       pop_lex(2);
    1683         [ #  # ]:          0 :       if (loop_break())
    1684                 :            :       {
    1685                 :          0 :         status = br_status;
    1686                 :          0 :         br_status = br_NONE;
    1687                 :          0 :         stop = gerepileuptoint(av2, gel(done,1));
    1688                 :            :       }
    1689                 :            :     }
    1690                 :          0 :     gel(a,1) = incloop(gel(a,1));
    1691         [ #  # ]:          0 :     if (!stop) avma = av2;
    1692                 :            :   }
    1693                 :          0 :   avma = av2;
    1694                 :          0 :   mt_queue_end(&pt);
    1695                 :          0 :   br_status = status;
    1696                 :          0 :   avma = av;
    1697                 :            : }
    1698                 :            : 
    1699                 :            : void
    1700                 :          0 : parforprime(GEN a, GEN b, GEN code, GEN code2)
    1701                 :            : {
    1702                 :          0 :   pari_sp av = avma, av2;
    1703                 :          0 :   long running, pending = 0;
    1704                 :          0 :   long status = br_NONE;
    1705                 :          0 :   GEN worker = snm_closure(is_entry("_parfor_worker"), mkvec(code));
    1706                 :          0 :   GEN done, stop = NULL;
    1707                 :            :   struct pari_mt pt;
    1708                 :            :   forprime_t T;
    1709                 :            : 
    1710         [ #  # ]:          0 :   if (!forprime_init(&T, a,b)) { avma = av; return; }
    1711                 :          0 :   mt_queue_start(&pt, worker);
    1712                 :          0 :   av2 = avma;
    1713 [ #  # ][ #  # ]:          0 :   while ((running = (!stop && forprime_next(&T))) || pending)
         [ #  # ][ #  # ]
    1714                 :            :   {
    1715         [ #  # ]:          0 :     mt_queue_submit(&pt, 0, running ? mkvec(T.pp): NULL);
    1716                 :          0 :     done = mt_queue_get(&pt, NULL, &pending);
    1717 [ #  # ][ #  # ]:          0 :     if (code2 && done && (!stop || cmpii(gel(done,1),stop) < 0))
         [ #  # ][ #  # ]
    1718                 :            :     {
    1719                 :          0 :       push_lex(gel(done,1), code2);
    1720                 :          0 :       push_lex(gel(done,2), NULL);
    1721                 :          0 :       closure_evalvoid(code2);
    1722                 :          0 :       pop_lex(2);
    1723         [ #  # ]:          0 :       if (loop_break())
    1724                 :            :       {
    1725                 :          0 :         status = br_status;
    1726                 :          0 :         br_status = br_NONE;
    1727                 :          0 :         stop = gerepileuptoint(av2, gel(done,1));
    1728                 :            :       }
    1729                 :            :     }
    1730         [ #  # ]:          0 :     if (!stop) avma = av2;
    1731                 :            :   }
    1732                 :          0 :   avma = av2;
    1733                 :          0 :   mt_queue_end(&pt);
    1734                 :          0 :   br_status = status;
    1735                 :          0 :   avma = av;
    1736                 :            : }
    1737                 :            : 
    1738                 :            : GEN
    1739                 :       7994 : closure_callgen2(GEN C, GEN x, GEN y)
    1740                 :            : {
    1741                 :       7994 :   long i, ar = closure_arity(C);
    1742                 :       7994 :   st_alloc(ar);
    1743                 :       7994 :   gel(st,sp++) = x;
    1744                 :       7994 :   gel(st,sp++) = y;
    1745         [ -  + ]:       7994 :   for(i=3; i<=ar; i++) gel(st,sp++) = NULL;
    1746                 :       7994 :   return closure_returnupto(C);
    1747                 :            : }
    1748                 :            : 
    1749                 :            : GEN
    1750                 :         14 : closure_callgenvec(GEN C, GEN args)
    1751                 :            : {
    1752                 :         14 :   long i, l = lg(args), ar = closure_arity(C);
    1753                 :         14 :   st_alloc(ar);
    1754         [ +  + ]:         28 :   for (i = 1; i < l;   i++) gel(st,sp++) = gel(args,i);
    1755         [ -  + ]:         14 :   for(      ; i <= ar; i++) gel(st,sp++) = NULL;
    1756                 :         14 :   return closure_returnupto(C);
    1757                 :            : }
    1758                 :            : 
    1759                 :            : GEN
    1760                 :          0 : closure_callgenall(GEN C, long n, ...)
    1761                 :            : {
    1762                 :            :   va_list ap;
    1763                 :          0 :   long i, ar = closure_arity(C);
    1764                 :          0 :   va_start(ap,n);
    1765                 :          0 :   st_alloc(ar);
    1766 [ #  # ][ #  # ]:          0 :   for (i = 1; i <=n;  i++) gel(st,sp++) = va_arg(ap, GEN);
                 [ #  # ]
    1767         [ #  # ]:          0 :   for(      ; i <=ar; i++) gel(st,sp++) = NULL;
    1768                 :          0 :   va_end(ap);
    1769                 :          0 :   return closure_returnupto(C);
    1770                 :            : }
    1771                 :            : 
    1772                 :            : GEN
    1773                 :   17215093 : gp_eval(void *E, GEN x)
    1774                 :            : {
    1775                 :   17215093 :   GEN code = (GEN)E;
    1776                 :   17215093 :   set_lex(-1,x);
    1777                 :   17215093 :   return closure_evalnobrk(code);
    1778                 :            : }
    1779                 :            : 
    1780                 :            : GEN
    1781                 :       4445 : gp_evalupto(void *E, GEN x)
    1782                 :            : {
    1783                 :       4445 :   pari_sp av = avma;
    1784                 :       4445 :   return copyupto(gp_eval(E,x), (GEN)av);
    1785                 :            : }
    1786                 :            : 
    1787                 :            : long
    1788                 :     166663 : gp_evalbool(void *E, GEN x)
    1789                 :            : {
    1790                 :     166663 :   pari_sp av = avma;
    1791                 :     166663 :   long res  = !gequal0(gp_eval(E,x));
    1792                 :     166663 :   avma = av; return res;
    1793                 :            : }
    1794                 :            : 
    1795                 :            : long
    1796                 :    3648680 : gp_evalvoid(void *E, GEN x)
    1797                 :            : {
    1798                 :    3648680 :   GEN code = (GEN)E;
    1799                 :    3648680 :   set_lex(-1,x);
    1800                 :    3648680 :   closure_evalvoid(code);
    1801                 :    3648680 :   return loop_break();
    1802                 :            : }
    1803                 :            : 
    1804                 :            : GEN
    1805                 :       1946 : gp_call(void *E, GEN x)
    1806                 :            : {
    1807                 :       1946 :   GEN code = (GEN)E;
    1808                 :       1946 :   return closure_callgen1(code, x);
    1809                 :            : }
    1810                 :            : 
    1811                 :            : long
    1812                 :        504 : gp_callbool(void *E, GEN x)
    1813                 :            : {
    1814                 :        504 :   pari_sp av = avma;
    1815                 :        504 :   GEN code = (GEN)E;
    1816                 :        504 :   long res  = !gequal0(closure_callgen1(code, x));
    1817                 :        504 :   avma = av; return res;
    1818                 :            : }
    1819                 :            : 
    1820                 :            : long
    1821                 :          0 : gp_callvoid(void *E, GEN x)
    1822                 :            : {
    1823                 :          0 :   GEN code = (GEN)E;
    1824                 :          0 :   closure_callvoid1(code, x);
    1825                 :          0 :   return loop_break();
    1826                 :            : }
    1827                 :            : 
    1828                 :            : INLINE const char *
    1829                 :          0 : disassemble_cast(long mode)
    1830                 :            : {
    1831   [ #  #  #  #  :          0 :   switch (mode)
                      # ]
    1832                 :            :   {
    1833                 :            :   case Gsmall:
    1834                 :          0 :     return "small";
    1835                 :            :   case Ggen:
    1836                 :          0 :     return "gen";
    1837                 :            :   case Gvar:
    1838                 :          0 :     return "var";
    1839                 :            :   case Gvoid:
    1840                 :          0 :     return "void";
    1841                 :            :   default:
    1842                 :          0 :     return "unknown";
    1843                 :            :   }
    1844                 :            : }
    1845                 :            : 
    1846                 :            : void
    1847                 :          0 : closure_disassemble(GEN C)
    1848                 :            : {
    1849                 :            :   const char * code;
    1850                 :            :   GEN oper;
    1851                 :            :   long i;
    1852         [ #  # ]:          0 :   if (typ(C)!=t_CLOSURE) pari_err_TYPE("disassemble",C);
    1853                 :          0 :   code=closure_codestr(C);
    1854                 :          0 :   oper=closure_get_oper(C);
    1855         [ #  # ]:          0 :   for(i=1;i<lg(oper);i++)
    1856                 :            :   {
    1857                 :          0 :     op_code opcode=(op_code) code[i];
    1858                 :          0 :     long operand=oper[i];
    1859                 :          0 :     pari_printf("%05ld\t",i);
    1860   [ #  #  #  #  :          0 :     switch(opcode)
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
          #  #  #  #  #  
             #  #  #  #  
                      # ]
    1861                 :            :     {
    1862                 :            :     case OCpushlong:
    1863                 :          0 :       pari_printf("pushlong\t%ld\n",operand);
    1864                 :          0 :       break;
    1865                 :            :     case OCpushgnil:
    1866                 :          0 :       pari_printf("pushgnil\n");
    1867                 :          0 :       break;
    1868                 :            :     case OCpushgen:
    1869                 :          0 :       pari_printf("pushgen\t\t%ld\n",operand);
    1870                 :          0 :       break;
    1871                 :            :     case OCpushreal:
    1872                 :          0 :       pari_printf("pushreal\t%ld\n",operand);
    1873                 :          0 :       break;
    1874                 :            :     case OCpushstoi:
    1875                 :          0 :       pari_printf("pushstoi\t%ld\n",operand);
    1876                 :          0 :       break;
    1877                 :            :     case OCpushvar:
    1878                 :            :       {
    1879                 :          0 :         entree *ep = (entree *)operand;
    1880                 :          0 :         pari_printf("pushvar\t%s\n",ep->name);
    1881                 :          0 :         break;
    1882                 :            :       }
    1883                 :            :     case OCpushdyn:
    1884                 :            :       {
    1885                 :          0 :         entree *ep = (entree *)operand;
    1886                 :          0 :         pari_printf("pushdyn\t\t%s\n",ep->name);
    1887                 :          0 :         break;
    1888                 :            :       }
    1889                 :            :     case OCpushlex:
    1890                 :          0 :       pari_printf("pushlex\t\t%ld\n",operand);
    1891                 :          0 :       break;
    1892                 :            :     case OCstoredyn:
    1893                 :            :       {
    1894                 :          0 :         entree *ep = (entree *)operand;
    1895                 :          0 :         pari_printf("storedyn\t%s\n",ep->name);
    1896                 :          0 :         break;
    1897                 :            :       }
    1898                 :            :     case OCstorelex:
    1899                 :          0 :       pari_printf("storelex\t%ld\n",operand);
    1900                 :          0 :       break;
    1901                 :            :     case OCstoreptr:
    1902                 :          0 :       pari_printf("storeptr\n");
    1903                 :          0 :       break;
    1904                 :            :     case OCsimpleptrdyn:
    1905                 :            :       {
    1906                 :          0 :         entree *ep = (entree *)operand;
    1907                 :          0 :         pari_printf("simpleptrdyn\t%s\n",ep->name);
    1908                 :          0 :         break;
    1909                 :            :       }
    1910                 :            :     case OCsimpleptrlex:
    1911                 :          0 :       pari_printf("simpleptrlex\t%ld\n",operand);
    1912                 :          0 :       break;
    1913                 :            :     case OCnewptrdyn:
    1914                 :            :       {
    1915                 :          0 :         entree *ep = (entree *)operand;
    1916                 :          0 :         pari_printf("newptrdyn\t%s\n",ep->name);
    1917                 :          0 :         break;
    1918                 :            :       }
    1919                 :            :     case OCnewptrlex:
    1920                 :          0 :       pari_printf("newptrlex\t%ld\n",operand);
    1921                 :          0 :       break;
    1922                 :            :     case OCpushptr:
    1923                 :          0 :       pari_printf("pushptr\n");
    1924                 :          0 :       break;
    1925                 :            :     case OCstackgen:
    1926                 :          0 :       pari_printf("stackgen\t%ld\n",operand);
    1927                 :          0 :       break;
    1928                 :            :     case OCendptr:
    1929                 :          0 :       pari_printf("endptr\t\t%ld\n",operand);
    1930                 :          0 :       break;
    1931                 :            :     case OCprecreal:
    1932                 :          0 :       pari_printf("precreal\n");
    1933                 :          0 :       break;
    1934                 :            :     case OCprecdl:
    1935                 :          0 :       pari_printf("precdl\n");
    1936                 :          0 :       break;
    1937                 :            :     case OCstoi:
    1938                 :          0 :       pari_printf("stoi\n");
    1939                 :          0 :       break;
    1940                 :            :     case OCutoi:
    1941                 :          0 :       pari_printf("utoi\n");
    1942                 :          0 :       break;
    1943                 :            :     case OCitos:
    1944                 :          0 :       pari_printf("itos\t\t%ld\n",operand);
    1945                 :          0 :       break;
    1946                 :            :     case OCitou:
    1947                 :          0 :       pari_printf("itou\t\t%ld\n",operand);
    1948                 :          0 :       break;
    1949                 :            :     case OCtostr:
    1950                 :          0 :       pari_printf("tostr\t\t%ld\n",operand);
    1951                 :          0 :       break;
    1952                 :            :     case OCvarn:
    1953                 :          0 :       pari_printf("varn\t\t%ld\n",operand);
    1954                 :          0 :       break;
    1955                 :            :     case OCcopy:
    1956                 :          0 :       pari_printf("copy\n");
    1957                 :          0 :       break;
    1958                 :            :     case OCcopyifclone:
    1959                 :          0 :       pari_printf("copyifclone\n");
    1960                 :          0 :       break;
    1961                 :            :     case OCcompo1:
    1962                 :          0 :       pari_printf("compo1\t\t%s\n",disassemble_cast(operand));
    1963                 :          0 :       break;
    1964                 :            :     case OCcompo1ptr:
    1965                 :          0 :       pari_printf("compo1ptr\n");
    1966                 :          0 :       break;
    1967                 :            :     case OCcompo2:
    1968                 :          0 :       pari_printf("compo2\t\t%s\n",disassemble_cast(operand));
    1969                 :          0 :       break;
    1970                 :            :     case OCcompo2ptr:
    1971                 :          0 :       pari_printf("compo2ptr\n");
    1972                 :          0 :       break;
    1973                 :            :     case OCcompoC:
    1974                 :          0 :       pari_printf("compoC\n");
    1975                 :          0 :       break;
    1976                 :            :     case OCcompoCptr:
    1977                 :          0 :       pari_printf("compoCptr\n");
    1978                 :          0 :       break;
    1979                 :            :     case OCcompoL:
    1980                 :          0 :       pari_printf("compoL\n");
    1981                 :          0 :       break;
    1982                 :            :     case OCcompoLptr:
    1983                 :          0 :       pari_printf("compoLptr\n");
    1984                 :          0 :       break;
    1985                 :            :     case OCcheckargs:
    1986                 :          0 :       pari_printf("checkargs\t0x%lx\n",operand);
    1987                 :          0 :       break;
    1988                 :            :     case OCcheckargs0:
    1989                 :          0 :       pari_printf("checkargs0\t0x%lx\n",operand);
    1990                 :          0 :       break;
    1991                 :            :     case OCcheckuserargs:
    1992                 :          0 :       pari_printf("checkuserargs\t%ld\n",operand);
    1993                 :          0 :       break;
    1994                 :            :     case OCdefaultlong:
    1995                 :          0 :       pari_printf("defaultlong\t%ld\n",operand);
    1996                 :          0 :       break;
    1997                 :            :     case OCdefaultulong:
    1998                 :          0 :       pari_printf("defaultulong\t%ld\n",operand);
    1999                 :          0 :       break;
    2000                 :            :     case OCdefaultgen:
    2001                 :          0 :       pari_printf("defaultgen\t%ld\n",operand);
    2002                 :          0 :       break;
    2003                 :            :     case OCgetargs:
    2004                 :          0 :       pari_printf("getargs\t\t%ld\n",operand);
    2005                 :          0 :       break;
    2006                 :            :     case OCdefaultarg:
    2007                 :          0 :       pari_printf("defaultarg\t%ld\n",operand);
    2008                 :          0 :       break;
    2009                 :            :     case OClocalvar:
    2010                 :            :       {
    2011                 :          0 :         entree *ep = (entree *)operand;
    2012                 :          0 :         pari_printf("localvar\t%s\n",ep->name);
    2013                 :          0 :         break;
    2014                 :            :       }
    2015                 :            :     case OClocalvar0:
    2016                 :            :       {
    2017                 :          0 :         entree *ep = (entree *)operand;
    2018                 :          0 :         pari_printf("localvar0\t%s\n",ep->name);
    2019                 :          0 :         break;
    2020                 :            :       }
    2021                 :            :     case OCcallgen:
    2022                 :            :       {
    2023                 :          0 :         entree *ep = (entree *)operand;
    2024                 :          0 :         pari_printf("callgen\t\t%s\n",ep->name);
    2025                 :          0 :         break;
    2026                 :            :       }
    2027                 :            :     case OCcallgen2:
    2028                 :            :       {
    2029                 :          0 :         entree *ep = (entree *)operand;
    2030                 :          0 :         pari_printf("callgen2\t%s\n",ep->name);
    2031                 :          0 :         break;
    2032                 :            :       }
    2033                 :            :     case OCcalllong:
    2034                 :            :       {
    2035                 :          0 :         entree *ep = (entree *)operand;
    2036                 :          0 :         pari_printf("calllong\t%s\n",ep->name);
    2037                 :          0 :         break;
    2038                 :            :       }
    2039                 :            :     case OCcallint:
    2040                 :            :       {
    2041                 :          0 :         entree *ep = (entree *)operand;
    2042                 :          0 :         pari_printf("callint\t\t%s\n",ep->name);
    2043                 :          0 :         break;
    2044                 :            :       }
    2045                 :            :     case OCcallvoid:
    2046                 :            :       {
    2047                 :          0 :         entree *ep = (entree *)operand;
    2048                 :          0 :         pari_printf("callvoid\t%s\n",ep->name);
    2049                 :          0 :         break;
    2050                 :            :       }
    2051                 :            :     case OCcalluser:
    2052                 :          0 :       pari_printf("calluser\t%ld\n",operand);
    2053                 :          0 :       break;
    2054                 :            :     case OCvec:
    2055                 :          0 :       pari_printf("vec\t\t%ld\n",operand);
    2056                 :          0 :       break;
    2057                 :            :     case OCcol:
    2058                 :          0 :       pari_printf("col\t\t%ld\n",operand);
    2059                 :          0 :       break;
    2060                 :            :     case OCmat:
    2061                 :          0 :       pari_printf("mat\t\t%ld\n",operand);
    2062                 :          0 :       break;
    2063                 :            :     case OCnewframe:
    2064                 :          0 :       pari_printf("newframe\t%ld\n",operand);
    2065                 :          0 :       break;
    2066                 :            :     case OCsaveframe:
    2067                 :          0 :       pari_printf("saveframe\t%ld\n", operand);
    2068                 :          0 :       break;
    2069                 :            :     case OCpop:
    2070                 :          0 :       pari_printf("pop\t\t%ld\n",operand);
    2071                 :          0 :       break;
    2072                 :            :     case OCdup:
    2073                 :          0 :       pari_printf("dup\t\t%ld\n",operand);
    2074                 :          0 :       break;
    2075                 :            :     case OCavma:
    2076                 :          0 :       pari_printf("avma\n",operand);
    2077                 :          0 :       break;
    2078                 :            :     case OCgerepile:
    2079                 :          0 :       pari_printf("gerepile\n",operand);
    2080                 :          0 :       break;
    2081                 :            :     case OCcowvardyn:
    2082                 :            :       {
    2083                 :          0 :         entree *ep = (entree *)operand;
    2084                 :          0 :         pari_printf("cowvardyn\t%s\n",ep->name);
    2085                 :          0 :         break;
    2086                 :            :       }
    2087                 :            :     case OCcowvarlex:
    2088                 :          0 :       pari_printf("cowvarlex\t%ld\n",operand);
    2089                 :          0 :       break;
    2090                 :            :     }
    2091                 :            :   }
    2092                 :          0 : }
    2093                 :            : 
    2094                 :            : static int
    2095                 :          0 : opcode_need_relink(op_code opcode)
    2096                 :            : {
    2097      [ #  #  # ]:          0 :   switch(opcode)
    2098                 :            :   {
    2099                 :            :   case OCpushlong:
    2100                 :            :   case OCpushgen:
    2101                 :            :   case OCpushgnil:
    2102                 :            :   case OCpushreal:
    2103                 :            :   case OCpushstoi:
    2104                 :            :   case OCpushlex:
    2105                 :            :   case OCstorelex:
    2106                 :            :   case OCstoreptr:
    2107                 :            :   case OCsimpleptrlex:
    2108                 :            :   case OCnewptrlex:
    2109                 :            :   case OCpushptr:
    2110                 :            :   case OCstackgen:
    2111                 :            :   case OCendptr:
    2112                 :            :   case OCprecreal:
    2113                 :            :   case OCprecdl:
    2114                 :            :   case OCstoi:
    2115                 :            :   case OCutoi:
    2116                 :            :   case OCitos:
    2117                 :            :   case OCitou:
    2118                 :            :   case OCtostr:
    2119                 :            :   case OCvarn:
    2120                 :            :   case OCcopy:
    2121                 :            :   case OCcopyifclone:
    2122                 :            :   case OCcompo1:
    2123                 :            :   case OCcompo1ptr:
    2124                 :            :   case OCcompo2:
    2125                 :            :   case OCcompo2ptr:
    2126                 :            :   case OCcompoC:
    2127                 :            :   case OCcompoCptr:
    2128                 :            :   case OCcompoL:
    2129                 :            :   case OCcompoLptr:
    2130                 :            :   case OCcheckargs:
    2131                 :            :   case OCcheckargs0:
    2132                 :            :   case OCcheckuserargs:
    2133                 :            :   case OCgetargs:
    2134                 :            :   case OCdefaultarg:
    2135                 :            :   case OCdefaultgen:
    2136                 :            :   case OCdefaultlong:
    2137                 :            :   case OCdefaultulong:
    2138                 :            :   case OCcalluser:
    2139                 :            :   case OCvec:
    2140                 :            :   case OCcol:
    2141                 :            :   case OCmat:
    2142                 :            :   case OCnewframe:
    2143                 :            :   case OCsaveframe:
    2144                 :            :   case OCdup:
    2145                 :            :   case OCpop:
    2146                 :            :   case OCavma:
    2147                 :            :   case OCgerepile:
    2148                 :            :   case OCcowvarlex:
    2149                 :          0 :     break;
    2150                 :            :   case OCpushvar:
    2151                 :            :   case OCpushdyn:
    2152                 :            :   case OCstoredyn:
    2153                 :            :   case OCsimpleptrdyn:
    2154                 :            :   case OCnewptrdyn:
    2155                 :            :   case OClocalvar:
    2156                 :            :   case OClocalvar0:
    2157                 :            :   case OCcallgen:
    2158                 :            :   case OCcallgen2:
    2159                 :            :   case OCcalllong:
    2160                 :            :   case OCcallint:
    2161                 :            :   case OCcallvoid:
    2162                 :            :   case OCcowvardyn:
    2163                 :          0 :     return 1;
    2164                 :            :   }
    2165                 :          0 :   return 0;
    2166                 :            : }
    2167                 :            : 
    2168                 :            : static void
    2169                 :          0 : closure_relink(GEN C, hashtable *table)
    2170                 :            : {
    2171                 :          0 :   const char *code = closure_codestr(C);
    2172                 :          0 :   GEN oper = closure_get_oper(C);
    2173                 :          0 :   GEN fram = gel(closure_get_dbg(C),3);
    2174                 :            :   long i, j;
    2175         [ #  # ]:          0 :   for(i=1;i<lg(oper);i++)
    2176 [ #  # ][ #  # ]:          0 :     if (oper[i] && opcode_need_relink((op_code)code[i]))
    2177                 :          0 :       oper[i] = (long) hash_search(table,(void*) oper[i])->val;
    2178         [ #  # ]:          0 :   for (i=1;i<lg(fram);i++)
    2179         [ #  # ]:          0 :     for (j=1;j<lg(gel(fram,i));j++)
    2180         [ #  # ]:          0 :       if (mael(fram,i,j))
    2181                 :          0 :         mael(fram,i,j) = (long) hash_search(table,(void*) mael(fram,i,j))->val;
    2182                 :          0 : }
    2183                 :            : 
    2184                 :            : void
    2185                 :          0 : gen_relink(GEN x, hashtable *table)
    2186                 :            : {
    2187                 :          0 :   long i, lx, tx = typ(x);
    2188   [ #  #  #  # ]:          0 :   switch(tx)
    2189                 :            :   {
    2190                 :            :     case t_CLOSURE:
    2191                 :          0 :       closure_relink(x, table);
    2192                 :          0 :       gen_relink(closure_get_data(x), table);
    2193         [ #  # ]:          0 :       if (lg(x)==8) gen_relink(closure_get_frame(x), table);
    2194                 :          0 :       break;
    2195                 :            :     case t_LIST:
    2196                 :          0 :       gen_relink(list_data(x), table);
    2197                 :          0 :       break;
    2198                 :            :     case t_VEC: case t_COL: case t_MAT: case t_ERROR:
    2199                 :          0 :       lx = lg(x);
    2200         [ #  # ]:          0 :       for (i=lontyp[tx]; i<lx; i++) gen_relink(gel(x,i), table);
    2201                 :            :   }
    2202                 :          0 : }
    2203                 :            : 
    2204                 :            : static void
    2205                 :          0 : closure_unlink(GEN C)
    2206                 :            : {
    2207                 :          0 :   const char *code = closure_codestr(C);
    2208                 :          0 :   GEN oper = closure_get_oper(C);
    2209                 :          0 :   GEN fram = gel(closure_get_dbg(C),3);
    2210                 :            :   long i, j;
    2211         [ #  # ]:          0 :   for(i=1;i<lg(oper);i++)
    2212 [ #  # ][ #  # ]:          0 :     if (oper[i] && opcode_need_relink((op_code) code[i]))
    2213                 :            :     {
    2214                 :          0 :       long n = pari_stack_new(&s_relocs);
    2215                 :          0 :       relocs[n] = (entree *) oper[i];
    2216                 :            :     }
    2217         [ #  # ]:          0 :   for (i=1;i<lg(fram);i++)
    2218         [ #  # ]:          0 :     for (j=1;j<lg(gel(fram,i));j++)
    2219         [ #  # ]:          0 :       if (mael(fram,i,j))
    2220                 :            :       {
    2221                 :          0 :         long n = pari_stack_new(&s_relocs);
    2222                 :          0 :         relocs[n] = (entree *) mael(fram,i,j);
    2223                 :            :       }
    2224                 :          0 : }
    2225                 :            : 
    2226                 :            : static void
    2227                 :          0 : gen_unlink(GEN x)
    2228                 :            : {
    2229                 :          0 :   long i, lx, tx = typ(x);
    2230   [ #  #  #  # ]:          0 :   switch(tx)
    2231                 :            :   {
    2232                 :            :     case t_CLOSURE:
    2233                 :          0 :       closure_unlink(x);
    2234                 :          0 :       gen_unlink(closure_get_data(x));
    2235         [ #  # ]:          0 :       if (lg(x)==8) gen_unlink(closure_get_frame(x));
    2236                 :          0 :       break;
    2237                 :            :     case t_LIST:
    2238                 :          0 :       gen_unlink(list_data(x));
    2239                 :          0 :       break;
    2240                 :            :     case t_VEC: case t_COL: case t_MAT: case t_ERROR:
    2241                 :          0 :       lx = lg(x);
    2242         [ #  # ]:          0 :       for (i = lontyp[tx]; i<lx; i++) gen_unlink(gel(x,i));
    2243                 :            :   }
    2244                 :          0 : }
    2245                 :            : 
    2246                 :            : GEN
    2247                 :          0 : copybin_unlink(GEN C)
    2248                 :            : {
    2249                 :          0 :   long i, l , n, nold = s_relocs.n;
    2250                 :            :   GEN v, w, V, res;
    2251         [ #  # ]:          0 :   if (C)
    2252                 :          0 :     gen_unlink(C);
    2253                 :            :   else
    2254                 :            :   { /* contents of all variables */
    2255                 :          0 :     long v, maxv = pari_var_next();
    2256         [ #  # ]:          0 :     for (v=0; v<maxv; v++)
    2257                 :            :     {
    2258                 :          0 :       entree *ep = varentries[v];
    2259 [ #  # ][ #  # ]:          0 :       if (!ep || !ep->value) continue;
    2260                 :          0 :       gen_unlink((GEN)ep->value);
    2261                 :            :     }
    2262                 :            :   }
    2263                 :          0 :   n = s_relocs.n-nold;
    2264                 :          0 :   v = cgetg(n+1, t_VECSMALL);
    2265         [ #  # ]:          0 :   for(i=0; i<n; i++)
    2266                 :          0 :     v[i+1] = (long) relocs[i];
    2267                 :          0 :   s_relocs.n = nold;
    2268                 :          0 :   w = vecsmall_uniq(v); l = lg(w);
    2269                 :          0 :   res = cgetg(3,t_VEC);
    2270                 :          0 :   V = cgetg(l, t_VEC);
    2271         [ #  # ]:          0 :   for(i=1; i<l; i++)
    2272                 :            :   {
    2273                 :          0 :     entree *ep = (entree*) w[i];
    2274                 :          0 :     gel(V,i) = strtoGENstr(ep->name);
    2275                 :            :   }
    2276                 :          0 :   gel(res,1) = vecsmall_copy(w);
    2277                 :          0 :   gel(res,2) = V;
    2278                 :          0 :   return res;
    2279                 :            : }
    2280                 :            : 
    2281                 :            : static ulong
    2282                 :          0 : hash_id(void *x) { return (ulong)x; }
    2283                 :            : static int
    2284                 :          0 : eq_id(void *x, void *y) { return x == y; }
    2285                 :            : 
    2286                 :            : /* e = t_VECSMALL of entree *ep [ addresses ],
    2287                 :            :  * names = t_VEC of strtoGENstr(ep.names),
    2288                 :            :  * Return hashtable : ep => is_entry(ep.name) */
    2289                 :            : hashtable *
    2290                 :          0 : hash_from_link(GEN e, GEN names, int use_stack)
    2291                 :            : {
    2292                 :          0 :   long i, l = lg(e);
    2293                 :          0 :   hashtable *h = hash_create(l-1, &hash_id, &eq_id, use_stack);
    2294         [ #  # ]:          0 :   if (lg(names) != l) pari_err_DIM("hash_from_link");
    2295         [ #  # ]:          0 :   for (i = 1; i < l; i++)
    2296                 :            :   {
    2297                 :          0 :     char *s = GSTR(gel(names,i));
    2298                 :          0 :     hash_insert(h, (void*)e[i], (void*)fetch_entry(s, strlen(s)));
    2299                 :            :   }
    2300                 :          0 :   return h;
    2301                 :            : }
    2302                 :            : 
    2303                 :            : void
    2304                 :          0 : bincopy_relink(GEN C, GEN V)
    2305                 :            : {
    2306                 :          0 :   pari_sp av = avma;
    2307                 :          0 :   hashtable *table = hash_from_link(gel(V,1),gel(V,2),1);
    2308                 :          0 :   gen_relink(C, table);
    2309                 :          0 :   avma = av;
    2310                 :          0 : }

Generated by: LCOV version 1.9