core_shell_microgelsmagnetic_modelticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change
on this file since ead25cb was
bdd08df,
checked in by butler, 7 years ago
|
corrects fractal core shell S(q) to be based on radius of full particle
(core radius + thickness) and tested. Fixes #792
|
-
Property mode set to
100644
|
File size:
836 bytes
|
Rev | Line | |
---|
[217590b] | 1 | static double |
---|
| 2 | form_volume(double radius, double thickness) |
---|
[7d4b2ae] | 3 | { |
---|
[3a48772] | 4 | return M_4PI_3 * cube(radius + thickness); |
---|
[7d4b2ae] | 5 | } |
---|
| 6 | |
---|
[217590b] | 7 | static double |
---|
| 8 | Iq(double q, |
---|
| 9 | double radius, |
---|
| 10 | double thickness, |
---|
| 11 | double core_sld, |
---|
| 12 | double shell_sld, |
---|
| 13 | double solvent_sld, |
---|
| 14 | double volfraction, |
---|
| 15 | double fractal_dim, |
---|
| 16 | double cor_length) |
---|
| 17 | { |
---|
[bdd08df] | 18 | //The radius for the building block of the core shell particle that is |
---|
| 19 | //needed by the Teixeira fractal S(q) is the radius of the whole particle. |
---|
| 20 | const double cs_radius = radius + thickness; |
---|
| 21 | const double sq = fractal_sq(q, cs_radius, fractal_dim, cor_length); |
---|
[6d96b66] | 22 | const double pq = core_shell_kernel(q, radius, thickness, |
---|
| 23 | core_sld, shell_sld, solvent_sld); |
---|
[7d4b2ae] | 24 | |
---|
[6d96b66] | 25 | // Note: core_shell_kernel already performs the 1e-4 unit conversion |
---|
| 26 | return volfraction * sq * pq; |
---|
[7d4b2ae] | 27 | } |
---|
| 28 | |
---|
Note: See
TracBrowser
for help on using the repository browser.