1 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
5 ! Authors: Paul N. Swarztrauber and Richard A. Valent
7 ! $Id: sintmi.f,v 1.2 2004/06/15 21:29:20 rodney Exp $
9 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
11 SUBROUTINE SINTMI (N, WSAVE, LENSAV, IER)
12 INTEGER N, LENSAV, IER
17 IF (LENSAV .LT. N/2 + N + INT(LOG(REAL(N))) +4) THEN
19 CALL XERFFT ('SINTMI', 3)
29 WSAVE(K) = 2.*SIN(K*DT)
31 LNSV = NP1 + INT(LOG(REAL(NP1))) +4
32 CALL RFFTMI (NP1, WSAVE(NS2+1), LNSV, IER1)
35 CALL XERFFT ('SINTMI',-5)