1 | /** |
---|
2 | * Scattering model for a cylinder |
---|
3 | * @author: Mathieu Doucet / UTK |
---|
4 | */ |
---|
5 | |
---|
6 | #include "cylinder.h" |
---|
7 | #include "disperse_cylinder.h" |
---|
8 | #include <math.h> |
---|
9 | #include "libCylinder.h" |
---|
10 | |
---|
11 | double c_disperser( double (*eval)(), double dp[], int n_pars, |
---|
12 | int *idList, double *sigmaList, int n_pts, double q, double phi ); |
---|
13 | |
---|
14 | /** |
---|
15 | * Function to evaluate 1D scattering function |
---|
16 | * @param pars: parameters of the cylinder |
---|
17 | * @param q: q-value |
---|
18 | * @return: function value |
---|
19 | */ |
---|
20 | double disperse_cylinder_analytical_1D(DispCylinderParameters *pars, double q) { |
---|
21 | double dp[5]; |
---|
22 | |
---|
23 | // Fill paramater array |
---|
24 | dp[0] = pars->scale; |
---|
25 | dp[1] = pars->radius; |
---|
26 | dp[2] = pars->length; |
---|
27 | dp[3] = pars->contrast; |
---|
28 | dp[4] = pars->background; |
---|
29 | |
---|
30 | return CylinderForm(dp, q); |
---|
31 | } |
---|
32 | |
---|
33 | double disp_cyl_run(double pars[], double q, double phi) { |
---|
34 | CylinderParameters danse_pars; |
---|
35 | danse_pars.scale = pars[0]; |
---|
36 | danse_pars.radius = pars[1]; |
---|
37 | danse_pars.length = pars[2]; |
---|
38 | danse_pars.contrast = pars[3]; |
---|
39 | danse_pars.background = pars[4]; |
---|
40 | danse_pars.cyl_theta = pars[5]; |
---|
41 | danse_pars.cyl_phi = pars[6]; |
---|
42 | |
---|
43 | return cylinder_analytical_2D(&danse_pars, q, phi); |
---|
44 | |
---|
45 | } |
---|
46 | |
---|
47 | /** |
---|
48 | * Function to evaluate 2D scattering function |
---|
49 | * @param pars: parameters of the cylinder |
---|
50 | * @param q: q-value |
---|
51 | * @return: function value |
---|
52 | */ |
---|
53 | double disperse_cylinder_analytical_2D(DispCylinderParameters *pars, double q, double phi) { |
---|
54 | CylinderParameters cyl_pars; |
---|
55 | double dp[6]; |
---|
56 | int paramList[3]; |
---|
57 | double sigmaList[3]; |
---|
58 | int npts; |
---|
59 | |
---|
60 | |
---|
61 | // Fill paramater struct |
---|
62 | dp[0] = pars->scale; |
---|
63 | dp[1] = pars->radius; |
---|
64 | dp[2] = pars->length; |
---|
65 | dp[3] = pars->contrast; |
---|
66 | dp[4] = pars->background; |
---|
67 | dp[6] = pars->cyl_phi; |
---|
68 | dp[5] = pars->cyl_theta; |
---|
69 | |
---|
70 | paramList[0] = 6; |
---|
71 | paramList[1] = 5; |
---|
72 | paramList[2] = 1; |
---|
73 | |
---|
74 | |
---|
75 | sigmaList[0] = pars->sigma_phi; |
---|
76 | sigmaList[1] = pars->sigma_theta; |
---|
77 | sigmaList[2] = pars->sigma_radius; |
---|
78 | |
---|
79 | npts = (int)(floor(pars->n_pts)); |
---|
80 | |
---|
81 | return c_disperser( &disp_cyl_run, &dp, 3, |
---|
82 | paramList, sigmaList, npts, q, phi ); |
---|
83 | |
---|
84 | |
---|
85 | } |
---|
86 | |
---|