問題タブ [gamma-function]

gammu-smsd を実行するとエラーが発生する

バージョン 1.31.90 の my gammu とバージョン 1.34.0 のデータベース構造。

gammu のバージョンをアップグレードしてみる

ただし、gammu バージョンはバージョン 1.31.90 のままです

この問題を解決するのを手伝ってもらえますか?

ありがとうございました

C ++ 10でtgamma関数を使用しようとしています.32ビットマシンで実行しています。を含めましたが、「identifier tgamma is undefined」というエラーが引き続き表示されます。何が理由でしょうか？

モンテカルロ法を使用して pi (pi = gamma(1/2)**2) を推定するコードを Python で記述しようとしています。sqrt(pi) は、関数 gamma(x=1/2) を囲む正方形の面積に、生成されたランダム ポイントの総数を掛けて、関数内で終了し、ポイントの総数で割ったものであるというコードを書く必要があります。生成されます。

私が話しているコードの種類を示すために、ここではガンマ関数を使用せずにそれを行いました (まだモンテカルロ法です)。ここで Pi は、生成されたポイントの総数で除算された、ユニット 1 の円で終了するポイントの数を 2x2 平方倍した面積として評価されます。ポイントが単位円内にあるための基準は、ピタゴラス (np.sqrt(x x+y y)<=1) です。

Gamma(1/2) を囲む領域をどのように決定するのか、関数内で終了するポイントの基準がここにあるのかどうかはわかりません。

誰にもアイデアがありますか？

ありがとう！

R の次の関数に問題があります。

次の値を挿入すると、次のようになります。

失敗し、NaN. これは、Excel で正確な関数を実装すると、数値ではない結果が得られるためです。J32 は のセルalpha、K32betaと L32 は のセルであるため、Excel での実装は簡単ですN。結果のセルの実装を以下に示します。

したがって、関数はゼロより大きいアルファとベータ、およびゼロ以上のnに対してのみ定義されているため、これは正しい答えを与えるようです。したがって、ここで何が起こっているのだろうか？パッケージ Rmpf を試して数値精度を上げましたが、何もしていないようです。

ありがとう

コンテキスト: Python と sympy を使用して、非常に単純なガンマ関数フラクタル プロットを作成しようとしています。最初は、それがどのように機能するかを理解するための非常に単純なバージョンです (counter=0 または 1 の値に基づく 2 つのマッピング カラー)。

基本的に、コード (以下) はガンマ関数を呼び出し、いくつかの複素数比較を行います。複素数 "nextcomplex=gamma(mycomplex)" が最初の "mycomplex" 複素数よりも "1+0i" に近いことを確認するだけです。 . フラクタルを作成するための最終的なアルゴリズムはそれよりも複雑ですが、基本的な計算はそれらと同じであるため、その単純なコードの速度を向上させる必要があります。

短い間隔では問題なく動作し、値をプロットできますが、大きな間隔では非常に遅くなります。現在実行していますが、1 時間以上実行されており、合計で test_limitn x test_limitm=1000x1000 要素に対して実行されています。(たとえば、100x100 までは問題なく、値をプロットして非常に基本的なフラクタルを確認できます)

私の質問は次のとおりです。コードを強化して高速化するにはどうすればよいですか? (たとえば、他の Python ライブラリ、または比較を行うのにはるかに適した関数が他にあるなど)。

どんなヒントでも大歓迎です、ありがとう！

この C ソース ファイルgamma()を指す関数の R ソース コードを調べていました。

コードのスニペットは次のとおりです。

これで、このコードが、などngamの変数の値を初期化していることがわかりました。しかし、 が何であるかはわかりません。xsmldxrelNOMORE_FOR_THREADS

それが定義されている場所を見つけることができなかったので、いつ#ifdef句に入るか、それ以外の場合は句に入っ#elseてIEEE標準に従って変数を定義するかについて混乱しています。

誰かが定義されているものNOMORE_FOR_THREADSや定義されている場所に光を当てることができれば、それは本当に役に立ちます.

助けてくれてありがとう。

私は問題を試していますが、最初にプロットする必要がありますr

n固定整数です。

このルーチンを作成する方法はありますか?

2 つの不完全なガンマ関数を除算するために必要な方程式に取り組んでいます。IG(a,b) と IG(2a,b) であるため、最初の項が異なるだけです。除算の閉じた形式または単純な近似はありますか?

を使用Gamma(x) = x^(a-1) * e^-xする場合、0 から無限大までの積分を見つけるにはどうすればよいですか?a = 3/2scipy.integrate.quad

これまでのところ、いくつかのコードを試しましたが、「quad: 最初の引数は呼び出し可能ではありません」と言われています。

これは私が使用したコードです：