Commit 8e38d797 authored by Yaroslav Gevorkov's avatar Yaroslav Gevorkov
Browse files

finer multi seed refining

parent 2cfffd64
...@@ -101,12 +101,12 @@ void PinkIndexer::refine(Lattice& indexedLattice, const Matrix3Xf& ucsDirections ...@@ -101,12 +101,12 @@ void PinkIndexer::refine(Lattice& indexedLattice, const Matrix3Xf& ucsDirections
break; break;
case RefinementType::firstFixedThenVariableLatticeParametersMultiSeed: case RefinementType::firstFixedThenVariableLatticeParametersMultiSeed:
{ {
#define REFINEMENT_TRIES 250 #define REFINEMENT_TRIES 750
int fittedNodesCount[REFINEMENT_TRIES]; int fittedNodesCount[REFINEMENT_TRIES];
double fittedNodesMeanDefects[REFINEMENT_TRIES]; double fittedNodesMeanDefects[REFINEMENT_TRIES];
Lattice fittedLattices[REFINEMENT_TRIES]; Lattice fittedLattices[REFINEMENT_TRIES];
float maxRelativeDeviation = 0.012; float maxRelativeDeviation = 0.0125;
Array<float, 1, 3> columnDeviationNorms = indexedLattice.getBasis().colwise().norm() * maxRelativeDeviation; Array<float, 1, 3> columnDeviationNorms = indexedLattice.getBasis().colwise().norm() * maxRelativeDeviation;
#pragma omp parallel for #pragma omp parallel for
for (int i = 0; i < REFINEMENT_TRIES; ++i) for (int i = 0; i < REFINEMENT_TRIES; ++i)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment