Changeset 610 in Sophya for trunk/SophyaLib/SkyT/nupower.cc

Timestamp:

Nov 22, 1999, 12:25:47 AM (26 years ago)

Author:

ansari

Message:

Corrections diverses - programme tgsky.cc (generation aleatoire de ciel)

et skymixer.cc completee - Reza 21/11/99

File:

• trunk/SophyaLib/SkyT/nupower.cc (modified) (5 diffs)

trunk/SophyaLib/SkyT/nupower.cc

@@ -1 +1 @@
 //--------------------------------------------------------------------------
 // File and Version Information:
-//      $Id: nupower.cc,v 1.2 1999-11-20 21:00:49 ansari Exp $
+//      $Id: nupower.cc,v 1.3 1999-11-21 23:25:45 ansari Exp $
 //
 // Description:
@@ -24 +24 @@
 // Constructor --
 //----------------
-PowerLawSpectra::PowerLawSpectra(double a, double b, double nu0, double numin, double numax)
+PowerLawSpectra::PowerLawSpectra(double a, double b, double nu0, double dnu, double numin, double numax)
   : RadSpectra(numin, numax)
 {
@@ -30 +30 @@
   _b = b;
   _nu0 = nu0;
+  _dnu = (dnu > 1.e-19) ? dnu : 1.;
 }
 
@@ -40 +41 @@
 PowerLawSpectra::flux(double nu) const
 {
-  if ((nu < _numin) || (nu > _numax)) return(0.);
-  else return(  _a * pow((nu-_nu0), _b) );
+  if ((nu< _numin) || (nu > _numax)) return(0.);
+  double tmp = (nu-_nu0)/_dnu;
+  if (tmp <= 0.)  return(0.); 
+  else return(  _a * pow(tmp, _b) );
 }
 
@@ -47 +50 @@
 PowerLawSpectra::Print(ostream& os) const
 {
-  os << "PowerLawSpectra::Print f = a (nu-nu0)^b " << _a << "(nu-" << _nu0 
-     << ") ^ " << _b << endl;
+  os << "PowerLawSpectra::Print f = a ((nu-nu0)/dnu)^b " << _a << "((nu-" << _nu0 
+     << ") / " << _dnu << ") ^ " << _b << endl;
   os << " - Fmin,Fmax= " << minFreq() << "," << maxFreq() << endl;
   os << "MeanFreq= " << meanFreq() << "  Emission= " << flux(meanFreq()) << endl;