@@ -222,8 +222,6 @@ ComputeJacobianTerms<TFixedImage, TTransform>::Compute(double & TrC, double & Tr
222
222
/* * Initialize band matrix. */
223
223
CovarianceMatrixType bandcov (numberOfParameters, bandcovsize, 0.0 );
224
224
225
- unsigned int samplenr = 0 ;
226
-
227
225
/* *
228
226
* TERM 1
229
227
*
@@ -238,10 +236,6 @@ ComputeJacobianTerms<TFixedImage, TTransform>::Compute(double & TrC, double & Tr
238
236
}
239
237
for (const auto & sample : *sampleContainer)
240
238
{
241
- /* * Print progress 0-50%.
242
- *progressObserver->UpdateAndPrintProgress( samplenr );*/
243
- ++samplenr;
244
-
245
239
/* * Read fixed coordinates and get Jacobian J_j. */
246
240
const FixedImagePointType & point = sample.m_ImageCoordinates ;
247
241
this ->m_Transform ->GetJacobian (point, jacj, jacind);
@@ -412,7 +406,6 @@ ComputeJacobianTerms<TFixedImage, TTransform>::Compute(double & TrC, double & Tr
412
406
JacobianType jacjcovjacj (outdim, outdim);
413
407
itk::Array<SizeValueType> jacindExpanded (numberOfParameters);
414
408
415
- samplenr = 0 ;
416
409
for (const auto & sample : *sampleContainer)
417
410
{
418
411
/* * Read fixed coordinates and get Jacobian. */
@@ -504,10 +497,6 @@ ComputeJacobianTerms<TFixedImage, TTransform>::Compute(double & TrC, double & Tr
504
497
/* * Max_j [JCJ_j]. */
505
498
maxJCJ = std::max (maxJCJ, JCJ_j);
506
499
507
- /* * Show progress 50-100%. */
508
- // progressObserver->UpdateAndPrintProgress( samplenr + nrofsamples );
509
- ++samplenr;
510
-
511
500
} // end loop over sample container
512
501
513
502
/* * Finalize progress information. */
0 commit comments