IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 39248 for trunk


Ignore:
Timestamp:
Dec 11, 2015, 9:39:10 AM (11 years ago)
Author:
eugene
Message:

use ohana_ namespace for gaussdev related functions

Location:
trunk/Ohana/src
Files:
18 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/addstar/src/mkcmf.c

    r39242 r39248  
    280280  }
    281281   
    282   gaussdev_init ();
     282  ohana_gaussdev_init ();
    283283
    284284  // load test stars from a file:
    … …  
    578578
    579579    if (ADDNOISE) {
    580       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    581       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    582       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     580      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     581      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     582      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    583583      flux = pow (10.0, -0.4*M[i]);
    584584      fSN = 1.0 / sqrt(flux);
    … …  
    624624
    625625    if (ADDNOISE) {
    626       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    627       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    628       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     626      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     627      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     628      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    629629      flux = pow (10.0, -0.4*M[i]);
    630630      fSN = 1.0 / sqrt(flux);
    … …  
    673673
    674674    if (ADDNOISE) {
    675       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    676       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    677       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     675      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     676      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     677      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    678678      flux = pow (10.0, -0.4*M[i]);
    679679      fSN = 1.0 / sqrt(flux);
    … …  
    724724
    725725    if (ADDNOISE) {
    726       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    727       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    728       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     726      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     727      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     728      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    729729      flux = pow (10.0, -0.4*M[i]);
    730730      fSN = 1.0 / sqrt(flux);
    … …  
    781781
    782782    if (ADDNOISE) {
    783       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    784       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    785       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     783      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     784      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     785      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    786786      flux = pow (10.0, -0.4*M[i]);
    787787      fSN = 1.0 / sqrt(flux);
    … …  
    842842
    843843    if (ADDNOISE) {
    844       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    845       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    846       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     844      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     845      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     846      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    847847      flux = pow (10.0, -0.4*M[i]);
    848848      fSN = 1.0 / sqrt(flux);
    … …  
    911911
    912912    if (ADDNOISE) {
    913       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    914       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    915       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     913      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     914      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     915      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    916916      flux = pow (10.0, -0.4*M[i]);
    917917      fSN = 1.0 / sqrt(flux);
    … …  
    10181018
    10191019    if (ADDNOISE) {
    1020       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    1021       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    1022       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     1020      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     1021      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     1022      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    10231023      flux = pow (10.0, -0.4*M[i]);
    10241024      fSN = 1.0 / sqrt(flux);
    … …  
    11221122
    11231123    if (ADDNOISE) {
    1124       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    1125       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    1126       M[i] += fSN * gaussdev_rnd(0.0, 1.0);
     1124      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     1125      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     1126      M[i] += fSN * ohana_gaussdev_rnd(0.0, 1.0);
    11271127      flux = pow(10.0, -0.4 * M[i]);
    11281128      fSN = 1.0 / sqrt(flux);
    … …  
    12481248
    12491249    if (ADDNOISE) {
    1250       X[i] += FX * fSN * gaussdev_rnd(0.0, 1.0);
    1251       Y[i] += FY * fSN * gaussdev_rnd(0.0, 1.0);
    1252       M[i] += fSN*gaussdev_rnd(0.0, 1.0);
     1250      X[i] += FX * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     1251      Y[i] += FY * fSN * ohana_gaussdev_rnd(0.0, 1.0);
     1252      M[i] += fSN*ohana_gaussdev_rnd(0.0, 1.0);
    12531253      flux = pow (10.0, -0.4*M[i]);
    12541254      fSN = 1.0 / sqrt(flux);

  • trunk/Ohana/src/fakeastro/src/fakeastro.c

    r39244 r39248  
    33int main (int argc, char **argv) {
    44
    5   gaussdev_init ();
     5  ohana_gaussdev_init ();
    66
    77  /* get configuration info, args */

  • trunk/Ohana/src/fakeastro/src/make_2mass_measures.c

    r39242 r39248  
    8080
    8181    // uR,uD in linear arcsec
    82     double dRsee = gaussdev_rnd (0.0, 1.0 / SN);
    83     double dDsee = gaussdev_rnd (0.0, 1.0 / SN);
     82    double dRsee = ohana_gaussdev_rnd (0.0, 1.0 / SN);
     83    double dDsee = ohana_gaussdev_rnd (0.0, 1.0 / SN);
    8484
    8585    double dRoff = (dRpm + dRsee) / 3600.0;

  • trunk/Ohana/src/fakeastro/src/make_fake_stars_catalog.c

    r39242 r39248  
    113113
    114114    // uR,uD in linear arcsec
    115     double dRsee = gaussdev_rnd (0.0, 1.0 / SN);
    116     double dDsee = gaussdev_rnd (0.0, 1.0 / SN);
     115    double dRsee = ohana_gaussdev_rnd (0.0, 1.0 / SN);
     116    double dDsee = ohana_gaussdev_rnd (0.0, 1.0 / SN);
    117117
    118118    double dRoff = (dRpm + dRsee) / 3600.0;

  • trunk/Ohana/src/fakeastro/src/make_fakeqsos.c

    r39242 r39248  
    5151
    5252    // Mr will be interpretted as m_r
    53     double Mr = gaussdev_rnd (18.0, 1.5);
     53    double Mr = ohana_gaussdev_rnd (18.0, 1.5);
    5454   
    5555    stars[i].R = R;

  • trunk/Ohana/src/fakeastro/src/make_fakestars.c

    r39242 r39248  
    4848      int inPatch = FALSE;
    4949      while (!inPatch) {
    50         z = gaussdev_rnd (0.0, FAKEASTRO_ZGAL);
     50        z = ohana_gaussdev_rnd (0.0, FAKEASTRO_ZGAL);
    5151        r = sqrt(drand48()) * FAKEASTRO_RGAL;
    5252        Lrad = drand48() * 2 * M_PI;
    … …  
    9898    double Mr;
    9999    if (bigPeak) {
    100       Mr = gaussdev_rnd (11.25, 1.0);
     100      Mr = ohana_gaussdev_rnd (11.25, 1.0);
    101101    } else {
    102       Mr = gaussdev_rnd (10.00, 2.5);
     102      Mr = ohana_gaussdev_rnd (10.00, 2.5);
    103103    }
    104104   

  • trunk/Ohana/src/libohana/include/ohana.h

    r39242 r39248  
    487487
    488488/* in gaussdev.c */
    489 double gaussian (double x, double mean, double sigma);
    490 void gaussdev_init (void);
    491 void gaussdev_free (void);
    492 double gaussdev_rnd (double mean, double sigma);
    493 
    494 # endif
     489double ohana_gaussian (double x, double mean, double sigma);
     490void ohana_gaussdev_init (void);
     491void ohana_gaussdev_free (void);
     492double ohana_gaussdev_rnd (double mean, double sigma);
     493
     494# endif

  • trunk/Ohana/src/libohana/src/gaussdev.c

    r39242 r39248  
    88
    99// return the value of a gaussian at position x with mean,sigma
    10 double gaussian (double x, double mean, double sigma) {
     10double ohana_gaussian (double x, double mean, double sigma) {
    1111
    1212  double f;
    … …  
    2121/* integrate a gaussian from -10 sigma to +10 sigma */
    2222// how long does this take with a reasonably high resolution?
    23 void gaussdev_init (void) {
     23void ohana_gaussdev_init (void) {
    2424 
    2525  int i;
    … …  
    4242 
    4343  for (i = 0, x = -10.0; (i < Ngaussint) && (x <= 10.0); x += dx)  {
    44     df = (3.0*gaussian(x    , mean, sigma) +
    45           9.0*gaussian(x+dx1, mean, sigma) +
    46           9.0*gaussian(x+dx2, mean, sigma) +
    47           3.0*gaussian(x+dx3, mean, sigma)) * (dx1/8.0);
     44    df = (3.0*ohana_gaussian(x    , mean, sigma) +
     45          9.0*ohana_gaussian(x+dx1, mean, sigma) +
     46          9.0*ohana_gaussian(x+dx2, mean, sigma) +
     47          3.0*ohana_gaussian(x+dx3, mean, sigma)) * (dx1/8.0);
    4848    val += df;
    4949    if (val > (i + 0.5) / (double) Ngaussint) {
    … …  
    5555}
    5656
    57 void gaussdev_free (void) {
     57void ohana_gaussdev_free (void) {
    5858  free (gaussint);
    5959}
    … …  
    6161// return a number drawn from a gaussian deviate with mean, sigma
    6262// must first call gaussdev_init()
    63 double gaussdev_rnd (double mean, double sigma) {
     63double ohana_gaussdev_rnd (double mean, double sigma) {
    6464 
    6565  int i;
    … …  
    7777// return the value of \int gaussian(x,mean,sigma)
    7878// double gaussian_integral (double x, double mean, double sigma) {
    79 double gaussian_integral (int i) {
     79double ohana_gaussian_integral (int i) {
    8080
    8181  // int i = (x - mean) / (sigma * Ngaussint);

  • trunk/Ohana/src/misc/src/fakedist.c

    r39242 r39248  
    3939  init_outmatrix ();
    4040 
    41   gaussdev_init ();
     41  ohana_gaussdev_init ();
    4242
    4343  read_datafiles ();
    … …  
    343343        if (noise < dMo) {
    344344          *Noise = noise;
    345           *mag = gaussdev_rnd (v, noise);
    346           *color = gaussdev_rnd (uv, 1.4*noise);
     345          *mag = ohana_gaussdev_rnd (v, noise);
     346          *color = ohana_gaussdev_rnd (uv, 1.4*noise);
    347347          return;
    348348        }
    … …  
    436436      Ngood = 1.0;
    437437      for (k = 0; k < NTRY; k++) {
    438         v = gaussdev_rnd (mag, noise);
    439         uv = gaussdev_rnd (color, 1.4*noise);
     438        v = ohana_gaussdev_rnd (mag, noise);
     439        uv = ohana_gaussdev_rnd (color, 1.4*noise);
    440440        x = (uv - UV0 - 0.7*A_V) / DUV;
    441441        y = (v - V0 - Dist - A_V) / DV;

  • trunk/Ohana/src/misc/src/fakepop.c

    r39242 r39248  
    3636  gfits_scan (&UV_h, "DEC_Y", "%lf", 1, &ldM);
    3737 
    38   gaussdev_init ();
     38  ohana_gaussdev_init ();
    3939
    4040  fprintf (stderr, "beginning main loop\n");
    … …  
    5757        v = gfits_get_matrix_value (&V_i, X, Y) + d + Av;
    5858        noise = dVo*sqrt(1.0 + pow (10.0, (0.4*(v - dVref))));
    59         V = gaussdev_rnd (v, noise);
    60         UV = gaussdev_rnd (uv, 1.4*noise);
     59        V = ohana_gaussdev_rnd (v, noise);
     60        UV = ohana_gaussdev_rnd (uv, 1.4*noise);
    6161        if (noise < dMo) {
    6262          fprintf (stdout, "%f %f   %f %f   %f  %f\n", V, noise, V+UV, noise, mass, age);

  • trunk/Ohana/src/misc/src/fakestars.c

    r39242 r39248  
    4141  offset = atof (argv[4]);
    4242
    43   gaussdev_init ();
     43  ohana_gaussdev_init ();
    4444
    4545  /* load header, open file */
    … …  
    9494    dMr = 0.001 * stars[i].dM;
    9595    dMo = hypot (dMs, dMr);
    96     dM  = gaussdev_rnd (0.0, dMo);
     96    dM  = ohana_gaussdev_rnd (0.0, dMo);
    9797    stars[i].M += dM + offset;
    9898  }

  • trunk/Ohana/src/misc/src/fitdist.c

    r39242 r39248  
    3939  fprintf (stderr, "read in parameters\n");
    4040
    41   gaussdev_init ();
     41  ohana_gaussdev_init ();
    4242
    4343  read_datafiles ();
    … …  
    384384      Ngood = 1.0;
    385385      for (k = 0; k < NTRY; k++) {
    386         v = gaussdev_rnd (Mag, Noise);
    387         uv = gaussdev_rnd (Color, 1.4*Noise);
     386        v = ohana_gaussdev_rnd (Mag, Noise);
     387        uv = ohana_gaussdev_rnd (Color, 1.4*Noise);
    388388        x = (uv - UV0 - 0.7*A_V) / DUV;
    389389        y = (v - V0 - Dist - A_V) / DV;

  • trunk/Ohana/src/misc/src/magtoage.c

    r39242 r39248  
    2121  long A, B;
    2222 
    23   gaussdev_init();
     23  ohana_gaussdev_init();
    2424
    2525  lAo =   0.0;
    … …  
    112112      }
    113113      else {
    114         v = gaussdev_rnd (V, dV);
    115         uv = gaussdev_rnd ((U-V), dUV);
     114        v = ohana_gaussdev_rnd (V, dV);
     115        uv = ohana_gaussdev_rnd ((U-V), dUV);
    116116      }
    117117      x = (uv - UV0 - 0.7*Av) / DUV;

  • trunk/Ohana/src/misc/src/magtomass.c

    r39242 r39248  
    4141  fprintf (stderr, "using mags in columns %d & %d\n", col1, col2);
    4242
    43   gaussdev_init ();
     43  ohana_gaussdev_init ();
    4444
    4545  lAo = 0.1;
    … …  
    136136        Ngood = 1.0;
    137137        for (k = 0; k < Ntry; k++) {
    138           v = gaussdev_rnd (V, dV);
    139           uv = gaussdev_rnd ((U-V), dUV);
     138          v = ohana_gaussdev_rnd (V, dV);
     139          uv = ohana_gaussdev_rnd ((U-V), dUV);
    140140          x = (uv - UV0 - 0.7*Av) / DUV;
    141141          y = (v - d - V0 - Av) / DV;

  • trunk/Ohana/src/misc/src/mkgauss.c

    r39242 r39248  
    1515  Npts = atoi (argv[3]);
    1616
    17   gaussdev_init ();
     17  ohana_gaussdev_init ();
    1818
    1919  for (i = 0; i < Npts; i++) {
    20     f = gaussdev_rnd (mean, sigma);
     20    f = ohana_gaussdev_rnd (mean, sigma);
    2121    fprintf (stdout, "%f\n", f);
    2222  }

  • trunk/Ohana/src/opihi/cmd.data/gaussdeviate.c

    r39245 r39248  
    1717  ResetVector (vec, OPIHI_FLT, Npts);
    1818
    19   gaussdev_init ();
     19  ohana_gaussdev_init ();
    2020  for (i = 0; i < Npts; i++) {
    21     vec[0].elements.Flt[i] = gaussdev_rnd (mean, sigma);
     21    vec[0].elements.Flt[i] = ohana_gaussdev_rnd (mean, sigma);
    2222  }
    2323  return (TRUE);
    … …  
    4848  ResetVector (vec, OPIHI_FLT, Npts);
    4949
    50   gaussdev_init ();
     50  ohana_gaussdev_init ();
    5151  for (i = 0; i < Npts; i++) {
    5252    vec[0].elements.Flt[i] = gaussian_int (i);

  • trunk/Ohana/src/opihi/mana/simsignal.c

    r39242 r39248  
    2727    ivalue = scale * 0.5 * (sin (i*2*M_PI/period) + 1) + 0.5;
    2828    /*
    29     dvalue = gaussdev_rnd (cvalue, sigma);
     29    dvalue = ohana_gaussdev_rnd (cvalue, sigma);
    3030    cvalue = (dvalue + range) / (2.0*range);
    3131    dvalue = MAX (0, MIN (0.99999, cvalue));
    … …  
    43438 bit = 2^8
    4444
    45   gaussdev_init ();
     45  ohana_gaussdev_init ();
    4646  sigma = 2.0 / SN;
    4747  range = 1 + 5*sigma;

  • trunk/Ohana/src/relastro/src/fitpm.c

    r39246 r39248  
    6565  }
    6666
    67   gaussdev_init ();
     67  ohana_gaussdev_init ();
    6868
    6969  FitStats *fitStats = FitStatsInit (N_POINTS + N_OUTLIERS, N_BOOTSTRAP);
    … …  
    154154  }
    155155  // return exit_status();
    156   gaussdev_free();
     156  ohana_gaussdev_free();
    157157
    158158  exit (0);
    … …  
    176176    ParFactor (&pR, &pD, Ro, Do, points[i].T);
    177177
    178     double dR = gaussdev_rnd(0.0, POS_ERROR);
    179     double dD = gaussdev_rnd(0.0, POS_ERROR);
     178    double dR = ohana_gaussdev_rnd(0.0, POS_ERROR);
     179    double dD = ohana_gaussdev_rnd(0.0, POS_ERROR);
    180180
    181181    points[i].R = Ro + (dR + plx*pR + uR*(points[i].T - MJD_REF_YRS))/3600/dcos(Do);
Note: See TracChangeset for help on using the changeset viewer.