Commit 36158b13 authored by Thomas Pitiot 's avatar Thomas Pitiot

working

parent 96c538dc
......@@ -606,9 +606,7 @@ void ParticleCell3D<PFP>::edgeState(const VEC3& current)
Geom::Plane3D<typename PFP::REAL> pl = Algo::Surface::Geometry::facePlane<PFP>(m,d,position);
VEC3 p3 = pl.normal()+this->m_position;
Geom::Plane3D<typename PFP::REAL> plOrtho(this->m_position,current,p3);
VEC3 e(position[m.phi1(d)]-position[d]);
Geom::intersectionPlaneRay(plOrtho,this->m_position,current-this->m_position,this->m_position);
Geom::intersectionPlaneRay(plOrtho,this->m_position,current-this->m_position,this->m_position);
edgeState(current);
return;
......@@ -639,8 +637,6 @@ void ParticleCell3D<PFP>::edgeState(const VEC3& current)
Geom::Plane3D<typename PFP::REAL> pl = Algo::Surface::Geometry::facePlane<PFP>(m,d,position);
VEC3 p3 = pl.normal()+this->m_position;
Geom::Plane3D<typename PFP::REAL> plOrtho(this->m_position,current,p3);
VEC3 e(position[m.phi1(d)]-position[d]);
Geom::intersectionPlaneRay(plOrtho,this->m_position,current-this->m_position,this->m_position);
edgeState(current);
return;
......@@ -667,7 +663,6 @@ void ParticleCell3D<PFP>::edgeState(const VEC3& current)
Geom::Plane3D<typename PFP::REAL> pl = Algo::Surface::Geometry::facePlane<PFP>(m,d,position);
VEC3 p3 = pl.normal()+this->m_position;
Geom::Plane3D<typename PFP::REAL> plOrtho(this->m_position,current,p3);
VEC3 e(position[m.phi1(d)]-position[d]);
Geom::intersectionPlaneRay(plOrtho,this->m_position,current-this->m_position,this->m_position);
......
......@@ -369,8 +369,6 @@ void ParticleCell3DMemo<PFP>::edgeState(const VEC3& current, CellMarkerMemo<MAP,
Geom::Plane3D<typename PFP::REAL> pl = Algo::Surface::Geometry::facePlane<PFP>(this->m,this->d,this->position);
VEC3 p3 = pl.normal()+this->m_position;
Geom::Plane3D<typename PFP::REAL> plOrtho(this->m_position,current,p3);
VEC3 e(this->position[this->m.phi1(this->d)]-this->position[this->d]);
Geom::intersectionPlaneRay(plOrtho,this->m_position,current-this->m_position,this->m_position);
edgeState(current, memo_cross);
return;
......@@ -397,9 +395,7 @@ void ParticleCell3DMemo<PFP>::edgeState(const VEC3& current, CellMarkerMemo<MAP,
Geom::Plane3D<typename PFP::REAL> pl = Algo::Surface::Geometry::facePlane<PFP>(this->m,this->d,this->position);
VEC3 p3 = pl.normal()+this->m_position;
Geom::Plane3D<typename PFP::REAL> plOrtho(this->m_position,current,p3);
VEC3 e(this->position[this->m.phi1(this->d)]-this->position[this->d]);
Geom::intersectionPlaneRay(plOrtho,this->m_position,current-this->m_position,this->m_position);
Geom::intersectionPlaneRay(plOrtho,this->m_position,current-this->m_position,this->m_position);
this->m_positionFace = this->m_position;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment