Nektar::NekMeshUtils::OptiFace Class Reference

#include <OptimiseFunctions.h>

Inheritance diagram for Nektar::NekMeshUtils::OptiFace:

Collaboration diagram for Nektar::NekMeshUtils::OptiFace:

Public Member Functions
	OptiFace (Array< OneD, Array< OneD, NekDouble > > a, std::map< std::pair< int, int >, NekDouble > w, std::set< std::pair< int, int > > sp, CADSurfSharedPtr su)
	~OptiFace ()
NekDouble	F (Array< OneD, NekDouble > xitst)
DNekMat	dF (Array< OneD, NekDouble > xitst)
Array< OneD, NekDouble >	Getxi ()
Array< OneD, NekDouble >	Getli ()
Array< OneD, NekDouble >	Getui ()
void	Update (Array< OneD, NekDouble > xinew)
Array< OneD, Array< OneD, NekDouble > >	GetSolution ()
Public Member Functions inherited from Nektar::NekMeshUtils::OptiObj
	OptiObj ()
virtual	~OptiObj ()

Private Attributes
CADSurfSharedPtr	s
std::map< std::pair< int, int > , NekDouble >	z
std::set< std::pair< int, int > >	spring
Array< OneD, Array< OneD, NekDouble > >	uv
int	ni
int	np
int	nq

Friends
class	MemoryManager< OptiFace >

Detailed Description

Definition at line 84 of file OptimiseFunctions.h.

Constructor & Destructor Documentation

Nektar::NekMeshUtils::OptiFace::OptiFace ( Array< OneD, Array< OneD, NekDouble > >  a, std::map< std::pair< int, int >, NekDouble >  w, std::set< std::pair< int, int > >  sp, CADSurfSharedPtr  su  )

inline

Definition at line 89 of file OptimiseFunctions.h.

References ni, np, nq, s, spring, uv, and z.

    {
        uv     = a;
        z      = w;
        spring = sp;
        s      = su;
        np     = uv.num_elements();
        nq     = 0.5 * (-1 + sqrt(1 + 8 * np));
        ni     = (nq - 3) * (nq - 2) / 2;
    };

Nektar::NekMeshUtils::OptiFace::~OptiFace ( )

inline

Definition at line 103 of file OptimiseFunctions.h.

{};

Member Function Documentation

DNekMat Nektar::NekMeshUtils::OptiFace::dF ( Array< OneD, NekDouble >  xitst )

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 417 of file OptimiseFunctions.cpp.

References Nektar::NekMeshUtils::Dot(), Nektar::iterator, and Nektar::NekMeshUtils::Take().

{
    Array<OneD, Array<OneD, NekDouble> > val = uv;
    for (int i = np - ni; i < np; i++)
    {
        val[i][0] = xitst[(i - np + ni) * 2 + 0];
        val[i][1] = xitst[(i - np + ni) * 2 + 1];
    }

    DNekMat ret(ni * 2, 1, 0.0);

    vector<Array<OneD, NekDouble> > r, dru, drv;

    for (int i = 0; i < val.num_elements(); i++)
    {
        Array<OneD, NekDouble> ri(3), du(3), dv(3);
        Array<OneD, NekDouble> d1 = s->D1(val[i]);
        for (int i = 0; i < 3; i++)
        {
            ri[i] = d1[i];
            du[i] = d1[i + 3];
            dv[i] = d1[i + 6];
        }
        r.push_back(ri);
        dru.push_back(du);
        drv.push_back(dv);
    }

    for (int i = 0; i < ni * 2; i++) // for each of the varibles
    {
        int var = floor(i / 2) + np - ni;
        int tp  = i % 2; // 0 is a u 1 is a v

        set<pair<int, int> >::iterator it;
        for (it = spring.begin(); it != spring.end();
             it++) // for each of the springs
        {
            Array<OneD, NekDouble> dr1, dr2;

            if ((*it).first == var)
            {
                if (tp == 0)
                {
                    dr1 = dru[(*it).first];
                }
                else
                {
                    dr1 = drv[(*it).first];
                }
            }
            else
            {
                dr1 = Array<OneD, NekDouble>(3, 0.0);
            }

            if ((*it).second == var)
            {
                if (tp == 0)
                {
                    dr2 = dru[(*it).second];
                }
                else
                {
                    dr2 = drv[(*it).second];
                }
            }
            else
            {
                dr2 = Array<OneD, NekDouble>(3, 0.0);
            }

            ret(i, 0) +=
                2.0 / z[(*it)] *
                Dot(Take(r[(*it).first], r[(*it).second]), Take(dr1, dr2));
        }
    }
    return ret;
}

