core_shell_microgelscostrafo411magnetic_modelticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change
on this file since 3a48772 was
3a48772,
checked in by Paul Kienzle <pkienzle@…>, 8 years ago
|
use predefined constants for fractions of pi
|
-
Property mode set to
100644
|
File size:
539 bytes
|
Rev | Line | |
---|
[55b283e8] | 1 | double form_volume(void); |
---|
| 2 | |
---|
| 3 | double Iq(double q, double radius2, double arms); |
---|
| 4 | |
---|
| 5 | static double _mass_fractal_kernel(double q, double radius2, double arms) |
---|
| 6 | { |
---|
| 7 | |
---|
[3a48772] | 8 | double u_2 = radius2 * q * q; |
---|
[55b283e8] | 9 | double v = u_2 * arms / (3.0 * arms - 2.0); |
---|
| 10 | |
---|
[3a48772] | 11 | double term1 = v + expm1(-v); |
---|
| 12 | double term2 = ((arms - 1.0)/2.0) * square(expm1(-v)); |
---|
[55b283e8] | 13 | |
---|
[3a48772] | 14 | return (2.0 * (term1 + term2)) / (arms * v * v); |
---|
[55b283e8] | 15 | |
---|
| 16 | } |
---|
| 17 | |
---|
| 18 | double form_volume(void) |
---|
| 19 | { |
---|
| 20 | return 1.0; |
---|
| 21 | } |
---|
| 22 | |
---|
| 23 | double Iq(double q, double radius2, double arms) |
---|
| 24 | { |
---|
| 25 | return _mass_fractal_kernel(q, radius2, arms); |
---|
| 26 | } |
---|
Note: See
TracBrowser
for help on using the repository browser.