test_thermal.c File Reference

#include <stdarg.h>
#include <stddef.h>
#include <setjmp.h>
#include <cmocka.h>
#include <math.h>
#include <stdio.h>
#include <stdint.h>
#include <stdlib.h>
#include <string.h>
#include "stub.h"
#include <libgadget/config.h>
#include <libgenic/thermal.h>

Include dependency graph for test_thermal.c:

Go to the source code of this file.

Functions
static void	test_mean_velocity (void **state)
static void	test_thermal_vel (void **state)
int	main (void)

Function Documentation

main()

int main

(

void

)

Definition at line 72 of file test_thermal.c.

               {
    const struct CMUnitTest tests[] = {
        cmocka_unit_test(test_mean_velocity),
        cmocka_unit_test(test_thermal_vel)
    };
    return cmocka_run_group_tests_mpi(tests, NULL, NULL);
}

References test_mean_velocity(), and test_thermal_vel().

Here is the call graph for this function:

test_mean_velocity()

static void test_mean_velocity ( void **  state )

static

Definition at line 18 of file test_thermal.c.

{
    /*Check has units of velocity*/
    assert_true(fabs(NU_V0(1, 1, 1e3) - 100*NU_V0(1, 1, 1e5)) < 1e-6);
    /*Check scales linearly with neutrino mass*/
    assert_true(fabs(10*NU_V0(1, 0.1, 1e5) - NU_V0(1, 1, 1e5)) < 1e-6);
    /*Check scales as z (gadget's cosmological velocity unit is accounted for outside)*/
    assert_true(fabs(0.5*NU_V0(0.5, 1, 1e5) -  NU_V0(1, 1, 1e5)) < 1e-6);
}

References NU_V0().

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

test_thermal_vel()

static void test_thermal_vel ( void **  state )

static

Definition at line 29 of file test_thermal.c.

{
    /*Seed table with velocity of 100 km/s*/
    struct thermalvel nu_vels;
    init_thermalvel(&nu_vels, 100, 5000/100, 0);
 
    /*Test getting the distribution*/
    assert_true(fabs(nu_vels.fermi_dirac_vel[0]) < 1e-6);
    assert_true(fabs(nu_vels.fermi_dirac_vel[LENGTH_FERMI_DIRAC_TABLE - 1] -  MAX_FERMI_DIRAC) < 1e-3);
 
    /*Number verified by mathematica*/
    int ii = 0;
    while(nu_vels.fermi_dirac_cumprob[ii] < 0.5) {
        ii++;
    }
    assert_true(fabs(nu_vels.fermi_dirac_vel[ii] - 2.839075) < 0.002);
    /*Check some statistical properties (max, min, mean)*/
    double mean = 0;
    double max = 0;
    double min = 1e10;
    int nsample;
    float Vel[3] = {0};
    int64_t MaxID = 100000;
    gsl_rng * g_rng = gsl_rng_alloc(gsl_rng_ranlxd1);
    for (nsample=0; nsample < MaxID; nsample++)
    {
        add_thermal_speeds(&nu_vels, g_rng, Vel);
        double v2 = sqrt(Vel[0]*Vel[0]+Vel[1]*Vel[1]+Vel[2]*Vel[2]);
        if(v2 > max)
            max = v2;
        if(v2 < min)
            min = v2;
        mean+=v2;
        memset(Vel, 0, 3*sizeof(float));
    }
    gsl_rng_free(g_rng);
    mean/=nsample;
    /*Mean should be roughly 3*zeta(4)/zeta(3)*7/8/(3/4)* m_vamp*/
    assert_true(fabs(mean - 3*pow(M_PI,4)/90./1.202057*(7./8)/(3/4.)*100) < 1);
    assert_true(min > 0);
    assert_true( max < MAX_FERMI_DIRAC*100);
}

References add_thermal_speeds(), thermalvel::fermi_dirac_cumprob, thermalvel::fermi_dirac_vel, init_thermalvel(), LENGTH_FERMI_DIRAC_TABLE, and MAX_FERMI_DIRAC.

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function: