From 73745a4e82f5f9c3fa07f49cac0439c21ba50829 Mon Sep 17 00:00:00 2001
From: nsha <nsha@lle.rochester.edu>
Date: Fri, 16 Feb 2024 16:47:00 -0500
Subject: [PATCH] add interfaces for sine transform

---
 src/fftpack.f90 | 18 ++++++++++++++++++
 1 file changed, 18 insertions(+)

diff --git a/src/fftpack.f90 b/src/fftpack.f90
index cea79f9..ed14799 100644
--- a/src/fftpack.f90
+++ b/src/fftpack.f90
@@ -20,6 +20,8 @@ module fftpack
     public :: dct_t1i, dct_t1
     public :: dct_t23i, dct_t2, dct_t3
 
+    public :: dsinti, dsint
+
     public :: rk
 
     interface
@@ -177,6 +179,22 @@ pure subroutine dcost(n, x, wsave)
             real(kind=rk), intent(in) :: wsave(*)
         end subroutine dcost
 
+
+        pure subroutine dsinti(n, wsave)
+            import rk
+            integer, intent(in) :: n
+            real(kind=rk), intent(out) :: wsave(*)
+        end subroutine dsinti
+
+
+        pure subroutine dsint(n, x, wsave)
+            import rk
+            integer, intent(in) :: n
+            real(kind=rk), intent(inout) :: x(*)
+            real(kind=rk), intent(in) :: wsave(*)
+        end subroutine dsint
+
+
         !> Version: experimental
         !>
         !> Integer frequency values involved in complex FFT.