[dennou-ruby:000071] Re: gtrb & multi-D array

[Numaguti Atusi <numa@xxxxxx>]

�¸��Ǥ���

 > ���ȡ����Ǥ�����
 > http://www.jus.or.jp/workshop/ruby/

�ʤ�ۤɤʤ�ۤɡ��Ǥ���

�Ȥ����ǡ��ޤȤ�Ρ�
  ���ԡ��ɤ򤢤����ٳ��ݤ��ʤ���Фʤ�ʤ�
        ʸ������̸����񤤤ƸƤӽФ��Ȥ����ƥ��֥������Ȥ����������˿�
        �ͤ򿨤���ˡ��ɬ��
�Ϥɤ��������Ȥʤ�Ǥ��礦�� 

���֤򿩤��Τϡ���������ä��ǡ����γ����Ǥ��Ф�����ƽ�����
���뤳�Ȥ�¿���Τǡ��������Ƥ�(�ݥ����Ϥ���)ľ�� fortran 
�ʤ� C �ʤ�Ǥ�����褦�ʥ��饹���äƤ����С����Τ褦��
�褯�Ȥ����ƽ����򤹤�롼�����񤤤Ƥ����Ф褤�櫓�Ǥ���͡�
�ǡ���Ʊ�Τλ�§�黻��ʿ�Ѥ�Ȥ롤��ɸ�Ѵ� etc. etc.

¿ʬ��ꤽ���ʤ��Ȥ�®��Ū������ˤʤ�Τϡ�Ruby �ǽ񤤤�
�����Υ᥽�åɤ���������ǡ��������ƤϤ᤿�����ʤɤ�
�����Ȥ��Ǥ��礦����

������󡤤����᥽�åɤ����Ƥ�ñ��ʱ黻�ʤ顤���α黻��
Numeric ���饹�Ǥʤ�����ҤΤ褦�����󥯥饹��Ŭ�Ѥ����
�褤�櫓�Ǥ��������Ƚ�Ǥʤɤ����äƤ���ȶ줷����

�Ĥޤꡤ
   def t2theta(p)
       self * (P0/p)**K
   end
���餤�δ�ñ�ʼ��򤢤ƤϤ��Τʤ�С�
����γ�����Ʊ�Τλ�§�黻��Ƽ�����ؿ����餤�� C �� fortran ��
�Ѱդ��Ƥ����Ф褤��

��������
   def t2theta2(p)
       if (p>=0) 
         self * (P0/p)**K
       else
         self * (-p/P0)**K
       end
   end
�Ȥ����Ƚ�Ǥ����äƤ���Ⱦ����줷���ʤ롤�Ǥ⡤�ʤ�Ȥ��Ϥʤ롥
���Ȥ��С������Ѵ����ơ�
   def t2theta2(p)
     s = sign(1.0,p)
     self * ( (P0/p)**K * s + (-p/P0)**K * (1.0-s) )
   end
����ϡ���äȤ��ä��褯
   def t2theta2(p)
       ( self * (P0/p)**K  ).if_ge(p,0.0) |
       ( self * (-p/P0)**K )
   end
�Ȥ�
( if_ge(p,0.0) �� p���б����Ǥ� >=0 �ʤ� self, ����ʳ��ʤ� nil���֤���
  |(x) �� self �����Ǥ� nil �ʤ� x���б����Ǥ򡤤���ʳ��ʤ� self ���֤���
  �ºݤˤ� C �Ȥ��ǽ񤫤줿�᥽�å� )
     
�Ǥ⡤���湽¤��ʣ�������񤷤��Ǥ��͡��³������롥
�ΤΥ٥��ȥ�׻�����ȤäƤ���ߤ���������

-- 
�¸�  �ء��̳�ƻ�����ر��ϵ�Ķ��ʳظ���� �絤���η��Ķ��ʳ��칶
tel:011-706-2365 fax:011-706-4865 ��060-0810 ���ڻ��̶���10����5����  
numa@xxxxxx