cSAXS/cdtools
0
0
Fork 0
mirror of https://github.com/cdtools-developers/cdtools.git synced 2026-10-12 19:10:26 +02:00
Code Issues Packages Projects Releases Wiki Activity

some horrible state but I need to stop developing on perlmutter

Browse Source
This commit is contained in:
abe
2022-10-17 16:32:59 -07:00
parent cb954a6067
commit e13b519e22
5 changed files with 24 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/cdtools/models/simple_ptycho.py
+10 -3
View File
@@ -10,6 +10,11 @@ import numpy as np
__all__ = ['SimplePtycho']
class complexWrapper(t.Tensor):
def __new__(base_tensor):
return t.view_as_complex(base_tensor)
class SimplePtycho(CDIModel):
"""A simple ptychography model for exploring ideas and extensions
@@ -44,11 +49,14 @@ class SimplePtycho(CDIModel):
# object
self.register_buffer('probe_norm', t.max(t.abs(probe_guess)))
self.probe_data = t.nn.Parameter(t.view_as_real(probe_guess / self.probe_norm))
self.probe_data = complexParameter(probe_guess/self.probe_norm)
#self.probe_data = t.nn.Parameter(t.view_as_real(probe_guess / self.probe_norm))
self.obj_data = t.nn.Parameter(t.view_as_real(obj_guess))
@property
def probe(self):
probe = t.view_as_complex(self.probe_data)
return t.view_as_complex(self.probe_data)
@property
@@ -106,7 +114,6 @@ class SimplePtycho(CDIModel):
def interaction(self, index, translations):
pix_trans = tools.interactions.translations_to_pixel(self.probe_basis,
translations,
surface_normal=self.surface_normal)
@@ -181,7 +188,7 @@ class SimplePtycho(CDIModel):
('Object Phase',
lambda self, fig: p.plot_phase(self.obj, fig=fig, basis=self.probe_basis))
]
def ePIE(self, iterations, dataset, beta = 1.0):
"""Runs an ePIE reconstruction as described in `Maiden et al. (2017) <https://www.osapublishing.org/optica/abstract.cfm?uri=optica-4-7-736>`_.
src/cdtools/tools/interactions/interactions.py
+2 -2
View File
@@ -43,7 +43,6 @@ def translations_to_pixel(basis, translations, surface_normal=t.Tensor([0.,0.,1.
pixel_translations : torch.Tensor
A Jx2 stack of translations in internal (i,j) pixel-space, or a single translation
"""
projection_1 = t.as_tensor(np.array([[1,0,0],
[0,1,0],
[0,0,0]]),
@@ -53,7 +52,8 @@ def translations_to_pixel(basis, translations, surface_normal=t.Tensor([0.,0.,1.
projection_2[2] = t.as_tensor(-surface_normal / surface_normal[2],
dtype=projection_2.dtype,
device=projection_2.device)
projection_2 = t.inverse(projection_2)
projection_2 = t.linalg.inv(projection_2)
basis_vectors_inv = t.pinverse(basis).to(device=translations.device,
dtype=translations.dtype)