Nektar::NekMeshUtils::linesource Struct Reference

#include <Octree.h>

Collaboration diagram for Nektar::NekMeshUtils::linesource:

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 51 of file Octree.h.

Constructor & Destructor Documentation

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

inline

Definition at line 55 of file Octree.h.

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

Member Function Documentation

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

inline

Definition at line 93 of file Octree.h.

Referenced by withinRange().

    {
        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]));
    }

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

inline

Definition at line 63 of file Octree.h.

References Length().

    {
        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;
        }
    }

Member Data Documentation

NekDouble Nektar::NekMeshUtils::linesource::delta

Definition at line 54 of file Octree.h.

NekDouble Nektar::NekMeshUtils::linesource::R

Definition at line 54 of file Octree.h.

Array<OneD, NekDouble> Nektar::NekMeshUtils::linesource::x1

Definition at line 53 of file Octree.h.

Array<OneD, NekDouble> Nektar::NekMeshUtils::linesource::x2

Definition at line 53 of file Octree.h.