summaryrefslogtreecommitdiff
path: root/fortran/navy_four_levels
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2011-08-23 13:05:52 -0400
committerEugeniy Mikhailov <evgmik@gmail.com>2011-08-23 13:05:52 -0400
commit9214d37cb22ac28fc64c152dd6b222481dbcb351 (patch)
tree92257ba555fbb05044060c543ec920cad3fac89e /fortran/navy_four_levels
parent9d93adcbc78004e9b24b35b08c26c00ea5b7d2af (diff)
downloadNresonances-9214d37cb22ac28fc64c152dd6b222481dbcb351.tar.gz
Nresonances-9214d37cb22ac28fc64c152dd6b222481dbcb351.zip
Added note about time execution estimate
Diffstat (limited to 'fortran/navy_four_levels')
-rw-r--r--fortran/navy_four_levels/README18
1 files changed, 17 insertions, 1 deletions
diff --git a/fortran/navy_four_levels/README b/fortran/navy_four_levels/README
index 42c7e18..0f30a3a 100644
--- a/fortran/navy_four_levels/README
+++ b/fortran/navy_four_levels/README
@@ -1,4 +1,20 @@
-* Eugeniy is grumpy on 07/10/2007:
+* Eugeniy is grumpy on 2011/07/10:
Gosh, this is July 2011 and we are modifying Fortran 77 version to be Fortran 95.
Is there a sense?
+
+* Eugeniy on 2011/08/23:
+Navy chose very bad algorithm.
+It seems to be a good idea to put normilezed dt=dz to avoid some unnesesary calculations
+on each step. But results are devastating since code counter intuitively runs longer for
+a shorter cell, and not surprisingly proportionally to time interval.
+
+Execution time proportional to the total number of grid nodes which can be calculated as
+
+Ngrid = tptz*zpts = (tmax-tmin)*Gamma_super / ((zmax*Gamma_super/(c*zpts)) *zpts
+
+Ngrid = (tmax-tmin)*c/zmax * zpts^2
+
+where 'c' is speed of light, 'zmax' is length of the cell, and 'zpts' is number
+of desired points along z axis.
+