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 16746-c2cb716) Lines: 781 1269 61.5 %
Date: 2014-08-31 Functions: 62 93 66.7 %
Legend: Lines: hit not hit | Branches: + taken - not taken # not executed Branches: 375 798 47.0 %

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

Generated by: LCOV version 1.9