Changeset 4047 in Sophya for trunk/Cosmo/SimLSS


Ignore:
Timestamp:
Jan 12, 2012, 10:37:47 PM (14 years ago)
Author:
cmv
Message:

details..., cmv 12/01/2012

Location:
trunk/Cosmo/SimLSS
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/Cosmo/SimLSS/Makefile

    r3803 r4047  
    1111
    1212MYEXTINC = ${EXTLIBDIR}/${MACH}/include
    13 MYLIB  = $(SOPHYAEXTSLBLIST) -L$(LIB) -lcmvgenfluc -lcmvsimbao -lfftw3_threads -lfftw3 -lm
    14 MYLIB4 = $(SOPHYAEXTSLBLIST) -L$(LIB) -lcmvgenfluc4 -lcmvsimbao -lfftw3f_threads -lfftw3f -lm
     13MYLIB  = -L$(LIB) -lcmvgenfluc  -lcmvsimbao $(SOPHYAEXTSLBLIST) -lfftw3_threads -lfftw3 -lm
     14MYLIB4 = -L$(LIB) -lcmvgenfluc4 -lcmvsimbao $(SOPHYAEXTSLBLIST) -lfftw3f_threads -lfftw3f -lm
    1515
    1616#--------------------------------------------------------------------------
    … …  
    8383        rm -f $(LIBG) $(LIBGOBJ)
    8484        rm -f $(LIBG4) $(LIBGOBJ4)
    85         rm -rf $(OBJ)/CmvBAO_cxxrep/
    8685        rm -f $(PROGSTEST) $(PROGSTESTOBJ)
    8786        rm -f $(SPROGS) $(SPROGSOBJ)

  • trunk/Cosmo/SimLSS/constcosmo.h

    r3193 r4047  
    1919static const double ProtonMass_Cst = 1.67262171e-24; // 1 "proton mass" in "g"
    2020static const double Deg2Rad_Cst = (M_PI/180.); // 1 deg into rad
     21// 4 PI sr = 180/Pi * 2 * 360 = 360^2/Pi  -> deg^2 * 4*(Pi/360)^2 = sr
     22static const double Deg2toSr_Cst = M_PI*M_PI/32400.; // 1 deg^2 into steradian: 4*(Pi/360)^2
    2123
    2224static const double T_CMB_Par = 2.725;  // temperature du CMB en "K"

  • trunk/Cosmo/SimLSS/cosmocalc.cc

    r3970 r4047  
    235235  printf("Otot=%g  Ocurv=%g\n",Otot(z),Ocurv(z));
    236236  if(z <= 0.) return;
    237   printf("Distance comoving: los=%g Mpc, transv=%g Mpc\n",Dloscom(z),Dtrcom(z));
     237  printf("Distance comoving: los=%g Mpc (dlos/dz=%g), transv=%g Mpc\n",Dloscom(z),DloscomDz(z),Dtrcom(z));
    238238  printf("Distance: angular=%g Mpc, lum=%g Mpc\n",Dang(z),Dlum(z));
    239239  printf("Volume comoving element: %g Mpc^3/sr/dz=1\n",dVol(z));
    … …  
    292292{
    293293  return _Dhubble * NInteg(z);
     294}
     295
     296double CosmoCalc::DloscomDz(double z)
     297// d(Dloscom)/dz : derivee de la distance LOS comobile par rapport a z
     298{
     299  return SpeedOfLight_Cst / H(z);
    294300}
    295301

  • trunk/Cosmo/SimLSS/cosmocalc.h

    r3952 r4047  
    5353  double Dtrcom(double z);
    5454  double Dloscom(double z);
     55  double DloscomDz(double z);
    5556  double Dang(double z);
    5657  double Dlum(double z);

  • trunk/Cosmo/SimLSS/genefluct3d.cc

    r4008 r4047  
     1//#include "sopnamsp.h"
    12#include "machdefs.h"
    23#include <iostream>
Note: See TracChangeset for help on using the changeset viewer.