From 0e546af43ab00a51cdb635dd4e835f03d1155cec Mon Sep 17 00:00:00 2001
From: kamischi <karl-michael.schindler@web.de>
Date: Fri, 5 Jul 2024 18:55:41 +0200
Subject: [PATCH] Update sicot.f90
tested and works.
---
source/f90/sicot.f90 | 7 +------
1 file changed, 1 insertion(+), 6 deletions(-)
diff --git a/source/f90/sicot.f90 b/source/f90/sicot.f90
index 32e8ed4..1e34dd2 100644
--- a/source/f90/sicot.f90
+++ b/source/f90/sicot.f90
@@ -86,14 +86,9 @@ subroutine sicot(f, m, h, x0)
enddo
j = n
- do i = 2, n
- j = j - 1
+ do j = n - 1, 1, -1
f(j) = f(j) + f(j + 1)
enddo
-! alternative, but not yet tested:
-! do j = n - 1, 1, -1
-! f(j) = f(j) + f(j + 1)
-! enddo
msign = -iabs(m)
call sintr(f, msign, h)
return
--
GitLab