source: sasview/sansmodels/src/c_models/DiamCyl.cpp @ 6e10cff

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload
Last change on this file since 6e10cff was 67424cd, checked in by Mathieu Doucet <doucetm@…>, 13 years ago

Reorganizing models in preparation of cpp cleanup

  • Property mode set to 100644
File size: 3.8 KB
Line 
1/**
2        This software was developed by the University of Tennessee as part of the
3        Distributed Data Analysis of Neutron Scattering Experiments (DANSE)
4        project funded by the US National Science Foundation.
5
6        If you use DANSE applications to do scientific research that leads to
7        publication, we ask that you acknowledge the use of the software with the
8        following sentence:
9
10        "This work benefited from DANSE software developed under NSF award DMR-0520547."
11
12        copyright 2008, University of Tennessee
13 */
14
15/**
16 * Scattering model classes
17 * The classes use the IGOR library found in
18 *   sansmodels/src/libigor
19 *
20 *      TODO: refactor so that we pull in the old sansmodels.c_extensions
21 */
22
23#include <math.h>
24#include "models.hh"
25#include "parameters.hh"
26#include <stdio.h>
27using namespace std;
28
29extern "C" {
30        #include "libStructureFactor.h"
31        #include "DiamCyl.h"
32}
33
34DiamCylFunc :: DiamCylFunc() {
35        radius      = Parameter(20.0, true);
36        radius.set_min(0.0);
37        length      = Parameter(400, true);
38        length.set_min(0.0);
39}
40
41/**
42 * Function to evaluate 1D scattering function
43 * The NIST IGOR library is used for the actual calculation.
44 * @param q: q-value
45 * @return: function value
46 */
47double DiamCylFunc :: operator()(double q) {
48        double dp[2];
49
50        // Fill parameter array for IGOR library
51        // Add the background after averaging
52        dp[0] = radius();
53        dp[1] = length();
54
55        // Get the dispersion points for the radius
56        vector<WeightPoint> weights_rad;
57        radius.get_weights(weights_rad);
58
59        // Get the dispersion points for the length
60        vector<WeightPoint> weights_len;
61        length.get_weights(weights_len);
62
63        // Perform the computation, with all weight points
64        double sum = 0.0;
65        double norm = 0.0;
66
67        // Loop over radius weight points
68        for(size_t i=0; i<weights_rad.size(); i++) {
69                dp[0] = weights_rad[i].value;
70                // Loop over length weight points
71                for(size_t j=0; j<weights_len.size(); j++) {
72                        dp[1] = weights_len[j].value;
73
74                        sum += weights_rad[i].weight
75                                                * weights_len[j].weight *DiamCyl(dp[1], dp[0]);
76                        norm += weights_rad[i].weight
77                                * weights_len[j].weight;
78                }
79        }
80        return sum/norm ;
81}
82
83/**
84 * Function to evaluate 2D scattering function
85 * @param q_x: value of Q along x
86 * @param q_y: value of Q along y
87 * @return: function value
88 */
89double DiamCylFunc :: operator()(double qx, double qy) {
90        double q = sqrt(qx*qx + qy*qy);
91        return (*this).operator()(q);
92}
93
94/**
95 * Function to evaluate 2D scattering function
96 * @param pars: parameters of the cylinder
97 * @param q: q-value
98 * @param phi: angle phi
99 * @return: function value
100 */
101double DiamCylFunc :: evaluate_rphi(double q, double phi) {
102        double qx = q*cos(phi);
103        double qy = q*sin(phi);
104        return (*this).operator()(qx, qy);
105}
106/**
107 * Function to calculate effective radius
108 * @return: effective radius value
109 */
110double DiamCylFunc :: calculate_ER() {
111//NOT implemented yet!!!
112  return 0.0;
113}
114// Testing code
115
116/**
117int main(void)
118{
119        DiamCylFunc c = DiamCylFunc();
120
121        printf("I(Qx=%g,Qy=%g) = %g\n", 0.001, 0.001, c(0.001, 0.001));
122        printf("I(Q=%g) = %g\n", 0.001, c(0.001));
123        c.radius.dispersion = new GaussianDispersion();
124        c.radius.dispersion->npts = 100;
125        c.radius.dispersion->width = 20;
126
127        //c.length.dispersion = GaussianDispersion();
128        //c.length.dispersion.npts = 20;
129        //c.length.dispersion.width = 65;
130
131        //printf("I(Q=%g) = %g\n", 0.001, c(0.001));
132        //printf("I(Q=%g) = %g\n", 0.001, c(0.001));
133        //printf("I(Qx=%g, Qy=%g) = %g\n", 0.001, 0.001, c(0.001, 0.001));
134        //printf("I(Q=%g,  Phi=%g) = %g\n", 0.00447, .7854, c.evaluate_rphi(sqrt(0.00002), .7854));
135
136
137
138        double i_avg = c(0.01, 0.01);
139        double i_1d = c(sqrt(0.0002));
140
141        printf("\nI(Qx=%g, Qy=%g) = %g\n", 0.01, 0.01, i_avg);
142        printf("I(Q=%g)         = %g\n", sqrt(0.0002), i_1d);
143        printf("ratio %g %g\n", i_avg/i_1d, i_1d/i_avg);
144
145
146        return 0;
147}
148 **/
Note: See TracBrowser for help on using the repository browser.