<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<title></title>
</head>
<body>
Hi,<br>
<br>
I would like to reask the same question as an old mail sent to this
mailing list last summer (see below). Apparently, it didn't receive
any answer. I notice that m_DT is linked to curvature-based evolution
terms and m_WaveDT to the other evolution terms. Is it because
curvature-based terms is more sensitive to numerical instabilities than
the others and could need a smaller dt?<br>
<br>
How did you choose the following values for m_WaveDT and m_DT? Is
there another justification than the one given by Sethian?<br>
template< class TImageType ><br>
double LevelSetFunction<TImageType>::m_WaveDT = 1.0/(2.0 *
ImageDimension);<br>
<br>
template < class TImageType ><br>
double LevelSetFunction<TImageType>::m_DT = 1.0/(2.0 *
ImageDimension);<br>
<br>
Thanks in advance for your help,<br>
<br>
Valerie<br>
<h1>[Insight-users] time step in itkLevelSetFunction</h1>
<b>Amadeus W. M.</b> <a
href="mailto:insight-users%40itk.org?Subject=%5BInsight-users%5D%20time%20step%20in%20itkLevelSetFunction&In-Reply-To="
title="[Insight-users] time step in itkLevelSetFunction">amadeus84 at
cablespeed.com </a><br>
<i>Sat Aug 6 13:15:37 EDT 2005</i>
<ul>
<li>Previous message: <a
href="http://public.kitware.com/pipermail/insight-users/2005-August/014327.html">[Insight-users]
Re: Level set: Time vs. Iteration </a></li>
<li>Next message: <a
href="http://public.kitware.com/pipermail/insight-users/2005-August/014330.html">[Insight-users]
dicom stack orientation information </a></li>
<li> <b>Messages sorted by:</b> <a
href="http://public.kitware.com/pipermail/insight-users/2005-August/date.html#14329">[
date ]</a> <a
href="http://public.kitware.com/pipermail/insight-users/2005-August/thread.html#14329">[
thread ]</a> <a
href="http://public.kitware.com/pipermail/insight-users/2005-August/subject.html#14329">[
subject ]</a> <a
href="http://public.kitware.com/pipermail/insight-users/2005-August/author.html#14329">[
author ]</a> </li>
</ul>
<hr><!--beginarticle-->
<pre>Can someone explain why the time step is computed the way it is in
the function ComputeGlobalTimeStep in itkLevelSetFunction.txx?
I mean the mathematical justification. What are m_WaveDT and m_DT?
In Level Set Methods, Sethian only mentions (page 67, section 6.4.1)
max F \Delta T <= \Delta x
A pointer to the appropriate theory would suffice.
Thank you!
</pre>
<br>
</body>
</html>