On Sun, Dec 30, 2012 at 3:41 AM, Happyman <bahtiyor_zohi...@mail.ru> wrote: > nums=32 > rows=120 > cols=150 > > for k in range(0,nums): > for i in range(0,rows): > for j in range(0,cols): > if float ( R[ k ] [ i ] [ j ] ) == > 0.0:
why the float() -- what data type is R? > else: > val11[ i ] [ j ], val22[ i > ][ j ] = integrate.quad( lambda x : F1(x)*F2(x) , 0 , pi) this is odd -- Do F1 and F2 depend on i,j, or k somehow? or are you somehow integerting over the k-dimension? In which case, I'm guessing that integration is you time killer anyway -- do some profiling to know for sure. -Chris -- Christopher Barker, Ph.D. Oceanographer Emergency Response Division NOAA/NOS/OR&R (206) 526-6959 voice 7600 Sand Point Way NE (206) 526-6329 fax Seattle, WA 98115 (206) 526-6317 main reception chris.bar...@noaa.gov _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
