Nektar::NekMeshUtils::OptiEdge Class Reference

#include <OptimiseFunctions.h>

Inheritance diagram for Nektar::NekMeshUtils::OptiEdge:

Collaboration diagram for Nektar::NekMeshUtils::OptiEdge:

Public Member Functions
	OptiEdge (Array< OneD, NekDouble > a, Array< OneD, NekDouble > dis, CADObjSharedPtr ob)
	~OptiEdge ()
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, NekDouble >	GetSolution ()
Public Member Functions inherited from Nektar::NekMeshUtils::OptiObj
	OptiObj ()
virtual	~OptiObj ()

Private Attributes
CADObjSharedPtr	o
Array< OneD, NekDouble >	z
Array< OneD, NekDouble >	all

Friends
class	MemoryManager< OptiEdge >

Detailed Description

Definition at line 49 of file OptimiseFunctions.h.

Constructor & Destructor Documentation

Nektar::NekMeshUtils::OptiEdge::OptiEdge ( Array< OneD, NekDouble >  a, Array< OneD, NekDouble >  dis, CADObjSharedPtr  ob  )

inline

Definition at line 54 of file OptimiseFunctions.h.

References all, o, and z.

    {
        all = a;
        z   = dis;
        o   = ob;
    };

Nektar::NekMeshUtils::OptiEdge::~OptiEdge ( )

inline

Definition at line 63 of file OptimiseFunctions.h.

{};

Member Function Documentation

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

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 236 of file OptimiseFunctions.cpp.

References Nektar::NekMeshUtils::curve, Nektar::NekMeshUtils::Dot(), Nektar::NekMeshUtils::surf, and Nektar::NekMeshUtils::Take().

{
    Array<OneD, NekDouble> val(all.num_elements());

    if (o->GetType() == curve)
    {
        val[0] = all[0];
        for (int i = 0; i < xitst.num_elements(); i++)
        {
            val[i + 1] = xitst[i];
        }
        val[all.num_elements() - 1] = all[all.num_elements() - 1];
    }
    else if (o->GetType() == surf)
    {
        val[0] = all[0];
        val[1] = all[1];
        for (int i = 0; i < xitst.num_elements(); i++)
        {
            val[i + 2] = xitst[i];
        }
        val[all.num_elements() - 2] = all[all.num_elements() - 2];
        val[all.num_elements() - 1] = all[all.num_elements() - 1];
    }

    DNekMat ret;

    if (o->GetType() == curve)
    {
        CADCurveSharedPtr c = boost::dynamic_pointer_cast<CADCurve>(o);
        vector<Array<OneD, NekDouble> > r;
        vector<Array<OneD, NekDouble> > dr;

        for (int i = 0; i < all.num_elements(); i++)
        {
            Array<OneD, NekDouble> ri(3), dri(3);
            Array<OneD, NekDouble> d2 = c->D2(val[i]);
            for (int j = 0; j < 3; j++)
            {
                ri[j]  = d2[j];
                dri[j] = d2[j + 3];
            }
            r.push_back(ri);
            dr.push_back(dri);
        }

        DNekMat J(all.num_elements() - 2, 1, 0.0);
        for (int i = 0; i < all.num_elements() - 2; i++)
        {
            J(i, 0) =
                2.0 / (z[i + 1] - z[i]) * Dot(dr[i + 1], Take(r[i + 1], r[i])) -
                2.0 / (z[i + 2] - z[i + 1]) *
                    Dot(dr[i + 1], Take(r[i + 2], r[i + 1]));
        }

        ret = J;
    }
    else if (o->GetType() == surf)
    {
        CADSurfSharedPtr s = boost::dynamic_pointer_cast<CADSurf>(o);
        // need to organise the all array
        Array<OneD, Array<OneD, NekDouble> > uv(val.num_elements() / 2);
        for (int i = 0; i < val.num_elements() / 2; i++)
        {
            uv[i]    = Array<OneD, NekDouble>(2);
            uv[i][0] = val[i * 2 + 0];
            uv[i][1] = val[i * 2 + 1];
        }

        vector<Array<OneD, NekDouble> > r;
        vector<Array<OneD, NekDouble> > dru, drv;
        for (int i = 0; i < uv.num_elements(); i++)
        {
            Array<OneD, NekDouble> ri(3), drui(3), drvi(3);
            Array<OneD, NekDouble> d2 = s->D1(uv[i]);
            for (int j = 0; j < 3; j++)
            {
                ri[j]   = d2[j];
                drui[j] = d2[j + 3];
                drvi[j] = d2[j + 6];
            }
            r.push_back(ri);
            dru.push_back(drui);
            drv.push_back(drvi);
        }

        DNekMat J(2 * (uv.num_elements() - 2), 1, 0.0);
        for (int i = 0; i < uv.num_elements() - 2; i++)
        {
            J(2 * i + 0, 0) = 2.0 / (z[i + 1] - z[i]) *
                                  Dot(dru[i + 1], Take(r[i + 1], r[i])) +
                              2.0 / (z[i + 2] - z[i + 1]) *
                                  Dot(dru[i + 1], Take(r[i + 1], r[i + 2]));

            J(2 * i + 1, 0) = 2.0 / (z[i + 1] - z[i]) *
                                  Dot(drv[i + 1], Take(r[i + 1], r[i])) +
                              2.0 / (z[i + 2] - z[i + 1]) *
                                  Dot(drv[i + 1], Take(r[i + 1], r[i + 2]));
        }

        ret = J;
    }

    return ret;
}

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

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 176 of file OptimiseFunctions.cpp.

