1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
! { dg-do compile }
! Testcase from PR 17713
module fit_functions
implicit none
contains
subroutine gauss( x, a, y, dy, ma )
double precision, intent(in) :: x
double precision, intent(in) :: a(:)
double precision, intent(out) :: y
double precision, intent(out) :: dy(:)
integer, intent(in) :: ma
end subroutine gauss
end module fit_functions
subroutine mrqcof( x, y, sig, ndata, a, ia, ma )
use fit_functions
implicit none
double precision, intent(in) :: x(:), y(:), sig(:)
integer, intent(in) :: ndata
double precision, intent(in) :: a(:)
integer, intent(in) :: ia(:), ma
integer i
double precision yan, dyda(ma)
do i = 1, ndata
call gauss( x(i), a, yan, dyda, ma )
end do
end subroutine mrqcof
|
