<!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: itkFEMLinearSystemWrapperDenseVNL.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>itkFEMLinearSystemWrapperDenseVNL.h</h1><a href="itkFEMLinearSystemWrapperDenseVNL_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: itkFEMLinearSystemWrapperDenseVNL.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-29 21:55:14 $</span> <a name="l00007"></a>00007 <span class="comment"> Version: $Revision: 1.8 $</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 __itkFEMLinearSystemWrapperDenseVNL_h</span> <a name="l00019"></a>00019 <span class="preprocessor"></span><span class="preprocessor">#define __itkFEMLinearSystemWrapperDenseVNL_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_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/algo/vnl_svd.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_1LinearSystemWrapperDenseVNL.html">00038</a> <span class="keyword">class </span><a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html" title="LinearSystemWrapper class that uses VNL numeric library functions to define a sparse...">LinearSystemWrapperDenseVNL</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_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#789c06d332fcca846debae70f7cbfe36">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_1LinearSystemWrapperDenseVNL.html#a2677215da334c14ca9a53de9fe2f100">00049</a> <span class="keyword">typedef</span> <a class="code" href="classvnl__matrix.html">vnl_matrix<Float></a> <a class="code" href="classvnl__matrix.html">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_1LinearSystemWrapperDenseVNL.html#ede308765a12877f436645711c18525a">00052</a> <span class="keyword">typedef</span> std::vector< MatrixRepresentation* > <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#ede308765a12877f436645711c18525a">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_1LinearSystemWrapperDenseVNL.html#9b46161e838a87a9ad61cacf9ca80bf6">00055</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#9b46161e838a87a9ad61cacf9ca80bf6">LinearSystemWrapperDenseVNL</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_1LinearSystemWrapperDenseVNL.html#eaec1520a75c93fe83d2b0ce6d6d1ff7">~LinearSystemWrapperDenseVNL</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_1LinearSystemWrapperDenseVNL.html#f5fc1f284fc734e184a8a7fd424bc247">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_1LinearSystemWrapperDenseVNL.html#f313a5bfb8d488a09274a1aee685b731">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_1LinearSystemWrapperDenseVNL.html#6d02b75625f38bded37a1a1bbe1565fd">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_1LinearSystemWrapperDenseVNL.html#5394930cb8a58f6875e7b2a25dab3a3b">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_1LinearSystemWrapperDenseVNL.html#7a324c98fbb021094fbda218aea0bfd1">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_1LinearSystemWrapperDenseVNL.html#b4874db4b7e967ecb1313b977b8a68e0">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_1LinearSystemWrapperDenseVNL.html#97e3d7edc72b0236c006c99bfcf2786b">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_1LinearSystemWrapperDenseVNL.html#f27a6f76d299e16117f3974b12d7eee8">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_1LinearSystemWrapperDenseVNL.html#3273c3e94c40762950fa7e1fbbf8c413">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_1LinearSystemWrapperDenseVNL.html#7a4bc38e40e10f66185b06dfb0a2081a">00068</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#7a4bc38e40e10f66185b06dfb0a2081a">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_1LinearSystemWrapperDenseVNL.html#dc898deb62b0785b903d891ee61a5820">00072</a> <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">Float</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#dc898deb62b0785b903d891ee61a5820">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_1LinearSystemWrapperDenseVNL.html#786652f6246b5ff373f6d12624a98fcf">00073</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#786652f6246b5ff373f6d12624a98fcf">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_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#193233bc0d0c2a51a03d40ff0eed32ed">00074</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#193233bc0d0c2a51a03d40ff0eed32ed">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_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#3f5e5e3f25a84ed6a632fd2437909226">00075</a> <span class="keyword">virtual</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">Float</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#3f5e5e3f25a84ed6a632fd2437909226">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_1LinearSystemWrapperDenseVNL.html#f1a8855e7baa2c37873b590d26d6bfcb">00076</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#f1a8855e7baa2c37873b590d26d6bfcb">SetVectorValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#05b5c9637d67f66f82d3754b211d203c">00077</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#05b5c9637d67f66f82d3754b211d203c">AddVectorValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">Float</a> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#3d55b539cd3fcaff8f267f06bcc30df8">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_1LinearSystemWrapperDenseVNL.html#81bc7a1dc04d75a236bfaa3d04e3a955">00079</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#81bc7a1dc04d75a236bfaa3d04e3a955">SetSolutionValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#4145cb7eee23cb70f9d32b844aa22dfa">00080</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#4145cb7eee23cb70f9d32b844aa22dfa">AddSolutionValue</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#a76c4cd4a7fbf9bcca2414593e8cc8e3">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_1LinearSystemWrapperDenseVNL.html#f56d7476a8a2b2d50f083bc1affaedb6">ScaleMatrix</a>(<a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">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_1LinearSystemWrapperDenseVNL.html#247d5b658e4fc6dadfcc2a21a0f90bfb">ScaleVector</a>(<a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">Float</a> scale, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex); <a name="l00086"></a>00086 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#23f14bf62b91fec304bc819af29bc1f4">ScaleSolution</a>(<a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#652eeb019082bd8753c94baff1187505">Float</a> scale, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> solutionIndex); <a name="l00087"></a>00087 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#dec1a32c81e198c56cb81eb88a1c5b8f">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="l00088"></a>00088 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#e390a45b0693cec7dd21cffacf84be5a">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="l00089"></a>00089 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#f0c4ea350d45ec10dd75edfbac615c01">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="l00090"></a>00090 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#046c9e9b11fb87e688f80038e8191dd4">CopySolution2Vector</a>(<span class="keywordtype">unsigned</span> solutionIndex, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> vectorIndex); <a name="l00091"></a>00091 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#fb407c203d99881a4d8b8cdf5a2a6b55">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="l00092"></a>00092 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#bf896ffe8ae29edd83801748b17b6f0a">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="l00093"></a>00093 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#ef56a583e429c703263eaa2ed8057bdf">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="l00094"></a>00094 <a name="l00095"></a>00095 <span class="keyword">private</span>: <a name="l00096"></a>00096 <a name="l00098"></a>00098 <span class="comment">//std::vector< vnl_sparse_matrix<Float>* > *m_Matrices;</span> <a name="l00099"></a>00099 <a class="code" href="classitk_1_1fem_1_1LinearSystemWrapperDenseVNL.html#ede308765a12877f436645711c18525a">MatrixHolder</a> *m_Matrices; <a name="l00100"></a>00100 <a name="l00102"></a>00102 std::vector< vnl_vector<Float>* > *m_Vectors; <a name="l00103"></a>00103 <a name="l00105"></a>00105 std::vector< vnl_vector<Float>* > *m_Solutions; <a name="l00106"></a>00106 <a name="l00107"></a>00107 }; <a name="l00108"></a>00108 <a name="l00109"></a>00109 }} <span class="comment">// end namespace itk::fem</span> <a name="l00110"></a>00110 <a name="l00111"></a>00111 <span class="preprocessor">#endif</span> <a name="l00112"></a>00112 <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>