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 - members.c (source / functions) Hit Total Coverage
Test: PARI/GP v2.8.0 lcov report (development 16746-c2cb716) Lines: 290 296 98.0 %
Date: 2014-08-31 Functions: 53 55 96.4 %
Legend: Lines: hit not hit | Branches: + taken - not taken # not executed Branches: 201 214 93.9 %

           Branch data     Line data    Source code
       1                 :            : /* Copyright (C) 2000-2003  The PARI group.
       2                 :            : 
       3                 :            : This file is part of the PARI/GP 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                 :            : 
      17                 :            : /********************************************************************/
      18                 :            : /**                                                                **/
      19                 :            : /**                          MEMBER FUNCTIONS                      **/
      20                 :            : /**                                                                **/
      21                 :            : /********************************************************************/
      22                 :            : INLINE int
      23                 :       2000 : is_ell5(GEN x) {
      24                 :            :   long lx;
      25         [ +  + ]:       2000 :   if (typ(x) != t_VEC) return 0;
      26                 :       1750 :   lx = lg(x);
      27 [ +  + ][ +  + ]:       2000 :   return lx == 17 || (lx == 6 && !is_vec_t(typ(gel(x,2))));
                 [ +  + ]
      28                 :            : }
      29                 :       2170 : INLINE int is_ell(GEN x) {
      30                 :       2170 :   long lx = lg(x);
      31 [ +  + ][ +  + ]:       2170 :   return (typ(x) == t_VEC && lx == 17);
      32                 :            : }
      33                 :            : 
      34                 :            : static void
      35                 :       4790 : member_err(const char *s, GEN y) { pari_err_TYPE(s,y); }
      36                 :            : 
      37                 :            : GEN
      38                 :        130 : member_e(GEN x)
      39                 :            : {
      40                 :        130 :   GEN y = get_prid(x);
      41         [ +  + ]:        130 :   if (!y) member_err("e",x);
      42                 :         10 :   return gel(y,3);
      43                 :            : }
      44                 :            : 
      45                 :            : GEN
      46                 :        205 : member_f(GEN x)
      47                 :            : {
      48                 :        205 :   GEN y = get_prid(x);
      49         [ +  + ]:        205 :   if (!y)
      50                 :            :   {
      51         [ +  + ]:        195 :     if (typ(x) == t_FFELT) return utoipos(FF_f(x));
      52                 :        105 :     member_err("f",x);
      53                 :            :   }
      54                 :        100 :   return gel(y,4);
      55                 :            : }
      56                 :            : 
      57                 :            : GEN
      58                 :      35315 : member_p(GEN x)
      59                 :            : {
      60                 :      35315 :   long t; (void)get_nf(x,&t);
      61   [ +  +  +  +  :      35315 :   switch(t)
                      + ]
      62                 :            :   {
      63                 :          5 :     case typ_GAL: return gal_get_p(x);
      64      [ +  +  + ]:        140 :     case typ_ELL: switch(ell_get_type(x))
      65                 :            :     {
      66                 :            :       case t_ELL_Fp:
      67                 :         85 :       case t_ELL_Fq: return ellff_get_p(x);
      68                 :          5 :       case t_ELL_Qp: return ellQp_get_p(x);
      69                 :         50 :       default: member_err("p",x);
      70                 :            :     }
      71                 :          5 :     case typ_MODPR: x = get_prid(x);
      72                 :      35055 :     case typ_PRID: return pr_get_p(x);
      73                 :            :   }
      74      [ -  +  + ]:        115 :   switch(typ(x)) {
      75                 :          0 :     case t_PADIC: return gel(x,2);
      76                 :         35 :     case t_FFELT: return FF_p_i(x);
      77                 :            :   }
      78                 :         80 :   member_err("p",x);
      79                 :      35185 :   return NULL;
      80                 :            : }
      81                 :            : 
      82                 :            : GEN
      83                 :        135 : member_bid(GEN x)
      84                 :            : {
      85                 :        135 :   long t; (void)get_nf(x,&t);
      86      [ +  +  + ]:        135 :   switch(t) {
      87                 :         10 :     case typ_BNR: return bnr_get_bid(x);
      88                 :         10 :     case typ_BID: return x;
      89                 :            :   }
      90                 :        115 :   member_err("bid",x);
      91                 :         20 :   return NULL;
      92                 :            : }
      93                 :            : 
      94                 :            : GEN
      95                 :        130 : member_bnf(GEN x)
      96                 :            : {
      97                 :        130 :   long t; GEN y = get_bnf(x,&t);
      98         [ +  + ]:        130 :   if (!y) member_err("bnf",x);
      99                 :         10 :   return y;
     100                 :            : }
     101                 :            : 
     102                 :            : GEN
     103                 :        140 : member_nf(GEN x)
     104                 :            : {
     105                 :        140 :   long t; GEN y = get_nf(x,&t);
     106         [ +  + ]:        140 :   if (!y) {
     107         [ +  + ]:        120 :     if (t == typ_RNF) return gel(x,10);
     108 [ +  + ][ -  + ]:        105 :     if (t == typ_ELL && ell_get_type(x) == t_ELL_NF) return ellnf_get_nf(x);
     109                 :        105 :     member_err("nf",x);
     110                 :            :   }
     111                 :         35 :   return y;
     112                 :            : }
     113                 :            : 
     114                 :            : /* integral basis */
     115                 :            : GEN
     116                 :        145 : member_zk(GEN x)
     117                 :            : {
     118                 :        145 :   long t; GEN y = get_nf(x,&t);
     119         [ +  + ]:        145 :   if (!y)
     120                 :            :   {
     121      [ +  +  + ]:        110 :     switch(t)
     122                 :            :     {
     123                 :            :       case typ_Q:
     124                 :          5 :         y = cgetg(3,t_VEC);
     125                 :          5 :         gel(y,1) = gen_1;
     126                 :          5 :         gel(y,2) = pol_x(varn(gel(x,1))); return y;
     127                 :            :       case typ_RNF:
     128                 :          5 :         return gel(x,7);
     129                 :            :     }
     130                 :        100 :     member_err("zk",x);
     131                 :            :   }
     132                 :         45 :   return nf_get_zk(y);
     133                 :            : }
     134                 :            : 
     135                 :            : GEN
     136                 :       1685 : member_disc(GEN x) /* discriminant */
     137                 :            : {
     138                 :       1685 :   long t; GEN y = get_nf(x,&t);
     139         [ +  + ]:       1685 :   if (!y)
     140                 :            :   {
     141   [ +  +  +  + ]:       1655 :     switch(t)
     142                 :            :     {
     143                 :          5 :       case typ_Q  : return quad_disc(x);
     144                 :       1565 :       case typ_ELL: return ell_get_disc(x);
     145                 :          5 :       case typ_RNF: return rnf_get_disc(x);
     146                 :            :     }
     147                 :         80 :     member_err("disc",x);
     148                 :            :   }
     149                 :       1605 :   return nf_get_disc(y);
     150                 :            : }
     151                 :            : 
     152                 :            : GEN
     153                 :        245 : member_pol(GEN x) /* polynomial */
     154                 :            : {
     155                 :        245 :   long t; GEN y = get_nf(x,&t);
     156         [ +  + ]:        245 :   if (!y)
     157                 :            :   {
     158   [ +  +  +  +  :        125 :     switch(t)
                      + ]
     159                 :            :     {
     160                 :          5 :       case typ_POL: return x;
     161                 :          5 :       case typ_Q  : return gel(x,1);
     162                 :          5 :       case typ_GAL: return gal_get_pol(x);
     163                 :         20 :       case typ_RNF: return rnf_get_pol(x);
     164                 :            :     }
     165         [ +  + ]:         90 :     if (typ(x)==t_POLMOD) return gel(x,2);
     166         [ +  + ]:         85 :     if (typ(x)==t_FFELT) return FF_to_FpXQ(x);
     167                 :         70 :     member_err("pol",x);
     168                 :            :   }
     169                 :        175 :   return nf_get_pol(y);
     170                 :            : }
     171                 :            : 
     172                 :            : GEN
     173                 :        160 : member_polabs(GEN x)
     174                 :            : {
     175                 :        160 :   long t; (void)get_nf(x,&t);
     176         [ +  + ]:        160 :   if (t != typ_RNF) member_err("pol",x);
     177                 :         35 :   return rnf_get_polabs(x);
     178                 :            : }
     179                 :            : 
     180                 :            : GEN
     181                 :        150 : member_mod(GEN x) /* modulus */
     182                 :            : {
     183                 :        150 :   long t; (void)get_nf(x,&t);
     184   [ +  +  +  + ]:        150 :   switch(t) {
     185                 :          5 :     case typ_GAL: return gal_get_mod(x);
     186                 :          5 :     case typ_BNR: return bnr_get_mod(x);
     187                 :         10 :     case typ_BID: return gel(x,1);
     188                 :            :   }
     189   [ +  -  +  + ]:        130 :   switch(typ(x))
     190                 :            :   {
     191                 :         15 :     case t_INTMOD: case t_POLMOD: case t_QUAD: break;
     192                 :          0 :     case t_PADIC: return gel(x,3);
     193                 :         35 :     case t_FFELT: return FF_mod(x);
     194                 :         80 :     default: member_err("mod",x);
     195                 :            :   }
     196                 :         70 :   return gel(x,1);
     197                 :            : }
     198                 :            : 
     199                 :            : GEN
     200                 :        390 : member_sign(GEN x) /* signature */
     201                 :            : {
     202                 :        390 :   long t; GEN y = get_nf(x,&t);
     203         [ +  + ]:        390 :   if (!y) member_err("sign",x);
     204                 :         60 :   return gel(y,2);
     205                 :            : }
     206                 :            : GEN
     207                 :        130 : member_r1(GEN x) { return gel(member_sign(x), 1); }
     208                 :            : GEN
     209                 :        130 : member_r2(GEN x) { return gel(member_sign(x), 2); }
     210                 :            : 
     211                 :            : GEN
     212                 :        135 : member_index(GEN x)
     213                 :            : {
     214                 :        135 :   long t; GEN y = get_nf(x,&t);
     215         [ +  + ]:        135 :   if (!y)
     216                 :            :   {
     217         [ +  + ]:        110 :     if (t == typ_RNF) return rnf_get_index(x);
     218                 :        105 :     member_err("index",x);
     219                 :            :   }
     220                 :         30 :   return nf_get_index(y);
     221                 :            : }
     222                 :            : 
     223                 :            : /* x assumed to be output by get_nf: ie a t_VEC with length 11 */
     224                 :            : static GEN
     225                 :        390 : nfmats(GEN x)
     226                 :            : {
     227                 :            :   GEN y;
     228         [ +  + ]:        390 :   if (!x) return NULL;
     229                 :         60 :   y = gel(x,5);
     230 [ +  - ][ -  + ]:         60 :   if (typ(y) == t_VEC && lg(y) < 8) return NULL;
     231                 :        390 :   return y;
     232                 :            : }
     233                 :            : 
     234                 :            : GEN
     235                 :        130 : member_t2(GEN x) /* T2 matrix */
     236                 :            : {
     237                 :        130 :   long t; GEN y = nfmats(get_nf(x,&t));
     238         [ +  + ]:        130 :   if (!y) member_err("t2",x);
     239                 :         20 :   return gram_matrix(gel(y,2));
     240                 :            : }
     241                 :            : 
     242                 :            : GEN
     243                 :        130 : member_diff(GEN x) /* different */
     244                 :            : {
     245                 :        130 :   long t; GEN y = nfmats(get_nf(x,&t));
     246         [ +  + ]:        130 :   if (!y) member_err("diff",x);
     247                 :         20 :   return gel(y,5);
     248                 :            : }
     249                 :            : 
     250                 :            : GEN
     251                 :        130 : member_codiff(GEN x) /* codifferent */
     252                 :            : {
     253                 :        130 :   long t; GEN T, D, DinvT, nf = get_nf(x,&t), y = nfmats(nf);
     254         [ +  + ]:        130 :   if (!y) member_err("codiff",x);
     255                 :         20 :   T = gel(y,4);
     256                 :         20 :   D = absi(nf_get_disc(nf));
     257                 :         20 :   DinvT = ZM_inv(T,D);
     258                 :         20 :   return gdiv(ZM_hnfmod(DinvT, D), D);
     259                 :            : }
     260                 :            : 
     261                 :            : GEN
     262                 :        145 : member_roots(GEN x) /* roots */
     263                 :            : {
     264                 :        145 :   long t; GEN y = get_nf(x,&t);
     265         [ +  + ]:        145 :   if (!y)
     266                 :            :   {
     267         [ +  + ]:        125 :     if (t == typ_GAL) return gal_get_roots(x);
     268         [ +  + ]:        120 :     if (t == typ_ELL)
     269      [ +  +  + ]:         35 :       switch(ell_get_type(x))
     270                 :            :       {
     271                 :         10 :         case t_ELL_Qp: return mkcol( ellQp_root(x, ellQp_get_prec(x)) );
     272                 :            :         case t_ELL_Q:
     273                 :         15 :         case t_ELL_Rg: return ellR_roots(x, ellR_get_prec(x));
     274                 :            :       }
     275                 :         95 :     member_err("roots",x);
     276                 :            :   }
     277                 :         50 :   return nf_get_roots(y);
     278                 :            : }
     279                 :            : 
     280                 :            : /* assume x output by get_bnf: ie a t_VEC with length 10 */
     281                 :            : static GEN
     282                 :         70 : check_RES(GEN x, const char *s)
     283                 :            : {
     284                 :         70 :   GEN y = gel(x,8);
     285 [ +  - ][ -  + ]:         70 :   if (typ(y) != t_VEC || lg(y) < 4) member_err(s,x);
     286                 :         70 :   return y;
     287                 :            : }
     288                 :            : 
     289                 :            : /* y = get_bnf(x, &t) */
     290                 :            : static GEN
     291                 :       8250 : _member_clgp(GEN x, GEN y, long t) /* class group (3-component row vector) */
     292                 :            : {
     293         [ +  + ]:       8250 :   if (!y)
     294                 :            :   {
     295      [ +  +  + ]:       8170 :     switch(t)
     296                 :            :     {
     297                 :       3105 :       case typ_QUA: return mkvec3(gel(x,1), gel(x,2), gel(x,3));
     298                 :         40 :       case typ_BID: return gel(x,2);
     299                 :            :     }
     300         [ +  + ]:       5025 :     if (typ(x)==t_VEC)
     301         [ +  + ]:       4830 :       switch(lg(x))
     302                 :            :       {
     303                 :            :         case 3: /* no gen */
     304                 :       4680 :         case 4: return x;
     305                 :            :       }
     306                 :        345 :     member_err("clgp",x);
     307                 :            :   }
     308         [ +  + ]:         80 :   if (t==typ_BNR) return gel(x,5);
     309                 :       7905 :   y = check_RES(y, "clgp"); return gel(y,1);
     310                 :            : }
     311                 :            : static GEN
     312                 :       8250 : _check_clgp(GEN x, GEN y, long t)
     313                 :       8250 : { GEN c = _member_clgp(x,y,t); checkabgrp(c); return c; }
     314                 :            : GEN
     315                 :        150 : member_clgp(GEN x)
     316                 :        150 : { long t; GEN y = get_bnf(x,&t); return _check_clgp(x,y,t); }
     317                 :            : 
     318                 :            : GEN
     319                 :        140 : member_reg(GEN x) /* regulator */
     320                 :            : {
     321                 :        140 :   long t; GEN y = get_bnf(x,&t);
     322         [ +  + ]:        140 :   if (!y)
     323                 :            :   {
     324         [ +  + ]:        120 :     if (t == typ_QUA) return gel(x,4);
     325                 :        115 :     member_err("reg",x);
     326                 :            :   }
     327         [ +  + ]:         20 :   if (t == typ_BNR) pari_err_IMPL("ray regulator");
     328                 :         15 :   y = check_RES(y, "reg");
     329                 :         20 :   return gel(y,2);
     330                 :            : }
     331                 :            : 
     332                 :            : GEN
     333                 :        150 : member_fu(GEN x) /* fundamental units */
     334                 :            : {
     335                 :        150 :   long t; GEN y = get_bnf(x,&t);
     336         [ +  + ]:        150 :   if (!y)
     337                 :            :   {
     338         [ +  + ]:        120 :     switch(t)
     339                 :            :     {
     340                 :            :       case typ_Q:
     341                 :          5 :         x = quad_disc(x);
     342         [ +  - ]:          5 :         return (signe(x)<0)? cgetg(1,t_VEC): quadunit(x);
     343                 :            :     }
     344                 :        115 :     member_err("fu",x);
     345                 :            :   }
     346         [ +  + ]:         30 :   if (t == typ_BNR) pari_err_IMPL("ray units");
     347                 :         30 :   return matbasistoalg(y, bnf_get_fu(y));
     348                 :            : }
     349                 :            : 
     350                 :            : /* torsion units. return [w,e] where w is the number of roots of 1, and e a
     351                 :            :  * polymod generator */
     352                 :            : GEN
     353                 :        130 : member_tu(GEN x)
     354                 :            : {
     355                 :        130 :   long t; GEN bnf = get_bnf(x,&t), res = cgetg(3,t_VEC);
     356         [ +  + ]:        130 :   if (!bnf)
     357                 :            :   {
     358                 :            :     GEN y;
     359         [ +  + ]:        120 :     if (t != typ_Q) member_err("tu",x);
     360                 :          5 :     y = quad_disc(x);
     361 [ +  - ][ -  + ]:          5 :     if (signe(y) > 0 || cmpiu(y,4) > 0) return mkvec2(gen_m1, gen_2);
     362                 :            : 
     363         [ +  - ]:          5 :     gel(res,1) = utoipos((itos(y) == -4)? 4: 6);
     364                 :          5 :     gel(res,2) = gcopy(x);
     365                 :            :   }
     366                 :            :   else
     367                 :            :   {
     368         [ +  + ]:         10 :     if (t == typ_BNR) pari_err_IMPL("ray torsion units");
     369                 :          5 :     gel(res,1) = utoipos( bnf_get_tuN(bnf) );
     370                 :          5 :     gel(res,2) = basistoalg(bnf, bnf_get_tuU(bnf));
     371                 :            :   }
     372                 :         10 :   return res;
     373                 :            : }
     374                 :            : 
     375                 :            : GEN
     376                 :          0 : member_futu(GEN x) /*  concatenation of fu and tu, w is lost */
     377                 :            : {
     378                 :          0 :   return shallowconcat(member_fu(x), gel(member_tu(x),2));
     379                 :            : }
     380                 :            : 
     381                 :            : GEN
     382                 :          0 : member_tufu(GEN x) /*  concatenation of tu and fu, w is lost */
     383                 :            : {
     384                 :          0 :   return shallowconcat(gel(member_tu(x),2), member_fu(x));
     385                 :            : }
     386                 :            : 
     387                 :            : /* structure of (Z_K/m)^*, where x is an idealstarinit (with or without gen)
     388                 :            :  * or a bnrinit (with or without gen) */
     389                 :            : GEN
     390                 :        130 : member_zkst(GEN x)
     391                 :            : {
     392                 :        130 :   long t; (void)get_nf(x,&t);
     393      [ +  +  + ]:        130 :   switch(t)
     394                 :            :   {
     395                 :         10 :     case typ_BID: return gel(x,2);
     396                 :            :     case typ_BNR: {
     397                 :          5 :       GEN bid = bnr_get_bid(x);
     398 [ +  - ][ +  - ]:          5 :       if (typ(bid) == t_VEC && lg(bid) > 2) return gel(bid,2);
     399                 :            :     }
     400                 :            :   }
     401                 :        115 :   member_err("zkst",x);
     402                 :         15 :   return NULL; /* not reached */
     403                 :            : }
     404                 :            : 
     405                 :            : GEN
     406                 :       2120 : member_no(GEN x) /* number of elements of a group (of type clgp) */
     407                 :            : {
     408                 :       2120 :   pari_sp av = avma;
     409                 :       2120 :   long t; GEN y = get_bnf(x,&t);
     410 [ +  + ][ +  + ]:       2120 :   if (t == typ_ELL) switch(ell_get_type(x))
     411                 :            :   {
     412                 :            :     case t_ELL_Fp:
     413                 :         10 :     case t_ELL_Fq: return ellcard(x, NULL);
     414                 :            :   }
     415                 :       2110 :   x = _check_clgp(x,y,t);
     416                 :       2025 :   avma = av; return gel(x,1);
     417                 :            : }
     418                 :            : 
     419                 :            : GEN
     420                 :       4430 : member_cyc(GEN x) /* cyclic decomposition (SNF) of a group (of type clgp) */
     421                 :            : {
     422                 :       4430 :   pari_sp av = avma;
     423                 :       4430 :   long t; GEN y = get_bnf(x,&t);
     424 [ +  + ][ +  + ]:       4430 :   if (t == typ_ELL) switch(ell_get_type(x))
     425                 :            :   {
     426                 :            :     case t_ELL_Fp:
     427                 :         10 :     case t_ELL_Fq: return ellgroup(x, NULL);
     428                 :            :   }
     429                 :       4420 :   x = _check_clgp(x,y,t);
     430                 :       4335 :   avma = av; return gel(x,2);
     431                 :            : }
     432                 :            : 
     433                 :            : /* SNF generators of a group (of type clgp), or generators of a prime
     434                 :            :  * ideal */
     435                 :            : GEN
     436                 :       1650 : member_gen(GEN x)
     437                 :            : {
     438                 :            :   pari_sp av;
     439                 :       1650 :   long t; GEN y = get_bnf(x,&t);
     440   [ +  +  +  +  :       1650 :   switch(t)
                      + ]
     441                 :            :   {
     442                 :          5 :     case typ_MODPR: x = get_prid(x);
     443                 :         25 :     case typ_PRID: return mkvec2(gel(x,1), gel(x,2));
     444                 :          5 :     case typ_GAL: return gal_get_gen(x);
     445                 :         50 :     case typ_ELL: return ellgenerators(x);
     446                 :            :   }
     447                 :       1570 :   av = avma;
     448                 :       1570 :   x = _check_clgp(x,y,t);
     449         [ +  + ]:       1500 :   if (lg(x)!=4) member_err("gen",x);
     450                 :       1565 :   avma = av; return gel(x,3);
     451                 :            : }
     452                 :            : GEN
     453                 :        240 : member_group(GEN x)
     454                 :            : {
     455                 :        240 :   long t; (void)get_nf(x,&t);
     456         [ +  + ]:        240 :   if (t == typ_GAL) return gal_get_group(x);
     457         [ +  + ]:        160 :   if (t == typ_ELL) return ellgroup0(x, NULL, 1);
     458                 :        105 :   member_err("group",x);
     459                 :        125 :   return NULL; /* not reached */
     460                 :            : }
     461                 :            : GEN
     462                 :        130 : member_orders(GEN x)
     463                 :            : {
     464                 :        130 :   long t; (void)get_nf(x,&t);
     465         [ +  + ]:        130 :   if (t == typ_GAL) return gal_get_orders(x);
     466                 :        125 :   member_err("orders",x);
     467                 :          5 :   return NULL; /* not reached */
     468                 :            : }
     469                 :            : 
     470                 :            : GEN
     471                 :        400 : member_a1(GEN x)
     472                 :            : {
     473         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a1",x);
     474                 :        295 :   return ell_get_a1(x);
     475                 :            : }
     476                 :            : 
     477                 :            : GEN
     478                 :        400 : member_a2(GEN x)
     479                 :            : {
     480         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a2",x);
     481                 :        295 :   return ell_get_a2(x);
     482                 :            : }
     483                 :            : 
     484                 :            : GEN
     485                 :        400 : member_a3(GEN x)
     486                 :            : {
     487         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a3",x);
     488                 :        295 :   return ell_get_a3(x);
     489                 :            : }
     490                 :            : 
     491                 :            : GEN
     492                 :        400 : member_a4(GEN x)
     493                 :            : {
     494         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a4",x);
     495                 :        295 :   return ell_get_a4(x);
     496                 :            : }
     497                 :            : 
     498                 :            : GEN
     499                 :        400 : member_a6(GEN x)
     500                 :            : {
     501         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a6",x);
     502                 :        295 :   return ell_get_a6(x);
     503                 :            : }
     504                 :            : 
     505                 :            : GEN
     506                 :        130 : member_b2(GEN x)
     507                 :            : {
     508         [ +  + ]:        130 :   if (!is_ell(x)) member_err("b2",x);
     509                 :         20 :   return ell_get_b2(x);
     510                 :            : }
     511                 :            : 
     512                 :            : GEN
     513                 :        130 : member_b4(GEN x)
     514                 :            : {
     515         [ +  + ]:        130 :   if (!is_ell(x)) member_err("b4",x);
     516                 :         20 :   return ell_get_b4(x);
     517                 :            : }
     518                 :            : 
     519                 :            : GEN
     520                 :        130 : member_b6(GEN x)
     521                 :            : {
     522         [ +  + ]:        130 :   if (!is_ell(x)) member_err("b6",x);
     523                 :         20 :   return ell_get_b6(x);
     524                 :            : }
     525                 :            : 
     526                 :            : GEN
     527                 :        130 : member_b8(GEN x)
     528                 :            : {
     529         [ +  + ]:        130 :   if (!is_ell(x)) member_err("b8",x);
     530                 :         20 :   return ell_get_b8(x);
     531                 :            : }
     532                 :            : 
     533                 :            : GEN
     534                 :        135 : member_c4(GEN x)
     535                 :            : {
     536         [ +  + ]:        135 :   if (!is_ell(x)) member_err("c4",x);
     537                 :         25 :   return ell_get_c4(x);
     538                 :            : }
     539                 :            : 
     540                 :            : GEN
     541                 :        135 : member_c6(GEN x)
     542                 :            : {
     543         [ +  + ]:        135 :   if (!is_ell(x)) member_err("c6",x);
     544                 :         25 :   return ell_get_c6(x);
     545                 :            : }
     546                 :            : 
     547                 :            : GEN
     548                 :        775 : member_j(GEN x)
     549                 :            : {
     550         [ +  + ]:        775 :   if (!is_ell(x)) member_err("j",x);
     551                 :        665 :   return ell_get_j(x);
     552                 :            : }
     553                 :            : 
     554                 :            : static int
     555                 :        125 : ell_is_complex(GEN x)
     556 [ +  + ][ +  + ]:        125 : { long t = ell_get_type(x); return t == t_ELL_Q || t == t_ELL_Rg; }
     557                 :            : 
     558                 :            : GEN
     559                 :        165 : member_omega(GEN x)
     560                 :            : {
     561         [ +  + ]:        165 :   if (!is_ell(x)) member_err("omega",x);
     562         [ +  + ]:         55 :   if (!ell_is_complex(x)) pari_err_TYPE("omega [not defined over C]",x);
     563                 :         30 :   return ellR_omega(x, ellR_get_prec(x));
     564                 :            : }
     565                 :            : 
     566                 :            : GEN
     567                 :        155 : member_eta(GEN x)
     568                 :            : {
     569         [ +  + ]:        155 :   if (!is_ell(x)) member_err("eta",x);
     570         [ +  + ]:         45 :   if (!ell_is_complex(x)) pari_err_TYPE("eta [not defined over C]",x);
     571                 :         25 :   return ellR_eta(x, ellR_get_prec(x));
     572                 :            : }
     573                 :            : 
     574                 :            : GEN
     575                 :        135 : member_area(GEN x)
     576                 :            : {
     577                 :            :   GEN w, w1, w2, a,b,c,d;
     578                 :            :   long prec;
     579         [ +  + ]:        135 :   if (!is_ell(x)) member_err("area",x);
     580         [ +  + ]:         25 :   if (!ell_is_complex(x)) pari_err_TYPE("area [not defined over C]",x);
     581                 :          5 :   prec = ellR_get_prec(x);
     582                 :          5 :   w = ellR_omega(x, prec);
     583                 :          5 :   w1 = gel(w,1); a = real_i(w1); b = imag_i(w1);
     584                 :          5 :   w2 = gel(w,2); c = real_i(w2); d = imag_i(w2);
     585                 :          5 :   return gabs(gsub(gmul(a,d),gmul(b,c)), prec);
     586                 :            : }
     587                 :            : 
     588                 :            : GEN
     589                 :        150 : member_tate(GEN x)
     590                 :            : {
     591                 :            :   long prec;
     592         [ +  + ]:        150 :   if (!is_ell(x)) member_err("tate",x);
     593         [ +  + ]:         40 :   if (ell_get_type(x) != t_ELL_Qp)
     594                 :         20 :     pari_err_TYPE("tate [not defined over Qp]",x);
     595                 :         20 :   prec = ellQp_get_prec(x);
     596                 :         20 :   return ellQp_Tate_uniformization(x, prec);
     597                 :            : }

Generated by: LCOV version 1.9