[BACK]Return to speed.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / tune

Diff for /OpenXM_contrib/gmp/tune/Attic/speed.c between version 1.1 and 1.1.1.2

version 1.1, 2000/09/09 14:13:19 version 1.1.1.2, 2003/08/25 16:06:38
Line 1 
Line 1 
 /* Speed measuring program. */  /* Speed measuring program.
   
 /*  Copyright 1999, 2000, 2001, 2002 Free Software Foundation, Inc.
 Copyright (C) 1999, 2000 Free Software Foundation, Inc.  
   
 This file is part of the GNU MP Library.  This file is part of the GNU MP Library.
   
Line 18  License for more details.
Line 17  License for more details.
 You should have received a copy of the GNU Lesser General Public License  You should have received a copy of the GNU Lesser General Public License
 along with the GNU MP Library; see the file COPYING.LIB.  If not, write to  along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
 the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,  the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 MA 02111-1307, USA.  MA 02111-1307, USA. */
 */  
   
 /* Usage message is in the code below, run with no arguments to print it.  /* Usage message is in the code below, run with no arguments to print it.
    See README for interesting applications.     See README for interesting applications.
Line 35  MA 02111-1307, USA.
Line 33  MA 02111-1307, USA.
   
 #include "config.h"  #include "config.h"
   
 #if HAVE_GETOPT_H  
 #include <getopt.h>  /* for getopt_long() */  
 #endif  
 #include <limits.h>  #include <limits.h>
 #include <stdio.h>  #include <stdio.h>
 #include <stdlib.h>  #include <stdlib.h>
 #include <string.h>  #include <string.h>
 #include <unistd.h> /* for getpid() */  
 #include <sys/time.h>  /* for struct timeval for sys/resource.h */  #if HAVE_UNISTD_H
   #include <unistd.h>  /* for getpid, R_OK */
   #endif
   
   #if TIME_WITH_SYS_TIME
   # include <sys/time.h>  /* for struct timeval */
   # include <time.h>
   #else
   # if HAVE_SYS_TIME_H
   #  include <sys/time.h>
   # else
   #  include <time.h>
   # endif
   #endif
   
   #if HAVE_SYS_RESOURCE_H
 #include <sys/resource.h>  /* for getrusage() */  #include <sys/resource.h>  /* for getrusage() */
   #endif
   
   
 #include "gmp.h"  #include "gmp.h"
 #include "gmp-impl.h"  #include "gmp-impl.h"
   #include "longlong.h"  /* for the benefit of speed-many.c */
   #include "tests.h"
 #include "speed.h"  #include "speed.h"
   
   
 #if !HAVE_DECL_OPTARG  #if !HAVE_DECL_OPTARG
 extern char *optarg;  extern char *optarg;
 extern int optind, opterr;  extern int optind, opterr;
Line 63  extern int optind, opterr;
Line 77  extern int optind, opterr;
 #ifdef SPEED_EXTRA_PROTOS  #ifdef SPEED_EXTRA_PROTOS
 SPEED_EXTRA_PROTOS  SPEED_EXTRA_PROTOS
 #endif  #endif
   #ifdef SPEED_EXTRA_PROTOS2
   SPEED_EXTRA_PROTOS2
   #endif
   
   
 #define MPN_FILL(ptr, size, n)                  \  #define MPN_FILL(ptr, size, n)          \
   do {                                          \    do {                                  \
     mp_size_t  i;                               \      mp_size_t __i;                      \
     for (i = 0; i < size; i++)                  \      ASSERT ((size) >= 0);               \
       ptr[i] = n;                               \      for (__i = 0; __i < (size); __i++)  \
         (ptr)[__i] = (n);                 \
   } while (0)    } while (0)
   
   
   #if BITS_PER_MP_LIMB == 32
   #define GMP_NUMB_0xAA  (CNST_LIMB(0xAAAAAAAA) & GMP_NUMB_MASK)
   #endif
   #if BITS_PER_MP_LIMB == 64
   #define GMP_NUMB_0xAA  (CNST_LIMB(0xAAAAAAAAAAAAAAAA) & GMP_NUMB_MASK)
   #endif
   
   
 #define CMP_ABSOLUTE     1  #define CMP_ABSOLUTE     1
 #define CMP_RATIO        2  #define CMP_RATIO        2
 #define CMP_DIFFERENCE   3  #define CMP_DIFFERENCE   3
Line 86  int  option_unit = UNIT_SECONDS;
Line 113  int  option_unit = UNIT_SECONDS;
 #define DATA_RANDOM   1  #define DATA_RANDOM   1
 #define DATA_RANDOM2  2  #define DATA_RANDOM2  2
 #define DATA_ZEROS    3  #define DATA_ZEROS    3
 #define DATA_FFS      4  #define DATA_AAS      4
 #define DATA_2FD      5  #define DATA_FFS      5
   #define DATA_2FD      6
 int  option_data = DATA_RANDOM;  int  option_data = DATA_RANDOM;
   
 int        option_square = 0;  int        option_square = 0;