References Nektar::NekMeshUtils::curve, Nektar::NekMeshUtils::surf, and Nektar::NekMeshUtils::Take().

{
    Array<OneD, NekDouble> val(all.num_elements());

    if (o->GetType() == curve)
    {
        val[0] = all[0];
        for (int i = 0; i < xitst.num_elements(); i++)
        {
            val[i + 1] = xitst[i];
        }
        val[all.num_elements() - 1] = all[all.num_elements() - 1];
    }
    else if (o->GetType() == surf)
    {
        val[0] = all[0];
        val[1] = all[1];
        for (int i = 0; i < xitst.num_elements(); i++)
        {
            val[i + 2] = xitst[i];
        }
        val[all.num_elements() - 2] = all[all.num_elements() - 2];
        val[all.num_elements() - 1] = all[all.num_elements() - 1];
    }

    NekDouble ret = 0.0;
    if (o->GetType() == curve)
    {
        CADCurveSharedPtr c = boost::dynamic_pointer_cast<CADCurve>(o);

        for (int i = 0; i < all.num_elements() - 1; i++)
        {
            Array<OneD, NekDouble> dis = Take(c->P(val[i + 1]), c->P(val[i]));
            NekDouble norm =
                dis[0] * dis[0] + dis[1] * dis[1] + dis[2] * dis[2];
            ret += norm / (z[i + 1] - z[i]);
        }
    }
    else if (o->GetType() == surf)
    {
        CADSurfSharedPtr s = boost::dynamic_pointer_cast<CADSurf>(o);
        // need to organise the val array
        Array<OneD, Array<OneD, NekDouble> > uv(val.num_elements() / 2);
        for (int i = 0; i < val.num_elements() / 2; i++)
        {
            uv[i]    = Array<OneD, NekDouble>(2);
            uv[i][0] = val[i * 2 + 0];
            uv[i][1] = val[i * 2 + 1];
        }
        for (int i = 0; i < uv.num_elements() - 1; i++)
        {
            Array<OneD, NekDouble> dis = Take(s->P(uv[i + 1]), s->P(uv[i]));
            NekDouble norm =
                dis[0] * dis[0] + dis[1] * dis[1] + dis[2] * dis[2];
            ret += norm / (z[i + 1] - z[i]);
        }
    }
    return ret;
}

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

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 102 of file OptimiseFunctions.cpp.

References ASSERTL0, Nektar::NekMeshUtils::curve, Nektar::NekMeshUtils::surf, and Nektar::NekMeshUtils::vert.

