94 ref=scaled/np.sum(w)  variable
106 assert_allclose(ref,result,rtol=1e-6,atol=1e-6)
114 ref=np.dot(s,w)/np.sum(w)  variable
115 print(ref)
120 assert_allclose(ref,res,2e-5)
126 ref=entropy(s)  variable
127 print(ref)
130 assert_allclose(ref,res,1e-6)
134 assert_allclose(ref,res,1e-10)
143 ref=entropy(sa,sb)  variable
144 print(ref)
147 assert_allclose(ref,res,1e-6)
151 assert_allclose(ref,res,1e-10)
157 ref=logsumexp(s)  variable
158 print(ref)
161 assert_allclose(ref,res,1e-6)
176 ref=np.log(np.dot(sa,sb))  variable
177 print(ref)
184 assert_allclose(ref,res,3e-6)
189 ref = np.exp(sa)  variable
190 print(ref)
194 assert_allclose(ref,res,1e-6)
198 assert_allclose(ref,res,1e-10)
204 ref = np.log(sa)  variable
205 print(ref)
209 assert_allclose(ref,res,2e-5,1e-5)
213 assert_allclose(ref,res,2e-9,1e-9)
218 ref=cholesky(a,lower=True)  variable
219 print(ref)
223 assert_allclose(ref,res,1e-6,1e-6)
227 assert_allclose(ref,res,1e-10,1e-10)
312 ref =a.dot(b)  variable
313 print(ref)
317 assert_allclose(ref,res,1e-10,1e-10)
323 ref = a - b  variable
324 print(ref)
328 assert_allclose(ref,res,1e-10,1e-10)
332 ref = np.abs(s)  variable
334 print(ref)
336 assert_allclose(ref,res,1e-10,1e-10)
341 ref = sa + sb  variable
343 print(ref)
345 assert_allclose(ref,res,1e-10,1e-10)
350 ref = sa - sb  variable
352 print(ref)
354 assert_allclose(ref,res,1e-10,1e-10)
359 ref = sa.dot(sb)  variable
361 print(ref)
363 assert_allclose(ref,res,1e-10,1e-10)
368 ref = sa * sb  variable
370 print(ref)
372 assert_allclose(ref,res,1e-10,1e-10)
376 ref = -sa  variable
378 print(ref)
380 assert_allclose(ref,res,1e-10,1e-10)
384 ref = sa + 0.1  variable
386 print(ref)
388 assert_allclose(ref,res,1e-10,1e-10)
392 ref = sa * 0.1  variable
394 print(ref)
396 assert_allclose(ref,res,1e-10,1e-10)
400 ref = np.mean(sa)  variable
402 print(ref)
404 assert_allclose(ref,res,1e-10,1e-10)
408 ref = np.sum(sa * sa)  variable
410 print(ref)
412 assert_allclose(ref,res,1e-10,1e-10)
416 ref = tstd(sa)  variable
418 print(ref)
420 assert_allclose(ref,res,1e-10,1e-10)
424 ref = tvar(sa)  variable
426 print(ref)
428 assert_allclose(ref,res,1e-10,1e-10)
432 ref = np.ones(nb)*4.0   variable
434 assert_allclose(ref,res,1e-10,1e-10)
439 ref = sa  variable
441 assert_allclose(ref,res,1e-10,1e-10)
446 ref=den//num  variable
449 print(ref)
465 ref = np.abs(ca)  variable
466 print(ref)
469 assert_allclose(ref,res,1e-6,1e-6)
472 assert_allclose(ref,res,1e-10,1e-10)
477 ref = np.abs(ca) * np.abs(ca)  variable
478 print(ref)
481 assert_allclose(ref,res,1e-6,1e-6)
484 assert_allclose(ref,res,1e-10,1e-10)
491 ref = imToReal1D(ca * cb)  variable
492 print(ref)
495 assert_allclose(ref,res,1e-6,1e-6)
499 assert_allclose(ref,res,1e-10,1e-10)