PolynomialPatch.h

Go to the documentation of this file.

/* 
 *  Copyright (c) 2015, Chris Rogers
 *  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 STFC 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 OWNER 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.
 */

#include "Fields/Interpolation/SquarePolynomialVector.h"
#include "Fields/Interpolation/VectorMap.h"
#include "Fields/Interpolation/Mesh.h"

#ifndef PolynomialPatch_hh
#define PolynomialPatch_hh

namespace interpolation {

class PolynomialPatch : public VectorMap {
  public:
    PolynomialPatch(Mesh* grid_points_,
                    Mesh* validity_region,
                    std::vector<SquarePolynomialVector*> polynomials_);

    PolynomialPatch();

    ~PolynomialPatch();

    VectorMap* clone() const; //copy function

    virtual void function(const double* point, double* value) const;

    inline unsigned int getPointDimension() const {return point_dimension_;}

    inline unsigned int getValueDimension() const {return value_dimension_;}

    SquarePolynomialVector* getPolynomialVector(const double* point) const;

    Mesh* getMesh() const {return validity_region_;}

    std::vector<SquarePolynomialVector*> getPolynomials() const {return points_;}

  private:
    Mesh* validity_region_;
    Mesh* grid_points_;
    std::vector<SquarePolynomialVector*> points_;
    unsigned int point_dimension_;
    unsigned int value_dimension_;

    PolynomialPatch(const PolynomialPatch& rhs);
    PolynomialPatch& operator=(const PolynomialPatch& rhs);
};

}

#endif // PolynomialPatch_hh