1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
/***************************************************************************
*
* Author: "Sjors H.W. Scheres"
* MRC Laboratory of Molecular Biology
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* This complete copyright notice must be included in any revised version of the
* source code. Additional authorship citations may be added, but existing
* author citations must be preserved.
***************************************************************************/
#include "src/ctffind_runner_mpi.h"
void CtffindRunnerMpi::read(int argc, char **argv)
{
// Define a new MpiNode
node = new MpiNode(argc, argv);
// First read in non-parallelisation-dependent variables
CtffindRunner::read(argc, argv);
// Don't put any output to screen for mpi slaves
verb = (node->isMaster()) ? 1 : 0;
// Possibly also read parallelisation-dependent variables here
// Print out MPI info
printMpiNodesMachineNames(*node);
}
void CtffindRunnerMpi::run()
{
if (!do_only_join_results)
{
// Each node does part of the work
long int my_first_micrograph, my_last_micrograph, my_nr_micrographs;
divide_equally(fn_micrographs.size(), node->size, node->rank, my_first_micrograph, my_last_micrograph);
my_nr_micrographs = my_last_micrograph - my_first_micrograph + 1;
int barstep;
if (verb > 0)
{
std::cout << " Estimating CTF parameters using Niko Grigorieff's CTFFIND ..." << std::endl;
init_progress_bar(my_nr_micrographs);
barstep = XMIPP_MAX(1, my_nr_micrographs / 60);
}
for (long int imic = my_first_micrograph; imic <= my_last_micrograph; imic++)
{
if (verb > 0 && imic % barstep == 0)
progress_bar(imic);
executeCtffind(fn_micrographs[imic]);
}
if (verb > 0)
progress_bar(my_nr_micrographs);
}
MPI_Barrier(MPI_COMM_WORLD);
// Only the master writes the joined result file
if (node->isMaster())
{
joinCtffindResults();
}
}
|
