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 16393-29b9383) Lines: 289 295 98.0 %
Date: 2014-04-24 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                 :       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                 :       2125 : INLINE int is_ell(GEN x) {
      30                 :       2125 :   long lx = lg(x);
      31 [ +  + ][ +  + ]:       2125 :   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                 :        330 : member_p(GEN x)
      59                 :            : {
      60                 :        330 :   long t; (void)get_nf(x,&t);
      61   [ +  +  +  +  :        330 :   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                 :         70 :     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                 :        200 :   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 :     member_err("nf",x);
     109                 :            :   }
     110                 :         35 :   return y;
     111                 :            : }
     112                 :            : 
     113                 :            : /* integral basis */
     114                 :            : GEN
     115                 :        145 : member_zk(GEN x)
     116                 :            : {
     117                 :        145 :   long t; GEN y = get_nf(x,&t);
     118         [ +  + ]:        145 :   if (!y)
     119                 :            :   {
     120      [ +  +  + ]:        110 :     switch(t)
     121                 :            :     {
     122                 :            :       case typ_Q:
     123                 :          5 :         y = cgetg(3,t_VEC);
     124                 :          5 :         gel(y,1) = gen_1;
     125                 :          5 :         gel(y,2) = pol_x(varn(gel(x,1))); return y;
     126                 :            :       case typ_RNF:
     127                 :          5 :         return gel(x,7);
     128                 :            :     }
     129                 :        100 :     member_err("zk",x);
     130                 :            :   }
     131                 :         45 :   return nf_get_zk(y);
     132                 :            : }
     133                 :            : 
     134                 :            : GEN
     135                 :       1560 : member_disc(GEN x) /* discriminant */
     136                 :            : {
     137                 :       1560 :   long t; GEN y = get_nf(x,&t);
     138         [ +  + ]:       1560 :   if (!y)
     139                 :            :   {
     140   [ +  +  +  + ]:       1530 :     switch(t)
     141                 :            :     {
     142                 :          5 :       case typ_Q  : return quad_disc(x);
     143                 :       1440 :       case typ_ELL: return ell_get_disc(x);
     144                 :          5 :       case typ_RNF: return rnf_get_disc(x);
     145                 :            :     }
     146                 :         80 :     member_err("disc",x);
     147                 :            :   }
     148                 :       1480 :   return nf_get_disc(y);
     149                 :            : }
     150                 :            : 
     151                 :            : GEN
     152                 :        245 : member_pol(GEN x) /* polynomial */
     153                 :            : {
     154                 :        245 :   long t; GEN y = get_nf(x,&t);
     155         [ +  + ]:        245 :   if (!y)
     156                 :            :   {
     157   [ +  +  +  +  :        125 :     switch(t)
                      + ]
     158                 :            :     {
     159                 :          5 :       case typ_POL: return x;
     160                 :          5 :       case typ_Q  : return gel(x,1);
     161                 :          5 :       case typ_GAL: return gal_get_pol(x);
     162                 :         20 :       case typ_RNF: return rnf_get_pol(x);
     163                 :            :     }
     164         [ +  + ]:         90 :     if (typ(x)==t_POLMOD) return gel(x,2);
     165         [ +  + ]:         85 :     if (typ(x)==t_FFELT) return FF_to_FpXQ(x);
     166                 :         70 :     member_err("pol",x);
     167                 :            :   }
     168                 :        175 :   return nf_get_pol(y);
     169                 :            : }
     170                 :            : 
     171                 :            : GEN
     172                 :        160 : member_polabs(GEN x)
     173                 :            : {
     174                 :        160 :   long t; (void)get_nf(x,&t);
     175         [ +  + ]:        160 :   if (t != typ_RNF) member_err("pol",x);
     176                 :         35 :   return rnf_get_polabs(x);
     177                 :            : }
     178                 :            : 
     179                 :            : GEN
     180                 :        150 : member_mod(GEN x) /* modulus */
     181                 :            : {
     182                 :        150 :   long t; (void)get_nf(x,&t);
     183   [ +  +  +  + ]:        150 :   switch(t) {
     184                 :          5 :     case typ_GAL: return gal_get_mod(x);
     185                 :          5 :     case typ_BNR: return bnr_get_mod(x);
     186                 :         10 :     case typ_BID: return gel(x,1);
     187                 :            :   }
     188   [ +  -  +  + ]:        130 :   switch(typ(x))
     189                 :            :   {
     190                 :         15 :     case t_INTMOD: case t_POLMOD: case t_QUAD: break;
     191                 :          0 :     case t_PADIC: return gel(x,3);
     192                 :         35 :     case t_FFELT: return FF_mod(x);
     193                 :         80 :     default: member_err("mod",x);
     194                 :            :   }
     195                 :         70 :   return gel(x,1);
     196                 :            : }
     197                 :            : 
     198                 :            : GEN
     199                 :        390 : member_sign(GEN x) /* signature */
     200                 :            : {
     201                 :        390 :   long t; GEN y = get_nf(x,&t);
     202         [ +  + ]:        390 :   if (!y) member_err("sign",x);
     203                 :         60 :   return gel(y,2);
     204                 :            : }
     205                 :            : GEN
     206                 :        130 : member_r1(GEN x) { return gel(member_sign(x), 1); }
     207                 :            : GEN
     208                 :        130 : member_r2(GEN x) { return gel(member_sign(x), 2); }
     209                 :            : 
     210                 :            : GEN
     211                 :        135 : member_index(GEN x)
     212                 :            : {
     213                 :        135 :   long t; GEN y = get_nf(x,&t);
     214         [ +  + ]:        135 :   if (!y)
     215                 :            :   {
     216         [ +  + ]:        110 :     if (t == typ_RNF) return rnf_get_index(x);
     217                 :        105 :     member_err("index",x);
     218                 :            :   }
     219                 :         30 :   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                 :        390 : nfmats(GEN x)
     225                 :            : {
     226                 :            :   GEN y;
     227         [ +  + ]:        390 :   if (!x) return NULL;
     228                 :         60 :   y = gel(x,5);
     229 [ +  - ][ -  + ]:         60 :   if (typ(y) == t_VEC && lg(y) < 8) return NULL;
     230                 :        390 :   return y;
     231                 :            : }
     232                 :            : 
     233                 :            : GEN
     234                 :        130 : member_t2(GEN x) /* T2 matrix */
     235                 :            : {
     236                 :        130 :   long t; GEN y = nfmats(get_nf(x,&t));
     237         [ +  + ]:        130 :   if (!y) member_err("t2",x);
     238                 :         20 :   return gram_matrix(gel(y,2));
     239                 :            : }
     240                 :            : 
     241                 :            : GEN
     242                 :        130 : member_diff(GEN x) /* different */
     243                 :            : {
     244                 :        130 :   long t; GEN y = nfmats(get_nf(x,&t));
     245         [ +  + ]:        130 :   if (!y) member_err("diff",x);
     246                 :         20 :   return gel(y,5);
     247                 :            : }
     248                 :            : 
     249                 :            : GEN
     250                 :        130 : member_codiff(GEN x) /* codifferent */
     251                 :            : {
     252                 :        130 :   long t; GEN T, D, DinvT, nf = get_nf(x,&t), y = nfmats(nf);
     253         [ +  + ]:        130 :   if (!y) member_err("codiff",x);
     254                 :         20 :   T = gel(y,4);
     255                 :         20 :   D = absi(nf_get_disc(nf));
     256                 :         20 :   DinvT = ZM_inv(T,D);
     257                 :         20 :   return gdiv(ZM_hnfmod(DinvT, D), D);
     258                 :            : }
     259                 :            : 
     260                 :            : GEN
     261                 :        145 : member_roots(GEN x) /* roots */
     262                 :            : {
     263                 :        145 :   long t; GEN y = get_nf(x,&t);
     264         [ +  + ]:        145 :   if (!y)
     265                 :            :   {
     266         [ +  + ]:        125 :     if (t == typ_GAL) return gal_get_roots(x);
     267         [ +  + ]:        120 :     if (t == typ_ELL)
     268      [ +  +  + ]:         35 :       switch(ell_get_type(x))
     269                 :            :       {
     270                 :         10 :         case t_ELL_Qp: return mkcol( ellQp_root(x, ellQp_get_prec(x)) );
     271                 :            :         case t_ELL_Q:
     272                 :         15 :         case t_ELL_Rg: return ellR_roots(x, ellR_get_prec(x));
     273                 :            :       }
     274                 :         95 :     member_err("roots",x);
     275                 :            :   }
     276                 :         50 :   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                 :         70 : check_RES(GEN x, const char *s)
     282                 :            : {
     283                 :         70 :   GEN y = gel(x,8);
     284 [ +  - ][ -  + ]:         70 :   if (typ(y) != t_VEC || lg(y) < 4) member_err(s,x);
     285                 :         70 :   return y;
     286                 :            : }
     287                 :            : 
     288                 :            : /* y = get_bnf(x, &t) */
     289                 :            : static GEN
     290                 :       8250 : _member_clgp(GEN x, GEN y, long t) /* class group (3-component row vector) */
     291                 :            : {
     292         [ +  + ]:       8250 :   if (!y)
     293                 :            :   {
     294      [ +  +  + ]:       8170 :     switch(t)
     295                 :            :     {
     296                 :       3105 :       case typ_QUA: return mkvec3(gel(x,1), gel(x,2), gel(x,3));
     297                 :         40 :       case typ_BID: return gel(x,2);
     298                 :            :     }
     299         [ +  + ]:       5025 :     if (typ(x)==t_VEC)
     300         [ +  + ]:       4830 :       switch(lg(x))
     301                 :            :       {
     302                 :            :         case 3: /* no gen */
     303                 :       4680 :         case 4: return x;
     304                 :            :       }
     305                 :        345 :     member_err("clgp",x);
     306                 :            :   }
     307         [ +  + ]:         80 :   if (t==typ_BNR) return gel(x,5);
     308                 :       7905 :   y = check_RES(y, "clgp"); return gel(y,1);
     309                 :            : }
     310                 :            : static GEN
     311                 :       8250 : _check_clgp(GEN x, GEN y, long t)
     312                 :       8250 : { GEN c = _member_clgp(x,y,t); checkabgrp(c); return c; }
     313                 :            : GEN
     314                 :        150 : member_clgp(GEN x)
     315                 :        150 : { long t; GEN y = get_bnf(x,&t); return _check_clgp(x,y,t); }
     316                 :            : 
     317                 :            : GEN
     318                 :        140 : member_reg(GEN x) /* regulator */
     319                 :            : {
     320                 :        140 :   long t; GEN y = get_bnf(x,&t);
     321         [ +  + ]:        140 :   if (!y)
     322                 :            :   {
     323         [ +  + ]:        120 :     if (t == typ_QUA) return gel(x,4);
     324                 :        115 :     member_err("reg",x);
     325                 :            :   }
     326         [ +  + ]:         20 :   if (t == typ_BNR) pari_err_IMPL("ray regulator");
     327                 :         15 :   y = check_RES(y, "reg");
     328                 :         20 :   return gel(y,2);
     329                 :            : }
     330                 :            : 
     331                 :            : GEN
     332                 :        145 : member_fu(GEN x) /* fundamental units */
     333                 :            : {
     334                 :        145 :   long t; GEN y = get_bnf(x,&t);
     335         [ +  + ]:        145 :   if (!y)
     336                 :            :   {
     337         [ +  + ]:        120 :     switch(t)
     338                 :            :     {
     339                 :            :       case typ_Q:
     340                 :          5 :         x = quad_disc(x);
     341         [ +  - ]:          5 :         return (signe(x)<0)? cgetg(1,t_VEC): quadunit(x);
     342                 :            :     }
     343                 :        115 :     member_err("fu",x);
     344                 :            :   }
     345         [ +  + ]:         25 :   if (t == typ_BNR) pari_err_IMPL("ray units");
     346                 :         25 :   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                 :        130 : member_tu(GEN x)
     353                 :            : {
     354                 :        130 :   long t; GEN bnf = get_bnf(x,&t), res = cgetg(3,t_VEC);
     355         [ +  + ]:        130 :   if (!bnf)
     356                 :            :   {
     357                 :            :     GEN y;
     358         [ +  + ]:        120 :     if (t != typ_Q) member_err("tu",x);
     359                 :          5 :     y = quad_disc(x);
     360 [ +  - ][ -  + ]:          5 :     if (signe(y) > 0 || cmpiu(y,4) > 0) return mkvec2(gen_m1, gen_2);
     361                 :            : 
     362         [ +  - ]:          5 :     gel(res,1) = utoipos((itos(y) == -4)? 4: 6);
     363                 :          5 :     gel(res,2) = gcopy(x);
     364                 :            :   }
     365                 :            :   else
     366                 :            :   {
     367         [ +  + ]:         10 :     if (t == typ_BNR) pari_err_IMPL("ray torsion units");
     368                 :          5 :     gel(res,1) = utoipos( bnf_get_tuN(bnf) );
     369                 :          5 :     gel(res,2) = basistoalg(bnf, bnf_get_tuU(bnf));
     370                 :            :   }
     371                 :         10 :   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                 :        130 : member_zkst(GEN x)
     390                 :            : {
     391                 :        130 :   long t; (void)get_nf(x,&t);
     392      [ +  +  + ]:        130 :   switch(t)
     393                 :            :   {
     394                 :         10 :     case typ_BID: return gel(x,2);
     395                 :            :     case typ_BNR: {
     396                 :          5 :       GEN bid = bnr_get_bid(x);
     397 [ +  - ][ +  - ]:          5 :       if (typ(bid) == t_VEC && lg(bid) > 2) return gel(bid,2);
     398                 :            :     }
     399                 :            :   }
     400                 :        115 :   member_err("zkst",x);
     401                 :         15 :   return NULL; /* not reached */
     402                 :            : }
     403                 :            : 
     404                 :            : GEN
     405                 :       2120 : member_no(GEN x) /* number of elements of a group (of type clgp) */
     406                 :            : {
     407                 :       2120 :   pari_sp av = avma;
     408                 :       2120 :   long t; GEN y = get_bnf(x,&t);
     409 [ +  + ][ +  + ]:       2120 :   if (t == typ_ELL) switch(ell_get_type(x))
     410                 :            :   {
     411                 :            :     case t_ELL_Fp:
     412                 :         10 :     case t_ELL_Fq: return ellcard(x, NULL);
     413                 :            :   }
     414                 :       2110 :   x = _check_clgp(x,y,t);
     415                 :       2025 :   avma = av; return gel(x,1);
     416                 :            : }
     417                 :            : 
     418                 :            : GEN
     419                 :       4430 : member_cyc(GEN x) /* cyclic decomposition (SNF) of a group (of type clgp) */
     420                 :            : {
     421                 :       4430 :   pari_sp av = avma;
     422                 :       4430 :   long t; GEN y = get_bnf(x,&t);
     423 [ +  + ][ +  + ]:       4430 :   if (t == typ_ELL) switch(ell_get_type(x))
     424                 :            :   {
     425                 :            :     case t_ELL_Fp:
     426                 :         10 :     case t_ELL_Fq: return ellgroup(x, NULL);
     427                 :            :   }
     428                 :       4420 :   x = _check_clgp(x,y,t);
     429                 :       4335 :   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                 :       1650 : member_gen(GEN x)
     436                 :            : {
     437                 :            :   pari_sp av;
     438                 :       1650 :   long t; GEN y = get_bnf(x,&t);
     439   [ +  +  +  +  :       1650 :   switch(t)
                      + ]
     440                 :            :   {
     441                 :          5 :     case typ_MODPR: x = get_prid(x);
     442                 :         25 :     case typ_PRID: return mkvec2(gel(x,1), gel(x,2));
     443                 :          5 :     case typ_GAL: return gal_get_gen(x);
     444                 :         50 :     case typ_ELL: return ellgenerators(x);
     445                 :            :   }
     446                 :       1570 :   av = avma;
     447                 :       1570 :   x = _check_clgp(x,y,t);
     448         [ +  + ]:       1500 :   if (lg(x)!=4) member_err("gen",x);
     449                 :       1565 :   avma = av; return gel(x,3);
     450                 :            : }
     451                 :            : GEN
     452                 :        240 : member_group(GEN x)
     453                 :            : {
     454                 :        240 :   long t; (void)get_nf(x,&t);
     455         [ +  + ]:        240 :   if (t == typ_GAL) return gal_get_group(x);
     456         [ +  + ]:        160 :   if (t == typ_ELL) return ellgroup0(x, NULL, 1);
     457                 :        105 :   member_err("group",x);
     458                 :        125 :   return NULL; /* not reached */
     459                 :            : }
     460                 :            : GEN
     461                 :        130 : member_orders(GEN x)
     462                 :            : {
     463                 :        130 :   long t; (void)get_nf(x,&t);
     464         [ +  + ]:        130 :   if (t == typ_GAL) return gal_get_orders(x);
     465                 :        125 :   member_err("orders",x);
     466                 :          5 :   return NULL; /* not reached */
     467                 :            : }
     468                 :            : 
     469                 :            : GEN
     470                 :        400 : member_a1(GEN x)
     471                 :            : {
     472         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a1",x);
     473                 :        295 :   return ell_get_a1(x);
     474                 :            : }
     475                 :            : 
     476                 :            : GEN
     477                 :        400 : member_a2(GEN x)
     478                 :            : {
     479         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a2",x);
     480                 :        295 :   return ell_get_a2(x);
     481                 :            : }
     482                 :            : 
     483                 :            : GEN
     484                 :        400 : member_a3(GEN x)
     485                 :            : {
     486         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a3",x);
     487                 :        295 :   return ell_get_a3(x);
     488                 :            : }
     489                 :            : 
     490                 :            : GEN
     491                 :        400 : member_a4(GEN x)
     492                 :            : {
     493         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a4",x);
     494                 :        295 :   return ell_get_a4(x);
     495                 :            : }
     496                 :            : 
     497                 :            : GEN
     498                 :        400 : member_a6(GEN x)
     499                 :            : {
     500         [ +  + ]:        400 :   if (!is_ell5(x)) member_err("a6",x);
     501                 :        295 :   return ell_get_a6(x);
     502                 :            : }
     503                 :            : 
     504                 :            : GEN
     505                 :        130 : member_b2(GEN x)
     506                 :            : {
     507         [ +  + ]:        130 :   if (!is_ell(x)) member_err("b2",x);
     508                 :         20 :   return ell_get_b2(x);
     509                 :            : }
     510                 :            : 
     511                 :            : GEN
     512                 :        130 : member_b4(GEN x)
     513                 :            : {
     514         [ +  + ]:        130 :   if (!is_ell(x)) member_err("b4",x);
     515                 :         20 :   return ell_get_b4(x);
     516                 :            : }
     517                 :            : 
     518                 :            : GEN
     519                 :        130 : member_b6(GEN x)
     520                 :            : {
     521         [ +  + ]:        130 :   if (!is_ell(x)) member_err("b6",x);
     522                 :         20 :   return ell_get_b6(x);
     523                 :            : }
     524                 :            : 
     525                 :            : GEN
     526                 :        130 : member_b8(GEN x)
     527                 :            : {
     528         [ +  + ]:        130 :   if (!is_ell(x)) member_err("b8",x);
     529                 :         20 :   return ell_get_b8(x);
     530                 :            : }
     531                 :            : 
     532                 :            : GEN
     533                 :        130 : member_c4(GEN x)
     534                 :            : {
     535         [ +  + ]:        130 :   if (!is_ell(x)) member_err("c4",x);
     536                 :         20 :   return ell_get_c4(x);
     537                 :            : }
     538                 :            : 
     539                 :            : GEN
     540                 :        130 : member_c6(GEN x)
     541                 :            : {
     542         [ +  + ]:        130 :   if (!is_ell(x)) member_err("c6",x);
     543                 :         20 :   return ell_get_c6(x);
     544                 :            : }
     545                 :            : 
     546                 :            : GEN
     547                 :        760 : member_j(GEN x)
     548                 :            : {
     549         [ +  + ]:        760 :   if (!is_ell(x)) member_err("j",x);
     550                 :        650 :   return ell_get_j(x);
     551                 :            : }
     552                 :            : 
     553                 :            : static int
     554                 :        105 : ell_is_complex(GEN x)
     555 [ +  + ][ +  + ]:        105 : { long t = ell_get_type(x); return t == t_ELL_Q || t == t_ELL_Rg; }
     556                 :            : 
     557                 :            : GEN
     558                 :        145 : member_omega(GEN x)
     559                 :            : {
     560         [ +  + ]:        145 :   if (!is_ell(x)) member_err("omega",x);
     561         [ +  + ]:         35 :   if (!ell_is_complex(x)) pari_err_TYPE("omega [not defined over C]",x);
     562                 :         10 :   return ellR_omega(x, ellR_get_prec(x));
     563                 :            : }
     564                 :            : 
     565                 :            : GEN
     566                 :        155 : member_eta(GEN x)
     567                 :            : {
     568         [ +  + ]:        155 :   if (!is_ell(x)) member_err("eta",x);
     569         [ +  + ]:         45 :   if (!ell_is_complex(x)) pari_err_TYPE("eta [not defined over C]",x);
     570                 :         25 :   return ellR_eta(x, ellR_get_prec(x));
     571                 :            : }
     572                 :            : 
     573                 :            : GEN
     574                 :        135 : member_area(GEN x)
     575                 :            : {
     576                 :            :   GEN w, w1, w2, a,b,c,d;
     577                 :            :   long prec;
     578         [ +  + ]:        135 :   if (!is_ell(x)) member_err("area",x);
     579         [ +  + ]:         25 :   if (!ell_is_complex(x)) pari_err_TYPE("area [not defined over C]",x);
     580                 :          5 :   prec = ellR_get_prec(x);
     581                 :          5 :   w = ellR_omega(x, prec);
     582                 :          5 :   w1 = gel(w,1); a = real_i(w1); b = imag_i(w1);
     583                 :          5 :   w2 = gel(w,2); c = real_i(w2); d = imag_i(w2);
     584                 :          5 :   return gabs(gsub(gmul(a,d),gmul(b,c)), prec);
     585                 :            : }
     586                 :            : 
     587                 :            : GEN
     588                 :        150 : member_tate(GEN x)
     589                 :            : {
     590                 :            :   long prec;
     591         [ +  + ]:        150 :   if (!is_ell(x)) member_err("tate",x);
     592         [ +  + ]:         40 :   if (ell_get_type(x) != t_ELL_Qp)
     593                 :         20 :     pari_err_TYPE("tate [not defined over Qp]",x);
     594                 :         20 :   prec = ellQp_get_prec(x);
     595                 :         20 :   return ellQp_Tate_uniformization(x, prec);
     596                 :            : }

Generated by: LCOV version 1.9