fortranのcombination_with_replacement関数ですか？

Question

私はy総コストに収まるx製品のすべての可能性を得るために、小さなモジュールをPythonで作成しました。モジュールはうまく動作しますが、遅くなります。 6つの製品を各製品の30回の反復まで計算するのに約6時間かかります。だから、私はFORTRANでスクリプトを書き直すことを考えていましたが、もっと速いスピードでミルクを吐くことができないかどうかを確認しました。残念なことに、私はFORTRANを初めて使っています。今ではほとんどのライブラリなどはありません。

pythonのitertools.combinations_with_replacement(pool, r)と同じモジュール/機能がFORTRANにありますか、同じことを達成するモジュールはありますか？

Answer 1

Do not do that。あなたは、アルゴリズムの最適化を行う必要があるときにマイクロ最適化をしようとしています。これが私があなたにitertoolsに基づいてshould not use the exponential-time solutionと主張した理由ですが、better, recursive solutionです。

Answer 2

FWIWのitertools documentationには、_と_置き換え（）の純粋なPythonがあります。それは簡潔であり、Fortranに翻訳するのが難しいはずがありません。

def combinations_with_replacement(iterable, r): 
    # combinations_with_replacement('ABC', 2) --> AA AB AC BB BC CC 
    pool = tuple(iterable) 
    n = len(pool) 
    if not n and r: 
     return 
    indices = [0] * r 
    yield tuple(pool[i] for i in indices) 
    while True: 
     for i in reversed(range(r)): 
      if indices[i] != n - 1: 
       break 
     else: 
      return 
     indices[i:] = [indices[i] + 1] * (r - i) 
     yield tuple(pool[i] for i in indices)