C++ (Cpp) tf1Plus Exemples

Exemple #1

Fichier : tensorField.C Projet : CFMS/foam-extend-foam-extend-3.2

// This is a nasty hack for flat geometries.  For serious SVD, please use hinv
// (Singular value decomposition and Householder transfromations)
// HJ, 24/Oct/2009
void inv(Field<tensor>& tf, const UList<tensor>& tf1)
{
    if (tf.empty())
    {
        return;
    }

    scalar scale = magSqr(tf1[0]);

    // Fixed terrible hack.  HJ, 20/Jan/2011
    boolList removeCmpts(3);
    removeCmpts[0] = magSqr(tf1[0].xx())/scale < SMALL;
    removeCmpts[1] = magSqr(tf1[0].yy())/scale < SMALL;
    removeCmpts[2] = magSqr(tf1[0].zz())/scale < SMALL;


    if (removeCmpts[0] || removeCmpts[1] || removeCmpts[2])
    {
        tensorField tf1Plus(tf1);

        if (removeCmpts[0])
        {
            tf1Plus += tensor(1,0,0,0,0,0,0,0,0);
        }

        if (removeCmpts[1])
        {
            tf1Plus += tensor(0,0,0,0,1,0,0,0,0);
        }

        if (removeCmpts[2])
        {
            tf1Plus += tensor(0,0,0,0,0,0,0,0,1);
        }

        TFOR_ALL_F_OP_FUNC_F(tensor, tf, =, inv, tensor, tf1Plus)

        if (removeCmpts[0])
        {
            tf -= tensor(1,0,0,0,0,0,0,0,0);
        }

        if (removeCmpts[1])
        {
            tf -= tensor(0,0,0,0,1,0,0,0,0);
        }

        if (removeCmpts[2])
        {
            tf -= tensor(0,0,0,0,0,0,0,0,1);
        }
    }

Exemple #2

Fichier : tensorField.C Projet : 000861/OpenFOAM-2.1.x

void inv(Field<tensor>& tf, const UList<tensor>& tf1)
{
    if (tf.empty())
    {
        return;
    }

    scalar scale = magSqr(tf1[0]);
    Vector<bool> removeCmpts
    (
        magSqr(tf1[0].xx())/scale < SMALL,
        magSqr(tf1[0].yy())/scale < SMALL,
        magSqr(tf1[0].zz())/scale < SMALL
    );

    if (removeCmpts.x() || removeCmpts.y() || removeCmpts.z())
    {
        tensorField tf1Plus(tf1);

        if (removeCmpts.x())
        {
            tf1Plus += tensor(1,0,0,0,0,0,0,0,0);
        }

        if (removeCmpts.y())
        {
            tf1Plus += tensor(0,0,0,0,1,0,0,0,0);
        }

        if (removeCmpts.z())
        {
            tf1Plus += tensor(0,0,0,0,0,0,0,0,1);
        }

        TFOR_ALL_F_OP_FUNC_F(tensor, tf, =, inv, tensor, tf1Plus)

        if (removeCmpts.x())
        {
            tf -= tensor(1,0,0,0,0,0,0,0,0);
        }

        if (removeCmpts.y())
        {
            tf -= tensor(0,0,0,0,1,0,0,0,0);
        }

        if (removeCmpts.z())
        {
            tf -= tensor(0,0,0,0,0,0,0,0,1);
        }
    }