OpenWalnut  1.4.0
WDataSetDTI.cpp
1 //---------------------------------------------------------------------------
2 //
3 // Project: OpenWalnut ( http://www.openwalnut.org )
4 //
5 // Copyright 2009 OpenWalnut Community, BSV@Uni-Leipzig and CNCF@MPI-CBS
6 // For more information see http://www.openwalnut.org/copying
7 //
8 // This file is part of OpenWalnut.
9 //
10 // OpenWalnut is free software: you can redistribute it and/or modify
11 // it under the terms of the GNU Lesser General Public License as published by
12 // the Free Software Foundation, either version 3 of the License, or
13 // (at your option) any later version.
14 //
15 // OpenWalnut is distributed in the hope that it will be useful,
16 // but WITHOUT ANY WARRANTY; without even the implied warranty of
17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 // GNU Lesser General Public License for more details.
19 //
20 // You should have received a copy of the GNU Lesser General Public License
21 // along with OpenWalnut. If not, see <http://www.gnu.org/licenses/>.
22 //
23 //---------------------------------------------------------------------------
24 
25 #include "WDataSetDTI.h"
26 
27 // prototype instance as singleton
28 boost::shared_ptr< WPrototyped > WDataSetDTI::m_prototype = boost::shared_ptr< WPrototyped >();
29 
30 
31 WDataSetDTI::WDataSetDTI( boost::shared_ptr< WValueSetBase > newValueSet, boost::shared_ptr< WGrid > newGrid )
32  : WDataSetSingle( newValueSet, newGrid )
33 {
34  WAssert( newValueSet, "No value set given." );
35  WAssert( newGrid, "No grid given." );
36  WAssert( newValueSet->size() == newGrid->size(), "Number of values unequal number of positions in grid." );
37  WAssert( newValueSet->order() == 1, "The value set does not contain vectors." );
38  WAssert( newValueSet->dimension() == 6, "The size of the vetors does not match symmetric matrices: must be equal to 6." );
39 }
40 
42 {
43 }
44 
45 WDataSetSingle::SPtr WDataSetDTI::clone( boost::shared_ptr< WValueSetBase > newValueSet ) const
46 {
47  return WDataSetSingle::SPtr( new WDataSetDTI( newValueSet, getGrid() ) );
48 }
49 
50 WDataSetSingle::SPtr WDataSetDTI::clone( boost::shared_ptr< WGrid > newGrid ) const
51 {
52  return WDataSetSingle::SPtr( new WDataSetDTI( getValueSet(), newGrid ) );
53 }
54 
56 {
58 }
59 
61 {
62  boost::shared_ptr< WValueSet< float > > values = boost::dynamic_pointer_cast< WValueSet< float > >( m_valueSet );
63  WAssert( values, "The value set of a WDataSetDTI must be a WValueSet< float >, nothing else!" );
64  return WTensorSym< 2, 3, float >( values->getWValue( index ) );
65 }
66