Skip to content
Snippets Groups Projects
Commit b351f450 authored by Gizem Özler's avatar Gizem Özler
Browse files

Update electrostatics.f90

parent efb25cc4
Branches
Tags
1 merge request!5Master
Hide whitespace changes
Inline Side-by-side
Showing
with 2 additions and 2 deletions
src/electrostatics.f90
+ 2
− 2
View file @ b351f450
......@@ -397,7 +397,7 @@
call random_number(random)
random_normal=cos(2._pr*pi*random(1))*sqrt(-2._pr*log(random(2)+1.e-12_pr)) ! Box-Muller method
dqp=sigma*random_normal+mu
dqp=sigma*random_normal+mu ! eqn. 1
q_elPast(n)=q_elNext(n)
q_elNext(n)=q_elNext(n)+dqp
......@@ -453,7 +453,7 @@
call random_number(random)
random_normal=cos(2._pr*pi*random(1))*sqrt(-2._pr*log(random(2)+1.e-12_pr)) ! Box-Muller method
dqp=sigma*random_normal+mu
dqp=sigma*random_normal+mu ! eqn. 1
q_elNext(n1)=q_elNext(n1)+dqp
q_elNext(n2)=q_elNext(n2)-dqp
term_prev(n1)=term_prev(n1) + c_rat_prev(n1)*Ac_rat*dqp/(mu+1.e-19_pr) ! eqn. 30
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment