Mercurial > repos > rv43 > tomo

changeset 29:59c5c0b53579 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/rolfverberg/galaxytools commit dbc4ff10f088b0217bab0ce5d6d5a7763619c779"
author rv43
date Mon, 18 Apr 2022 21:38:30 +0000
parents c7211931fc99
children c4a90d5b1ed2
files tomo.py
diffstat 1 files changed, 4 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/tomo.py	Mon Apr 18 21:26:01 2022 +0000
+++ b/tomo.py	Mon Apr 18 21:38:30 2022 +0000
@@ -1984,15 +1984,17 @@
                 continue
             print(f'lower_row = {lower_row} upper_row = {upper_row} self.tomo_stacks[i].shape[0] = {self.tomo_stacks[i].shape[0]}')
             assert(0 <= lower_row < upper_row < self.tomo_stacks[i].shape[0])
-            print(f'OK4 {i} c')
             center_offsets = [lower_center_offset-lower_row*center_slope,
                     upper_center_offset+(self.tomo_stacks[i].shape[0]-1-upper_row)*center_slope]
             t0 = time()
+            num_core = 1
+            print(f'OK4 {i} c')
             self.tomo_recon_stacks[i]= self._reconstructOneTomoStack(self.tomo_stacks[i],
                     thetas, center_offsets=center_offsets, sigma=0.1, num_core=num_core,
                     algorithm='gridrec', run_secondary_sirt=True, secondary_iter=25)
+            print(f'OK4 {i} d')
+            exit('Done')
             logging.info(f'Reconstruction of stack {index} took {time()-t0:.2f} seconds!')
-            print(f'OK4 {i} d')
             if self.galaxy_flag:
                 x_slice = int(self.tomo_stacks[i].shape[0]/2) 
                 title = f'{basetitle} {index} xslice{x_slice}'