Print

Print


Hi Matteo,
Apparently, the movie particles have not been extracted using the exact 
same coordinate files as for the "average particles". This is crucial. 
Make sure you follow exactly this page: 
http://www2.mrc-lmb.cam.ac.uk/relion/index.php/Process_movies
HTH, S


On 02/27/2014 08:37 AM, Matteo Allegretti wrote:
> Dear Sjors,
> I am encountering an error when I try to launch autorefine "continue old
> run" with movies. Actually the image number in the error reported below
> change...
> What would you suggest?
>
> Another question: can the grouping procedure be coupled with the movie-
> refinement?
>
> Thank you so much,
> Matteo
>
>   Expectation iteration 21
> 10.93/50.82 min ............~~(,_,">master encountered error: readMRC:
> Image number 73 exceeds stack size 72
> File: ./src/rwMRC.h line: 189
> +++ RELION: command line arguments (with defaults for optional ones
> between parantheses) +++
> ====== Continue options =====
>                           --continue : _optimiser.star file of the
> iteration after which to continue
> ====== General options =====
>                        --o (OLD_ctX) : Output rootname
>                         --iter (OLD) : Maximum number of iterations to perform
>                   --tau2_fudge (OLD) : Regularisation parameter (values
> higher than 1 give more weight to the data)
>              --flatten_solvent (OLD) : Switch on masking on the references?
>                 --solvent_mask (OLD) : User-provided mask for the references
>                --solvent_mask2 (OLD) : User-provided secondary mask
>                          --tau (OLD) : STAR file with input tau2-spectrum
> (to be kept constant)
>            --particle_diameter (OLD) : Diameter of the circular mask that
> will be applied to the experimental images (in Angstroms)
>         --join_random_halves (false) : Join previously split random halves
> again (typically to perform a final reconstruction).
> ====== Re-align movie frames =====
>            --realign_movie_frames () : Input STAR file with the movie frames
>                --nr_frames_prior (5) : Number of movie frames to calculate
> running-average priors
>       --movie_frames_running_avg (3) : Number of movie frames in each
> running average
> ====== Orientations =====
>                 --oversampling (OLD) : Adaptive oversampling order to
> speed-up calculations (0=no oversampling, 1=2x, 2=4x, etc)
>                --healpix_order (OLD) : Healpix order for the angular
> sampling rate on the sphere (before oversampling):
> hp2=15deg, hp3=7.5deg, etc
>                     --psi_step (OLD) : Angular sampling (before
> oversampling) for the in-plane angle (default=10deg for
> 2D, hp sampling for 3D)
>                 --offset_range (OLD) : Search range for origin offsets (in
> pixels)
>                  --offset_step (OLD) : Sampling rate for origin offsets (in
> pixels)
>     --auto_local_healpix_order (OLD) : Minimum healpix order (before
> oversampling) from which auto-refine procedure will use local searches
>                    --sigma_ang (OLD) : Stddev on all three Euler angles for
> local angular searches (of +/- 3 stddev)
>                    --sigma_rot (OLD) : Stddev on the first Euler angle for
> local angular searches (of +/- 3 stddev)
>                   --sigma_tilt (OLD) : Stddev on the first Euler angle for
> local angular searches (of +/- 3 stddev)
>                    --sigma_psi (OLD) : Stddev on the in-plane angle for
> local angular searches (of +/- 3 stddev)
>                    --sigma_off (OLD) : Stddev. on the translations
> ====== Corrections =====
>                        --scale (OLD) : Switch on intensity-scale
> corrections on image groups
>                         --norm (OLD) : Switch on normalisation-error
> correction
> ====== Computation =====
>                              --j (1) : Number of threads to run in parallel
> (only useful on multi-core machines)
>                         --pool (OLD) : Number of images to be processed
> together
>                  --cpus_per_node (8) : Maximum number of models for which
> the FT can be set simultaneously (decrease in case of
> memory problems in this step)
>                           --verb (1) : Verbosity (1=normal, 0=silent)
> ====== Expert options =====
>           --strict_highres_exp (OLD) : Resolution limit (in Angstrom) to
> restrict probability calculations in the expectation step
> --------------------------------------------------------------------------
> MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
> with errorcode 1.
>
> NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
> You may or may not see output from other processes, depending on
> exactly when Open MPI kills them.
>
>
> --
>
> Matteo Allegretti
> PhD candidate
>
> Max Planck Institut für Biophysik
> Max-von-Laue-Straße 3
> D-60438 Frankfurt am Main
> Office: +49-69-6303-3015

-- 
Sjors Scheres
MRC Laboratory of Molecular Biology
Francis Crick Avenue, Cambridge Biomedical Campus
Cambridge CB2 0QH, U.K.
tel: +44 (0)1223 267061
http://www2.mrc-lmb.cam.ac.uk/groups/scheres