- const TransformJacobianType & jacobian = transform
- ->GetJacobian( fixedImagePoint );
+ TransformJacobianType jacobian;
+ transform->ComputeJacobianWithRespectToParameters(fixedImagePoint, jacobian);
for(unsigned int par=0; par<this->m_NumberOfParameters; par++) {
RealType sumF = itk::NumericTraits< RealType >::Zero;
for(unsigned int par=0; par<this->m_NumberOfParameters; par++) {
RealType sumF = itk::NumericTraits< RealType >::Zero;
for( unsigned int threadID = 0; threadID<this->m_NumberOfThreads; threadID++ ) {
memset( m_ThreaderDerivativeF[threadID].data_block(),
0,
for( unsigned int threadID = 0; threadID<this->m_NumberOfThreads; threadID++ ) {
memset( m_ThreaderDerivativeF[threadID].data_block(),
0,