qmidpoint.hh

Go to the documentation of this file.

/*!
 *
 * Copyright (C) 2015 Technical University of Liberec.  All rights reserved.
 * 
 * This program is free software; you can redistribute it and/or modify it under
 * the terms of the GNU General Public License version 3 as published by the
 * Free Software Foundation. (http://www.gnu.org/licenses/gpl-3.0.en.html)
 * 
 * This program is distributed in the hope that it will be useful, but WITHOUT
 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
 * FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
 *
 * 
 * @file    qmidpoint.hh
 * @brief   Midpoint rule qudrature.
 * @author  Pavel Exner
 */

#ifndef QMIDPOINT_HH_
#define QMIDPOINT_HH_

#include "system/global_defs.h"
#include "quadrature/quadrature.hh"
#include "mesh/point.hh"


/** @brief Class representing midpoint rule, with uniformly distributed points of the same weight.
 */
class QMidpoint : public Quadrature<1> {
public:
    /// Empty constructor
    QMidpoint(const unsigned int n_quadrature_points){
        
        double qweight = 1.0/n_quadrature_points;
        this->weights.resize(n_quadrature_points,qweight);
        this->quadrature_points.resize(n_quadrature_points);
        for(unsigned int q=0; q < n_quadrature_points; q++)
            this->set_point(q, arma::vec({0.5*qweight + q*qweight}));
    }
};

#endif // QMIDPOINT_HH_