26 #ifndef _PROG_COMMONLINES_HH 27 #define _PROG_COMMONLINES_HH 99 void produceSideInfo();
102 void processBlock(
int i,
int j);
105 void getAndPrepareBlock(
int i,
106 std::vector<
MultidimArray<std::complex<double> > > &blockRTFs,
113 void qualifyCommonLines();
116 void solveForShifts();
143 #define POW2(x) (x*x)
double lpf
Low pass filter.
String outputStyle
Output style.
void setIndex(int i, double theta)
void quaternionCommonLines(const DMatrix &quaternions, CommonLineInfo &clInfo)
void setImages(int k1, int k2)
constexpr signed int SMALL_TRIANGLE
void commonlineMatrixCheat(const DMatrix &quaternions, size_t nRays, DMatrix &clMatrix, DMatrix &clCorr)
double hpf
High pass filter.
double angj
Angle of the best common line in image j.
void randomQuaternions(int k, DMatrix &qArray)
void anglesRotationMatrix(size_t nRays, int i, int j, DMatrix &U)
double angi
Angle of the best common line in image i.
ql0001_ & k(htemp+1),(cvec+1),(atemp+1),(bj+1),(bl+1),(bu+1),(x+1),(clamda+1), &iout, infoqp, &zero,(w+1), &lenw,(iw+1), &leniw, &glob_grd.epsmac
std::vector< CommonLine > CLmatrix
int Nthr
Number of threads.
int tripletRotationMatrix(const DMatrix &clMatrix, size_t nRays, int k1, int k2, int k3, DMatrix &R)
void computeSyncMatrix(const DMatrix &clMatrix, size_t nRays, DMatrix &sMatrix)
int Nmpi
Number of processors.
void quaternionToMatrix(const DVector &q, DMatrix &rotMatrix)
double stepAng
Angular sampling.
double outlierFraction
Outlier fraction.
double distanceij
Distance between both common lines.
FileName fn_sel
input file
void saveMatrix(const char *fn, DMatrix &array)
FileName fn_out
output file
CommonLineInfo(int n, int k1=0, int k2=0)
void rotationsFromSyncMatrix(const DMatrix &sMatrix)
double percentile
Percentile (good common lines have very high percentiles)
bool scaleDistance
Scale output measure.