diff options
Diffstat (limited to 'src/refine.c')
-rw-r--r-- | src/refine.c | 7 |
1 files changed, 5 insertions, 2 deletions
diff --git a/src/refine.c b/src/refine.c index d0fc8cc..b3ac84a 100644 --- a/src/refine.c +++ b/src/refine.c @@ -212,7 +212,7 @@ static void refine_simplex_contract(SimplexVertex *s, int v, int v_best) { s[v].dcx = s[v_best].dcx + 0.5 * (s[v].dcx - s[v_best].dcx); s[v].dcy = s[v_best].dcy + 0.5 * (s[v].dcy - s[v_best].dcy); s[v].dcz = s[v_best].dcz + 0.5 * (s[v].dcz - s[v_best].dcz); - + } static double refine_iteration(SimplexVertex *s, Deviation *d, int nf) { @@ -326,7 +326,7 @@ double refine_do_cell(ControlContext *ctx) { Deviation *d; double delta; int i, nf, f, it, maxiter; - const double tol = 0.0001e9; /* Stopping condition */ + const double tol = 0.00001e9; /* Stopping condition */ //const double tol = 0.001; /* For testing */ if ( !ctx->cell_lattice ) { @@ -342,6 +342,7 @@ double refine_do_cell(ControlContext *ctx) { /* Determine the size of the 'deviation table' */ nf = 0; for ( i=0; i<ctx->images->n_images; i++ ) { + int j; if ( !ctx->images->images[i].rflist ) { @@ -352,6 +353,8 @@ double refine_do_cell(ControlContext *ctx) { if ( ctx->images->images[i].rflist->features[j].partner != NULL ) nf++; } + printf("%i features from image %i\n", nf, i); + } if ( REFINE_DEBUG ) printf("RF: There are %i partnered features in total\n", nf); |