IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Dec 9, 2015, 3:42:26 PM (11 years ago)
Author:
eugene
Message:

move multple copies of gaussian deviates to libohana

Location:
trunk/Ohana/src/opihi
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/Ohana/src/opihi/cmd.data/gaussdeviate.c

    r34088 r39242  
    1717  ResetVector (vec, OPIHI_FLT, Npts);
    1818
    19   gauss_init (2048);
     19  gaussdev_init ();
    2020  for (i = 0; i < Npts; i++) {
    21     vec[0].elements.Flt[i] = rnd_gauss (mean, sigma);
     21    vec[0].elements.Flt[i] = gaussdev_rnd (mean, sigma);
    2222  }
    2323  return (TRUE);
     
    2828   
    2929}
    30 
    31 double int_gauss (int i);
    3230
    3331int gaussintegral (int argc, char **argv) {
     
    4644  ResetVector (vec, OPIHI_FLT, Npts);
    4745
    48   gauss_init (Npts);
     46  gaussdev_init ();
    4947  for (i = 0; i < Npts; i++) {
    50     vec[0].elements.Flt[i] = int_gauss (i);
     48    vec[0].elements.Flt[i] = gaussian_int (i);
    5149  }
    5250  return (TRUE);
  • trunk/Ohana/src/opihi/lib.data/Makefile

    r36679 r39242  
    3030$(SDIR)/starfuncs.$(ARCH).o             \
    3131$(SDIR)/hermitian.$(ARCH).o             \
    32 $(SDIR)/gaussian.$(ARCH).o              \
    3332$(SDIR)/graphtools.$(ARCH).o            \
    3433$(SDIR)/queues.$(ARCH).o                \
     
    3635$(SDIR)/open_kapa.$(ARCH).o             \
    3736$(SDIR)/style_args.$(ARCH).o
     37
     38# moved to libohana
     39# $(SDIR)/gaussian.$(ARCH).o            \
    3840
    3941# fix malloc
  • trunk/Ohana/src/opihi/mana/simsignal.c

    r7917 r39242  
    66  float *buf;
    77  double cvalue, dvalue, sigma, SN, period;
    8   double rnd_gauss ();
    98
    109  if (argc != 5) {
     
    2827    ivalue = scale * 0.5 * (sin (i*2*M_PI/period) + 1) + 0.5;
    2928    /*
    30     dvalue = rnd_gauss (cvalue, sigma);
     29    dvalue = gaussdev_rnd (cvalue, sigma);
    3130    cvalue = (dvalue + range) / (2.0*range);
    3231    dvalue = MAX (0, MIN (0.99999, cvalue));
     
    44438 bit = 2^8
    4544
    46   gauss_init (2*scale);
     45  gaussdev_init ();
    4746  sigma = 2.0 / SN;
    4847  range = 1 + 5*sigma;
Note: See TracChangeset for help on using the changeset viewer.