Changeset 37481 for branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/fakeastro_galaxy.c

Timestamp:

Oct 10, 2014, 3:56:43 PM (12 years ago)

Author:

eugene

Message:

adding fake qsos

File:

• branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/fakeastro_galaxy.c (modified) (1 diff)

branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/fakeastro_galaxy.c

@@ -36 +36 @@
     FakeAstro_Stars *stars = make_fakestars (Nstars);
 
-    FakeAstro_Stars *qso_icrf = make_fakeqsos (Nstars);
+    // only generate the QSOs 1x:
+    if (n == 0) {
+      FakeAstro_Stars *qso_icrf = make_fakeqsos (FAKEASTRO_NQSO_ICRF, TRUE); // quasars with known positions
+      REALLOCATE (stars, FakeAstro_Stars, Nstars + FAKEASTRO_NQSO_ICRF);
+      for (i = 0; i < FAKEASTRO_NQSO_ICRF; i++) {
+        stars[Nstars+i] = qso_icrf[i];
+      }
+      Nstars += FAKEASTRO_NQSO_ICRF;
+      FakeAstro_Stars *qso_zero = make_fakeqsos (FAKEASTRO_NQSO_ZERO, FALSE); // quasars with unknown positions, but zero velocity
+      REALLOCATE (stars, FakeAstro_Stars, Nstars + FAKEASTRO_NQSO_ZERO);
+      for (i = 0; i < FAKEASTRO_NQSO_ZERO; i++) {
+        stars[Nstars+i] = qso_zero[i];
+      }
+      Nstars += FAKEASTRO_NQSO_ZERO;
+    }
 
     MARKTIME ("generate %d fake stars in %f sec\n", Nstars, dtime);