InterpolatedProfile.h

Go to the documentation of this file.

//- ****************************************************************************
//-
//- Copyright 2009 National Technology & Engineering Solutions of Sandia, LLC
//- (NTESS). Under the terms of Contract DE-NA0003525 with NTESS, the U.S.
//- Government retains certain rights in this software.
//-
//- BSD Open Source License
//- All rights reserved.
//-
//- Redistribution and use in source and binary forms, with or without
//- modification, are permitted provided that the following conditions are met:
//-
//-   1. Redistributions of source code must retain the above copyright notice,
//-      this list of conditions and the following disclaimer.
//-
//-   2. Redistributions in binary form must reproduce the above copyright
//-      notice, this list of conditions and the following disclaimer in the
//-      documentation and/or other materials provided with the distribution.
//-
//-   3. Neither the name of the copyright holder nor the names of its
//-      contributors may be used to endorse or promote products derived from
//-      this software without specific prior written permission.
//-
//- THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
//- AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
//- IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
//- ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
//- LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
//- CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
//- SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
//- INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
//- CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
//- ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
//- POSSIBILITY OF SUCH DAMAGE.
//-
//- ****************************************************************************
 
#ifndef InterpolatedProfile_H
#define InterpolatedProfile_H
 
// **** _SYSTEM INCLUDES_ ******************************************************
 
#include <vector>
 
using namespace std;
 
// **** _LOCAL INCLUDES_ *******************************************************
 
#include "SLBMGlobals.h"
#include "GridProfile.h"
 
// **** _BEGIN SLBM NAMESPACE_ **************************************************
 
namespace slbm {
 
class SLBM_EXP_IMP InterpolatedProfile
{
public:
 
    InterpolatedProfile();
 
    InterpolatedProfile(Grid& grid, Location& location);
 
    InterpolatedProfile(const InterpolatedProfile& other);
 
    virtual ~InterpolatedProfile();
 
    InterpolatedProfile& operator=(const InterpolatedProfile& other);
 
    bool operator==(const InterpolatedProfile& other);
 
    bool operator!=(const InterpolatedProfile& other) {return ! (*this == other);};
 
    int getNCoefficients() { return nodes.size(); }
 
    vector<GridProfile*>& getNodes() { return nodes; }
 
    GridProfile* getNode(const int& i) { return nodes[i]; }
 
    vector<int>& getNodeIds() { return nodeIds; }
 
    int getNodeId(const int& i) { return nodeIds[i]; }
 
    vector<double>& getCoefficients() { return coefficients; };
 
    double getCoefficient(const int& i) { return coefficients[i]; }
 
     void getNodeIds(int* nodeIds, int& size);
 
    void getCoefficients(double* coeff, int& size);
 
    void getWeights(int* nodeIds, double* coeff, int&size);
 
    void interpRadius(const int& k, double& radius);
 
    void interpDepth(const int& k, double& depth);
 
    void interpVelocity(const int& type, const int& k, double& velocity);
 
    void interpGradient(const int& type, double& gradient);
 
    bool isActiveProfile();
 
    virtual size_t memSize();
 
    static int getClassCount();
 
protected:
 
    static int interpolatedProfileClassCount;
 
    vector<GridProfile*> nodes;
 
    vector<int> nodeIds;
 
    vector<double> coefficients;
 
};
 
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
//  INLINE FUNCTIONS
//
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
inline void InterpolatedProfile::getWeights(int* nodeids, double* coeff, int& size)
{
    size = nodes.size();
    for (int i=0; i<(int)nodes.size(); i++)
    {
        nodeids[i] = nodes[i]->getNodeId();
        coeff[i] = coefficients[i];
    }
}
 
inline void InterpolatedProfile::getNodeIds(int* nodeids, int& size)
{
    size = nodes.size();
    for (int i=0; i<(int)nodes.size(); i++)
        nodeids[i] = nodes[i]->getNodeId();
}
 
inline void InterpolatedProfile::getCoefficients(double* coeff, int& size)
{
    size = nodes.size();
    for (int i=0; i<(int)nodes.size(); i++)
        coeff[i] = coefficients[i];
}
 
inline void InterpolatedProfile::interpRadius(const int& k, double& radius)
{
    radius = 0.;
    for (int i=0; i<(int)nodes.size(); i++) if (nodes[i])
        radius += nodes[i]->getInterfaceRadius(k) * coefficients[i];
}
 
inline void InterpolatedProfile::interpDepth(const int& k, double& depth)
{
    depth = 0.;
    for (int i=0; i<(int)nodes.size(); i++) if (nodes[i])
        depth += nodes[i]->getInterfaceDepth(k) * coefficients[i];
}
 
inline void InterpolatedProfile::interpVelocity(const int& type, const int& k, double& velocity)
{
    velocity = 0.;
    for (int i=0; i<(int)nodes.size(); i++) if (nodes[i])
        velocity += nodes[i]->getVelocity(type,k) * coefficients[i];
}
 
inline void InterpolatedProfile::interpGradient(const int& type, double& gradient)
{
    gradient = 0.;
    for (int i=0; i<(int)nodes.size(); i++) if (nodes[i])
        gradient += nodes[i]->getMantleGradient(type) * coefficients[i];
}
 
inline bool InterpolatedProfile::isActiveProfile()
{
    for (int i=0; i<(int)nodes.size(); i++)
        if (nodes[i]->getActiveNodeId() < 0)
            return false;
    return true;
}
 
} // end slbm namespace
 
#endif // InterpolatedProfile.h