{
    Array<OneD, NekDouble> li;
    Array<OneD, NekDouble> bnds;
    switch (o->GetType())
    {
        case curve:
            li   = Array<OneD, NekDouble>(all.num_elements() - 2);
            bnds = boost::dynamic_pointer_cast<CADCurve>(o)->Bounds();
            for (int i = 1; i < all.num_elements() - 1; i++)
            {
                li[i - 1] = bnds[0];
            }
            break;

        case surf:
            li   = Array<OneD, NekDouble>(all.num_elements() - 4);
            bnds = boost::dynamic_pointer_cast<CADSurf>(o)->GetBounds();
            for (int i = 2; i < all.num_elements() - 2; i++)
            {
                if (i % 2 == 0)
                {
                    li[i - 2] = bnds[0];
                }
                else
                {
                    li[i - 2] = bnds[2];
                }
            }
            break;

        case vert:
            ASSERTL0(false, "Should not be able to pass vert");
    }
    return li;
}

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

inline

Definition at line 72 of file OptimiseFunctions.h.

References all.

    {
        return all;
    };

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

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 139 of file OptimiseFunctions.cpp.

References ASSERTL0, Nektar::NekMeshUtils::curve, Nektar::NekMeshUtils::surf, and Nektar::NekMeshUtils::vert.

{
    Array<OneD, NekDouble> ui;
    Array<OneD, NekDouble> bnds;
    switch (o->GetType())
    {
        case curve:
            ui   = Array<OneD, NekDouble>(all.num_elements() - 2);
            bnds = boost::dynamic_pointer_cast<CADCurve>(o)->Bounds();
            for (int i = 1; i < all.num_elements() - 1; i++)
            {
                ui[i - 1] = bnds[1];
            }
            break;

        case surf:
            ui   = Array<OneD, NekDouble>(all.num_elements() - 4);
            bnds = boost::dynamic_pointer_cast<CADSurf>(o)->GetBounds();
            for (int i = 2; i < all.num_elements() - 2; i++)
            {
                if (i % 2 == 0)
                {
                    ui[i - 2] = bnds[1];
                }
                else
                {
                    ui[i - 2] = bnds[3];
                }
            }
            break;

        case vert:
            ASSERTL0(false, "Should not be able to pass vert");
    }
    return ui;
}

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

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 74 of file OptimiseFunctions.cpp.

References ASSERTL0, Nektar::NekMeshUtils::curve, Nektar::NekMeshUtils::surf, and Nektar::NekMeshUtils::vert.

{
    Array<OneD, NekDouble> xi;

    switch (o->GetType())
    {
        case curve:
            xi = Array<OneD, NekDouble>(all.num_elements() - 2);
            for (int i = 1; i < all.num_elements() - 1; i++)
            {
                xi[i - 1] = all[i];
            }
            break;

        case surf:
            xi = Array<OneD, NekDouble>(all.num_elements() - 4);
            for (int i = 2; i < all.num_elements() - 2; i++)
            {
                xi[i - 2] = all[i];
            }
            break;

        case vert:
            ASSERTL0(false, "Should not be able to pass vert");
    }
    return xi;
}

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

virtual

Reimplemented from Nektar::NekMeshUtils::OptiObj.

Definition at line 342 of file OptimiseFunctions.cpp.

References Nektar::NekMeshUtils::curve, and Nektar::NekMeshUtils::surf.

{
    if (o->GetType() == curve)
    {
        for (int i = 0; i < xinew.num_elements(); i++)
        {
            all[i + 1] = xinew[i];
        }
    }
    else if (o->GetType() == surf)
    {
        for (int i = 0; i < xinew.num_elements(); i++)
        {
            all[i + 2] = xinew[i];
        }
    }
}

Friends And Related Function Documentation

friend class MemoryManager< OptiEdge >

friend

Definition at line 52 of file OptimiseFunctions.h.

Member Data Documentation

Array<OneD, NekDouble> Nektar::NekMeshUtils::OptiEdge::all

private

Definition at line 80 of file OptimiseFunctions.h.

Referenced by GetSolution(), and OptiEdge().

CADObjSharedPtr Nektar::NekMeshUtils::OptiEdge::o

private

Definition at line 75 of file OptimiseFunctions.h.

Referenced by OptiEdge().

Array<OneD, NekDouble> Nektar::NekMeshUtils::OptiEdge::z

private

Definition at line 79 of file OptimiseFunctions.h.

Referenced by OptiEdge().