126 |
*/ |
*/ |
127 |
#ifdef TEST |
#ifdef TEST |
128 |
|
|
129 |
#include "ideal_impl.h" |
#include "../ideal_impl.h" |
|
|
|
130 |
#include "../test.h" |
#include "../test.h" |
131 |
#include <math.h> |
#include <math.h> |
132 |
#include <stdio.h> |
#include <stdio.h> |
133 |
#include <assert.h> |
#include <assert.h> |
134 |
|
#include <stdlib.h> |
135 |
|
|
136 |
double phi0(double tau, double del){ |
double phi0(double tau, double del){ |
137 |
return log(del) - log(tau) - 12.76953 - 0.007841630*tau + 3.5*log(tau) - 1.934819e-4/tau - 1.247742e-5/(tau*tau) + 6.678326e-8/(tau*tau*tau) + 1.012941*log(1 - exp(-26.65788*tau)); |
return log(del) - log(tau) - 12.76953 - 0.007841630*tau + 3.5*log(tau) - 1.934819e-4/tau - 1.247742e-5/(tau*tau) + 6.678326e-8/(tau*tau*tau) + 1.012941*log(1 - exp(-26.65788*tau)); |