[dennou-ruby:000898] Re: numfft

[Takeshi Horinouchi <horinout@xxxxxx>]

��Ƿ��Ǥ���

�ⶶ(����)�͡�

> ���� NArray �ѥå������򥤥󥹥ȡ��뤷��NArray �Υ������򸫤ʤ��� 
> fft �ץ�������Τ��Τ��Ƚ񤭴����Ƥ������Ǥ���
> 10����ޤ��餤�� fftpack�Τߴʰ��ǣ����꡼���Ǥ����餤���ʡ�
> �ȻפäƤ���ޤ���
> # ���������ɥ�����Ȥ�񤭻Ϥ�ʤ��㡣

���äȡ���Ȥν��֤ˤĤ��Ƥ⤳����δ�˾������äƤʤ��ä��褦��
���ߤޤ���fftpack �ϸ�󤷤� fftw �򤪴ꤤ���ޤ������Τۤ�����
�Υץ�����ब���ͤ˽�������������Ǥ��ʥץ�����ߥ󥰾�ã�Τ���
�ϡ������ץ�������ޤͤ�Ǥ��ˡ�

�Ǥ�äơ�fftpack �Τۤ��ϼ¤� ruby �ץ�����ߥ󥰤Τߤ�OK�Ȥ���
���Ȥ˵��Ť��ޤ���������ʤ�ǽ餫����äƤ���ȸ��������Ȥ�����
���礦����������Լ�ݤǿ���������ޤ��󡣼¤ϡ�fftpack �� DCL 
�˴ޤޤ�Ƥޤ������ʲ��Υɥ�����Ȥ�ߤ�ĺ����Ф�ʬ����Τ褦
�ˡ�DCL�����äƤ�FFT�ϡ�FFTPACK�������Ƥ��ޤ���
http://www.gfd-dennou.org/arch/dcl/dcl-f77doc/rc1/math2/node3.html
�Ǥ�äơ�����ϴ��� f2c ����� ruby �����ľ�ܤˤϸƤ٤�褦��
�ʤäƤޤ���http://www.gfd-dennou.org/arch/ruby/ �ޤ���
http://www.gfd-dennou.org/arch/ruby/products/rbdcl/ �򻲾Ȥ��Ƥ�
�������ˡ��Ǥ����顢�켡������˴ؤ��Ƥϴ��˻Ȥ���Ϥ��Ǥ���������
�ޤ�������ư��뤫�ɤ�����ǧ���Ƥޤ��󡣤ޤ���fftw �ǹԤ���
�褦��¿�����ؤ��б��������ޤ��Ǥ����������ä� FFTPACK �˴ؤ���
���ꤤ�������Τϡ��ʲ��Τ褦�ˤʤ�ޤ��ʷ����֤��ޤ�����fftw����
�ˤ��ơ����Τ��ȤǼ���Ȥ�Τ��ɤ��Ǥ��ˡ�

  (1) rbdcl (����� Ruby/DCL �������˶᡹��̾�����Ȼפ��ޤ���
      �����դ��Ʋ��������פ���� DCL ��ľ�ܤΥ�åѡ��Ǥ���)
      ��� math2/fftlib �γƴؿ���ȤäƤߤơ������̤��ư��
      �򤹤뤫�ɤ����Τ���롣--- ���餯����פ������Ȼפ��ޤ�
      ���������äƤߤ�ɬ�פ�����ޤ����ʤ��������Ѱ����ϡ��ؿ���
      ����ͤȤʤäƤ�Ϥ��ʤΤ����դ��Ƥ�������������ʳ��ΰ�����
      ���֤�ޤ��Ʊ���Ǥ����Ǥ����顢�ʲ��Τ褦���б��ȤʤäƤ�
      �Ϥ��Ǥ���

	Fotrtran:  
		call sinqi(n,wsave) 
		call sinqf(n,r,wsave)
        Ruby:
		wsave = DCL.sinqi(n)
		DCL.sinqf(n,r,wsave)    # ���� r (input��output) �Ͼ��

      �����γ�ǧ�ϡ����󥹥ȡ��뤹��Ƚ���� rbdcl �� src/ �ǥ�
      �쥯�ȥ꡼�β��γƥѥå������Υ������ե�����򸫤Ƥ���������

  (2) ¿�����б�������� ruby �ΤߤΥץ�����ߥ󥰤Ȥʤ�ޤ���
      ���κݡ�������������Ǽ�ưŪ�˹Ԥ��٤��Ǥ���̵�̤�ʤ�����
      �� N ���Ѳ������Ȥ��Τ߽������Ƥ֤Ȥ������פ�Ǥ��ޤ��ˡ�
      FFTPACK ��Ȥ��Τϡ������˽񤤤��褦�� fftw �ˤϤʤ� SIN �Ѵ���
      COSINE �Ѵ��Υ��ݡ��Ȥ���Ū�ʤΤǡ�SIN��COS �ǻϤޤ�̾����
      �ؿ��ˤĤ��Ƥ��ꤤ���ޤ����������Ѵ�������������ʤ�����
      �ʤΤǡ����Ѵ��κݤˡ�FFTPACK �ν��Ϥ� 4N, 2(N-1) �ʤɤ���
      �����Τ���Υե��������dz�ä����������Ʋ���������Ψ�塢��
      ��� div! ��ȤäƤȤ������Ȥˤʤ�ޤ��͡�
      �ʤ��������ޤǤ��� SIN/COS �Ѵ��ʳ��⥵�ݡ��Ȥ���Τϴ�
      ñ�Ǥ��礦���顢;�Ϥ�����ФĤ��Ǥ˻Ĥ���ä�ĺ�����鹬
      ���Ǥ���

�ʾ塢������δ�˾�����Ƚ����λ��ĺ���ޤ��������ֻ����٤��ʤ�
�ޤ������Ȥ��ͤӤ��ޤ���

��Ƿ�� ��                    horinout@xxxxxx
�������������Ȳʳظ��楻�󥿡�     611-0011 �����Ըޥ���