17 template<
unsigned int N,
unsigned int M>
19 { local_bcoords_A_.zeros();
20 local_bcoords_B_.zeros();
28 template<
unsigned int N,
unsigned int M>
34 template<
unsigned int N,
unsigned int M>
36 const arma::vec::fixed<M+1> &lcB,
39 : local_bcoords_A_(lcA), local_bcoords_B_(lcB), dim_A_(dim_A), dim_B_(dim_B)
43 template<
unsigned int N,
unsigned int M>
73 default: idx_B_ = idx_B;
81 template<
unsigned int N,
unsigned int M>
100 template<
unsigned int N,
unsigned int M>
111 template<
unsigned int N,
unsigned int M>
116 arma::vec::fixed< 3 > c;
118 for(
unsigned int i=0; i<N+1; i++)
119 c += local_bcoords_A_[i]**ele.
node(i);
132 template<
unsigned int N,
unsigned int M>
136 dim_A_ == other.
dim_A() &&
137 dim_B_ == other.
dim_B() &&
138 idx_A_ == other.
idx_A() &&
139 idx_B_ == other.
idx_B();
146 os <<
"Local coords on element A(id=" << s.
idx_A_ <<
", dim=" << s.
dim_A_ <<
")" << endl;
148 os <<
"Local coords on element B(id=" << s.
idx_B_ <<
", dim=" << s.
dim_B_ <<
")" << endl;
150 os <<
"Result: " << int(s.
result_) << endl;
NodeAccessor< 3 > node(unsigned int ni) const
Internal auxiliary class represents an intersection point of simplex<N> and simplex<M>.
unsigned int idx_B() const
Returns the index of Simplex<M>.
bool topology_equal(const IntersectionPointAux< N, M > &other) const
Returns true, if other intersection point has the same topology.
void set_coordinates(const arma::vec::fixed< N+1 > &lcA, const arma::vec::fixed< M+1 > &lcB)
Setter for coordinates.
unsigned int dim_A_
Dimension of the object A of intersection. Equal N, by default.
IntersectionResult result() const
Result: 0 - negative sign, 1 - positive sign, 2 - degenerate (zero for all sides),...
void set_result(IntersectionResult result)
Setter orientation flag.
void set_topology(unsigned int idx_A, unsigned int dim_A, unsigned int idx_B, unsigned int dim_B)
Setter for topology data. See description of idx_A_ and dim_A_.
unsigned int idx_A_
Local indices of element objects that intersects.
unsigned int dim_A() const
Returns dimension of object A.
unsigned int dim_B_
Dimension of the object B of intersection. Equal M, by default.
unsigned int idx_A() const
Returns the index of Simplex<N>.
arma::vec::fixed< N+1 > local_bcoords_A_
Barycentric coordinates of an IP on simplex<N>.
IntersectionPointAux()
Default constructor.
unsigned int dim_B() const
Returns dimension of object B.
void clear()
Resets the object to default values.
arma::vec::fixed< M+1 > local_bcoords_B_
Barycentric coordinates of an IP on simplex<M>.
IntersectionResult result_
Intersection result according to Plucker products.
const arma::vec::fixed< N+1 > & local_bcoords_A() const
Returns barycentric coordinates in the Simplex<N>.
const arma::vec::fixed< M+1 > & local_bcoords_B() const
Returns barycentric coordinates in the Simplex<M>.
arma::vec::fixed< 3 > coords(ElementAccessor< 3 > ele) const
Computes real coordinates of IP, given the element ele in which IP lies.
IntersectionPointAux< M, N > switch_objects() const
Switches the object A and B.
static const IdxVector<(InDim >OutDim ? InDim+1 :dim-InDim) > interact(TInteraction< OutDim, InDim > interaction, bool inv=false)
template ostream & operator<<< 1, 3 >(ostream &os, const IntersectionPointAux< 1, 3 > &s)
template ostream & operator<<< 3, 2 >(ostream &os, const IntersectionPointAux< 3, 2 > &s)
template ostream & operator<<< 2, 1 >(ostream &os, const IntersectionPointAux< 2, 1 > &s)
template ostream & operator<<< 3, 1 >(ostream &os, const IntersectionPointAux< 3, 1 > &s)
template ostream & operator<<< 1, 2 >(ostream &os, const IntersectionPointAux< 1, 2 > &s)
ostream & operator<<(ostream &os, const IntersectionPointAux< N, M > &s)
template ostream & operator<<< 2, 3 >(ostream &os, const IntersectionPointAux< 2, 3 > &s)
template ostream & operator<<< 2, 2 >(ostream &os, const IntersectionPointAux< 2, 2 > &s)
Internal class representing intersection point.
Class RefElement defines numbering of vertices, sides, calculation of normal vectors etc.