Line 107  struct speed_params  sp;
Line 135  struct speed_params  sp;
   
 #define COLUMN_WIDTH  13  /* for the free-form output */  #define COLUMN_WIDTH  13  /* for the free-form output */
   
 #define FLAG_R            (1<<0)  #define FLAG_R            (1<<0)  /* require ".r" */
 #define FLAG_R_OPTIONAL   (1<<1)  #define FLAG_R_OPTIONAL   (1<<1)  /* optional ".r" */
 #define FLAG_RSIZE        (1<<2)  #define FLAG_RSIZE        (1<<2)
   #define FLAG_NODATA       (1<<3)  /* don't alloc xp, yp */
   
 const struct routine_t {  const struct routine_t {
   /* constants */    /* constants */
Line 123  const struct routine_t {
Line 152  const struct routine_t {
   { "noop_wxs",          speed_noop_wxs             },    { "noop_wxs",          speed_noop_wxs             },
   { "noop_wxys",         speed_noop_wxys            },    { "noop_wxys",         speed_noop_wxys            },
   
   { "mpn_add_n",         speed_mpn_add_n            },    { "mpn_add_n",         speed_mpn_add_n,     FLAG_R_OPTIONAL },
   { "mpn_sub_n",         speed_mpn_sub_n            },    { "mpn_sub_n",         speed_mpn_sub_n,     FLAG_R_OPTIONAL },
   { "mpn_add_n_self",    speed_mpn_add_n_self       },  
   { "mpn_add_n_inplace", speed_mpn_add_n_inplace    },  
   
   { "mpn_addmul_1",      speed_mpn_addmul_1,  FLAG_R },    { "mpn_addmul_1",      speed_mpn_addmul_1,  FLAG_R },
   { "mpn_submul_1",      speed_mpn_submul_1,  FLAG_R },    { "mpn_submul_1",      speed_mpn_submul_1,  FLAG_R },
   { "mpn_mul_1",         speed_mpn_mul_1,     FLAG_R },    { "mpn_mul_1",         speed_mpn_mul_1,     FLAG_R },
     { "mpn_mul_1_inplace", speed_mpn_mul_1_inplace, FLAG_R },
   #if HAVE_NATIVE_mpn_mul_2
     { "mpn_mul_2",         speed_mpn_mul_2,     FLAG_R },
   #endif
   
   { "mpn_divrem_1",      speed_mpn_divrem_1,  FLAG_R },    { "mpn_divrem_1",      speed_mpn_divrem_1,  FLAG_R },
   { "mpn_divrem_1f",     speed_mpn_divrem_1f, FLAG_R },    { "mpn_divrem_1f",     speed_mpn_divrem_1f, FLAG_R },
Line 142  const struct routine_t {
Line 173  const struct routine_t {
 #if HAVE_NATIVE_mpn_mod_1c  #if HAVE_NATIVE_mpn_mod_1c
   { "mpn_mod_1c",        speed_mpn_mod_1c,    FLAG_R },    { "mpn_mod_1c",        speed_mpn_mod_1c,    FLAG_R },
 #endif  #endif
     { "mpn_preinv_divrem_1",  speed_mpn_preinv_divrem_1,  FLAG_R },
     { "mpn_preinv_divrem_1f", speed_mpn_preinv_divrem_1f, FLAG_R },
     { "mpn_preinv_mod_1",  speed_mpn_preinv_mod_1, FLAG_R },
   
     { "mpn_divrem_1_div",  speed_mpn_divrem_1_div,  FLAG_R },
     { "mpn_divrem_1_inv",  speed_mpn_divrem_1_inv,  FLAG_R },
     { "mpn_divrem_1f_div", speed_mpn_divrem_1f_div, FLAG_R },
     { "mpn_divrem_1f_inv", speed_mpn_divrem_1f_inv, FLAG_R },
     { "mpn_mod_1_div",     speed_mpn_mod_1_div,     FLAG_R },
     { "mpn_mod_1_inv",     speed_mpn_mod_1_inv,     FLAG_R },
   
   { "mpn_divrem_2",      speed_mpn_divrem_2,        },    { "mpn_divrem_2",      speed_mpn_divrem_2,        },
   { "mpn_divexact_by3",  speed_mpn_divexact_by3     },    { "mpn_divrem_2_div",  speed_mpn_divrem_2_div,    },
     { "mpn_divrem_2_inv",  speed_mpn_divrem_2_inv,    },
   
   { "mpn_bz_divrem_n",   speed_mpn_bz_divrem_n      },    { "mpn_divexact_1",    speed_mpn_divexact_1,    FLAG_R },
   { "mpn_bz_divrem_sb",  speed_mpn_bz_divrem_sb     },    { "mpn_divexact_by3",  speed_mpn_divexact_by3          },
   { "mpn_bz_tdiv_qr",    speed_mpn_bz_tdiv_qr       },  
   
   #if HAVE_NATIVE_mpn_modexact_1c_odd
     { "mpn_modexact_1_odd",  speed_mpn_modexact_1_odd,  FLAG_R },
   #endif
     { "mpn_modexact_1c_odd", speed_mpn_modexact_1c_odd, FLAG_R },
   
     { "mpn_dc_tdiv_qr",       speed_mpn_dc_tdiv_qr       },
     { "mpn_dc_divrem_n",      speed_mpn_dc_divrem_n      },
     { "mpn_dc_divrem_sb",     speed_mpn_dc_divrem_sb     },
     { "mpn_dc_divrem_sb_div", speed_mpn_dc_divrem_sb_div },
     { "mpn_dc_divrem_sb_inv", speed_mpn_dc_divrem_sb_inv },
   
     { "mpn_sb_divrem_m3",     speed_mpn_sb_divrem_m3     },
     { "mpn_sb_divrem_m3_div", speed_mpn_sb_divrem_m3_div },
     { "mpn_sb_divrem_m3_inv", speed_mpn_sb_divrem_m3_inv },
   
   { "mpn_lshift",        speed_mpn_lshift, FLAG_R   },    { "mpn_lshift",        speed_mpn_lshift, FLAG_R   },
   { "mpn_rshift",        speed_mpn_rshift, FLAG_R   },    { "mpn_rshift",        speed_mpn_rshift, FLAG_R   },
   
   { "mpn_and_n",         speed_mpn_and_n            },    { "mpn_and_n",         speed_mpn_and_n,  FLAG_R_OPTIONAL },
   { "mpn_andn_n",        speed_mpn_andn_n           },    { "mpn_andn_n",        speed_mpn_andn_n, FLAG_R_OPTIONAL },
   { "mpn_nand_n",        speed_mpn_nand_n           },    { "mpn_nand_n",        speed_mpn_nand_n, FLAG_R_OPTIONAL },
   { "mpn_ior_n",         speed_mpn_ior_n            },    { "mpn_ior_n",         speed_mpn_ior_n,  FLAG_R_OPTIONAL },
   { "mpn_iorn_n",        speed_mpn_iorn_n           },    { "mpn_iorn_n",        speed_mpn_iorn_n, FLAG_R_OPTIONAL },
   { "mpn_nior_n",        speed_mpn_nior_n           },    { "mpn_nior_n",        speed_mpn_nior_n, FLAG_R_OPTIONAL },
   { "mpn_xor_n",         speed_mpn_xor_n            },    { "mpn_xor_n",         speed_mpn_xor_n,  FLAG_R_OPTIONAL },
   { "mpn_xnor_n",        speed_mpn_xnor_n           },    { "mpn_xnor_n",        speed_mpn_xnor_n, FLAG_R_OPTIONAL },
     { "mpn_com_n",         speed_mpn_com_n            },
   
   { "mpn_popcount",      speed_mpn_popcount         },    { "mpn_popcount",      speed_mpn_popcount         },
   { "mpn_hamdist",       speed_mpn_hamdist          },    { "mpn_hamdist",       speed_mpn_hamdist          },
   
   { "mpn_gcdext",        speed_mpn_gcdext           },    { "mpn_gcd_1",         speed_mpn_gcd_1,  FLAG_R_OPTIONAL },
   { "mpn_gcd",           speed_mpn_gcd              },    { "mpn_gcd_1N",        speed_mpn_gcd_1N, FLAG_R_OPTIONAL },
   { "mpn_gcd_1",         speed_mpn_gcd_1, FLAG_R_OPTIONAL },  
   
     { "mpn_gcd",           speed_mpn_gcd                    },
     { "mpn_gcd_binary",    speed_mpn_gcd_binary             },
     { "find_a",            speed_find_a,        FLAG_NODATA },
   #if HAVE_NATIVE_mpn_gcd_finda
     { "mpn_gcd_finda",     speed_mpn_gcd_finda, FLAG_NODATA },
   #endif
   
     { "mpn_gcdext",            speed_mpn_gcdext            },
     { "mpn_gcdext_single",     speed_mpn_gcdext_single     },
     { "mpn_gcdext_double",     speed_mpn_gcdext_double     },
     { "mpn_gcdext_one_single", speed_mpn_gcdext_one_single },
     { "mpn_gcdext_one_double", speed_mpn_gcdext_one_double },
   
     { "mpz_jacobi",        speed_mpz_jacobi           },
   { "mpn_jacobi_base",   speed_mpn_jacobi_base      },    { "mpn_jacobi_base",   speed_mpn_jacobi_base      },
     { "mpn_jacobi_base_1", speed_mpn_jacobi_base_1    },
     { "mpn_jacobi_base_2", speed_mpn_jacobi_base_2    },
     { "mpn_jacobi_base_3", speed_mpn_jacobi_base_3    },
   
   { "mpn_mul_basecase",  speed_mpn_mul_basecase, FLAG_R_OPTIONAL },    { "mpn_mul_basecase",  speed_mpn_mul_basecase, FLAG_R_OPTIONAL },
   { "mpn_sqr_basecase",  speed_mpn_sqr_basecase     },    { "mpn_sqr_basecase",  speed_mpn_sqr_basecase     },
   #if HAVE_NATIVE_mpn_sqr_diagonal
     { "mpn_sqr_diagonal",  speed_mpn_sqr_diagonal     },
   #endif
   
   { "mpn_mul_n",         speed_mpn_mul_n            },    { "mpn_mul_n",         speed_mpn_mul_n            },
   { "mpn_sqr_n",         speed_mpn_sqr_n            },    { "mpn_sqr_n",         speed_mpn_sqr_n            },
Line 181  const struct routine_t {
Line 257  const struct routine_t {
   { "mpn_kara_sqr_n",    speed_mpn_kara_sqr_n       },    { "mpn_kara_sqr_n",    speed_mpn_kara_sqr_n       },
   { "mpn_toom3_mul_n",   speed_mpn_toom3_mul_n      },    { "mpn_toom3_mul_n",   speed_mpn_toom3_mul_n      },
   { "mpn_toom3_sqr_n",   speed_mpn_toom3_sqr_n      },    { "mpn_toom3_sqr_n",   speed_mpn_toom3_sqr_n      },
   { "mpn_mul_fft_full",     speed_mpn_mul_fft_full     },    { "mpn_mul_fft_full",      speed_mpn_mul_fft_full      },
   { "mpn_mul_fft_full_sqr", speed_mpn_mul_fft_full_sqr },    { "mpn_mul_fft_full_sqr",  speed_mpn_mul_fft_full_sqr  },
   
   { "mpn_mul_fft",       speed_mpn_mul_fft,     FLAG_R_OPTIONAL },    { "mpn_mul_fft",       speed_mpn_mul_fft,     FLAG_R_OPTIONAL },
   { "mpn_mul_fft_sqr",   speed_mpn_mul_fft_sqr, FLAG_R_OPTIONAL },    { "mpn_mul_fft_sqr",   speed_mpn_mul_fft_sqr, FLAG_R_OPTIONAL },
   
     { "mpn_toom3_mul_n_mpn",   speed_mpn_toom3_mul_n_mpn   },
     { "mpn_toom3_mul_n_open",  speed_mpn_toom3_mul_n_open  },
     { "mpn_toom3_sqr_n_mpn",   speed_mpn_toom3_sqr_n_mpn   },
     { "mpn_toom3_sqr_n_open",  speed_mpn_toom3_sqr_n_open  },
   
     { "mpn_get_str",       speed_mpn_get_str,  FLAG_R_OPTIONAL },
   
     { "mpn_set_str",          speed_mpn_set_str,          FLAG_R_OPTIONAL },
     { "mpn_set_str_basecase", speed_mpn_set_str_basecase, FLAG_R_OPTIONAL },
     { "mpn_set_str_subquad",  speed_mpn_set_str_subquad,  FLAG_R_OPTIONAL },
   
     { "mpn_sqrtrem",       speed_mpn_sqrtrem          },
   
     { "mpn_fib2_ui",       speed_mpn_fib2_ui,    FLAG_NODATA },
     { "mpz_fib_ui",        speed_mpz_fib_ui,     FLAG_NODATA },
     { "mpz_fib2_ui",       speed_mpz_fib2_ui,    FLAG_NODATA },
     { "mpz_lucnum_ui",     speed_mpz_lucnum_ui,  FLAG_NODATA },
     { "mpz_lucnum2_ui",    speed_mpz_lucnum2_ui, FLAG_NODATA },
   
   { "mpz_add",           speed_mpz_add              },    { "mpz_add",           speed_mpz_add              },
   { "mpz_bin_uiui",      speed_mpz_bin_uiui, FLAG_R_OPTIONAL },    { "mpz_bin_uiui",      speed_mpz_bin_uiui, FLAG_NODATA | FLAG_R_OPTIONAL },
   { "mpz_fac_ui",        speed_mpz_fac_ui           },    { "mpz_fac_ui",        speed_mpz_fac_ui,   FLAG_NODATA   },
   { "mpz_fib_ui",        speed_mpz_fib_ui           },  
   { "mpz_powm",          speed_mpz_powm             },    { "mpz_powm",          speed_mpz_powm             },
     { "mpz_powm_mod",      speed_mpz_powm_mod         },
     { "mpz_powm_redc",     speed_mpz_powm_redc        },
     { "mpz_powm_ui",       speed_mpz_powm_ui          },
   
     { "mpz_mod",           speed_mpz_mod              },
     { "redc",              speed_redc                 },
   
   { "MPN_COPY",          speed_MPN_COPY             },    { "MPN_COPY",          speed_MPN_COPY             },
   { "MPN_COPY_INCR",     speed_MPN_COPY_INCR        },    { "MPN_COPY_INCR",     speed_MPN_COPY_INCR        },
   { "MPN_COPY_DECR",     speed_MPN_COPY_DECR        },    { "MPN_COPY_DECR",     speed_MPN_COPY_DECR        },
   { "memcpy",            speed_memcpy               },    { "memcpy",            speed_memcpy               },
   #if HAVE_NATIVE_mpn_copyi
     { "mpn_copyi",         speed_mpn_copyi            },
   #endif
   #if HAVE_NATIVE_mpn_copyd
     { "mpn_copyd",         speed_mpn_copyd            },
   #endif
   
   { "modlimb_invert",    speed_modlimb_invert       },    { "MPN_ZERO",          speed_MPN_ZERO             },
   
   { "malloc_free",                 speed_malloc_free                 },    { "modlimb_invert",       speed_modlimb_invert,       FLAG_NODATA },
   { "malloc_realloc_free",         speed_malloc_realloc_free         },    { "modlimb_invert_mul1",  speed_modlimb_invert_mul1,  FLAG_NODATA },
   { "mp_allocate_free",            speed_mp_allocate_free            },    { "modlimb_invert_loop",  speed_modlimb_invert_loop,  FLAG_NODATA },
   { "mp_allocate_reallocate_free", speed_mp_allocate_reallocate_free },    { "modlimb_invert_cond",  speed_modlimb_invert_cond,  FLAG_NODATA },
   { "mpz_init_clear",              speed_mpz_init_clear              },    { "modlimb_invert_arith", speed_modlimb_invert_arith, FLAG_NODATA },
   { "mpq_init_clear",              speed_mpq_init_clear              },  
   { "mpf_init_clear",              speed_mpf_init_clear              },  
   { "mpz_init_realloc_clear",      speed_mpz_init_realloc_clear      },  
   
     { "malloc_free",                  speed_malloc_free                  },
     { "malloc_realloc_free",          speed_malloc_realloc_free          },
     { "gmp_allocate_free",            speed_gmp_allocate_free            },
     { "gmp_allocate_reallocate_free", speed_gmp_allocate_reallocate_free },
     { "mpz_init_clear",               speed_mpz_init_clear               },
     { "mpq_init_clear",               speed_mpq_init_clear               },
     { "mpf_init_clear",               speed_mpf_init_clear               },
     { "mpz_init_realloc_clear",       speed_mpz_init_realloc_clear       },
   
   { "umul_ppmm",         speed_umul_ppmm,     FLAG_R_OPTIONAL },    { "umul_ppmm",         speed_umul_ppmm,     FLAG_R_OPTIONAL },
 #if HAVE_NATIVE_mpn_umul_ppmm  #if HAVE_NATIVE_mpn_umul_ppmm
   { "mpn_umul_ppmm",     speed_mpn_umul_ppmm, FLAG_R_OPTIONAL },    { "mpn_umul_ppmm",     speed_mpn_umul_ppmm, FLAG_R_OPTIONAL },
 #endif  #endif
   
     { "count_leading_zeros",  speed_count_leading_zeros,  FLAG_NODATA | FLAG_R_OPTIONAL },
     { "count_trailing_zeros", speed_count_trailing_zeros, FLAG_NODATA | FLAG_R_OPTIONAL },
   
   { "udiv_qrnnd",             speed_udiv_qrnnd,             FLAG_R_OPTIONAL },    { "udiv_qrnnd",             speed_udiv_qrnnd,             FLAG_R_OPTIONAL },
   { "udiv_qrnnd_preinv",      speed_udiv_qrnnd_preinv,      FLAG_R_OPTIONAL },    { "udiv_qrnnd_preinv",      speed_udiv_qrnnd_preinv,      FLAG_R_OPTIONAL },
   { "udiv_qrnnd_preinv2norm", speed_udiv_qrnnd_preinv2norm, FLAG_R_OPTIONAL },    { "udiv_qrnnd_preinv2norm", speed_udiv_qrnnd_preinv2norm, FLAG_R_OPTIONAL },
     { "udiv_qrnnd_c",           speed_udiv_qrnnd_c,           FLAG_R_OPTIONAL },
 #if HAVE_NATIVE_mpn_udiv_qrnnd  #if HAVE_NATIVE_mpn_udiv_qrnnd
   { "mpn_udiv_qrnnd",         speed_mpn_udiv_qrnnd,         FLAG_R_OPTIONAL },    { "mpn_udiv_qrnnd",         speed_mpn_udiv_qrnnd,         FLAG_R_OPTIONAL },
 #endif  #endif
     { "invert_limb",            speed_invert_limb,            FLAG_R_OPTIONAL },
   
     { "operator_div",           speed_operator_div,           FLAG_R_OPTIONAL },
     { "operator_mod",           speed_operator_mod,           FLAG_R_OPTIONAL },
   
 #ifdef SPEED_EXTRA_ROUTINES  #ifdef SPEED_EXTRA_ROUTINES
   SPEED_EXTRA_ROUTINES    SPEED_EXTRA_ROUTINES
 #endif  #endif
   #ifdef SPEED_EXTRA_ROUTINES2
     SPEED_EXTRA_ROUTINES2
   #endif
 };  };
   
   
 struct choice_t {  struct choice_t {
   const struct routine_t  *p;    const struct routine_t  *p;
   int                     r;    mp_limb_t               r;
     double                  scale;
   double                  time;    double                  time;
   int                     no_time;    int                     no_time;
   double                  prev_time;    double                  prev_time;
Line 252  data_fill (mp_ptr ptr, mp_size_t size)
Line 376  data_fill (mp_ptr ptr, mp_size_t size)
   case DATA_ZEROS:    case DATA_ZEROS:
     MPN_ZERO (ptr, size);      MPN_ZERO (ptr, size);
     break;      break;
     case DATA_AAS:
       MPN_FILL (ptr, size, GMP_NUMB_0xAA);
       break;
   case DATA_FFS:    case DATA_FFS:
     MPN_FILL (ptr, size, MP_LIMB_T_MAX);      MPN_FILL (ptr, size, GMP_NUMB_MAX);
     break;      break;
   case DATA_2FD:    case DATA_2FD:
     MPN_FILL (ptr, size, MP_LIMB_T_MAX);      MPN_FILL (ptr, size, GMP_NUMB_MAX);
     ptr[0] -= 2;      ptr[0] -= 2;
     break;      break;
   default:    default:
Line 277  void
Line 404  void
 run_one (FILE *fp, struct speed_params *s, mp_size_t prev_size)  run_one (FILE *fp, struct speed_params *s, mp_size_t prev_size)
 {  {
   const char  *first_open_fastest, *first_open_notfastest, *first_close;    const char  *first_open_fastest, *first_open_notfastest, *first_close;
   int         i, fastest;    int         i, fastest, want_data;
   double      fastest_time;    double      fastest_time;
   TMP_DECL (marker);    TMP_DECL (marker);
   
   TMP_MARK (marker);    TMP_MARK (marker);
   sp.xp = SPEED_TMP_ALLOC_LIMBS (s->size, s->align_xp);  
   sp.yp = SPEED_TMP_ALLOC_LIMBS (s->size, s->align_yp);  
   
   data_fill (s->xp, s->size);    /* allocate data, unless all routines are NODATA */
   data_fill (s->yp, s->size);    want_data = 0;
     for (i = 0; i < num_choices; i++)
       want_data |= ((choice[i].p->flag & FLAG_NODATA) == 0);
   
     if (want_data)
       {
         sp.xp = SPEED_TMP_ALLOC_LIMBS (s->size, s->align_xp);
         sp.yp = SPEED_TMP_ALLOC_LIMBS (s->size, s->align_yp);
   
         data_fill (s->xp, s->size);
         data_fill (s->yp, s->size);
       }
     else
       {
         sp.xp = NULL;
         sp.yp = NULL;
       }
   
   if (prev_size == -1 && option_cmp == CMP_DIFFPREV)    if (prev_size == -1 && option_cmp == CMP_DIFFPREV)
     {      {
       first_open_fastest = "(#";        first_open_fastest = "(#";
Line 308  run_one (FILE *fp, struct speed_params *s, mp_size_t p
Line 449  run_one (FILE *fp, struct speed_params *s, mp_size_t p
       s->r = choice[i].r;        s->r = choice[i].r;
       choice[i].time = speed_measure (choice[i].p->fun, s);        choice[i].time = speed_measure (choice[i].p->fun, s);
       choice[i].no_time = (choice[i].time == -1.0);        choice[i].no_time = (choice[i].time == -1.0);
         choice[i].time *= choice[i].scale;
   
   
       /* Apply the effect of CMP_DIFFPREV, but the new choice[i].prev_time        /* Apply the effect of CMP_DIFFPREV, but the new choice[i].prev_time
          is before any differences.  */           is before any differences.  */
       {        {
Line 510  fclose_written (FILE *fp, const char *filename)
Line 651  fclose_written (FILE *fp, const char *filename)
   
   
 void  void
 run_gnuplot (void)  run_gnuplot (int argc, char *argv[])
 {  {
   char  *plot_filename;    char  *plot_filename;
   char  *data_filename;    char  *data_filename;
   FILE  *fp;    FILE  *fp;
   int   i;    int   i;
   
   plot_filename = (char *) (*_mp_allocate_func)    plot_filename = (char *) (*__gmp_allocate_func)
     (strlen (option_gnuplot_basename) + 20);      (strlen (option_gnuplot_basename) + 20);
   data_filename = (char *) (*_mp_allocate_func)    data_filename = (char *) (*__gmp_allocate_func)
     (strlen (option_gnuplot_basename) + 20);      (strlen (option_gnuplot_basename) + 20);
   
   sprintf (plot_filename, "%s.gnuplot", option_gnuplot_basename);    sprintf (plot_filename, "%s.gnuplot", option_gnuplot_basename);
Line 527  run_gnuplot (void)
Line 668  run_gnuplot (void)
   
   fp = fopen_for_write (plot_filename);    fp = fopen_for_write (plot_filename);
   
     fprintf (fp, "# Generated with:\n");
     fprintf (fp, "#");
     for (i = 0; i < argc; i++)
       fprintf (fp, " %s", argv[i]);
     fprintf (fp, "\n");
     fprintf (fp, "\n");
   
   /* Putting the key at the top left is usually good, and you can change it    /* Putting the key at the top left is usually good, and you can change it
      interactively if it's not. */       interactively if it's not. */
   fprintf (fp, "set key left\n");    fprintf (fp, "set key left\n");
Line 559  run_gnuplot (void)
Line 707  run_gnuplot (void)
 }  }
   
   
 /* Return a long with n many one bits (starting from the least significant) */  /* Return a limb with n many one bits (starting from the least significant) */
 #define LONG_ONES(n) \  
   ((n) == BITS_PER_LONGINT ? -1L : (n) == 0 ? 0L : (1L << (n)) - 1)  
   
 long  #define LIMB_ONES(n) \
     ((n) == BITS_PER_MP_LIMB ? MP_LIMB_T_MAX      \
       : (n) == 0 ? CNST_LIMB(0)                   \
       : (CNST_LIMB(1) << (n)) - 1)
   
   mp_limb_t
 r_string (const char *s)  r_string (const char *s)
 {  {
   const char  *s_orig = s;    const char  *s_orig = s;
   long  n;    long        n;
   
     if (strcmp (s, "aas") == 0)
       return GMP_NUMB_0xAA;
   
     {
       mpz_t      z;
       mp_limb_t  l;
       int        set, siz;
   
       mpz_init (z);
       set = mpz_set_str (z, s, 0);
       siz = SIZ(z);
       l = (siz == 0 ? 0 : siz > 0 ? PTR(z)[0] : -PTR(z)[0]);
       mpz_clear (z);
       if (set == 0)
         {
           if (siz > 1 || siz < -1)
             printf ("Warning, r parameter %s truncated to %d bits\n",
                     s_orig, BITS_PER_MP_LIMB);
           return l;
         }
     }
   
   if (s[0] == '0' && (s[1] == 'x' || s[1] == 'X'))    if (s[0] == '0' && (s[1] == 'x' || s[1] == 'X'))
     n = strtoul (s+2, (char **) &s, 16);      n = strtoul (s+2, (char **) &s, 16);
   else    else
Line 577  r_string (const char *s)
Line 750  r_string (const char *s)
   if (strcmp (s, "bits") == 0)    if (strcmp (s, "bits") == 0)
     {      {
       mp_limb_t  l;        mp_limb_t  l;
       if (n > BITS_PER_LONGINT)        if (n > BITS_PER_MP_LIMB)
         {          {
           fprintf (stderr, "%ld bit parameter invalid (max %d bits)\n",            fprintf (stderr, "%ld bit parameter invalid (max %d bits)\n",
                    n, BITS_PER_LONGINT);                     n, BITS_PER_MP_LIMB);
           exit (1);            exit (1);
         }          }
       mpn_random (&l, 1);        mpn_random (&l, 1);
       return (l | (1 << (n-1))) & LONG_ONES(n);        return (l | (1 << (n-1))) & LIMB_ONES(n);
     }      }
   else  if (strcmp (s, "ones") == 0)    else  if (strcmp (s, "ones") == 0)
     {      {
       if (n > BITS_PER_LONGINT)        if (n > BITS_PER_MP_LIMB)
         {          {
           fprintf (stderr, "%ld bit parameter invalid (max %d bits)\n",            fprintf (stderr, "%ld bit parameter invalid (max %d bits)\n",
                    n, BITS_PER_LONGINT);                     n, BITS_PER_MP_LIMB);
           exit (1);            exit (1);
         }          }
       return LONG_ONES (n);        return LIMB_ONES (n);
     }      }
   else if (*s != '\0')    else if (*s != '\0')
     {      {
Line 607  r_string (const char *s)
Line 780  r_string (const char *s)
   
   
 void  void
 routine_find (struct choice_t *c, const char *s)  routine_find (struct choice_t *c, const char *s_orig)
 {  {
     const char  *s;
   int     i;    int     i;
   size_t  nlen;    size_t  nlen;
   
     c->name = s_orig;
     s = strchr (s_orig, '*');
     if (s != NULL)
       {
         c->scale = atof(s_orig);
         s++;
       }
     else
       {
         c->scale = 1.0;
         s = s_orig;
       }
   
   for (i = 0; i < numberof (routine); i++)    for (i = 0; i < numberof (routine); i++)
     {      {
       nlen = strlen (routine[i].name);        nlen = strlen (routine[i].name);
Line 624  routine_find (struct choice_t *c, const char *s)
Line 811  routine_find (struct choice_t *c, const char *s)
   
           if (! (routine[i].flag & (FLAG_R|FLAG_R_OPTIONAL)))            if (! (routine[i].flag & (FLAG_R|FLAG_R_OPTIONAL)))
             {              {
               fprintf (stderr, "Choice %s bad: doesn't take a \".<r>\" paramater\n", s);                fprintf (stderr,
                          "Choice %s bad: doesn't take a \".<r>\" parameter\n",
                          s_orig);
               exit (1);                exit (1);
             }              }
   
           c->p = &routine[i];            c->p = &routine[i];
           c->r = r_string (s + nlen + 1);            c->r = r_string (s + nlen + 1);
           c->name = s;  
           return;            return;
         }          }
   
Line 640  routine_find (struct choice_t *c, const char *s)
Line 828  routine_find (struct choice_t *c, const char *s)
   
           if (routine[i].flag & FLAG_R)            if (routine[i].flag & FLAG_R)
             {              {
               fprintf (stderr, "Choice %s bad: needs a \".<r>\" paramater\n", s);                fprintf (stderr,
                          "Choice %s bad: needs a \".<r>\" parameter\n",
                          s_orig);
               exit (1);                exit (1);
             }              }
   
           c->p = &routine[i];            c->p = &routine[i];
           c->r = 0;            c->r = 0;
           c->name = s;  
           return;            return;
         }          }
     }      }
   
   fprintf (stderr, "Choice %s unrecognised\n", s);    fprintf (stderr, "Choice %s unrecognised\n", s_orig);
   exit (1);    exit (1);
 }  }
   
Line 661  usage (void)
Line 850  usage (void)
 {  {
   int  i;    int  i;
   
     speed_time_init ();
   
   printf ("\    printf ("\
 Usage: speed [-options] -s size <routine>...\n\  Usage: speed [-options] -s size <routine>...\n\
 Measure the speed of some routines.\n\  Measure the speed of some routines.\n\
Line 678  Times are in seconds, accuracy is shown.\n\
Line 869  Times are in seconds, accuracy is shown.\n\
    -C         show times in cycles per limb\n\     -C         show times in cycles per limb\n\
    -u         print resource usage (memory) at end\n\     -u         print resource usage (memory) at end\n\
    -P name    output plot files \"name.gnuplot\" and \"name.data\"\n\     -P name    output plot files \"name.gnuplot\" and \"name.data\"\n\
    -a <type>  use given data: random(default), random2, zeros, ffs\n\     -a <type>  use given data: random(default), random2, zeros, aas, ffs, 2fd\n\
    -x, -y, -w, -W <align>  specify data alignments, sources and dests\n\     -x, -y, -w, -W <align>  specify data alignments, sources and dests\n\
    -o addrs   print addresses of data blocks\n\     -o addrs   print addresses of data blocks\n\
 \n\  \n\
Line 710  The available routines are as follows.\n\
Line 901  The available routines are as follows.\n\
   printf ("\n\    printf ("\n\
 Routines with a \".r\" need an extra parameter, for example mpn_lshift.6\n\  Routines with a \".r\" need an extra parameter, for example mpn_lshift.6\n\
 r should be in decimal, or use 0xN for hexadecimal.\n\  r should be in decimal, or use 0xN for hexadecimal.\n\
 Special forms for r are Nbits for a random N bit number, and Nones for N one\n\  
 bits.\n\  
 \n\  \n\
   Special forms for r are \"<N>bits\" for a random N bit number, \"<N>ones\" for\n\
   N one bits, or \"aas\" for 0xAA..AA.\n\
   \n\
 Times for sizes out of the range accepted by a routine are shown as 0.\n\  Times for sizes out of the range accepted by a routine are shown as 0.\n\
 The fastest routine at each size is marked with a # (free form output only).\n\  The fastest routine at each size is marked with a # (free form output only).\n\
 \n\  \n\
Line 726  Quickplot home page http://www.kachinatech.com/~quickp
Line 918  Quickplot home page http://www.kachinatech.com/~quickp
 int  int
 main (int argc, char *argv[])  main (int argc, char *argv[])
 {  {
 #define OPTION_ALIGN     1000  
 #if HAVE_GETOPT_LONG  
   static const struct option  longopts[] = {  
     { "align",           required_argument, NULL, OPTION_ALIGN    },  
     { "align-x",         required_argument, NULL, 'x' },  
     { "align-y",         required_argument, NULL, 'y' },  
     { "align-w",         required_argument, NULL, 'w' },  
     { "align-w2",        required_argument, NULL, 'W' },  
     { "data",            required_argument, NULL, 'a' },  
     { "cycles",          no_argument,       NULL, 'c' },  
     { "cycles-per-limb", no_argument,       NULL, 'C' },  
     { "diff",            no_argument,       NULL, 'd' },  
     { "diff-prev",       no_argument,       NULL, 'D' },  
     { "difference",      no_argument,       NULL, 'd' },  
     { "difference-prev", no_argument,       NULL, 'D' },  
     { "factor",          required_argument, NULL, 'f' },  
     { "plot",            no_argument,       NULL, 'P' },  
     { "precision",       required_argument, NULL, 'p' },  
     { "ratio",           no_argument,       NULL, 'r' },  
     { "randomize",       no_argument,       NULL, 'R' },  
     { "sizes",           required_argument, NULL, 's' },  
     { "step",            required_argument, NULL, 't' },  
     { "resources",       required_argument, NULL, 'u' },  
     { "uncached",        no_argument,       NULL, 'z' },  
     { NULL }  
   };  
 #endif  
   
   int  i;    int  i;
   int  opt;    int  opt;
   
   /* Unbuffered so output goes straight out when directed to a pipe or file    /* Unbuffered so output goes straight out when directed to a pipe or file
      and isn't lost if you kill the program half way.  */       and isn't lost on killing the program half way.  */
   setbuf (stdout, NULL);    setbuf (stdout, NULL);
   
 #define OPTSTRING  "a:CcDdEFf:o:p:P:rRs:t:ux:y:w:W:z"    for (;;)
 #if HAVE_GETOPT_LONG      {
   while ((opt = getopt_long(argc, argv, OPTSTRING, longopts, NULL))        opt = getopt(argc, argv, "a:CcDdEFf:o:p:P:rRs:t:ux:y:w:W:z");
          != EOF)        if (opt == EOF)
 #else          break;
     while ((opt = getopt(argc, argv, OPTSTRING)) != EOF)  
 #endif  
       {  
         switch (opt) {  
         case 'a':  
           if (strcmp (optarg, "random") == 0)       option_data = DATA_RANDOM;  
           else if (strcmp (optarg, "random2") == 0) option_data = DATA_RANDOM2;  
           else if (strcmp (optarg, "zeros") == 0)   option_data = DATA_ZEROS;  
           else if (strcmp (optarg, "ffs") == 0)     option_data = DATA_FFS;  
           else if (strcmp (optarg, "2fd") == 0)     option_data = DATA_2FD;  
           else  
             {  
               fprintf (stderr, "unrecognised data option: %s\n", optarg);  
               exit (1);  
             }  
           break;  
         case 'C':  
           if (option_unit  != UNIT_SECONDS) goto bad_unit;  
           option_unit = UNIT_CYCLESPERLIMB;  
           break;  
         case 'c':  
           if (option_unit != UNIT_SECONDS)  
             {  
             bad_unit:  
               fprintf (stderr, "cannot use more than one of -c, -C\n");  
               exit (1);  
             }  
           option_unit = UNIT_CYCLES;  
           break;  
         case 'D':  
           if (option_cmp != CMP_ABSOLUTE) goto bad_cmp;  
           option_cmp = CMP_DIFFPREV;  
           break;  
         case 'd':  
           if (option_cmp != CMP_ABSOLUTE)  
             {  
             bad_cmp:  
               fprintf (stderr, "cannot use more than one of -d, -D, -r\n");  
               exit (1);  
             }  
           option_cmp = CMP_DIFFERENCE;  
           break;  
         case 'E':  
           option_square = 1;  
           break;  
         case 'F':  
           option_square = 2;  
           break;  
         case 'f':  
           option_factor = atof (optarg);  
           if (option_factor <= 1.0)  
             {  
               fprintf (stderr, "-f factor must be > 1.0\n");  
               exit (1);  
             }  
           break;  
         case 'o':  
           speed_option_set (optarg);  
           break;  
         case 'P':  
           option_gnuplot = 1;  
           option_gnuplot_basename = optarg;  
           break;  
         case 'p':  
           speed_precision = atoi (optarg);  
           break;  
         case 'R':  
           option_seed = time (NULL);  
           break;  
         case 'r':  
           if (option_cmp != CMP_ABSOLUTE)  
             goto bad_cmp;  
           option_cmp = CMP_RATIO;  
           break;  
         case 's':  
           {  
             char  *s;  
             for (s = strtok (optarg, ","); s != NULL; s = strtok (NULL, ","))  
               {  
                 if (size_num == size_allocnum)  
                   {  
                     size_array = (struct size_array_t *)  
                       _mp_allocate_or_reallocate  
                       (size_array,  
                        size_allocnum * sizeof(size_array[0]),  
                        (size_allocnum+10) * sizeof(size_array[0]));  
                     size_allocnum += 10;  
                   }  
                 if (sscanf (s, "%ld-%ld",  
                             &size_array[size_num].start,  
                             &size_array[size_num].end) != 2)  
                   {  
                     size_array[size_num].start = size_array[size_num].end  
                       = atol (s);  
                   }  
   
                 if (size_array[size_num].start < 1        switch (opt) {
                     || size_array[size_num].end < 1        case 'a':
                     || size_array[size_num].start > size_array[size_num].end)          if (strcmp (optarg, "random") == 0)       option_data = DATA_RANDOM;
                   {          else if (strcmp (optarg, "random2") == 0) option_data = DATA_RANDOM2;
                     fprintf (stderr, "invalid size parameter: %s\n", s);          else if (strcmp (optarg, "zeros") == 0)   option_data = DATA_ZEROS;
                     exit (1);          else if (strcmp (optarg, "aas") == 0)     option_data = DATA_AAS;
                   }          else if (strcmp (optarg, "ffs") == 0)     option_data = DATA_FFS;
           else if (strcmp (optarg, "2fd") == 0)     option_data = DATA_2FD;
                 size_num++;          else
               }            {
               fprintf (stderr, "unrecognised data option: %s\n", optarg);
               exit (1);
           }            }
           break;          break;
         case 't':        case 'C':
           option_step = atol (optarg);          if (option_unit  != UNIT_SECONDS) goto bad_unit;
           if (option_step < 1)          option_unit = UNIT_CYCLESPERLIMB;
           break;
         case 'c':
           if (option_unit != UNIT_SECONDS)
             {
             bad_unit:
               fprintf (stderr, "cannot use more than one of -c, -C\n");
               exit (1);
             }
           option_unit = UNIT_CYCLES;
           break;
         case 'D':
           if (option_cmp != CMP_ABSOLUTE) goto bad_cmp;
           option_cmp = CMP_DIFFPREV;
           break;
         case 'd':
           if (option_cmp != CMP_ABSOLUTE)
             {
             bad_cmp:
               fprintf (stderr, "cannot use more than one of -d, -D, -r\n");
               exit (1);
             }
           option_cmp = CMP_DIFFERENCE;
           break;
         case 'E':
           option_square = 1;
           break;
         case 'F':
           option_square = 2;
           break;
         case 'f':
           option_factor = atof (optarg);
           if (option_factor <= 1.0)
             {
               fprintf (stderr, "-f factor must be > 1.0\n");
               exit (1);
             }
           break;
         case 'o':
           speed_option_set (optarg);
           break;
         case 'P':
           option_gnuplot = 1;
           option_gnuplot_basename = optarg;
           break;
         case 'p':
           speed_precision = atoi (optarg);
           break;
         case 'R':
           option_seed = time (NULL);
           break;
         case 'r':
           if (option_cmp != CMP_ABSOLUTE)
             goto bad_cmp;
           option_cmp = CMP_RATIO;
           break;
         case 's':
           {
             char  *s;
             for (s = strtok (optarg, ","); s != NULL; s = strtok (NULL, ","))
             {              {
               fprintf (stderr, "-t step must be >= 1\n");                if (size_num == size_allocnum)
               exit (1);                  {
                     size_array = (struct size_array_t *)
                       __gmp_allocate_or_reallocate
                       (size_array,
                        size_allocnum * sizeof(size_array[0]),
                        (size_allocnum+10) * sizeof(size_array[0]));
                     size_allocnum += 10;
                   }
                 if (sscanf (s, "%ld-%ld",
                             &size_array[size_num].start,
                             &size_array[size_num].end) != 2)
                   {
                     size_array[size_num].start = size_array[size_num].end
                       = atol (s);
                   }
   
                 if (size_array[size_num].start < 0
                     || size_array[size_num].end < 0
                     || size_array[size_num].start > size_array[size_num].end)
                   {
                     fprintf (stderr, "invalid size parameter: %s\n", s);
                     exit (1);
                   }
   
                 size_num++;
             }              }
           break;  
         case 'u':  
           option_resource_usage = 1;  
           break;  
         case 'z':  
           sp.cache = 1;  
           break;  
         case OPTION_ALIGN:  
           abort();  
         case 'x':  
           sp.align_xp = atol (optarg);  
           break;  
         case 'y':  
           sp.align_yp = atol (optarg);  
           break;  
         case 'w':  
           sp.align_wp = atol (optarg);  
           break;  
         case 'W':  
           sp.align_wp2 = atol (optarg);  
           break;  
         case '?':  
           exit(1);  
         }          }
           break;
         case 't':
           option_step = atol (optarg);
           if (option_step < 1)
             {
               fprintf (stderr, "-t step must be >= 1\n");
               exit (1);
             }
           break;
         case 'u':
           option_resource_usage = 1;
           break;
         case 'z':
           sp.cache = 1;
           break;
         case 'x':
           sp.align_xp = atol (optarg);
           break;
         case 'y':
           sp.align_yp = atol (optarg);
           break;
         case 'w':
           sp.align_wp = atol (optarg);
           break;
         case 'W':
           sp.align_wp2 = atol (optarg);
           break;
         case '?':
           exit(1);
       }        }
       }
   
   if (optind >= argc)    if (optind >= argc)
     {      {
Line 919  main (int argc, char *argv[])
Line 1080  main (int argc, char *argv[])
       exit (1);        exit (1);
     }      }
   
   srand (option_seed);    gmp_randseed_ui (RANDS, option_seed);
   srandom (option_seed);  
   srand48 (option_seed);  
   
   choice = (struct choice_t *) (*_mp_allocate_func)    choice = (struct choice_t *) (*__gmp_allocate_func)
     ((argc - optind) * sizeof(choice[0]));      ((argc - optind) * sizeof(choice[0]));
   for ( ; optind < argc; optind++)    for ( ; optind < argc; optind++)
     {      {
Line 940  main (int argc, char *argv[])
Line 1099  main (int argc, char *argv[])
     }      }
   
   speed_time_init ();    speed_time_init ();
     if (option_unit == UNIT_CYCLES || option_unit == UNIT_CYCLESPERLIMB)
       speed_cycletime_need_cycles ();
     else
       speed_cycletime_need_seconds ();
   
   if ((option_unit == UNIT_CYCLES || option_unit == UNIT_CYCLESPERLIMB)  
       && speed_cycletime == 1.0)  
     {  
       fprintf (stderr, "Times in cycles requested, but CPU frequency unknown.\n");  
       fprintf (stderr, "Use environment variable GMP_CPU_FREQUENCY in Hertz, eg. 450e6\n");  
       exit (1);  
     }  
   
   if (option_gnuplot)    if (option_gnuplot)
     {      {
       run_gnuplot ();        run_gnuplot (argc, argv);
     }      }
   else    else
     {      {
Line 960  main (int argc, char *argv[])
Line 1115  main (int argc, char *argv[])
       else        else
         printf ("overhead %.2f cycles",          printf ("overhead %.2f cycles",
                 speed_measure (speed_noop, NULL) / speed_cycletime);                  speed_measure (speed_noop, NULL) / speed_cycletime);
       printf (", precision %d units of %.2e secs, cycle %.1e\n",        printf (", precision %d units of %.2e secs",
               speed_precision, speed_unittime, speed_cycletime);                speed_precision, speed_unittime);
   
         if (speed_cycletime == 1.0 || speed_cycletime == 0.0)
           printf (", CPU freq unknown\n");
         else
           printf (", CPU freq %.2f MHz\n", 1e-6/speed_cycletime);
   
       printf ("       ");        printf ("       ");
       for (i = 0; i < num_choices; i++)        for (i = 0; i < num_choices; i++)
Line 973  main (int argc, char *argv[])
Line 1133  main (int argc, char *argv[])
   
   if (option_resource_usage)    if (option_resource_usage)
     {      {
 #if defined(linux)  #if HAVE_GETRUSAGE
       /* This is Linux kernel specific. */        {
       char  buf[128];          /* This doesn't give data sizes on linux 2.0.x, only utime. */
       sprintf (buf, "cat /proc/%d/status", getpid());          struct rusage  r;
       system (buf);          if (getrusage (RUSAGE_SELF, &r) != 0)
             perror ("getrusage");
           else
             printf ("getrusage(): utime %ld.%06ld data %ld stack %ld maxresident %ld\n",
                     r.ru_utime.tv_sec, r.ru_utime.tv_usec,
                     r.ru_idrss, r.ru_isrss, r.ru_ixrss);
         }
 #else  #else
       /* This doesn't give data sizes on Linux 2.0.36, only utime. */        printf ("getrusage() not available\n");
       struct rusage  r;  
       if (getrusage (RUSAGE_SELF, &r) != 0)  
         perror ("getrusage");  
       else  
         printf ("utime %ld.%06ld data %ld stack %ld maxresident %ld\n",  
                 r.ru_utime.tv_sec, r.ru_utime.tv_usec,  
                 r.ru_idrss, r.ru_isrss, r.ru_ixrss);  
 #endif  #endif
   
         /* Linux kernel. */
         {
           char  buf[128];
           sprintf (buf, "/proc/%d/status", getpid());
           if (access (buf, R_OK) == 0)
             {
               sprintf (buf, "cat /proc/%d/status", getpid());
               system (buf);
             }
   
         }
     }      }
   
   return 0;    return 0;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.1.2
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>