Perform the 1D integration.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(integration_class_1d), | intent(inout) | :: | me | |||
real(kind=wp), | intent(out) | :: | ans | |||
integer, | intent(out) | :: | ierr | |||
real(kind=wp), | intent(out) | :: | err |
subroutine integrate_1d (me, ans, ierr, err)
implicit none
class(integration_class_1d),intent(inout) :: me
real(wp),intent(out) :: ans
integer,intent(out) :: ierr
real(wp),intent(out) :: err
!call the low-level routine:
call me%dgauss_generic(me%a, me%b, me%tol, ans, ierr, err)
end subroutine integrate_1d