1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
|
/* const/test.c
*
* Copyright (C) 2003 Brian Gough
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or (at
* your option) any later version.
*
* This program is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
#include <config.h>
#include <stdlib.h>
#include <math.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_const.h>
#include <gsl/gsl_test.h>
#include <gsl/gsl_ieee_utils.h>
int
main (void)
{
gsl_ieee_env_setup ();
/* Basic check to make sure the header files are functioning */
{
double c = GSL_CONST_MKS_SPEED_OF_LIGHT;
double eps = GSL_CONST_MKS_VACUUM_PERMITTIVITY;
double mu = GSL_CONST_MKS_VACUUM_PERMEABILITY;
gsl_test_rel (c, 1.0/sqrt(eps*mu), 1e-6, "speed of light (mks)");
}
{
double ly = GSL_CONST_CGS_LIGHT_YEAR;
double c = GSL_CONST_CGS_SPEED_OF_LIGHT;
double y = 365.2425 * GSL_CONST_CGS_DAY;
gsl_test_rel (ly, c * y, 1e-6, "light year (cgs)");
}
{
double c = GSL_CONST_MKSA_SPEED_OF_LIGHT;
double eps = GSL_CONST_MKSA_VACUUM_PERMITTIVITY;
double mu = GSL_CONST_MKSA_VACUUM_PERMEABILITY;
gsl_test_rel (c, 1.0/sqrt(eps*mu), 1e-6, "speed of light (mksa)");
}
{
double ly = GSL_CONST_CGSM_LIGHT_YEAR;
double c = GSL_CONST_CGSM_SPEED_OF_LIGHT;
double y = 365.2425 * GSL_CONST_CGSM_DAY;
gsl_test_rel (ly, c * y, 1e-6, "light year (cgsm)");
}
{
double micro = GSL_CONST_NUM_MICRO;
double mega = GSL_CONST_NUM_MEGA;
double kilo = GSL_CONST_NUM_KILO;
gsl_test_rel (mega/kilo, 1/(micro*kilo), 1e-10, "kilo (mega/kilo, 1/(micro*kilo))");
}
{
double d = GSL_CONST_MKSA_DEBYE;
double c = GSL_CONST_MKSA_SPEED_OF_LIGHT;
double desu = d * c * 1000.0;
gsl_test_rel (desu, 1e-18, 1e-10, "debye (esu)");
}
{
double k = GSL_CONST_MKSA_BOLTZMANN;
double c = GSL_CONST_MKSA_SPEED_OF_LIGHT;
double h = GSL_CONST_MKSA_PLANCKS_CONSTANT_H;
double s = 2 * pow(M_PI, 5.0) * pow(k, 4.0) / (15 * pow(c, 2.0) * pow(h, 3.0));
double sigma = GSL_CONST_MKSA_STEFAN_BOLTZMANN_CONSTANT;
gsl_test_rel(s, sigma, 1e-10, "stefan boltzmann constant");
}
exit (gsl_test_summary ());
}
|
