cmbfast

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Anjan Ananda Sen
Posts: 13
Joined: March 30 2005
Affiliation: Center For Theoretical Physics, Jamia Millia Islamia, new delhi

cmbfast

Post by Anjan Ananda Sen » February 25 2009

Hi
I am trying to use the cmbfast with a file for the EOS. The EOS file I have have a running from 0.001 till 1 with uniform spacing. The total number of points is 200.

But while running cmbfast it is giving this problem:
computation stopped in dverk with the followingvalues -
ind = -3 tol = 1.061757D-17 x = 1.293044897493724D+04
n = 41 hmin = 1.000000D-08 xend = 1.303557128226869D+04
nw = 451 hmax = 2.000000D+00 prev xend = 1.290538910108542D+04
no of successful steps = 11144.
no of successive failures = 3.
no of function evals = 119743.
the components of y are

0.800999999999954D+00 -0.397046046650373D-06 -0.999996864413950D+00 0.179376777178244D-02 0.000000000000000D+00
0.179376776887656D-02 0.124489226448000D-16 0.238945873044038D-02 0.507752614395165D-09 -0.949709657993321D-03
-0.186344686744157D-04 -0.280098096395631D-06 -0.342631954883440D-08 -0.353701111065008D-10 -0.315903730160624D-12
-0.248632786879691D-14 -0.174904184160498D-16 -0.111219747672685D-18 -0.645215581984291D-21 -0.344045882871472D-23
-0.867046041114655D-07 -0.218905737320890D-08 -0.173887173109574D-07 -0.938976083632132D-09 -0.265646032158638D-10
-0.514712306451214D-12 -0.761548340398389D-14 -0.913150244596468D-16 -0.921423321015099D-18 -0.803092955548685D-20
-0.616269681196601D-22 -0.422514579680380D-24 -0.241556892268637D-26 0.238945479564562D-02 0.508213362135174D-09
-0.950385300510256D-03 -0.186723972585293D-04 -0.281201520128645D-06 -0.344821155016188D-08 -0.357009819144811D-10
-0.318417020989132D-12


Can somebody help me in this regard?

anjan sen

Ben Gold
Posts: 81
Joined: September 25 2004
Affiliation: University of Minnesota
Contact:

cmbfast

Post by Ben Gold » February 25 2009

What EOS are you trying to integrate and do you have dark energy perturbations turned on? Approaching w=-1 with perturbations on can cause problems, as the sound speed will diverge.

Anjan Ananda Sen
Posts: 13
Joined: March 30 2005
Affiliation: Center For Theoretical Physics, Jamia Millia Islamia, new delhi

cmbfast

Post by Anjan Ananda Sen » February 25 2009

the eos.dat file is perfectly fine

