Fortran90からC＃​​への変換の問題

Question

Fortran90コードをC＃に変換しています。私はFortran77についていくつかの知識を持っていますが、Fortran90には慣れていません。私は翻訳する方法がわからない次のコード行を実行しました。私はこのことを考えてい

C1 = real(product((/(-1,i1=1,m-1)/))*product((/(i1,i1=2,m)/))) 

のように変換する必要があります。

int product1 = -1; int product2 = 1; 
for (int i1 = 1 ; i1 <= (m-1); i1++) 
{ 
    product1 *= -1; 
} 
for (int i2 = 2, i2 <= m; i2++) 
{ 
    product2 *= i2; 
} 
float C1 = (float)(product1 * product2); 

私の不確実性が暗黙のは、配列を初期化するためのループの構築を行うが存在するという事実から生じます。すなわち

A = (/2*I, I = 1,5/) 

が、私は、問題のFortranの文のように使用単語「製品」を見たことがありません。私は知っているベクトルまたは行列の乗算PRODUCTと呼ばれる本質的な関数がありますが、 "製品"は私が働いているコードの配列ではないintrisic関数の構文PRODUCTは非常に明確に私の文はこの関数を使用していません。

洞察または助けがあれば幸いです。ありがとうございました。

Answer 1

あなたが部品を分解し、それらを印刷する場合は、これらは、単に簡潔なベクトル化用語を使用して作成された用語ですがわかります：

! given: (/(expr, start, end)/) 
! 
! (/(-1, i1=1, m-1)/) = vector, -1 repeated m-1 times 
! 
! (/(i1, i1=2, m)/) = vector, 2..m 
! 
! Both are vectors with m-1 terms 

注意すべき他の事はproduct()は3つの引数を必要としないことです。 2番目の引数（使用する次元）と3番目の引数（配列マスク）は必要ありません。

この時点で、最初の製品が実際-1m-1であり、第二製品がm!であることが明らかとなります。

ので、適切な（必ずしも効率的ではない）翻訳は次のようになります。

// product((/(-1,i1=1,m-1)/)) => -1^m-1 
double i = (m % 2 == 0 ? -1 : 1); 

// product((/(i1,i1=2,m)/)) => m! 
double mfact = i; 
for (int jj = 2; jj < m; ++jj) 
{ 
    mfact *= jj; 
} // C1 = mfact; 

、簡潔F90と "精神" が近いが、ほとんど効果的な：

double i = (m % 2 == 0 ? -1 : 1); 
double C1 = Enumerable.Range(2, m) 
         .Aggregate(i, (x, y) => x * y);