! { dg-do compile }
|
! { dg-do compile }
|
! Testcase from PR 17713
|
! Testcase from PR 17713
|
module fit_functions
|
module fit_functions
|
implicit none
|
implicit none
|
contains
|
contains
|
subroutine gauss( x, a, y, dy, ma )
|
subroutine gauss( x, a, y, dy, ma )
|
double precision, intent(in) :: x
|
double precision, intent(in) :: x
|
double precision, intent(in) :: a(:)
|
double precision, intent(in) :: a(:)
|
double precision, intent(out) :: y
|
double precision, intent(out) :: y
|
double precision, intent(out) :: dy(:)
|
double precision, intent(out) :: dy(:)
|
integer, intent(in) :: ma
|
integer, intent(in) :: ma
|
end subroutine gauss
|
end subroutine gauss
|
end module fit_functions
|
end module fit_functions
|
|
|
subroutine mrqcof( x, y, sig, ndata, a, ia, ma )
|
subroutine mrqcof( x, y, sig, ndata, a, ia, ma )
|
use fit_functions
|
use fit_functions
|
|
|
implicit none
|
implicit none
|
double precision, intent(in) :: x(:), y(:), sig(:)
|
double precision, intent(in) :: x(:), y(:), sig(:)
|
integer, intent(in) :: ndata
|
integer, intent(in) :: ndata
|
double precision, intent(in) :: a(:)
|
double precision, intent(in) :: a(:)
|
integer, intent(in) :: ia(:), ma
|
integer, intent(in) :: ia(:), ma
|
|
|
integer i
|
integer i
|
double precision yan, dyda(ma)
|
double precision yan, dyda(ma)
|
|
|
do i = 1, ndata
|
do i = 1, ndata
|
call gauss( x(i), a, yan, dyda, ma )
|
call gauss( x(i), a, yan, dyda, ma )
|
end do
|
end do
|
end subroutine mrqcof
|
end subroutine mrqcof
|
|
|
! { dg-final { cleanup-modules "fit_functions" } }
|
! { dg-final { cleanup-modules "fit_functions" } }
|
|
|
