Hello Brian,<br><br>thank you for your answer. I am still confuse about this. In the header file of the Crank-Nicolson solver (itkFEMSolverCrankNicolson.h) one can find the next:<br><br>" * FEM Solver for time dependent problems; uses Crank-Nicolson implicit discretization scheme.<br>
*It solves the following problem:<br> *<br> * ( M + alpha*dt* K )*U_t=(M - (1.- alpha)*dt* K)* U_{t-1} + dt*(alpha*f_{n+1} + (1-alpha)*f_n) "<br><br>This problem formulation differs from the one obtained for the heat equation (for example). So, I guess that it comes from an specific PDE by applying Crank-Nicolson discretization scheme to it. My questions then are:<br>
1. From which PDE does it come from? Which one is the problem that it is trying to solve?<br>2. If the PDE is different depending on the element that one chooses, should not one get different problem formulations to solve by choosing different elements (PDE)? And then is it the solver generic or valid for only one PDE?<br>
<br>Thank you!<br><br>Cristina<br><br><div class="gmail_quote">On Fri, Dec 10, 2010 at 8:59 PM, brian avants <span dir="ltr"><<a href="mailto:stnava@gmail.com">stnava@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hi Cristina<br>
<br>
The specific PDE depends upon the specific elements you choose. For<br>
instance, if you set up a triangular mesh via the<br>
2DC0LinearTriangularMembrane element, then you will use a membrane<br>
energy penalty term.<br>
<br>
The Crank-Nicolson solver uses the Crank-Nicolson approach to<br>
accumulate the solution to the registration problem over time.<br>
<br>
Wikipedia provides a reasonable general explanation of the C-N<br>
approach and gives an example for the heat equation;<br>
<br>
<a href="http://en.wikipedia.org/wiki/Crank%E2%80%93Nicolson_method" target="_blank">http://en.wikipedia.org/wiki/Crank%E2%80%93Nicolson_method</a><br>
<br>
B.<br>
<div><div></div><div class="h5"><br>
<br>
On Fri, Dec 10, 2010 at 10:07 AM, Cristina Oyarzun<br>
<<a href="mailto:coyarzunlaura@googlemail.com">coyarzunlaura@googlemail.com</a>> wrote:<br>
> Hello,<br>
><br>
> I am using itk FEM registration filter. This filter makes use of<br>
> itkFEMSolverCrankNicolson. In the header file of the solver one finds the<br>
> equation that the solver tries to solve. However, one can only find there<br>
> its Crank-Nicolson formulation but not the original PDE that leads to that<br>
> equation.<br>
><br>
> I tried to get the discretized equation by using the equation of motion that<br>
> is described in a book and Crank-Nicolson scheme but did not arrive to the<br>
> same solution. Could someone tell me which on is the original PDE that is<br>
> solved in itkFEMSolverCrankNicolson?<br>
><br>
> Thank you!!<br>
><br>
> Greetings,<br>
><br>
> Cristina<br>
><br>
</div></div>> _____________________________________<br>
> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at<br>
> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Kitware offers ITK Training Courses, for more information visit:<br>
> <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
><br>
> Please keep messages on-topic and check the ITK FAQ at:<br>
> <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
><br>
><br>
<font color="#888888"><br>
<br>
<br>
--<br>
ß®∫∆π<br>
</font></blockquote></div><br>