Line data    Source code

       1              : 
       2              : 
       3              : namespace ippl {
       4              :     template <typename T, unsigned NumNodes1D, typename ElementType>
       5          252 :     MidpointQuadrature<T, NumNodes1D, ElementType>::MidpointQuadrature(
       6              :         const ElementType& ref_element)
       7          252 :         : Quadrature<T, NumNodes1D, ElementType>(ref_element) {
       8          252 :         this->degree_m = 1;
       9              : 
      10          252 :         this->a_m = 0.0;
      11          252 :         this->b_m = 1.0;
      12              : 
      13          252 :         computeNodesAndWeights();
      14          252 :     }
      15              : 
      16              :     template <typename T, unsigned NumNodes1D, typename ElementType>
      17          252 :     void MidpointQuadrature<T, NumNodes1D, ElementType>::computeNodesAndWeights() {
      18          252 :         const T segment_length = (this->b_m - this->a_m) / NumNodes1D;
      19              : 
      20          252 :         this->weights_m = Vector<T, NumNodes1D>(segment_length);
      21              : 
      22          252 :         this->integration_nodes_m = Vector<T, NumNodes1D>();
      23          576 :         for (unsigned i = 0; i < NumNodes1D; ++i) {
      24          324 :             this->integration_nodes_m[i] = 0.5 * segment_length + i * segment_length + this->a_m;
      25              :         }
      26          252 :     }
      27              : 
      28              : }  // namespace ippl