;
; PURPOSE: benchmark for SPIDER back projection and projection of volume
; USES:    BP RP - Back Projection - 3D, iterative, with constraints, ||
;
; timebprp.dat creates all necessary files
;
; USAGE: time $spib/spider4 spi @timebprp

; Turn off  un-necessary output
MD
VB OFF

; Use all available processors
MD 
SET MP
0             ; All available processors

; Create 75 x 75 x75 volume with sphere inside.
MO 3
jnkvol001        ; Volume     (output)
75,75,75      ; Size             
SP            ; Sphere
N
1
9,3           ; Radii
30,29,28      ; Center
0,0,0         ; No more spheres

; Create angles doc file for Eulerian angles --> 5088 dir.
VO EA X53
2
0.0,0.0
0.0,0.0
jnk_angles   ; Angles doc file   (output)

; X53 is now 5088

; Create selection doc file
DO I=1,X53
   SD X0,X0
   jnk_select ; Selection file    (output)
ENDDO

SD E
jnk_select   ; Selection file    

; Create empty inline stack for 5088 images.
MS           ; Make stack
_8@          ; Inline stack      (output)
75,75,1      ; Size             
x53          ; Number of images

; Compute projection of volume using Eulerian angles.
PJ 3Q        ; Projection
jnkvol001    ; Volume            (input)
35           ; Radius
jnk_select   ; Selection file    (input)
jnk_angles   ; Angles doc file   (input)
_8@****      ; Projection stack  (output)

; Run iterative backprojection now.
BP RP        ; Back projection
_8@****      ; Image file template (input)
jnk_select   ; Selection file      (input)
23           ; Radius
jnk_angles   ; Angles doc file     (input)
*            ; No symmetries
jnk_vol002   ; Volume              (output)
0.2E-5,0.0   ; Lambda & correction limit
10,0         ; Iteration limit, Mode
.5,.5        ; Min & max
.5           ; Smoothing constant

EN