BoundaryCondition override is not respected for rtkBackwardDifferenceDivergenceImageFilter.hxx about rtk HOT 1 CLOSED

//Author \Hans J. Johnson
// Test program to find bug in periodic boundary conditions.

#include <rtkForwardDifferenceGradientImageFilter.h>
#include <rtkBackwardDifferenceDivergenceImageFilter.h>

const size_t IM2DSIZE=3;

template
void DumpImage(IMTYPE out)
{
typename IMTYPE::ObjectType::IndexType idx;
std::cout << "==========================================" << std::endl;
for(size_t k = 0; k < IM2DSIZE; ++k)
{
idx[2] = k;
for(size_t j = 0; j < IM2DSIZE; ++j)
{
idx[1]=j;
for(size_t i = 0; i< IM2DSIZE; ++i)
{
idx[0]=i;
std::cout.width(16);
std::cout.fill(' ');
std::cout << std::right << out->GetPixel(idx);
}
std::cout << std::endl;
}
std::cout << ". - . - . - . - . - . - . - . - . - . - . " << std::endl;
}
std::cout << "==========================================" << std::endl;
}

int main( int argc, char * argv[])
{
typedef itk::Image<float,3> Image2D;
Image2D::Pointer im2 = Image2D::New();
Image2D::RegionType region;
Image2D::RegionType::SizeType size;
size.Fill(IM2DSIZE);
region.SetSize(size);
im2->SetRegions(region);
im2->Allocate(true);
Image2D::IndexType idx;
idx.Fill(IM2DSIZE/2);
idx.Fill(0);
im2->SetPixel(idx,1.0F);
//typedef itk::GradientImageFilter GradientType;
typedef rtk::ForwardDifferenceGradientImageFilter GradientType;
GradientType::Pointer grad = GradientType::New();
grad->SetInput(im2);
typedef itk::PeriodicBoundaryCondition FloatBoundaryType;
grad->OverrideBoundaryCondition( new FloatBoundaryType );
grad->Update();
GradientType::OutputImageType::Pointer gradIm = grad->GetOutput();
//typedef itk::DivergenceImageFilterGradientType::OutputImageType DivergenceType;
typedef rtk::BackwardDifferenceDivergenceImageFilterGradientType::OutputImageType DivergenceType;
DivergenceType::Pointer div = DivergenceType::New();
div->SetInput(gradIm);
typedef itk::PeriodicBoundaryConditionGradientType::OutputImageType CVBoundaryType;
div->OverrideBoundaryCondition(new CVBoundaryType);
div->Update();
DivergenceType::OutputImageType::Pointer divIm = div->GetOutput();
std::cout << "Intensity Image" << std::endl;
DumpImage(im2);
DumpImage(gradIm);
DumpImage(divIm);

return 0;
}

from rtk.