LobattoQuadrature.hpp

Go to the documentation of this file.

//==============================================================================
//
//                       Seismo Virtual Laboratory
//             Module for Serial and Parallel Analysis of seismic 
//         wave propagation and soil-structure interaction simulation
//         Copyright (C) 2018-2021, The California Institute of Technology
//                         All Rights Reserved.
//
// Commercial use of this program without express permission of the California
// Institute of Technology, is strictly  prohibited. See  file "COPYRIGHT"  in
// main  directory  for  information on  usage  and  redistribution, and for a
// DISCLAIMER OF ALL WARRANTIES.
//
//==============================================================================
//
// Written by:
//   Danilo S. Kusanovic (dkusanov@caltech.edu)
//   Elnaz E. Seylabi    (elnaze@unr.edu)
//
// Supervised by:
//   Domniki M. Asimaki  (domniki@caltech.edu)
//
// References : 
//  [1] 
//
// Description:
//------------------------------------------------------------------------------

#ifndef _LOBATTOQUADRATURE_HPP_
#define _LOBATTOQUADRATURE_HPP_

#include <string>
#include <Eigen/Dense>

#include "QuadratureRule.hpp"

class LobattoQuadrature : public QuadratureRule {

    public:
        LobattoQuadrature(std::string name, unsigned int nQp);

        ~LobattoQuadrature();

        unsigned int GetNumberOfQuadraturePoints();

        void GetQuadraturePoints(std::string name, Eigen::VectorXd &wi, Eigen::MatrixXd &xi);

    private:
        unsigned int nQuadraturePoints;

        unsigned int nOrderQuadrature;

        void SetLineQuadraturePoints(Eigen::VectorXd &wi, Eigen::MatrixXd &xi);

        void SetTriaQuadraturePoints(Eigen::VectorXd &wi, Eigen::MatrixXd &xi);

        void SetQuadQuadraturePoints(Eigen::VectorXd &wi, Eigen::MatrixXd &xi);

        void SetTetraQuadraturePoints(Eigen::VectorXd &wi, Eigen::MatrixXd &xi);

        void SetHexaQuadraturePoints(Eigen::VectorXd &wi, Eigen::MatrixXd &xi);
};

#endif