core_shell_microgelsmagnetic_modelticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change
on this file since d277229 was
d277229,
checked in by grethevj, 6 years ago
|
Models updated to include choices for effective interaction radii
|
-
Property mode set to
100644
|
File size:
1.5 KB
|
Line | |
---|
1 | |
---|
2 | static double |
---|
3 | f_constant(double q, double r, double sld) |
---|
4 | { |
---|
5 | const double bes = sas_3j1x_x(q * r); |
---|
6 | const double vol = M_4PI_3 * cube(r); |
---|
7 | return sld * vol * bes; |
---|
8 | } |
---|
9 | |
---|
10 | static double |
---|
11 | form_volume(double core_radius, double fp_n, double thickness[]) |
---|
12 | { |
---|
13 | double r = core_radius; |
---|
14 | int n = (int)(fp_n+0.5); |
---|
15 | for (int i=0; i < n; i++) { |
---|
16 | r += thickness[i]; |
---|
17 | } |
---|
18 | return M_4PI_3 * cube(r); |
---|
19 | } |
---|
20 | |
---|
21 | static double |
---|
22 | outer_radius(double core_radius, double fp_n, double thickness[]) |
---|
23 | { |
---|
24 | double r = core_radius; |
---|
25 | int n = (int)(fp_n+0.5); |
---|
26 | for (int i=0; i < n; i++) { |
---|
27 | r += thickness[i]; |
---|
28 | } |
---|
29 | return r; |
---|
30 | } |
---|
31 | |
---|
32 | static double |
---|
33 | effective_radius(int mode, double core_radius, double fp_n, double thickness[]) |
---|
34 | { |
---|
35 | if (mode == 1) { |
---|
36 | double r = core_radius; |
---|
37 | int n = (int)(fp_n+0.5); |
---|
38 | for (int i=0; i < n; i++) { |
---|
39 | r += thickness[i]; |
---|
40 | } |
---|
41 | return r; |
---|
42 | //return outer_radius(core_radius,fp_n,thickness); |
---|
43 | } else { |
---|
44 | return core_radius; |
---|
45 | } |
---|
46 | } |
---|
47 | |
---|
48 | static void |
---|
49 | Fq(double q, double *F1, double *F2, double core_sld, double core_radius, |
---|
50 | double solvent_sld, double fp_n, double sld[], double thickness[]) |
---|
51 | { |
---|
52 | const int n = (int)(fp_n+0.5); |
---|
53 | double f, r, last_sld; |
---|
54 | r = core_radius; |
---|
55 | last_sld = core_sld; |
---|
56 | f = 0.; |
---|
57 | for (int i=0; i<n; i++) { |
---|
58 | f += M_4PI_3 * cube(r) * (sld[i] - last_sld) * sas_3j1x_x(q*r); |
---|
59 | last_sld = sld[i]; |
---|
60 | r += thickness[i]; |
---|
61 | } |
---|
62 | f += M_4PI_3 * cube(r) * (solvent_sld - last_sld) * sas_3j1x_x(q*r); |
---|
63 | *F1 = 1e-2 * f; |
---|
64 | *F2 = 1e-4 * f * f; |
---|
65 | } |
---|
Note: See
TracBrowser
for help on using the repository browser.