Changeset 39248 for trunk/Ohana/src/addstar
- Timestamp:
- Dec 11, 2015, 9:39:10 AM (11 years ago)
- File:
-
- 1 edited
-
trunk/Ohana/src/addstar/src/mkcmf.c (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Ohana/src/addstar/src/mkcmf.c
r39242 r39248 280 280 } 281 281 282 gaussdev_init ();282 ohana_gaussdev_init (); 283 283 284 284 // load test stars from a file: … … 578 578 579 579 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); 583 583 flux = pow (10.0, -0.4*M[i]); 584 584 fSN = 1.0 / sqrt(flux); … … 624 624 625 625 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); 629 629 flux = pow (10.0, -0.4*M[i]); 630 630 fSN = 1.0 / sqrt(flux); … … 673 673 674 674 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); 678 678 flux = pow (10.0, -0.4*M[i]); 679 679 fSN = 1.0 / sqrt(flux); … … 724 724 725 725 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); 729 729 flux = pow (10.0, -0.4*M[i]); 730 730 fSN = 1.0 / sqrt(flux); … … 781 781 782 782 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); 786 786 flux = pow (10.0, -0.4*M[i]); 787 787 fSN = 1.0 / sqrt(flux); … … 842 842 843 843 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); 847 847 flux = pow (10.0, -0.4*M[i]); 848 848 fSN = 1.0 / sqrt(flux); … … 911 911 912 912 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); 916 916 flux = pow (10.0, -0.4*M[i]); 917 917 fSN = 1.0 / sqrt(flux); … … 1018 1018 1019 1019 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); 1023 1023 flux = pow (10.0, -0.4*M[i]); 1024 1024 fSN = 1.0 / sqrt(flux); … … 1122 1122 1123 1123 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); 1127 1127 flux = pow(10.0, -0.4 * M[i]); 1128 1128 fSN = 1.0 / sqrt(flux); … … 1248 1248 1249 1249 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); 1253 1253 flux = pow (10.0, -0.4*M[i]); 1254 1254 fSN = 1.0 / sqrt(flux);
Note:
See TracChangeset
for help on using the changeset viewer.