NekDouble Nektar::NekMeshUtils::OptiFace::F ( Array< OneD, NekDouble >  xitst )

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 395 of file OptimiseFunctions.cpp.

References Nektar::iterator, and Nektar::NekMeshUtils::Take().

{
    Array<OneD, Array<OneD, NekDouble> > val = uv;
    for (int i = np - ni; i < np; i++)
    {
        val[i][0] = xitst[(i - np + ni) * 2 + 0];
        val[i][1] = xitst[(i - np + ni) * 2 + 1];
    }

    NekDouble ret = 0.0;
    set<pair<int, int> >::iterator it;
    for (it = spring.begin(); it != spring.end(); it++)
    {
        Array<OneD, NekDouble> dis =
            Take(s->P(uv[(*it).first]), s->P(uv[(*it).second]));
        NekDouble norm = dis[0] * dis[0] + dis[1] * dis[1] + dis[2] * dis[2];
        ret += norm / z[(*it)];
    }

    return ret;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::OptiFace::Getli ( )

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 371 of file OptimiseFunctions.cpp.

{
    Array<OneD, NekDouble> li(ni * 2);
    Array<OneD, NekDouble> bnds = s->GetBounds();
    for (int i = np - ni; i < np; i++)
    {
        li[2 * (i - np + ni) + 0] = bnds[0];
        li[2 * (i - np + ni) + 1] = bnds[2];
    }
    return li;
}

Array<OneD, Array<OneD, NekDouble> > Nektar::NekMeshUtils::OptiFace::GetSolution ( )

inline

Definition at line 112 of file OptimiseFunctions.h.

References uv.

    {
        return uv;
    };

Array< OneD, NekDouble > Nektar::NekMeshUtils::OptiFace::Getui ( )

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 383 of file OptimiseFunctions.cpp.

{
    Array<OneD, NekDouble> ui(ni * 2);
    Array<OneD, NekDouble> bnds = s->GetBounds();
    for (int i = np - ni; i < np; i++)
    {
        ui[2 * (i - np + ni) + 0] = bnds[1];
        ui[2 * (i - np + ni) + 1] = bnds[3];
    }
    return ui;
}

Array< OneD, NekDouble > Nektar::NekMeshUtils::OptiFace::Getxi ( )

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 360 of file OptimiseFunctions.cpp.

{
    Array<OneD, NekDouble> ret(ni * 2);
    for (int i = np - ni; i < np; i++)
    {
        ret[2 * (i - np + ni) + 0] = uv[i][0];
        ret[2 * (i - np + ni) + 1] = uv[i][1];
    }
    return ret;
}

void Nektar::NekMeshUtils::OptiFace::Update ( Array< OneD, NekDouble >  xinew )

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 496 of file OptimiseFunctions.cpp.

{
    for (int i = np - ni; i < np; i++)
    {
        uv[i][0] = xinew[2 * (i - np + ni) + 0];
        uv[i][1] = xinew[2 * (i - np + ni) + 1];
    }
}

Friends And Related Function Documentation

friend class MemoryManager< OptiFace >

friend

Definition at line 87 of file OptimiseFunctions.h.

Member Data Documentation

int Nektar::NekMeshUtils::OptiFace::ni

private

Definition at line 122 of file OptimiseFunctions.h.

Referenced by OptiFace().

int Nektar::NekMeshUtils::OptiFace::np

private

Definition at line 122 of file OptimiseFunctions.h.

Referenced by OptiFace().

int Nektar::NekMeshUtils::OptiFace::nq

private

Definition at line 122 of file OptimiseFunctions.h.

Referenced by OptiFace().

CADSurfSharedPtr Nektar::NekMeshUtils::OptiFace::s

private

Definition at line 115 of file OptimiseFunctions.h.

Referenced by OptiFace().

std::set<std::pair<int, int> > Nektar::NekMeshUtils::OptiFace::spring

private

Definition at line 120 of file OptimiseFunctions.h.

Referenced by OptiFace().

Array<OneD, Array<OneD, NekDouble> > Nektar::NekMeshUtils::OptiFace::uv

private

Definition at line 121 of file OptimiseFunctions.h.

Referenced by GetSolution(), and OptiFace().

std::map<std::pair<int, int>, NekDouble> Nektar::NekMeshUtils::OptiFace::z

private

Definition at line 119 of file OptimiseFunctions.h.

Referenced by OptiFace().