c(x, eta)=exp(-beta*delta_E)が計算されるようにしたい。

Question

cの値が1.00となってしまいます。
問題点がなかなか見つかりません。
c(x, eta) =1.00となり、elseが機能していない。なぜかがわからない。

コード
#include <stdio.h>
#include <math.h>
#include <stdlib.h> //for rand() srand()
#define L 4 //total site number
#define N 2 //total particle numer

double swap(int a, int b)
{
    int temp;

    a = temp;
    a = b;
    b = temp;
    return 0;
}

double H(int eta[]) //definition of hamiltonian
{
    int i;
    double Hami = 0.0;

    eta[0] = eta[L];
    for (i=0;i<L;i++) {
        Hami = Hami - eta[i]*eta[i+1];
    }
    if (Hami == 0)
    {
        return 1;
    }
    return Hami;
}

double c(int x, int eta[]) // definition of exchangable probability
{
    double delta_E, beta, temp_H;
    double swap(int, int);
    double H(int eta[]);

    beta = 1.0;
    temp_H = H(eta);
    swap(eta[x], eta[x+1]);
    delta_E = H(eta) - temp_H;  
    if (delta_E < 0) {
        return (1);
    } else
    {
        return (exp(-beta*delta_E));
    }
    
}

int main()
{
    int eta[L+1] = {0};
    int x = 2;

    eta[1] = 1;
    eta[2] = 1;
    printf("%f\t%f\n", H(eta), c(x, eta));
}

Comments

[y_waiwai]

で、しつもんはなんでしょうか

[Okuhira007]

elseが実行されません。どうすればいいでしょうか？

[episteme]

double swap(int a, int b)
{
int temp;
a = temp;
a = b;
b = temp;
}

swapしていないし、戻り値も返していません。
コンパイル・エラーにならないんですか?

[Okuhira007]

コンパイルエラーになりませんでした。

[mjk]

https://paiza.io/projects/48zh1QazUKDKwa4WFNZ81g
プログラム実行途中のデバッグなどは試してみましたか？
プログラム実行途中にprintfで各変数などを出力してみましたが、コードの通り実行されているように思います。
何を求めているプログラムなのかや、期待する結果や途中の値など手元で計算出来る例など提示すると回答が得られやすいと思いました。
H(eta) -1.000000
temp_H -1.000000
delta_E 0.000000
delta_E < 0 : false
exp(-beta*delta_E) 1.000000
H(eta) -1.000000 c(x, eta) 1.000000

[Okuhira007]

ありがとうございます。

Answer 1

答えが合うかどうかは知らないが、
double swap(int *a, int *b)
{
int temp;

temp = *a;
*a = *b;
*b = temp;
return 0;

}
と、
double swap(int *, int *);
と、
swap(&eta[x], &eta[x+1]);
あたりは変える必要がありますね。

Comments

[Okuhira007]

ありがとうございます。直りました！