Skip to content
Snippets Groups Projects
Commit f25a7cd3 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Correctly initialise everything in the non-cosmological case. Also extend the...

Correctly initialise everything in the non-cosmological case. Also extend the range that is probed in testCosmology.c
parent 61e1fa58
No related branches found
No related tags found
1 merge request!979Cosmological integration
...@@ -568,6 +568,8 @@ void cosmology_init_no_cosmo(struct cosmology *c) { ...@@ -568,6 +568,8 @@ void cosmology_init_no_cosmo(struct cosmology *c) {
c->a_begin = 1.; c->a_begin = 1.;
c->a_end = 1.; c->a_end = 1.;
c->log_a_begin = 0.;
c->log_a_end = 0.;
c->log_a_table_begin = 0.; c->log_a_table_begin = 0.;
c->log_a_table_end = 0.; c->log_a_table_end = 0.;
......
...@@ -72,7 +72,7 @@ int main(int argc, char *argv[]) { ...@@ -72,7 +72,7 @@ int main(int argc, char *argv[]) {
message("Start checking computation..."); message("Start checking computation...");
for (int i = 0; i < N_CHECK; i++) { for (int i = 0; i < N_CHECK; i++) {
double a = 0.1 + 0.9 * i / (N_CHECK - 1.); double a = 0.01 + 0.99 * i / (N_CHECK - 1.);
/* Compute a(t(a)) and check if same results */ /* Compute a(t(a)) and check if same results */
double time = cosmology_get_time_since_big_bang(&cosmo, a); double time = cosmology_get_time_since_big_bang(&cosmo, a);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment