[OpenWalnut] tensor line integration

[Mario Hlawitschka]

Hi Francisco,

If you had looked at the code as I told you, you would have found a file
called 
src/modules/eigenSystem/WMEigenSystem.h
which is a good starting point for knowing how to compute the eigensystem.
There is a module called WMDeterministicFTMori that does the line integration as it is done for diffusion tensor data. Technically, the code provides a solution to your problem.
Depending on your data and your treatment of singularities, you may want to modify minor things.

Regards,

	Mario


Am Donnerstag, den 30.06.2011, 17:00 -0400 schrieb Francisco J.
Beron-Vera: 
> Hi,
> 
> Let x denotes position on D \in R^2 and T(x) is a (symmetric) 2-tensor 
> field on D. Then
> 
> (1) T(x) e(x) = \lambda(x) e(x)
> 
> where (e(x),\lambda(x)) is an eigenvector field/eigenvalue field pair. 
> Then consider the ODE
> 
> (2) x'(s) = e(x(s)), x(0) = x_0,
> 
> where s is the natural parameterization (arclength).
> 
> I need to compute the integral curves of (1) - tensor lines. I have a 
> table of numbers specifying the components of e on a regular lattice. 
> How can do this computation using OpenWalnut?
> 
> Thank you in advance.
> 
> Francisco
>