added transfer function animations to example notebook;fix for rare r…

…ange bug; comment for exporting time evolution
cmbant · May 15, 2023 · c5b3f6a · c5b3f6a
1 parent 52ed676
commit c5b3f6a
Show file tree

Hide file tree

Showing 9 changed files with 897 additions and 105 deletions.
diff --git a/camb/__init__.py b/camb/__init__.py
@@ -7,7 +7,7 @@
 __author__ = "Antony Lewis"
 __contact__ = "antony at cosmologist dot info"
 __url__ = "https://camb.readthedocs.io"
-__version__ = "1.4.0"
+__version__ = "1.4.1"
 
 from . import baseconfig
 

diff --git a/camb/cambdll.dll b/camb/cambdll.dll
diff --git a/camb/results.py b/camb/results.py
@@ -501,6 +501,10 @@ def get_time_evolution(self, q, eta, vars=model.evolve_names, lAccuracyBoost=4,
         """
         Get the mode evolution as a function of conformal time for some k values.
 
+        Note that gravitational potentials (e.g. Weyl) are not integrated in the code and are
+        calculated as derived parameters; they may be numerically unstable far outside the horizon.
+        (use the series expansion result if needed far outside the horizon)
+
         :param q: wavenumber values to calculate (or array of k values)
         :param eta: array of requested conformal times to output
         :param vars: list of variable names or sympy symbolic expressions to output (using camb.symbolic)

diff --git a/docs/CAMBdemo.html b/docs/CAMBdemo.html
diff --git a/docs/CAMBdemo.ipynb b/docs/CAMBdemo.ipynb
diff --git a/docs/ScalEqs.ipynb b/docs/ScalEqs.ipynb
@@ -2127,7 +2127,7 @@
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "Python 3",
+   "display_name": "Python 3 (ipykernel)",
    "language": "python",
    "name": "python3"
   },
@@ -2141,9 +2141,9 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.7.3"
+   "version": "3.11.0"
   }
  },
  "nbformat": 4,
- "nbformat_minor": 1
+ "nbformat_minor": 4
 }
diff --git a/fortran/camb_python.f90 b/fortran/camb_python.f90
@@ -655,7 +655,7 @@ function CAMB_TimeEvolution(this,nq, q, ntimes, times, noutputs, outputs, &
 
     global_error_flag = 0
     outputs = 0
-    taustart = GetTauStart(maxval(q))
+    taustart = min(times(1),GetTauStart(maxval(q)))
     if (.not. this%ThermoData%HasTHermoData .or. taustart < this%ThermoData%tauminn) call this%ThermoData%Init(this,taustart)
     !$OMP PARALLEL DO DEFAUlT(SHARED),SCHEDUlE(DYNAMIC), PRIVATE(EV, q_ix)
     do q_ix= 1, nq

diff --git a/fortran/config.f90 b/fortran/config.f90
@@ -3,7 +3,7 @@ module config
     use constants, only: const_twopi
     implicit none
 
-    character(LEN=*), parameter :: version = '1.4.0'
+    character(LEN=*), parameter :: version = '1.4.1'
 
     integer :: FeedbackLevel = 0 !if >0 print out useful information about the model
 

diff --git a/forutils b/forutils