#include <Octree.h>

Public Member Functions
	linesource (Array< OneD, NekDouble > p1, Array< OneD, NekDouble > p2, NekDouble r, NekDouble d)

bool	withinRange (Array< OneD, NekDouble > p)

NekDouble	Length ()

Public Attributes
Array< OneD, NekDouble >	x1

Array< OneD, NekDouble >	x2

NekDouble	R

NekDouble	delta

Detailed Description

Definition at line 50 of file Octree.h.

Constructor & Destructor Documentation

◆ linesource()

Nektar::NekMeshUtils::linesource::linesource	(	Array< OneD, NekDouble >	p1,
		Array< OneD, NekDouble >	p2,
		NekDouble	r,
		NekDouble	d
	)

inline

Definition at line 54 of file Octree.h.

         : x1(p1), x2(p2), R(r), delta(d)
     {
     }

Member Function Documentation

◆ Length()

NekDouble Nektar::NekMeshUtils::linesource::Length ( )

inline

Definition at line 92 of file Octree.h.

     {
         return sqrt((x1[0] - x2[0]) * (x1[0] - x2[0]) +
                     (x1[1] - x2[1]) * (x1[1] - x2[1]) +
                     (x1[2] - x2[2]) * (x1[2] - x2[2]));
     }

References x1, and x2.

Referenced by withinRange().

◆ withinRange()

bool Nektar::NekMeshUtils::linesource::withinRange ( Array< OneD, NekDouble > p )

inline

Definition at line 62 of file Octree.h.

     {
         Array<OneD, NekDouble> Le(3), Re(3), s(3);
         for (int i = 0; i < 3; i++)
         {
             Le[i] = p[i] - x1[i];
             Re[i] = p[i] - x2[i];
             s[i]  = x2[i] - x1[i];
         }
         Array<OneD, NekDouble> dev(3);
         dev[0] = Le[1] * Re[2] - Re[1] * Le[2];
         dev[1] = Le[2] * Re[0] - Re[2] * Le[0];
         dev[2] = Le[0] * Re[1] - Re[0] * Le[1];
  
         NekDouble dist =
             sqrt(dev[0] * dev[0] + dev[1] * dev[1] + dev[2] * dev[2]) / Length();
  
         NekDouble t = -1.0 * ((x1[0] - p[0]) * s[0] + (x1[1] - p[1]) * s[1] +
                               (x1[2] - p[2]) * s[2]) / Length() / Length();
  
         if (dist < R && !(t > 1) && !(t < 0))
         {
             return true;
         }
         else
         {
             return false;
         }
     }