Commit a556ef19 authored by Reimar H Leike's avatar Reimar H Leike
Browse files

saved LOS-response from potentially nasty overflow

parent ada9e2d2
...@@ -206,7 +206,7 @@ class LOSResponse(LinearOperator): ...@@ -206,7 +206,7 @@ class LOSResponse(LinearOperator):
tmp += (fullidx[j]//boxsz)*fct tmp += (fullidx[j]//boxsz)*fct
fct *= self._local_shape[j] fct *= self._local_shape[j]
tmp += cnt/float(nlos) tmp += cnt/float(nlos)
tmp += iarr[ofs:ofs+nval]/float(nlos*npix) tmp += iarr[ofs:ofs+nval]/(float(nlos)*float(npix))
pri[ofs:ofs+nval] = tmp pri[ofs:ofs+nval] = tmp
ofs += nval ofs += nval
cnt += 1 cnt += 1
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment