Closed prkkumar closed 1 year ago
This PR implements a random initialization of Polarization for 2D problems. Prior to this we were initializing P using a deterministic approach.
Old way:
double tmp = (i%3 + k%4)/5.; pOld_z(i,j,k) = (-1.0 + 2.0*tmp)*0.002;
New way:
pOld_z(i,j,k) = (-1.0 + 2.0*rng[i + k*n_cell[2]])*0.002;
where rng is a 2D array filled with uniform random numbers between 0 and 1.
rng
0
1
This PR implements a random initialization of Polarization for 2D problems. Prior to this we were initializing P using a deterministic approach.
Old way:
New way:
where
rng
is a 2D array filled with uniform random numbers between0
and1
.