0.001 -1
0.005 -1.
0.01 -1.
0.015 -0.999999
+0.02 -0.999998
0.025 -0.999996
0.03 -0.999994
0.035 -0.99999
0.04 -0.999986
0.045 -0.999979
0.05 -0.999972
0.055 -0.999962
0.06 -0.999951
0.065 -0.999938
0.07 -0.999922
0.075 -0.999905
0.08 -0.999884
0.085 -0.999861
0.09 -0.999835
0.095 -0.999806
0.1 -0.999774
0.105 -0.999738
0.11 -0.999699
0.115 -0.999656
0.12 -0.99961
0.125 -0.999559
0.13 -0.999504
0.135 -0.999445
0.14 -0.999381
0.145 -0.999312
0.15 -0.999239
0.155 -0.999161
0.16 -0.999078
0.165 -0.998989
0.17 -0.998895
0.175 -0.998795
0.18 -0.99869
0.185 -0.998579
0.19 -0.998461
0.195 -0.998338
0.2 -0.998208
0.205 -0.998072
0.21 -0.997929
0.215 -0.99778
0.22 -0.997624
0.225 -0.997461
0.23 -0.997291
0.235 -0.997113
0.24 -0.996928
0.245 -0.996736
0.25 -0.996537
0.255 -0.99633
0.26 -0.996115
0.265 -0.995892
0.27 -0.995661
0.275 -0.995423
0.28 -0.995176
0.285 -0.994921
0.29 -0.994658
0.295 -0.994387
0.3 -0.994107
0.305 -0.993818
0.31 -0.993522
0.315 -0.993216
0.32 -0.992902
0.325 -0.99258
0.33 -0.992248
0.335 -0.991908
0.34 -0.99156
0.345 -0.991202
0.35 -0.990835
0.355 -0.99046
0.36 -0.990076
0.365 -0.989683
0.37 -0.989281
0.375 -0.98887
0.38 -0.988451
0.385 -0.988023
0.39 -0.987585
0.395 -0.987139
0.4 -0.986684
0.405 -0.986221
0.41 -0.985749
0.415 -0.985268
0.42 -0.984778
0.425 -0.98428
0.43 -0.983774
0.435 -0.983259
0.44 -0.982735
0.445 -0.982203
0.45 -0.981663
0.455 -0.981115
0.46 -0.980559
0.465 -0.979995
0.47 -0.979422
0.475 -0.978842
0.48 -0.978255
0.485 -0.977659
0.49 -0.977056
0.495 -0.976446
0.5 -0.975829
0.505 -0.975204
0.51 -0.974572
0.515 -0.973933
0.52 -0.973287
0.525 -0.972635
0.53 -0.971976
0.535 -0.971311
0.54 -0.970639
0.545 -0.969961
0.55 -0.969277
0.555 -0.968587
0.56 -0.967892
0.565 -0.967191
0.57 -0.966484
0.575 -0.965772
0.58 -0.965055
0.585 -0.964332
0.59 -0.963605
0.595 -0.962873
0.6 -0.962137
0.605 -0.961395
0.61 -0.96065
0.615 -0.959901
0.62 -0.959147
0.625 -0.95839
0.63 -0.957629
0.635 -0.956864
0.64 -0.956096
0.645 -0.955325
0.65 -0.954551
0.655 -0.953773
0.66 -0.952993
0.665 -0.95221
0.67 -0.951425
0.675 -0.950637
0.68 -0.949847
0.685 -0.949055
0.69 -0.948261
0.695 -0.947466
0.7 -0.946668
0.705 -0.945869
0.71 -0.945069
0.715 -0.944268
0.72 -0.943465
0.725 -0.942661
0.73 -0.941857
0.735 -0.941052
0.74 -0.940246
0.745 -0.939439
0.75 -0.938633
0.755 -0.937826
0.76 -0.937019
0.765 -0.936212
0.77 -0.935405
0.775 -0.934599
0.78 -0.933793
0.785 -0.932987
0.79 -0.932182
0.795 -0.931377
0.8 -0.930573
0.805 -0.92977
0.81 -0.928969
0.815 -0.928168
0.82 -0.927368
0.825 -0.92657
0.83 -0.925773
0.835 -0.924977
0.84 -0.924183
0.845 -0.923391
0.85 -0.9226
0.855 -0.921811
0.86 -0.921024
0.865 -0.920239
0.87 -0.919456
0.875 -0.918675
0.88 -0.917896
0.885 -0.91712
0.89 -0.916345
0.895 -0.915573
0.9 -0.914804
0.905 -0.914037
0.91 -0.913272
0.915 -0.91251
0.92 -0.911751
0.925 -0.910994
0.93 -0.91024
0.935 -0.909489
0.94 -0.908741
0.945 -0.907996
0.95 -0.907253
0.955 -0.906514
0.96 -0.905778
0.965 -0.905044
0.97 -0.904314
0.975 -0.903587
0.98 -0.902863
0.985 -0.902143
0.99 -0.901425
0.995 -0.900711
1. -0.9

I am not also using the perturbation for dark energy

Post Reply