poly_diff Subroutine

  1. pf_quadrature.f90
  2. pf_mod_quadrature
  3. poly_diff

public subroutine poly_diff(p, n)

Subroutine to differentiate polynomial (in place)

Arguments

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

Called by

proc~~poly_diff~~CalledByGraph proc~poly_diff poly_diff proc~sdc_qnodes sdc_qnodes proc~sdc_qnodes->proc~poly_diff proc~pf_quadrature pf_quadrature proc~pf_quadrature->proc~sdc_qnodes 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_diff

Source Code

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

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

    pp = 0.0_pfqp

    do j = 1, n
       pp(j-1) = j * p(j)
    end do

    p = pp
  end subroutine poly_diff