Changeset 3662 in Sophya for trunk/Cosmo/SimLSS
- Timestamp:
- Oct 23, 2009, 10:29:14 PM (16 years ago)
- Location:
- trunk/Cosmo/SimLSS
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Cosmo/SimLSS/genefluct3d.cc
r3617 r3662 1435 1435 // nlobe ~ 1/Dtrcom^2 1436 1436 // Mlim ~ sqrt(DNu) * (Dlum)^2 / Dtrcom 1437 if(cosmo_ == NULL || redsh_ref_<0. | loscom2zred_.size()<1) {1437 if(cosmo_ == NULL || redsh_ref_<0. || loscom2zred_.size()<1) { 1438 1438 const char *bla = "GeneFluct3D::AddNoise2Real_Error: set Observator and Cosmology first"; 1439 1439 cout<<bla<<endl; throw ParmError(bla); -
trunk/Cosmo/SimLSS/genefluct3d.h
r3619 r3662 155 155 double compute_power_carte(void); 156 156 void check_array_alloc(void); 157 inline double pixelfilter(double x) 157 inline double pixelfilter(double x) // ATTENTION: seulement pour x>0 158 158 {return (x<0.025) ? 1.-x*x/6.*(1.-x*x/20.): sin(x)/x;} 159 159 -
trunk/Cosmo/SimLSS/pkspectrum.cc
r3381 r3662 545 545 // pour x<0.01 |DL(x)-G(X)^2|<2.0e-13 546 546 if(typfilter_ == GAUSSIAN) 547 if(x<0.01) return 1.-x2*(1.-x2/2.); else return exp(-x2);547 {if(x<0.01) return 1.-x2*(1.-x2/2.); else return exp(-x2);} 548 548 549 549 // Filtre top-hat T(x) = 3*(sin(x)-x*cos(x))/x^3
Note:
See TracChangeset
for help on using the changeset viewer.