[Insight-users] finite differences solvers in itk

[Luis Ibanez]

Hi Moti,


You may want to look at the code of the LinearSystemWrapper
in the directory:

         Insight/Code/Numerics/FEM/


The relevant files are:

             itkFEMItpackSparseMatrix.h
             itkFEMLinearSystemWrapper.cxx
             itkFEMLinearSystemWrapper.h
             itkFEMLinearSystemWrapperDenseVNL.cxx
             itkFEMLinearSystemWrapperDenseVNL.h
             itkFEMLinearSystemWrapperItpack.cxx
             itkFEMLinearSystemWrapperItpack.h
             itkFEMLinearSystemWrapperVNL.cxx
             itkFEMLinearSystemWrapperVNL.h
             itkFEMLinearSystemWrappers.h
             itkFEMSolution.h
             itkFEMSolver.cxx
             itkFEMSolver.h
             itkFEMSolverCrankNicolson.h
             itkFEMSolverHyperbolic.cxx
             itkFEMSolverHyperbolic.h


You may also find interested the following classes
in directory

                Insight/Code/Common


             itkDenseFiniteDifferenceImageFilter.h
             itkFiniteDifferenceImageFilter.h
             itkFiniteDifferenceSparseImageFilter.h
             itkFiniteDifferenceSparseImageFilter.txx
             itkFiniteDifferenceSparseImageFunction.h
             itkFiniteDifferenceSparseImageFunction.txx
             itkSparseImage.h
             itkSparseImage.txx



     Regards,


        Luis


------------------
Moti Freiman wrote:
> Hello all,
> In my application i need to solve a very big but sparse linear system
> constructed from a 3D image (volume).
> This system represented a solution to a pde using finite differences scheme.
> I want to solve it using iterative methods such as the Gauss-Seidel
> method, with multi grid.
> 
> I did not found any tool dedicated to such a task in the itk library.
> Does any one know about such tools?
> In any case i want to try to write a general framework for such tasks,
> so users will be able to define the equation easily, and the library
> will solve it. Does anyone try do do so? or have any helpfully ideas?
> Thanks,
> Moti Freiman
> 
> 
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
> 
>