<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> <html><head><meta name="robots" content="noindex"> <meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1"> <title>ITK: itkFEMLinearSystemWrapperVNL.h Source File</title> <link href="DoxygenStyle.css" rel="stylesheet" type="text/css"> </head><body bgcolor="#ffffff"> <!-- Section customized for INSIGHT : Tue Jul 17 01:02:45 2001 --> <center> <a href="index.html" class="qindex">Main Page</a> <a href="modules.html" class="qindex">Groups</a> <a href="namespaces.html" class="qindex">Namespace List</a> <a href="hierarchy.html" class="qindex">Class Hierarchy</a> <a href="classes.html" class="qindex">Alphabetical List</a> <a href="annotated.html" class="qindex">Compound List</a> <a href="files.html" class="qindex">File List</a> <a href="namespacemembers.html" class="qindex">Namespace Members</a> <a href="functions.html" class="qindex">Compound Members</a> <a href="globals.html" class="qindex">File Members</a> <a href="pages.html" class="qindex">Concepts</a></center> <!-- Generated by Doxygen 1.5.9 --> <div class="navpath"><a class="el" href="dir_02e23d841c198393a6e961e9a0678fff.html">Code</a> » <a class="el" href="dir_2c33022d19e8077162952ea9ed57656d.html">Numerics</a> » <a class="el" href="dir_316e5ab5c8a9cb56a9dce1c9fc454a23.html">FEM</a> </div> <div class="contents"> <h1>itkFEMLinearSystemWrapperVNL.h</h1><a href="itkFEMLinearSystemWrapperVNL_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*=========================================================================</span> <a name="l00002"></a>00002 <span class="comment"></span> <a name="l00003"></a>00003 <span class="comment"> Program: Insight Segmentation & Registration Toolkit</span> <a name="l00004"></a>00004 <span class="comment"> Module: $RCSfile: itkFEMLinearSystemWrapperVNL.h,v $</span> <a name="l00005"></a>00005 <span class="comment"> Language: C++</span> <a name="l00006"></a>00006 <span class="comment"> Date: $Date: 2009-01-30 21:10:12 $</span> <a name="l00007"></a>00007 <span class="comment"> Version: $Revision: 1.15 $</span> <a name="l00008"></a>00008 <span class="comment"></span> <a name="l00009"></a>00009 <span class="comment"> Copyright (c) Insight Software Consortium. All rights reserved.</span> <a name="l00010"></a>00010 <span class="comment"> See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.</span> <a name="l00011"></a>00011 <span class="comment"></span> <a name="l00012"></a>00012 <span class="comment"> This software is distributed WITHOUT ANY WARRANTY; without even </span> <a name="l00013"></a>00013 <span class="comment"> the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR </span> <a name="l00014"></a>00014 <span class="comment"> PURPOSE. See the above copyright notices for more information.</span> <a name="l00015"></a>00015 <span class="comment"></span> <a name="l00016"></a>00016 <span class="comment">=========================================================================*/</span> <a name="l00017"></a>00017 <a name="l00018"></a>00018 <span class="preprocessor">#ifndef __itkFEMLinearSystemWrapperVNL_h</span> <a name="l00019"></a>00019 <span class="preprocessor"></span><span class="preprocessor">#define __itkFEMLinearSystemWrapperVNL_h</span> <a name="l00020"></a>00020 <span class="preprocessor"></span><span class="preprocessor">#include "<a class="code" href="itkFEMLinearSystemWrapper_8h.html">itkFEMLinearSystemWrapper.h</a>"</span> <a name="l00021"></a>00021 <span class="preprocessor">#include "vnl/vnl_sparse_matrix.h"</span> <a name="l00022"></a>00022 <span class="preprocessor">#include "vnl/vnl_vector.h"</span> <a name="l00023"></a>00023 <span class="preprocessor">#include <vnl/vnl_sparse_matrix_linear_system.h></span> <a name="l00024"></a>00024 <span class="preprocessor">#include <vnl/algo/vnl_lsqr.h></span> <a name="l00025"></a>00025 <span class="preprocessor">#include <vector></span> <a name="l00026"></a>00026 <a name="l00027"></a>00027 <a name="l00028"></a>00028 <span class="keyword">namespace </span>itk { <a name="l00029"></a>00029 <span class="keyword">namespace </span>fem { <a name="l00030"></a>00030 <a name="l00031"></a>00031 <a name="l00038"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html">00038</a> <span class="keyword">class </span><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html" title="LinearSystemWrapper class that uses VNL numeric library functions to define a sparse...">LinearSystemWrapperVNL</a> : <span class="keyword">public</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapper.html" title="Defines all functions required by Solver class to allocate, assemble and solve a...">LinearSystemWrapper</a> <a name="l00039"></a>00039 { <a name="l00040"></a>00040 <span class="keyword">public</span>: <a name="l00041"></a>00041 <a name="l00042"></a>00042 <span class="comment">/* values stored in matrices & vectors */</span> <a name="l00043"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">00043</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1fem_1_1Solution.html#c95fdd34e9d556c92091685c952d5397">LinearSystemWrapper::Float</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a>; <a name="l00044"></a>00044 <a name="l00045"></a>00045 <span class="comment">/* superclass */</span> <a name="l00046"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#026cbe1b4f8116a938665aadc7601b6d">00046</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapper.html" title="Defines all functions required by Solver class to allocate, assemble and solve a...">LinearSystemWrapper</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapper.html" title="Defines all functions required by Solver class to allocate, assemble and solve a...">SuperClass</a>; <a name="l00047"></a>00047 <a name="l00048"></a>00048 <span class="comment">/* matrix typedef */</span> <a name="l00049"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#7f384e406e7e194cbddc367d82ccadba">00049</a> <span class="keyword">typedef</span> vnl_sparse_matrix<Float> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#7f384e406e7e194cbddc367d82ccadba">MatrixRepresentation</a>; <a name="l00050"></a>00050 <a name="l00051"></a>00051 <span class="comment">/* matrix holder typedef */</span> <a name="l00052"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#f446a0d92254dc34e9cb6ffc1cffda73">00052</a> <span class="keyword">typedef</span> std::vector< MatrixRepresentation* > <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#f446a0d92254dc34e9cb6ffc1cffda73">MatrixHolder</a>; <a name="l00053"></a>00053 <a name="l00054"></a>00054 <span class="comment">/* constructor & destructor */</span> <a name="l00055"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#9acc28d385d7691733d53e41517ab0df">00055</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#9acc28d385d7691733d53e41517ab0df">LinearSystemWrapperVNL</a>() : <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapper.html" title="Defines all functions required by Solver class to allocate, assemble and solve a...">LinearSystemWrapper</a>(), m_Matrices(0), m_Vectors(0), m_Solutions(0) {} <a name="l00056"></a>00056 <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#b57df863d671fc54a1f69b5eaf82e847">~LinearSystemWrapperVNL</a>(); <a name="l00057"></a>00057 <a name="l00058"></a>00058 <span class="comment">/* memory management routines */</span> <a name="l00059"></a>00059 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#53fffb7a918a445bacbe72ad8f6d804c">InitializeMatrix</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex); <a name="l00060"></a>00060 <span class="keyword">virtual</span> <span class="keywordtype">bool</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#790b3f9290f70fd5a69dfb050642e204">IsMatrixInitialized</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex); <a name="l00061"></a>00061 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#fc5bbe044d67f7515bb33bc767b39229">DestroyMatrix</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex); <a name="l00062"></a>00062 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#b8159a2856f01af4c9925560a44381c5">InitializeVector</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex); <a name="l00063"></a>00063 <span class="keyword">virtual</span> <span class="keywordtype">bool</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#b9f031c9abde9d667cf69594e4720704">IsVectorInitialized</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex); <a name="l00064"></a>00064 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#ab4300ecbcf8bfe5aa11f7a0d2ce2278">DestroyVector</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex); <a name="l00065"></a>00065 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#ccea125d272634011314ba16f7b7fa4a">InitializeSolution</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex); <a name="l00066"></a>00066 <span class="keyword">virtual</span> <span class="keywordtype">bool</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#c86ab2eebaee17dd02398fdf9cab4071">IsSolutionInitialized</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex); <a name="l00067"></a>00067 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#40c5825537bda3a641c5d3cbe2daefe5">DestroySolution</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex); <a name="l00068"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#f8c404dc2205c769cf97459d045504e1">00068</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#f8c404dc2205c769cf97459d045504e1">SetMaximumNonZeroValuesInMatrix</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span>, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span>) {} <a name="l00069"></a>00069 <a name="l00070"></a>00070 <a name="l00071"></a>00071 <span class="comment">/* assembly & solving routines */</span> <a name="l00072"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#db3ed9376038ae0f69433a5121acc801">00072</a> <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#db3ed9376038ae0f69433a5121acc801">GetMatrixValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> j, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex)<span class="keyword"> const </span>{ <span class="keywordflow">return</span> (*((*m_Matrices)[matrixIndex]))(i,j); } <a name="l00073"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#d5f0bee12782dfe3da34753d1c2ef470">00073</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#d5f0bee12782dfe3da34753d1c2ef470">SetMatrixValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> j, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> value, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex) { (*((*m_Matrices)[matrixIndex]))(i,j) = value; } <a name="l00074"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#05ca6952807e46cc2df806afe9d6515e">00074</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#05ca6952807e46cc2df806afe9d6515e">AddMatrixValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> j, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> value, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex) { (*((*m_Matrices)[matrixIndex]))(i,j) += value; } <a name="l00075"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#93ae4c400af31fb88fbf13a3dc7f50de">00075</a> <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#93ae4c400af31fb88fbf13a3dc7f50de">GetVectorValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex)<span class="keyword"> const </span>{ <span class="keywordflow">return</span> (* ( (*m_Vectors)[vectorIndex] ) )[i]; } <a name="l00076"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#68c49564b489c1632d753753a80cc271">00076</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#68c49564b489c1632d753753a80cc271">SetVectorValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> value, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex) { (*((*m_Vectors)[vectorIndex]))(i) = value; } <a name="l00077"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#ec17d9bad284db14c9603451c626d200">00077</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#ec17d9bad284db14c9603451c626d200">AddVectorValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> value, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex) { (*((*m_Vectors)[vectorIndex]))(i) += value; } <a name="l00078"></a>00078 <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#3e56a6f3a5330b054e20435b9c3e7ef7">GetSolutionValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex) <span class="keyword">const</span>; <a name="l00079"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#97d5c9e3e7b33858094185666744eae3">00079</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#97d5c9e3e7b33858094185666744eae3">SetSolutionValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> value, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex) { (*((*m_Solutions)[solutionIndex]))(i) = value; } <a name="l00080"></a><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#a0a830264cafe7db59b9c2a2acc1f236">00080</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#a0a830264cafe7db59b9c2a2acc1f236">AddSolutionValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> value, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex) { (*((*m_Solutions)[solutionIndex]))(i) += value; } <a name="l00081"></a>00081 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#efbcb6cdc1c083e7cf63b3b843f4c281">Solve</a>(<span class="keywordtype">void</span>); <a name="l00082"></a>00082 <a name="l00083"></a>00083 <span class="comment">/* matrix & vector manipulation routines */</span> <a name="l00084"></a>00084 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#5e1a3fda1ab7f7a688b75978eee81bf7">ScaleMatrix</a>(<a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#0a78d2f37882a6b0a796753ab78e44b3">Float</a> scale, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex); <a name="l00085"></a>00085 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#aea7f17e881f23576f0aefeee1d2cdfb">SwapMatrices</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex1, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex2); <a name="l00086"></a>00086 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#44538aada0039dc52e4b4c7d64256a92">SwapVectors</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex1, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex2); <a name="l00087"></a>00087 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#aa5e6eba6296e060f68fdc0debc88d21">SwapSolutions</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex1, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex2); <a name="l00088"></a>00088 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#ba1f7771e7fc9ad1d38454d48602d9be">CopySolution2Vector</a>(<span class="keywordtype">unsigned</span> solutionIndex, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex); <a name="l00089"></a>00089 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#490c2ad96515fe83c6e7006b8b250de5">CopyVector2Solution</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex); <a name="l00090"></a>00090 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#3b1fcf36e654a3092715cc946f0516e3">MultiplyMatrixMatrix</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> resultMatrixIndex, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> leftMatrixIndex, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> rightMatrixIndex); <a name="l00091"></a>00091 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#b291ad2c9e9d9acee09e902492bc49a3">MultiplyMatrixVector</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> resultVectorIndex, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> matrixIndex, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex); <a name="l00092"></a>00092 <a name="l00093"></a>00093 <span class="keyword">private</span>: <a name="l00094"></a>00094 <a name="l00096"></a>00096 <span class="comment">//std::vector< vnl_sparse_matrix<Float>* > *m_Matrices;</span> <a name="l00097"></a>00097 <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperVNL.html#f446a0d92254dc34e9cb6ffc1cffda73">MatrixHolder</a> *m_Matrices; <a name="l00098"></a>00098 <a name="l00100"></a>00100 std::vector< vnl_vector<Float>* > *m_Vectors; <a name="l00101"></a>00101 <a name="l00103"></a>00103 std::vector< vnl_vector<Float>* > *m_Solutions; <a name="l00104"></a>00104 <a name="l00105"></a>00105 }; <a name="l00106"></a>00106 <a name="l00107"></a>00107 }} <span class="comment">// end namespace itk::fem</span> <a name="l00108"></a>00108 <a name="l00109"></a>00109 <span class="preprocessor">#endif</span> <a name="l00110"></a>00110 <span class="preprocessor"></span> </pre></div></div> <hr><address><small> Generated at Thu May 7 22:41:57 2009 for ITK by <a href="http://www.stack.nl/~dimitri/doxygen/index.html"> <img src="http://www.stack.nl/~dimitri/doxygen/doxygen.png" alt="doxygen" align="middle" border=0 width=110 height=53> </a> 1.5.9 written by <a href="mailto:dimitri@stack.nl">Dimitri van Heesch</a>, © 1997-2000</small></address> </body> </html>