core_shell_microgelsmagnetic_modelticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change
on this file since 7e0b281 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
|
Line | |
---|
1 | static double |
---|
2 | form_volume(double radius, double thickness) |
---|
3 | { |
---|
4 | return M_4PI_3 * cube(radius + thickness); |
---|
5 | } |
---|
6 | |
---|
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 | { |
---|
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); |
---|
22 | const double pq = core_shell_kernel(q, radius, thickness, |
---|
23 | core_sld, shell_sld, solvent_sld); |
---|
24 | |
---|
25 | // Note: core_shell_kernel already performs the 1e-4 unit conversion |
---|
26 | return volfraction * sq * pq; |
---|
27 | } |
---|
28 | |
---|
Note: See
TracBrowser
for help on using the repository browser.