CosmoCoffee

[Ido Ben-Dayan]

I'm trying to run the new, July 2011 version of cosmomc.

I've compiled it successfully including wmap7 data.
Trying to run the params.ini file gives the following memory bug:

Random seeds: 19770, 11114 rand_inst: 0
Computing tensors: F
Doing CMB lensing: T
lmax = 2100
Number of C_ls = 4
Varying 7 parameters ( 0 fast)
starting Monte-Carlo
Initialising BBN Helium data...
reading WMAP7 data
Initializing WMAP likelihood, version v4.1
*** glibc detected *** ./cosmomc: malloc(): memory corruption: 0x0b153260 ***
======= Backtrace: =========
/lib/libc.so.6[0xb6d7850b]
/lib/libc.so.6[0xb6d7b400]
/lib/libc.so.6(__libc_malloc+0x6a)[0xb6d7d0ba]
/usr/lib/libgfortran.so.3[0xb6ec617b]
/usr/lib/libgfortran.so.3[0xb6f50662]
/usr/lib/libgfortran.so.3[0xb6f4f44c]
./cosmomc[0x80d05dc]
./cosmomc[0x806ad06]
./cosmomc[0x80a114e]
./cosmomc[0x80aaea2]
./cosmomc[0x80b13d0]
./cosmomc[0x823c1fb]
/lib/libc.so.6(__libc_start_main+0xfe)[0xb6d21ace]
./cosmomc[0x804cde1]
======= Memory map: ========
08048000−08275000 r-xp 00000000 08:03 22268 /home/idobd/bin/cosmomcJuly2011/cosmomc/cosmomc
08275000−08276000 r–p 0022d000 08:03 22268 /home/idobd/bin/cosmomcJuly2011/cosmomc/cosmomc
08276000−08279000 rw-p 0022e000 08:03 22268 /home/idobd/bin/cosmomcJuly2011/cosmomc/cosmomc
08279000−0b224000 rw-p 00000000 00:00 0 [heap]
9de88000-a69b2000 rw-p 00000000 00:00 0
a6e8e000-a7136000 rw-p 00000000 00:00 0
a72b3000-b5de4000 rw-p 00000000 00:00 0
b6200000-b6221000 rw-p 00000000 00:00 0
b6221000-b6300000 —p 00000000 00:00 0
b635f000-b6c62000 rw-p 00000000 00:00 0
b6c62000-b6d09000 r-xp 00000000 08:02 127916 /usr/lib/libblas.so.3.2.1
b6d09000-b6d0a000 r–p 000a6000 08:02 127916 /usr/lib/libblas.so.3.2.1
b6d0a000-b6d0b000 rw-p 000a7000 08:02 127916 /usr/lib/libblas.so.3.2.1
b6d0b000-b6e64000 r-xp 00000000 08:02 31793 /lib/libc−2.10.1.so
b6e64000-b6e65000 —p 00159000 08:02 31793 /lib/libc−2.10.1.so
b6e65000-b6e67000 r–p 00159000 08:02 31793 /lib/libc−2.10.1.so
b6e67000-b6e68000 rw-p 0015b000 08:02 31793 /lib/libc−2.10.1.so
b6e68000-b6e6c000 rw-p 00000000 00:00 0
b6e6c000-b6e8a000 r-xp 00000000 08:02 1088 /lib/libgcc_s.so.1
b6e8a000-b6e8b000 r–p 0001d000 08:02 1088 /lib/libgcc_s.so.1
b6e8b000-b6e8c000 rw-p 0001e000 08:02 1088 /lib/libgcc_s.so.1
b6e8c000-b6eb3000 r-xp 00000000 08:02 31806 /lib/libm−2.10.1.so
b6eb3000-b6eb4000 r–p 00027000 08:02 31806 /lib/libm−2.10.1.so
b6eb4000-b6eb5000 rw-p 00028000 08:02 31806 /lib/libm−2.10.1.so
b6eb5000-b6f7b000 r-xp 00000000 08:02 126055 /usr/lib/libgfortran.so.3.0.0
b6f7b000-b6f7c000 r–p 000c5000 08:02 126055 /usr/lib/libgfortran.so.3.0.0
b6f7c000-b6f7d000 rw-p 000c6000 08:02 126055 /usr/lib/libgfortran.so.3.0.0
b6f7d000-b6f7e000 rw-p 00000000 00:00 0
b6f7e000-b7729000 r-xp 00000000 08:02 129694 /usr/lib/liblapack.so.3.2.1
b7729000-b772a000 r–p 007aa000 08:02 129694 /usr/lib/liblapack.so.3.2.1
b772a000-b772d000 rw-p 007ab000 08:02 129694 /usr/lib/liblapack.so.3.2.1
b772d000-b783a000 rw-p 00000000 00:00 0
b784f000-b7850000 rw-p 00000000 00:00 0
b7850000-b786e000 r-xp 00000000 08:02 370 /lib/ld−2.10.1.so
b786e000-b786f000 r–p 0001d000 08:02 370 /lib/ld−2.10.1.so
b786f000-b7870000 rw-p 0001e000 08:02 370 /lib/ld−2.10.1.so
bf82b000-bf882000 rw-p 00000000 00:00 0 [stack]
ffffe000-fffff000 r-xp 00000000 00:00 0 [vdso]
Aborted


Can anybody help please?
Ido.

[Ido Ben-Dayan]

a more accurate information.
We think the bug initiates at the following subroutine:

subroutine ReadSZTemplate(aset, aname, ascale)
Type (CMBdataset) :: aset
real, intent(in) :: ascale
character(LEN=*), intent(IN) :: aname
integer l, unit
real sz
allocate(aset%sz_template(2:lmax))
aset%sz_template = 0
aset%has_sz_template = .true.
call OpenTxtFile(aname, tmp_file_unit)
do
read(tmp_file_unit,*,end=2) l, sz
if (l>=2 .and. l<=lmax) aset%sz_template(l) = ascale * sz/(l*(l+1)/twopi)
end do

Can someone help?
Ido.