Swendsen–Wang algortihm for Monte Caro Simulation in Potts Model

2 次查看（过去 30 天）

MINA 2016-2-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/269299-swendsen-wang-algortihm-for-monte-caro-simulation-in-potts-model

编辑： Walter Roberson 2016-2-22

I have implemented this algorithm in matlab by I don't get the result which I am supposed to get. I would appreciate it if someone can tell me where I am making a mistake/s.

W: is k- mutual nearest neighbors matrix

J: interaction strength matrix

n: number of data points

q: # spins

T: temperature

[ S ] = init_S_gen( n,q );   % this generates random q spin values for n data point  
for i=1:15    
      T=0.005*i;
  [ chi(i),Cluster,S ] = Monte_Carlo_sim(W,J,M,n,q,T,S);
  end
function [ chi,Cluster,S ] = Monte_Carlo_sim(W,J,M,n,q,T,S)  
for ii=1:M  
      P=[];  
      [P]=frozen(W,J,S,T);  
      Cluster=[];  
      [ Cluster ] = SW_cluster( P );  
      %%Physical quantities  
      m_S(ii)=gen_magnetization(S,q,n);  
      [S]=  update_S(Cluster,q);      
end  
chi=n/T*(nanmean(m_S.^2)-nanmean(m_S)^2);  
end
function [ Cluster ] = SW_cluster( P )  
%This function finds the SW_clusters  
n=length(P);  
d = rand(n,1); % The diagonal values  
t = triu(bsxfun(@min,d,d.').*rand(n),1); % The upper trianglar random values  
rand_N = diag(d)+t+t.'; % Put them together in a symmetric matrix  
frozen_acc=rand_N<P;  
PP=frozen_acc;  
[nComponents,~,~,Cluster] = networkComponents(PP);   % This function finds the connected components of the graph based on its adjacency matrix.
end
function [m_S]=gen_magnetization(S,q,n)  
%This function computes the m_S, magnetization associated with a spin configuration S  
spin_vec=1:q;  
spin_mat=repmat(spin_vec,[n 1]);  
S_mat=repmat(S',[1 q]);  
N_S=sum((spin_mat== S_mat),1);  
N_max_S=max(N_S);  
m_S=(q*N_max_S-n)/((q-1)*n);  
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

John D'Errico 2016-2-22

编辑：Walter Roberson 2016-2-22

Why do you think there is a mistake? GIVE US A HINT! Do you expect us to:

Research what you want to do, so find someplace that describes the "Swendsen–Wang algortihm for Monte Caro Simulation in Potts Model"
Read your code, trying to understand what you did.
Run your code, debugging any and all bugs in the code.
Provide fixes for the problems that we found

I'm sorry, but I don't have a few hours to waste. When you have a problem, make it easy for someone to help you. Otherwise, the crystal ball gets so foggy.

请先登录，再进行评论。

请先登录，再回答此问题。