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 16375-9f41ae0) Lines: 289 295 98.0 %
Date: 2014-04-19 Functions: 53 55 96.4 %
Legend: Lines: hit not hit | Branches: + taken - not taken # not executed Branches: 198 210 94.3 %

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

Generated by: LCOV version 1.9