[Insight-users] MPI Fast Marching

[Luis Ibanez]

Hi Kevin,

This is great !

I hope you are planning on submitting it to the Insight Journal.


Regarding your point about the float comparison, I agree with
you in that an exact comparison is unrealistic.  You may want
to consider the use of the following "epsilon" values:

              vcl_numeric_limits<double>::epsilon()
              vcl_numeric_limits<float>::epsilon()

as tolerances for the comparison.


Another option is to use the value:

            itk::NumericTraits< float >::min()


     Luis


-------------------------------------------------------------------
On Fri, Oct 16, 2009 at 11:51 AM, Kevin H. Hobbs <hobbsk at ohiou.edu> wrote:
> On Mon, 2009-10-12 at 18:52 -0400, Luis Ibanez wrote:
>> Hi Kevin,
>>
>> Thanks for bringing this document to our attention.
>
> I have attached a working but incomplete MPI Fast Marching program based
> on Tugurlan's Dissertation.
>
> The communication of boundary values and iteration are done.
>
> I still have to generate new trial and alive points after each
> iteration. Right now I just use the original seeds and the received
> boundary points with values less than the minimum boundary point from
> the local piece.
>
> I still have to test for convergance, which I believe will happen when
> none of the MPI processes receives any boundary values less than the
> minimum boundary point from the local piece.
>
> I still have to add an extractor to cut the fast marching image to the
> final dimensions.
>