MATLAB Answers

How to add a force to Langevin equation with parfor?

조회 수: 98(최근 30일)
Erez
Erez 24 Dec 2020
댓글: Erez 30 Dec 2020
The code below iterates over a simple Langevin equation, representing Brownian motion of 10^5 particles .
tic
dt=0.1;
Particles=10^5;
x=zeros(1,Particles);
MaxTime=10000 ;
parfor ii=1:MaxTime
Gaussian=sqrt(pi*dt)*normrnd(0,1,[1,Particles]);
x=x+Gaussian;
end
toc
The problem is that I don't know how to add a force term to the equation, in a way which still allows me to use parfor.
I want to do something like this:
tic
dt=0.1;
Particles=10^5;
x=zeros(1,Particles);
Force=zeros(1,Particles);
MaxTime=10000 ;
parfor ii=1:MaxTime
Gaussian=sqrt(pi*dt)*normrnd(0,1,[1,Particles]);
x=x+Force*dt+Gaussian;
Force=-x ;
end
toc
This gives the following error message: "The parfor cannot run due to the way variable 'x' is used".
With a standard 'for' loop of course this is not a problem.
It seems that, even though I know that Matlab is naturally paralellized, when one uses vector & matrices correctly, still for the Brownian motion
parfor seems to be 4 times faster than the standard 'for' loop.
How to do this WITH the force term?

채택된 답변

Raunak Gupta
Raunak Gupta 30 Dec 2020
Hi Erez,
From the code I think parallelization is not possible in this case because if you see the Force variable is updated in every iteration and x is also dependent on the value of Force that is calculated in previous iteration. So, there is a dependency on the sequence of execution of loop. Also, if x is termed as a reduction variable it has to be initialized within parfor loop, which will not be correct as per the logic in code.
You may refer to the following documentations to get information about parfor loops.
  댓글 수: 3
Erez
Erez 30 Dec 2020
Ok, thank you. Yes, it solves most of the mystery. Matlab always uses as many cores as it can at the same time though, to parralelize jobs with vectors. So even without parfor I always see the calculation running on 4-6 cores. But I guess that using parfor cuts the code into blocks in a better way somehow, which works faster (related to Matlab's internal algorithm). So you are right, there's not much more to do beyond it in this case.

댓글을 달려면 로그인하십시오.

추가 답변(0개)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by