******************************************************************
  
  subroutine pda_rffti(n,wsave)
  
  ****************************************************************
  
  subroutine pda_rffti initializes the array wsave which is used in
  both pda_rfftf and pda_rfftb. the prime factorization of n together with
  a tabulation of the trigonometric functions are computed and
  stored in wsave.
  
  input parameter
  
  n       the length of the sequence to be transformed.
  
  output parameter
  
  wsave   a work array which must be dimensioned at least 2*n+15.
          the same work array can be used for both pda_rfftf and pda_rfftb
          as long as n remains unchanged. different wsave arrays
          are required for different values of n. the contents of
          wsave must not be changed between calls of pda_rfftf or pda_rfftb.