|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
The v2 and v3 were design to be the same except difference of light fields
time profiles in initial conditions. But looks like navy guys during edits
drift both versions from correct realization.
After bug fixes and mergers of insignificant details.
I went peace by peace comparison and removed insignificant differences
* print time elapsed
* print the same other statements to output
* same spaces in the statements so they are not shown as diffs
* same number of t point
* z frames number calculated the same
* L is the same
* density is the same
* remove trailing semicolons to make lines the same
* Omcpeak=0.1 in both files
Then I look at crucial differences
* k1=1-(W32+W31+W34)*dt -> W34 taken in account
* j2=W21*dt
* j8=W41*dt added to v3 to take in account decay of rho44 to rho11
* W42=1 in v2 it used to be zero
* Omc is constant in both cases
* rho14 in v3 had product of two rhos, now it the sum as in v2
* rho24 in v3 was calculated differently form v2: no conjg(Omc) was used
* the fix is taken from v2 and seems to match our xmds code
The only difference remained is output file names.
|