/usr/include/psurface/ContactMapping.h is in libpsurface-dev 2.0.0-2.
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 | #ifndef CONTACT_MAPPING_HH
#define CONTACT_MAPPING_HH
#include <vector>
#include <iostream>
// Check for VC9 / VS2008 with installed feature pack.
#if defined(_MSC_VER) && (_MSC_VER>=1500)
#if defined(_CPPLIB_VER) && _CPPLIB_VER>=505
#include <array>
#else
#error Please install the Visual Studio 2008 SP1 for TR1 support.
#endif
#else
#include <tr1/array>
#endif
#include "StaticVector.h"
#include "PSurface.h"
#include "IntersectionPrimitive.h"
#include "IntersectionPrimitiveCollector.h"
namespace psurface {
template <int dimworld, class ctype>
struct DirectionFunction;
template <int dim, class ctype>
class ContactMapping {};
template <class ctype>
class ContactMapping<2,ctype>
{
public:
void build(const std::vector<std::tr1::array<ctype,2> >& coords1, ///< The vertices of the first surface as \f$x_0 ,y_0 ,z_0, x_1, y_1, z_1 ...\f$
const std::vector<std::tr1::array<int,2> >& tri1, ///< The triangles of the first surface
const std::vector<std::tr1::array<ctype,2> >& coords2, ///< The vertices of the second surface
const std::vector<std::tr1::array<int,2> >& tri2,
const DirectionFunction<2,ctype>* domainDirection = NULL,
const DirectionFunction<2,ctype>* targetDirection = NULL
);
void getOverlaps(std::vector<IntersectionPrimitive<1,ctype> >& overlaps)
{
IntersectionPrimitiveCollector<ctype>::collect(&psurface_, overlaps);
}
// /////////////////////////////////////////////
private:
void computeDiscreteDomainDirections(const DirectionFunction<2,ctype>* direction,
std::vector<StaticVector<ctype,2> >& normals);
void computeDiscreteTargetDirections(const std::vector<std::tr1::array<int,2> >& elements,
const DirectionFunction<2,ctype>* direction,
std::vector<StaticVector<ctype,2> >& normals);
PSurface<1,ctype> psurface_;
};
template <class ctype>
class ContactMapping<3,ctype>
{
public:
void build(const std::vector<std::tr1::array<ctype,3> >& coords1, ///< The vertices of the first surface as \f$x_0 ,y_0 ,z_0, x_1, y_1, z_1 ...\f$
const std::vector<std::tr1::array<int,3> >& tri1, ///< The triangles of the first surface
const std::vector<std::tr1::array<ctype,3> >& coords2, ///< The vertices of the second surface
const std::vector<std::tr1::array<int,3> >& tri2, ///< The triangles of the second surface
const DirectionFunction<3,ctype>* domainDirection = NULL,
const DirectionFunction<3,ctype>* targetDirection = NULL
);
void getOverlaps(std::vector<IntersectionPrimitive<2,ctype> >& overlaps) {
IntersectionPrimitiveCollector<ctype>::collect(&psurface_, overlaps);
}
private:
PSurface<2,ctype> psurface_;
};
} // namespace psurface
#endif
|