Skip to content

Commit fc99e3c

Browse files
Copilotivan-pi
andauthored
chore: address review feedback on nonautonomous API changes
Agent-Logs-Url: https://github.com/ivan-pi/stiff3/sessions/67f692d9-43f0-4337-94d3-193b5edcdf71 Co-authored-by: ivan-pi <21085643+ivan-pi@users.noreply.github.com>
1 parent a5bd3c1 commit fc99e3c

2 files changed

Lines changed: 3 additions & 3 deletions

File tree

src/stiff3_solver.f90

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -306,8 +306,8 @@ subroutine fun_aug(naug,yi,fi)
306306
integer, intent(in) :: naug
307307
real(wp), intent(in) :: yi(naug)
308308
real(wp), intent(inout) :: fi(naug)
309-
call fun(n,yi(1),yi(2:),fi(2:))
310309
fi(1) = 1.0_wp
310+
call fun(n,yi(1),yi(2:),fi(2:))
311311
end subroutine
312312

313313
subroutine jac_aug(naug,yi,dfi)

test/nonautonomous.f90

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -4,7 +4,7 @@ program nonautonomous
44

55
implicit none
66

7-
integer, parameter :: n = 1, nout = 100
7+
integer, parameter :: n = 1, nprint = 100
88
real(wp), parameter :: atol = 5.0e-6_wp
99
real(wp) :: y(n), w(n), x0, x1, h0, eps, y_exact
1010

@@ -15,7 +15,7 @@ program nonautonomous
1515
h0 = 1.0e-3_wp
1616
eps = 1.0e-8_wp
1717

18-
call stiff3(n,fun,jac,dfdx,out,nout,x0,x1,h0,eps,w,y)
18+
call stiff3(n,fun,jac,dfdx,out,nprint,x0,x1,h0,eps,w,y)
1919

2020
y_exact = 0.5_wp*(sin(x1) - cos(x1) + exp(-x1))
2121
if (abs(y(1)-y_exact) > atol) then

0 commit comments

Comments
 (0)