Skip to content

TotalVariation proximal in place  #1650

@MargaretDuff

Description

@MargaretDuff

TotalVariation proximal is not calculated in place correctly. This was spotted when trying to run SIRT with a TV constraint in #1514.

data = dataexample.SIMPLE_PHANTOM_2D.get(size=(128,128))
plt.figure()
plt.imshow(data.as_array())
plt.show()
out=data.geometry.allocate('random')
TotalVariation().proximal(data, tau=1, out=out)
plt.figure()
plt.imshow(out.as_array())
plt.show()

gives
image
image

data = dataexample.SIMPLE_PHANTOM_2D.get(size=(128,128))
plt.figure()
plt.imshow(data.as_array())
plt.show()
TotalVariation().proximal(data, tau=1, out=data)
plt.figure()
plt.imshow(data.as_array())
plt.show()

gives
image
image

Metadata

Metadata

Assignees

Labels

bugSomething isn't working

Type

No type

Projects

Status

Todo

Relationships

None yet

Development

No branches or pull requests

Issue actions