core_shell_microgelscostrafo411magnetic_modelticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change
on this file since ab2aea8 was
c5ac2b2,
checked in by Paul Kienzle <pkienzle@…>, 8 years ago
|
working core multi shell
|
-
Property mode set to
100644
|
File size:
1.0 KB
|
Rev | Line | |
---|
[f7930be] | 1 | |
---|
| 2 | static double |
---|
| 3 | f_constant(double q, double r, double sld) |
---|
| 4 | { |
---|
| 5 | const double bes = sph_j1c(q * r); |
---|
[c5ac2b2] | 6 | const double vol = M_4PI_3 * cube(r); |
---|
[f7930be] | 7 | return sld * vol * bes; |
---|
| 8 | } |
---|
| 9 | |
---|
[c5ac2b2] | 10 | double |
---|
| 11 | form_volume(double core_radius, double n, double thickness[]); |
---|
| 12 | double |
---|
[f7930be] | 13 | form_volume(double core_radius, double n, double thickness[]) |
---|
| 14 | { |
---|
| 15 | double r = core_radius; |
---|
[c5ac2b2] | 16 | for (int i=0; i < n; i++) { |
---|
[f7930be] | 17 | r += thickness[i]; |
---|
| 18 | } |
---|
[c5ac2b2] | 19 | return M_4PI_3 * cube(r); |
---|
[f7930be] | 20 | } |
---|
| 21 | |
---|
[c5ac2b2] | 22 | double |
---|
| 23 | Iq(double q, double core_sld, double core_radius, |
---|
| 24 | double solvent_sld, double num_shells, double sld[], double thickness[]); |
---|
| 25 | double |
---|
| 26 | Iq(double q, double core_sld, double core_radius, |
---|
| 27 | double solvent_sld, double num_shells, double sld[], double thickness[]) |
---|
[f7930be] | 28 | { |
---|
[c5ac2b2] | 29 | const int n = (int)ceil(num_shells); |
---|
| 30 | double f, r, last_sld; |
---|
[f7930be] | 31 | r = core_radius; |
---|
[c5ac2b2] | 32 | last_sld = core_sld; |
---|
| 33 | f = 0.; |
---|
| 34 | for (int i=0; i<n; i++) { |
---|
| 35 | f += M_4PI_3 * cube(r) * (sld[i] - last_sld) * sph_j1c(q*r); |
---|
| 36 | last_sld = sld[i]; |
---|
[f7930be] | 37 | r += thickness[i]; |
---|
| 38 | } |
---|
[c5ac2b2] | 39 | f += M_4PI_3 * cube(r) * (solvent_sld - last_sld) * sph_j1c(q*r); |
---|
| 40 | return f * f * 1.0e-4; |
---|
[f7930be] | 41 | } |
---|
Note: See
TracBrowser
for help on using the repository browser.