From 09ddc3a9f3f52a142527fccab41a6981955b8beb Mon Sep 17 00:00:00 2001
From: kamischi <karl-michael.schindler@web.de>
Date: Sun, 25 Dec 2022 16:11:14 +0100
Subject: [PATCH] Update fun.f90
suggestion for optimization
---
source/f90/fun.f90 | 3 +++
1 file changed, 3 insertions(+)
diff --git a/source/f90/fun.f90 b/source/f90/fun.f90
index d01e398..73df08d 100644
--- a/source/f90/fun.f90
+++ b/source/f90/fun.f90
@@ -22,5 +22,8 @@ double precision function fun(phi)
fun = dsqrt((1.0d0 - elleps + elleps * sinphi**2) * &
(1.0d0 - sinpsi * sinphi) * &
(1.0d0 + sinpsi * sinphi))
+! optimized but untested version
+! sinphi2 = dsin(phi)**2
+! fun = dsqrt((1.0d0 + elleps*(sinphi2 - 1.0d0)) * (1.0d0 - sinpsi**2 * sinphi2))
return
end function fun
--
GitLab