/usr/include/sofa/component/topology/SimpleTesselatedTetraTopologicalMapping.h is in libsofa1-dev 1.0~beta4-6.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 | /******************************************************************************
* SOFA, Simulation Open-Framework Architecture, version 1.0 beta 4 *
* (c) 2006-2009 MGH, INRIA, USTL, UJF, CNRS *
* *
* This library is free software; you can redistribute it and/or modify it *
* under the terms of the GNU Lesser General Public License as published by *
* the Free Software Foundation; either version 2.1 of the License, or (at *
* your option) any later version. *
* *
* This library is distributed in the hope that it will be useful, but WITHOUT *
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or *
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License *
* for more details. *
* *
* You should have received a copy of the GNU Lesser General Public License *
* along with this library; if not, write to the Free Software Foundation, *
* Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. *
*******************************************************************************
* SOFA :: Modules *
* *
* Authors: The SOFA Team and external contributors (see Authors.txt) *
* *
* Contact information: contact@sofa-framework.org *
******************************************************************************/
#ifndef SOFA_COMPONENT_TOPOLOGY_SIMPLETESELATEDTETRATOPOLOGICALMAPPING_H
#define SOFA_COMPONENT_TOPOLOGY_SIMPLETESELATEDTETRATOPOLOGICALMAPPING_H
#include <sofa/core/componentmodel/topology/TopologicalMapping.h>
#include <sofa/defaulttype/Vec.h>
#include <map>
#include <set>
#include <sofa/core/BaseMapping.h>
#include <sofa/component/topology/PointData.h>
#include <sofa/component/topology/EdgeData.h>
#include <sofa/component/topology/TetrahedronData.h>
namespace sofa
{
namespace component
{
namespace topology
{
using namespace sofa::defaulttype;
using namespace sofa::component::topology;
using namespace sofa::core::componentmodel::topology;
using namespace sofa::core;
/**
* This class, called SimpleTesselatedTetraTopologicalMapping, is a specific implementation of the interface TopologicalMapping where :
*
* INPUT TOPOLOGY = TetrahedronSetTopology
* OUTPUT TOPOLOGY = Set of TetrahedronSetTopologies, as the boundary of the INPUT TOPOLOGY
*
* Each tetrahedron in the input Topology will be divided in eight tetrahedrom in the output topology
*
* SimpleTesselatedTetraTopologicalMapping class is templated by the pair (INPUT TOPOLOGY, OUTPUT TOPOLOGY)
*
*/
class SimpleTesselatedTetraTopologicalMapping : public TopologicalMapping
{
public:
/** \brief Constructor.
*
* @param from the topology issuing TopologyChange objects (the "source").
* @param to the topology for which the TopologyChange objects must be translated (the "target").
*/
SimpleTesselatedTetraTopologicalMapping ( In* from=NULL, Out* to=NULL );
/** \brief Destructor.
*
* Does nothing.
*/
virtual ~SimpleTesselatedTetraTopologicalMapping(){};
/** \brief Initializes the target BaseTopology from the source BaseTopology.
*/
virtual void init();
/** \brief Translates the TopologyChange objects from the source to the target.
*
* Translates each of the TopologyChange objects waiting in the source list so that they have a meaning and
* reflect the effects of the first topology changes on the second topology.
*
*/
virtual void updateTopologicalMappingTopDown();
/** \brief Translates the TopologyChange objects from the target to the source.
*
* Translates each of the TopologyChange objects waiting in the source list so that they have a meaning and
* reflect the effects of the second topology changes on the first topology.
*
*/
virtual void updateTopologicalMappingBottomUp();
/// Return true if this mapping is able to propagate topological changes from input to output model
virtual bool propagateFromInputToOutputModel() { return true; }
/// Return true if this mapping is able to propagate topological changes from output to input model
virtual bool propagateFromOutputToInputModel() { return true; }
/// Pre-construction check method called by ObjectFactory.
///
/// This implementation read the object1 and object2 attributes and check
/// if they are compatible with the input and output topology types of this
/// mapping.
template<class T>
static bool canCreate ( T*& obj, core::objectmodel::BaseContext* context, core::objectmodel::BaseObjectDescription* arg )
{
if ( arg->findObject ( arg->getAttribute ( "object1","../.." ) ) == NULL )
context->serr << "Cannot create "<<className ( obj ) <<" as object1 is missing." << context->sendl;
if ( arg->findObject ( arg->getAttribute ( "object2",".." ) ) == NULL )
context->serr << "Cannot create "<<className ( obj ) <<" as object2 is missing." << context->sendl;
if ( arg->findObject ( arg->getAttribute ( "object1","../.." ) ) == NULL || arg->findObject ( arg->getAttribute ( "object2",".." ) ) == NULL )
return false;
BaseMeshTopology* topoIn;
BaseMeshTopology* topoOut;
( dynamic_cast<sofa::core::objectmodel::BaseObject*> ( arg->findObject ( arg->getAttribute ( "object1","../.." ) ) ) )->getContext()->get ( topoIn );
( dynamic_cast<sofa::core::objectmodel::BaseObject*> ( arg->findObject ( arg->getAttribute ( "object2",".." ) ) ) )->getContext()->get ( topoOut );
if ( dynamic_cast<In*> ( topoIn ) == NULL )
return false;
if ( dynamic_cast<Out*> ( topoOut ) == NULL )
return false;
return BaseMapping::canCreate ( obj, context, arg );
}
/// Construction method called by ObjectFactory.
///
/// This implementation read the object1 and object2 attributes to
/// find the input and output topologies of this mapping.
template<class T>
static void create ( T*& obj, core::objectmodel::BaseContext* context, core::objectmodel::BaseObjectDescription* arg )
{
BaseMeshTopology* topoIn=NULL;
BaseMeshTopology* topoOut=NULL;
if ( arg )
{
if ( arg->findObject ( arg->getAttribute ( "object1","../.." ) ) != NULL )
( dynamic_cast<sofa::core::objectmodel::BaseObject*> ( arg->findObject ( arg->getAttribute ( "object1","../.." ) ) ) )->getContext()->get ( topoIn );
if ( arg->findObject ( arg->getAttribute ( "object2",".." ) ) != NULL )
( dynamic_cast<sofa::core::objectmodel::BaseObject*> ( arg->findObject ( arg->getAttribute ( "object2",".." ) ) ) )->getContext()->get ( topoOut );
}
obj = new T (
( arg?dynamic_cast<In*> ( topoIn ) :NULL ),
( arg?dynamic_cast<Out*> ( topoOut ) :NULL ) );
if ( context ) context->addObject ( obj );
if ( ( arg ) && ( arg->getAttribute ( "object1" ) ) )
{
obj->object1.setValue ( arg->getAttribute ( "object1" ) );
arg->removeAttribute ( "object1" );
}
if ( ( arg ) && ( arg->getAttribute ( "object2" ) ) )
{
obj->object2.setValue ( arg->getAttribute ( "object2" ) );
arg->removeAttribute ( "object2" );
}
if ( arg ) obj->parse ( arg );
}
const PointData<int>& getPointMappedFromPoint() const { return pointMappedFromPoint; }
const EdgeData<int>& getPointMappedFromEdge() const { return pointMappedFromEdge; }
const PointData<int>& getPointSource() const { return pointSource; }
protected:
Data< std::string > object1;
Data< std::string > object2;
PointData<int> pointMappedFromPoint; ///< Each point of the input topology is mapped to the same point.
EdgeData<int> pointMappedFromEdge; ///< Each edge of the input topology is mapped to his midpoint.
TetrahedronData< fixed_array<int, 8> > tetrasMappedFromTetra; ///< Each Tetrahedron of the input topology is mapped to the 8 tetrahedrons in which it can be divided.
PointData<int> pointSource; ///< Which input topology element map to a given point in the output topology : 0 -> none, > 0 -> point index + 1, < 0 , - edge index -1
TetrahedronData<int> tetraSource; ///<Which tetra from the input topology map to a given tetra in the output topology (-1 if none)
DataPtr< helper::vector<int> > d_pointMappedFromPoint; ///< Show d_pointMappedFromPoint in the gui for debug
DataPtr< helper::vector<int> > d_pointMappedFromEdge; ///< Show d_pointMappedFromEdge in the gui for debug
DataPtr< helper::vector<int> > d_pointSource; ///< show d_pointSource in the gui for debug
void swapOutputPoints(int i1, int i2);
void removeOutputPoints( const sofa::helper::vector<unsigned int>& tab );
void renumberOutputPoints( const sofa::helper::vector<unsigned int>& tab );
void swapOutputTetras(int i1, int i2);
void removeOutputTetras( const sofa::helper::vector<unsigned int>& tab );
void setPointSource(int i, int source)
{
helper::vector<int> *pointSourceData = pointSource.beginEdit();
helper::vector<int> *pointMappedFromPointData = pointMappedFromPoint.beginEdit();
helper::vector<int>& pointMappedFromEdgeData = *(pointMappedFromEdge.beginEdit());
if (i != -1)
(*pointSourceData)[i] = source;
if (source > 0)
{
(*pointMappedFromPointData)[source-1] = i;
}
else if (source < 0)
{
pointMappedFromEdgeData[-source-1] = i;
}
pointSource.endEdit();
pointMappedFromPoint.endEdit();
pointMappedFromEdge.endEdit();
}
std::set<unsigned int> tetrasToRemove;
void swapInputPoints(int i1, int i2);
void removeInputPoints( const sofa::helper::vector<unsigned int>& tab );
void renumberInputPoints( const sofa::helper::vector<unsigned int>& tab );
void swapInputEdges(int i1, int i2);
void removeInputEdges( const sofa::helper::vector<unsigned int>& tab );
void swapInputTetras(int i1, int i2);
void removeInputTetras( const sofa::helper::vector<unsigned int>& tab );
};
} // namespace topology
} // namespace component
} // namespace sofa
#endif // SOFA_COMPONENT_TOPOLOGY_TETRA2TRIANGLETOPOLOGICALMAPPING_H
|