Changeset 3768 in Sophya for trunk/Cosmo/SimLSS/pkspectrum.cc

Timestamp:

May 3, 2010, 4:08:34 PM (15 years ago)

Author:

cmv

Message:

• refonte du code pour creer uniquement des conditions initiales
• introduction du tirage des vitesse LOS pour les redshift-distortion

cmv 03/05/2010

File:

• trunk/Cosmo/SimLSS/pkspectrum.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/pkspectrum.cc

@@ -399 +399 @@
 }
 
+double GrowthFactor::DsDz(double z,double dzinc)
+// y-y0 = a*(x-x0)^2 + b*(x-x0)
+// dy = a*dx^2 + b*dx   avec  dx = x-x0
+// dy'(dx) = 2*a*dx + b -> pour x=x0 on a dy'(0) = b
+{
+ if(z<0. || dzinc<=0.) {
+   cout<<"GrowthFactor::DsDz: z<0 or dzinc<=0. !"<<endl;
+   throw ParmError("GrowthFactor::DsDz: z<0 or dzinc<=0. !");
+ }
+
+ double z1, z2;
+ if(z>dzinc/2.) {
+   // cas ou z est suffisamment loin de zero
+   // resolution avec 2 points encadrant x0=z
+   z1 = z - dzinc; if(z1<0.) z1 = 0.;
+   z2 = z + dzinc;
+ } else {
+   // cas ou z est proche de zero
+   // resolution avec 2 points au dessus, x0=z1
+   z1 = z + dzinc;
+   z2 = z + 2.*dzinc;
+ }
+
+ double gz  = (*this)(z);
+ double dgz1 = (*this)(z1) - gz;
+ double dgz2 = (*this)(z2) - gz;
+
+ z1 -= z;
+ z2 -= z;
+ return (dgz2*z1*z1 - dgz1*z2*z2)/(z1*z2*(z1-z2));
+}
+
 void GrowthFactor::SetParTo(double OmegaMatter0,double OmegaLambda0)
 {