[Insight-developers] FEM BoundaryCondition : Is it functional ?

[Luis Ibanez]

Brian, Jim, Vincent,

Julie and Chuck are making great progress
on evaluating the Numerical Libraries in ITK.


They have found an strange behavior in the
FEM Testing related to the use of Boundary
Conditions.

In particular in:

ITK/Modules/Numerics/FEM/include/
              itkFEMSolverCrankNicolson.hxx

in line 226-227:

  /* step over all types of BCs */
  this->ApplyBC();  // BUG  -- are BCs applied appropriately to the problem?


The surprising thing is that when the
line above is commented out, the tests
still pass.

We are wondering is the call to ApplyBC()
is needed at all, or whether the tests are
not performing a strict enough evaluation
of the output.

Could you please comment on this issue ?


   Thanks


        Luis