poly_int Subroutine

  1. pf_quadrature.f90
  2. pf_mod_quadrature
  3. poly_int

public subroutine poly_int(p, n)

Subroutine to integrate polynomial (in place)

Arguments

Type IntentOptional AttributesName
real(kind=pfqp), intent(inout) :: p(0:n)
integer, intent(in), value:: n

Called by

proc~~poly_int~~CalledByGraph proc~poly_int poly_int proc~sdc_qmats sdc_qmats proc~sdc_qmats->proc~poly_int proc~pf_quadrature pf_quadrature proc~pf_quadrature->proc~sdc_qmats proc~pf_init_sdcmats pf_init_sdcmats proc~pf_init_sdcmats->proc~pf_quadrature proc~pf_level_setup pf_level_setup proc~pf_level_setup->proc~pf_init_sdcmats proc~pf_pfasst_setup pf_pfasst_setup proc~pf_pfasst_setup->proc~pf_level_setup

Contents

Source Code

poly_int

Source Code

  subroutine poly_int(p, n) 
    integer, intent(in),   value :: n
    real(pfqp),       intent(inout) :: p(0:n)

    integer  :: j
    real(pfqp) :: pp(0:n)

    pp = 0.0_pfqp

    do j = 0, n-1
       pp(j+1) = p(j) / (j+1)
    end do

    p = pp
  end subroutine poly_int