# mass = mass[mass>0.0]
# # print the error and show the figure
# print '\nGlobal error = %.10e' % np.abs(np.sum(fields['m'])-1.0)
# hlp.scatterEllipses(mass, cm[:,1:], cov=sxx[:,1:,1:], title='Example 5: Summing over voxels')
# In[ ]:
reload(psi)
reload(hlp)
# Sum fields over a region
inds = np.where(fields['m'] > 0.00000)
# print inds
fsum = psi.sumVoxels(fields, inds)
cm = fsum['x'] / fsum['m']
sxx = fsum['xx'] / fsum['m']
# print fsum['vv']/fsum['m']
print cm
print sxx
hlp.scatterEllipses(mass, [cm[1:]],
cov=[sxx[1:, 1:]],
nstd=1,
title='Example 5: Summing over voxels')
# print np.max(fields['v'])
# # print the error and show the figure # print '\nGlobal error = %.10e' % np.abs(np.sum(fields['m'])-1.0) # hlp.scatterEllipses(mass, cm[:,1:], cov=sxx[:,1:,1:], title='Example 5: Summing over voxels') # In[ ]: reload(psi) reload(hlp) # Sum fields over a region inds = np.where(fields['m'] > 0.00000) # print inds fsum = psi.sumVoxels(fields, inds) cm = fsum['x']/fsum['m'] sxx = fsum['xx']/fsum['m'] # print fsum['vv']/fsum['m'] print cm print sxx hlp.scatterEllipses(mass, [cm[1:]], cov=[sxx[1:,1:]], nstd=1, title='Example 5: Summing over voxels') # print np.max(fields['v']) # In[ ]:
